[Degrib] Besoin d'aide

[MArcel456]

Bonjour à toutes et à tous,

Je travaille actuellement dans une boite internationale d'informatique, et je dois réaliser un programme assez complexe permettant le téléchargement, la lecture et l'exploitation des données Grib pour prévoir les phénomènes de crue, d'innondation, etc...

Le programme en lui-même est terminé et a été codé en Python. Il permet entre autre, de simuler les données de pluie dans un logiciel (HEC-HMS pour ceux qui connaissent) et d'en resortir automatiquement les prévisions.

Là où je bute, c'est au niveau de Degrib.exe; Je lui fourni les paramètres que je souhaite en option à l'aide de la commande "Find" sous Windows, en lui indiquant les coordonnées de mon projet. A la sortie, j'obtiens un fichier txt contenant plusieurs valeurs, pour une même date de prévision, et ceci après avoir sélectionné l'altitude souhaitée.

Le fichier Grib que je souhaite extraire a été téléchargé depuis le serveur Nomads de la NOAA. Au passage, quelqu'un peut-il me renseigner sur l'heure a indiqué dans le lien du téléchargement? Ce sont les chiifres du style YYYYMMDD00 ou YYYYMMDD06 ou encore YYYYMMDD12. Je sais que ce sont les heures d'emmission des fichiers, mais je souhaitais connaitre le temps de décalage entre eux et moi (Sachant que je suis en France).

De plus, je souhaite récupérer les données de pluie du Grib, quels paramètres dois-je utiliser avec Degrib? Actuellement, j'utilise ACPCP et APCP pour une altitude de 850mb. Or je souhaite la pluie au raz du sol (Je crois que c'est ce qu'il me faut).

Voici la commande que j'utilise pour lancer Degrib:

Degrib.exe gfs.t06z.mastergrb2f72 -P -pnt 48.45,-4.25 -Unit m | find "APCP" | find "850" >> apcp.txt

Malheureusement, le fichier texte contient des données assez étrange (pour moi) du style:

Pour la température:

TMP, [K], 2010012100, 2010012100, 71.23

TMP, [K], 2010012100, 2010012100, 12.01

TMP, [K], 2010012100, 2010012100, 14.30

TMP, [K], 2010012100, 2010012100, -2.74

... et ainsi de suite pas mal de fois.

Merci d'avance à toute personne qui me proposera son aide.

Cordialement.

[jackT]

Malheureusement, le fichier texte contient des données assez étrange (pour moi) ...

Le fichier contient plusieurs valeurs de température , probablement pour des altitudes différentes (niveaux de pression)...

Quelle commande as-tu lancé pour avoir ces données de température ? Celle-ci ?

Degrib.exe gfs.t06z.mastergrb2f72 -P -pnt 48.45,-4.25 -Unit m | find "TMP"  >> tmp.txt

[MArcel456]

Bonjour,

Tout d'abord merci de ta réponse.

Je viens de lancer la première commande (celle du premier message) avec le gfs de ce matin. J'obtiens le résultat suivant:

Degrib.exe gfs.t06z.mastergrb2f72 -P -pnt 48.45,-4.25 -Unit m | find "TMP" | find "850" >> apcp.txt

TMP, [K], 201001220600, 201001250600, -59.850

TMP, [K], 201001220600, 201001250600, -33.850

TMP, [K], 201001220600, 201001250600, 2.850

TMP, [K], 201001220600, 201001250600, 7.850

TMP, [K], 201001220600, 201001250600, -30.850

TMP, [K], 201001220600, 201001250600, 4.850

Avec la tienne, j'obtiens:

Degrib.exe gfs.t06z.mastergrb2f72 -P -pnt 48.45,-4.25 -Unit m | find "TMP"  >> tmp.txt

TMP, [K], 201001220600, 201001250600, -74.950

TMP, [K], 201001220600, 201001250600, -75.450

TMP, [K], 201001220600, 201001250600, -70.250

TMP, [K], 201001220600, 201001250600, -66.550

TMP, [K], 201001220600, 201001250600, -62.550

TMP, [K], 201001220600, 201001250600, -60.450

TMP, [K], 201001220600, 201001250600, -57.550

TMP, [K], 201001220600, 201001250600, -63.650

TMP, [K], 201001220600, 201001250600, -59.850

TMP, [K], 201001220600, 201001250600, -49.050

TMP, [K], 201001220600, 201001250600, -41.350

Et ainsi de suite, pas mal de fois...

J'imagine qu'il s'agit des différents données de températeur pour chaque altitude.

Ce que je souhaite, c'est récupérer les données de pluie (APCP en particulier) sur plusieurs jours (3 jours au moins).

Il faudrait donc que Degrib me ressorte les données au niveau du sol (Surface). Je lui indique ce paramètre directement depuis le serveur nomads lorsque j'envois ma requète.

Cordialement.

[jackT]

Le fichier texte obtenu par une commande degrib telle que

Degrib.exe gfs.t06z.mastergrb2f72 -P -pnt 48.45,-4.25 -Unit m

contient les valeurs de toutes les variables, avec leur nom 5par exemple TMP , APCP...) la date et heure d'initialisation, la date et heure de la prévision, et la valeur du paramètre.Par contre, il n'y a aucune indication du niveau d'altitude (en pression) dans ce fichier. Pour pouvoir y arriver, il faut générer un fichier d'index

Par exemple, pour les précipitations de surface, avec la commande

Degrib.exe gfs.t06z.mastergrb2f72 -I | find "APCP" | find "0-SFC" > apcpindex.txt

le fichier apcpindex.txt contiendra une ligne avec un format comparable à

227.0, 1258597, 2, APCP03="03 hr Total precipitation [kg/(m^2)]", 0-SFC, 01/23/2010 00:00, 01/23/2010 03:00, 3.00

Dans cet exemple, la première valeur de cette ligne (227) indique la position du paramètre APCP en surface dans le fichier de données.

Donc en récupérant toutes les données avec

Degrib.exe gfs.t06z.mastergrb2f72 -P -pnt 48.45,-4.25 -Unit m  > gfs.txt

la 228ème ligne (227 + 1 car la première ligne du fichier contient les entêtes) du fichier gfs.txt contiendra la prévision de pluie en surface :

APCP03, [kg/(m^2)], 201001230000, 201001230300, 0.000

ici 0,000 mmVoilà pour la théorie. Maintenant, si tu n'es intéréssé que par les valeurs de prévisions de précipitations en surface, tu peux utiliser directement

Degrib.exe gfs.t06z.mastergrb2f72 -P -pnt 48.45,-4.25 -Unit m | find "APCP03" >> apcp.txt

qui devrait te donner qu'une seule ligne

APCP03, [kg/(m^2)], 201001230000, 201001230300, 0.000

Dans les données GFS, la prévision pour ACPC03 n'est calculée qu'en surface et il n'y a pas d'autres valeurs pour d'autres niveaux d'altitude

[MArcel456]

Merci beaucoup ça marche!

Bonne journée.

[fiflot]

a supp

Modifié 11 février 2018 par fiflot