Les données modèles Arome au pas de 0.1° et Arpège 0.2° sont maintenant disponible sur le site de MétéoFrance données publiques.
Il ne s'agit pas ici de refaire un nouveau Metéociel ou Weather4D, mais d'étudier le format GRIB2 pour en extraire les données utiles.
Le but de cette étude est de récupérer les donnes de surface de températures du modèle arome, de les importer dans QGIS et de mettre en valeur les températures de 0 à 3 degrés sur les Alpes, 48h avant l'arrivée d'une perturbation. Tout ça pour représenter l'évolution de la limite pluie neige qui est particulièrement fluctuante sur les Alpes du Sud par flux de Sud-Ouest. La représentation graphique se fera sur des cartes dynamiques 3D.
Installation
Il existe de nombreux extracteurs de fichier Grib2 , le plus simple sur le site de Météo France permet de visualiser les cartes modèles de la France . Il y a aussi Zygrib . Mais ce que je recherche ce sont les données brutes de modèle .
Après recherche, je me suis arrêté sur WGRIB de Climate Prediction Center qui n'a pas présenté de difficulés pour l'installation sur Ubuntu 17.10. Il a été necessaire d'installer les compilateurs FORTRAN et C par les commandes
- sudo apt install GFORTRAN
- sudo apt install GCC
La page pour la compilation se trouve ici et le fichier compressé est à télécharger à l'adresse
ftp://ftp.cpc.ncep.noaa.gov/wd51we/wgrib2/wgrib2.tg
Le fichier décompressé est placé chez moi dans le répertoire grib2 il ne reste plus qu'a lancer la commande make dans un terminal
- cd grib2
- make
Utilisation
Le lancement de l'extracteur se fait en ligne de commande en se placant dans le répertoire d'installatuion
- cd grib2
- ./wgrib2/wgrib2 [suivi du fichier.grib2 et de la liste des options ]
La description des commandes est sur la page http://www.cpc.ncep.noaa.gov/products/wesley/wgrib2/index.html
et les options http://www.cpc.ncep.noaa.gov/products/wesley/wgrib2/short_cmd_list.html
Quelques lignes de commandes utilisées. pour extraire les températures
- Commande pour lister les metadonnées du fichier AROME_0.01_HP1_02H_201802190300.grib2 téléchargé sur le site météo france.
l'option -v donne une description plus détaillée
knobuntu@knobuntu:~/grib2$ ./wgrib2/wgrib2 AROME_0.01_HP1_02H_201802190300.grib2 -v - Commande de réduction du fichier grib2 pour ne garder que les températures
knobuntu@knobuntu:~/grib2$ ./wgrib2/wgrib2 AROME_0.01_HP1_02H_201802190300.grib2 -match "TMP" -grib tempe.grib2 - Option pour limiter l'emprise de la carte etre -10 et +10 en longitude et 40 à 50 de latitude pour les Alpes
knobuntu@knobuntu:~/grib2$ ./wgrib2/wgrib2 tempe.grib2 -set_grib_type C2 -small _grib -10:10 40:50 temperegion.grib - Transformation au format netCDF: Qgis a une extension crayfish mais dans la version 2.14 cette extension ne marche pas . J'enb ai trouvé une autre qui lit les fichiers netCDF (netCDF browser).
knobuntu@knobuntu:~/grib2$ ./wgrib2/wgrib2 temperegion.grib2 -netcdf tempehh.nc
Il ne reste plus qu'à automatiser tout ça....
Importation dans QGIS
On importe le fichier tempehh.nc dans QGis par l'extension netCDF browser. On initialise qGIs avec la projection Lambert 2154 en activant la projection à la volée.
La suite consiste à définir le style de couleur, polygoniser le raster, isoler les températures 0 à 3 degré , importer un raster relif alti.dem sur le site de la NASA, et les cartes openstreetmap puis exporter au format html avec QGIS2THREEJS.
voili voilou... je vais mettre tout ça au point...
Les premiers résultats
construction de 3 classes de températures d'une carte arome prévue à échéance de 3h
- En gris les températures négatives
- en orange températures 0 à 3 degré. C''est dans cettte couche que l'on va trouver la limite pluie neige
- en vert : températures de plus de 3 degrés dans cette zone tout tombe sous forme de pluie.
Le fond de carte est une couche d'élévation pour la représentation 3D et la détermination de l'altitude.
On va voir par la suite les très fortes variations de la limite pluie neige en altitude et les variations des zones pluvieuses à l'approche d'une perturbation de Sud-Ouest sur les Alpes du Sud. Le modèle AROME est de loin le plus précis pour la prévision de l'iso 0 ° . Une autre approche consiste à extraire la température potentielle 1 degré mais je ne sais pas si elle est disponible.
