Déchiffrage du vent dans fichiers Grib

Bonjour à tous et bonne année 2025 !

J'ai commencé à déchiffrer les fichiers Grib (pour les intégrer dans mon programme de nav) avec wgrib.exe. Après quelques errements, tout va bien, je peux obtenir pression atmosphérique, température,... MAIS : je cale sur le vent, et c'est surtout cela que je veux :-)

Le problème : quelle que soit la source de mon fichier Grib (Zygrib, Xygrib, etc...) et son type (1 ou 2), je ne trouve pour le vent que les données UGRD et VGRD, qui semblent données en m/s. Comment en déduire la vitesse et la direction du vent ? J'y comprend d'autant moins que certaines valeurs sont négatives, ce qui, pour une "vitesse" en m/s est incompréhensible...

Il y a bien un forum (opengribs.org[...] ) mais impossible de s'y inscrire, ça bloque.

Donc, merci aux "sachants" d'HO d'éclairer ma lanterne !

Bon vent à tous,

Michel de Nieuport

L'équipage
#grib#méteo
1180michel1180michel
3j
robertoroberto
3j3j

Ce sont les deux composantes N-S et W-E du vent, il faut les additionner vectoriellement. Pour la vitesse un petit Pythagore, pour l angle un peu de trigonometrie :)
Je ne me souviens plus des conventions mais positif/negatif c est vers N / vers S, ou viceversa; idem W/E

1180michel:Waouw, quelle réponse fulgurante !Merci, Roberto, c'est bien clair et je vais m'y atteler...·le 08 jan. 18:16
Ajouter un commentairePublier
Jeff76 (VCC)Jeff76 (VCC)
3j

Bonjour, bonne année.
Comment et sur quel(s) site(s) récupères-tu les fichiers grib ? Openskiron ? Directement sur Xygrib ? Et sur quel logiciel cherches-tu à les lire ?

1180michel:Je récupère les fichiers Grib avec Zygrib (mais j'ai fait aussi avec Xygrib et directement sur un site météo que j'ai oublié :-) ).Et je les déchiffre avec Wgrib ou Wgrib2.Mais Roberto a tiré plus vite que son ombre, et répondu à ma question !Merci.·le 08 jan. 18:19
Ajouter un commentairePublier
robertoroberto
3j

Suggestion additionnelle qu il ne faudrait pas faire: prends les valeurs de vent d un grib donne, modifie-les d un petit chouia qu on pourrait considerer representatif d'une 'erreur de prevision', un deux trois noeuds en plus ou moins en vitesse, 10 degres en plus ou moins en direction, ou decale les heures.
Fais les routages avec le fichier original, puis avec les 5-10++ fichiers modifies.
Si tu as le temps, essaye la meme chose avec des gribs d'une situation synoptique 'forte' (ex passage d un front), puis dans une situation plus indeterminee (marais barometrique, dorsales fuyantes etc).
Si tu as vraiment beaucoup de temps, apres avoir utilise des polaires 'officielles', modifie aussi celles-ci d quelques dixiemes de noeuds de vitesse.
Compare entre eux tous ces 10-qques dizaines de routages, tous faits avec des donnees 'plausibles'. Enjoy le resultat du long travail, comme on dit :)

1180michel:Merci, Roberto, mais je n'en suis pas encore au routage ! Afficher les vents sur la carte est mon premier objectif, après on verra !Bon vent :-)·le 10 jan. 20:51
Ajouter un commentairePublier
LithopsLithops
2h

As tu essayé d'utiliser gdal ? C'est une bibliothèque très riche, qui te permet en un premier temps de voir la structure des fichiers.
Par exemple, gdalinfo :

gdalinfo Channel12kmWRFWAM250110-12.grb

te donne, entre autres :
Band 426 Block=148x1 Type=Float64, ColorInterp=Undefined
Description = 10[m] HTGL (Specified height level above ground)
Min=-10.670 Max=16.470
Minimum=-10.670, Maximum=16.470, Mean=0.251, StdDev=2.960
Metadata:
GRIBCOMMENT=u-component of wind [m/s]
GRIBELEMENT=UGRD
GRIBFORECASTSECONDS=432000
GRIBREFTIME=1736488800
GRIBSHORTNAME=10-HTGL
GRIBUNIT=[m/s]
GRIBVALIDTIME=1736920800
STATISTICSMAXIMUM=16.47
STATISTICSMEAN=0.25063895417156
STATISTICSMINIMUM=-10.67
STATISTICSSTDDEV=2.9601575390137
STATISTICSVALIDPERCENT=100
Band 427 Block=148x1 Type=Float64, ColorInterp=Undefined
Description = 10[m] HTGL (Specified height level above ground)
Min=-13.060 Max=16.300
Minimum=-13.060, Maximum=16.300, Mean=2.338, StdDev=4.010
Metadata:
GRIBCOMMENT=v-component of wind [m/s]
GRIBELEMENT=VGRD
GRIBFORECASTSECONDS=432000
GRIBREFTIME=1736488800
GRIBSHORTNAME=10-HTGL
GRIBUNIT=[m/s]
GRIBVALIDTIME=1736920800
STATISTICSMAXIMUM=16.3
STATISTICSMEAN=2.3379670975323
STATISTICSMINIMUM=-13.06
STATISTICSSTDDEV=4.0100926463621
STATISTICSVALIDPERCENT=100

Et là te vois que les données de vent sont les composantes u et v, tu vois les unités, ...Et toutes les metadata du fichier. Ensuite il y a d'autres bibliothèques mais il faudra que tu code un peu.

Tags
#voilerie#First#VHF#Cartographie#Voyages#OpenCPN#Dufour#assurances#Gib'Sea#Kelt#Océanis#La Rochelle#TDM#Bavaria#Feeling#Cambuse#accastillage#Manoeuvres#Sun Odyssey#méteo

Phare du monde

  • 4.5 (81)

2022