Prévision GFS

[Meteo-Normandie]

Ah ca me fait plaisir de voir que ca parle un peu technique ici /emoticons/smile@2x.png 2x" width="20" height="20">

get_grib c'est le script qui ne sert qu'a telecharger un petit morceau du gros grib c ca?

Alors attention il ne faut pas se servir uniquement de ce truc car défois ca plante (du coté NOAA).

J'ai un double systeme, si petit fichier donne des infos corrects (cad qu'il est pas vide et qu'il me donne pas du -250°C à 2m du sol par exemple) utiliser petit fichier, sinon allé chercher le gros (et là on mets plus de temps à generer...)

[wydook66]

Bonjour.

Si je souhaite savoir la température à 12h00 pour paris demain, je doit faire comment ?

Je vais sur la ligne GFS 0.5x0.5 degree (5 days), puis la suite je suis perdu /emoticons/biggrin@2x.png 2x" width="20" height="20"> je comprend pas /emoticons/smile@2x.png 2x" width="20" height="20">

Merci

[jeangab]

pour l'instant, je fait un bête file_get_contents sur l'url nomads en php pour la recup des grib, mais on peut utiliser get_gfs.pl , qui est recommandé par les gars du ncep pour éviter les surcharges serveurs (et les plantages assez fréquent dont tu parle meteo-normandie).

en m'inspirant de ce script, je vais voir pour en faire un en php, il suffit d'envoyer les bon content-range pour récupérer que les variable/niveau interessants, voir passer ca ensuite par ggrib pour garder la partie souhaitée en coordonnée.

du coté de cnvgrib, j'ai pas réussi non plus a le compiler , il me manque apparemment une lib params.mod, je n'arrive pas a la trouver (j'ai essayer avec gfortran et g95)

au passage, snowman49, tu utilise quel paramètre pour le flux solaire, j'ai essayé avec DSWRF, mais ca ne m'a pas l'air concluant.

[Meteo-Normandie]

Je recommence pas mes script avec get_gfs sinon j'aurais jamais fini get_grib fait presque pareil wydook t mal barré si tu n'y arrives pas déja à ce niveau là...

Ah compiler les trucs de la NOAA si il y a un truc qui me donne des crises d'urticaire c'est bien ca

[wydook66]

Je recommence pas mes script avec get_gfs sinon j'aurais jamais fini get_grib fait presque pareil wydook t mal barré si tu n'y arrives pas déja à ce niveau là...

Ah compiler les trucs de la NOAA si il y a un truc qui me donne des crises d'urticaire c'est bien ca

Je sait. j'ai jamais utiliser ca et ca fait déjà bien longtemps que je cherche.

Je recherche de l'aide avec les moyens que j'ai , mais j'avance pas

Donc si y en a qui veulent bien m'aider, c'est soit ici soit par MP. Merci beaucoup.

[snowman43]

Avec VFLX jeangab.

Pour la légende de tes cartes, tu as choisis toutes les couleurs ou tu utilises un script gs de grads pour les choisir ?

C'est le plus ch*** dans Grads, faire les légendes.

[wydook66]

Je pense que j'ai avancé un peu.

Je prend la ligne dods et en y mettent mes coordonnées et tout le reste, je me retrouve sur une page du style pour le tmp2

tmp2m, [1][1][1]

[0][0], 286.9

time, [1]

732706.5

lat, [1]

42.0

lon, [1]

60.0

Le problème, je voit pas la température ? ou alors c'est pas en degré ? comment je fait pour l'avoir en degré ?

[snowman43]

Faut mieu utiliser wgrib (wgrib2 pour GFS 0.5°), DODS ca lag et c'est pas à jour pendant 15j parfois !

La température est donnée en Kelvin.

[wydook66]

Faut mieu utiliser wgrib (wgrib2 pour GFS 0.5°), DODS ca lag et c'est pas à jour pendant 15j parfois !

La température est donnée en Kelvin.

J'ai installer cygwin et wgrib, ensuite on doit faire quoi ?

[jeangab]

merci pour la variable,

pour les couleurs, je fixe un degradé avec toute une suite de 'set rgb #coul r v b', puis 'set clevs' avec les niveaux , puis un 'set ccols ' pour les assigner a un niveau,

par ex pour les precipitation:

'set rgb 16 226 222 250'

'set rgb 17 209 204 250'

'set rgb 18 184 176 251'

'set rgb 19 168 157 251'

'set rgb 20 145 130 252'

'set rgb 21 127 111 252'

'set rgb 22 103 84 253'

'set rgb 23 87 65 254'

'set rgb 24 63 37 254'

'set rgb 25 35 1 250'

'set rgb 26 58 1 225'

'set rgb 27 79 1 201'

'set rgb 28 79 1 201'

'set rgb 29 122 1 152'

'set rgb 30 148 1 122'

'set rgb 31 166 0 102'

'set rgb 32 206 0 56'

'set rgb 33 222 0 36'

'set clevs 0.1 0.2 0.5 1 2 5 10 15 20 25 30 35 40 45 50 55 60'

'set ccols 0 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33'

le tout dans un script precipitation.gs

la barre vient du script cbarn.gs, refaite un peu avec photoshop

ensuite, un script php recoupe les images genérée par grads pour ne garder que la carte, et inscrit la date et l'url du site.

je prefere cette solution, car les fontes utilisé par grads ne sont pas des plus terrible /emoticons/biggrin@2x.png 2x" width="20" height="20">

[CFR]

Intéressante, cette discussion /emoticons/biggrin@2x.png 2x" width="20" height="20">

Dans le cadre de mon stage, je récupère une partie des données brutes de GFS sur Nomads également, sur cette page précisément :

http://nomads6.ncdc.noaa.gov/cgi-bin/ftp2u_gfs0.5.sh

Pour ma part, j'utilise NCL, un langage interprêté très pratique pour gérer les grib et afficher les données.

J'avais quelques questions, si certains peuvent répondre, quelqu'un connaît il la signification de toutes les variables, du genre HPBL ou DSWRF, etc... J'ai essayé de chercher une page web qui explicitait ça mais je n'ai pas trouvé.

D'autre part, y-a-t-il un moyen de downloader automatiquement certains fichiers, toutes les 6 heures pour chaque nouveau run, sur un serveur dédié ?

Merci d'avance si certains peuvent répondre. /emoticons/biggrin@2x.png 2x" width="20" height="20">

[snowman43]

Re Bonsoir,

J'ai installer cygwin et wgrib, ensuite on doit faire quoi ?

Tu dl les fichier grib avec wget , wget http.... puis tu fais

wgrib gfs.t00z.pgrbf111 -lon 40.5 5.5 > out.txt

Personnellement j'extrait des données TXT pour l'Europe de l'Ouest de GFS 0.5° (grib2) tous les jours.

Pour JeanGab tu fais comme moi donc, pas très pratique je trouve.

Pour CFR, tu fais avec wgrib

wgrib gfs.t00z.pgrbf111 -v -s > gfs_111.idx

Ca te donnera un fichier de ce genre :

195:25076883:d=2007012900:PRES Pressure [Pa]:surface:123 hour fcst

196:25318164:d=2007012900:PWAT Precipitable water [kg/m^2]:NCEP level type 200 0:123 hour fcst

197:25417912:d=2007012900:RH Relative humidity [%]:NCEP level type 200 0:123 hour fcst

198:25472454:d=2007012900:HGT Geopotential height [gpm]:tropopause:123 hour fcst

199:25810511:d=2007012900:TMP Temp. [K]:tropopause:123 hour fcst

200:25908986:d=2007012900:PRES Pressure [Pa]:tropopause:123 hour fcst

201:26187233:d=2007012900:VWSH Vertical speed shear [1/s]:tropopause:123 hour fcst

202:26357371:d=2007012900:LFTX Surface lifted index [K]:surface:123 hour fcst

203:26454294:d=2007012900:CAPE Convective available potential energy [J/kg]:surface:123 hour fcst

204:26550290:d=2007012900:CIN Convective inhibition [J/kg]:surface:123 hour fcst

205:26624983:d=2007012900:4LFTX Best (4-layer) lifted index [K]:surface:123 hour fcst

206:26724617:d=2007012900:CAPE Convective available potential energy [J/kg]:180-0 mb above ground:123 hour fcst

207:26823010:d=2007012900:CIN Convective inhibition [J/kg]:180-0 mb above ground:123 hour fcst

208:26896704:d=2007012900:HGT Geopotential height [gpm]:max wind:123 hour fcst

209:27331360:d=2007012900:TMP Temp. [K]:max wind:123 hour fcst

210:27512531:d=2007012900:PRES Pressure [Pa]:max wind:123 hour fcst

211:27894591:d=2007012900:HGT Geopotential height [gpm]:surface:123 hour fcst

212:28107033:d=2007012900:PRMSL Pressure reduced to MSL [Pa]:mean sea level:123 hour fcst

213:28197915:d=2007012900:RH Relative humidity [%]:0.44-1 sigma layer:123 hour fcst

214:28273748:d=2007012900:RH Relative humidity [%]:0.72-0.94 sigma layer:123 hour fcst

En l'occurence pour ce qui t'interesse :

242:31677636:d=2007012900:DLWRF Downward longwave radiation flux [W/m^2]:surface:120-123 hour fcst

243:31800530:d=2007012900:ULWRF Upward long wave radiation flux [W/m^2]:surface:120-123 hour fcst

244:31882399:d=2007012900:ULWRF Upward long wave radiation flux [W/m^2]:top of atmosphere:120-123 hour fcst

245:32000208:d=2007012900:USWRF Upward short wave radiation flux [W/m^2]:top of atmosphere:120-123 hour fcst

246:32084623:d=2007012900:USWRF Upward short wave radiation flux [W/m^2]:surface:120-123 hour fcst

247:32143683:d=2007012900:DSWRF Downward shortwave radiation flux [W/m^2]:surface:120-123 hour fcst

272:35619768:d=2007012900:HPBL Planetary boundary layer height [m]:surface:123 hour fcst

D'autre part, y-a-t-il un moyen de downloader automatiquement certains fichiers, toutes les 6 heures pour chaque nouveau run, sur un serveur dédié ?

Oui, il suffit d'envoyer la requête à ftp2u avec curl ou de télécharger les fichiers avec wget.

Pour NCL je suis très interessé par cette librairie, tu aurais un exemple de script que tu utilises pour afficher les données de gfs ?

Merci d'avance.

Bonne soirée.

/emoticons/biggrin@2x.png 2x" width="20" height="20">

[Meteo-Normandie]

Ma fonction qui telecharge un grib en php

path2_provi => dossier destination provisoire des grib

tab_cor contient les variables voulues

path2_getinv et getgrib c'est les chemin vers les script perl

http_server c'est le server de la noaa

path2_inventory c'est le dossier pr les idx

dern_run => je me souviens plus trop la forme exacte, le chemin exact jusqu'au run je crois

forecast=3,6,9,12 etc... c'est run+x heure

d'apres le get_gfs on pourrait ameliorer ca en resserant sur les coord

function download_partial($dern_run,$forecast){		global $path2_provi;		global $tab_cor;		global $path2_getinv;		global $path2_getgrib;		global $http_server;		global $path2_inventory;		if(strlen($forecast)=="1") $forecast="0".$forecast;		$run=substr($dern_run,-2,2);		foreach($tab_cor as $key=>$value){				if($chaine_vars=="") $chaine_vars="$value";				else  $chaine_vars="$value|$chaine_vars";		}		if($forecast<"181"){				$file_search="gfs.t".$run."z.pgrb2f".$forecast;				$file=str_replace("pub/","",$dern_run)."/".$file_search;				$file_inventory="$http_server/$file.idx";				$file_inventory=trim($file_inventory);				if(!fopen($file_inventory,"r")){						add_log("Impossible d'ouvrir [$file_inventory]");						$chaine="wget $http_server/$file --output-document=$path2_provi/$file_search";				}				else {						exec("wget $file_inventory --output-document=$path2_inventory/gfs_master.idx");						$chaine="$path2_getinv $file_inventory | egrep '($chaine_vars)' | $path2_getgrib $http_server/$file $path2_provi/$file_search";				}		}		else {				$file_search="gfs.t".$run."z.pgrbf".$forecast.".grib2";				$file=str_replace("pub/","",$dern_run)."/".$file_search;				$file_inventory=str_replace("grib2","idx",$file);				$chaine="wget $http_server/$file --output-document=$path2_provi/$file_search";		}		//$chaine="$path2_getinv $http_server/$file_inventory | egrep '($chaine_vars)' | $path2_getgrib $http_server/$file $path2_provi/$file_search";		echo "\n==== \n $chaine \n  ==== \n";		add_log("Execution de $chaine");		exec("$chaine");		echo "\n==== \n FIN DL \n  ==== \n";		return $file_search;}

ca retourne le fichier grib enregistré en local, c'est que ca doit servir après /emoticons/biggrin@2x.png 2x" width="20" height="20">

Avant il y a des fonctions pour recuperer ce dernier run, regarder si il a deja été traité etc...

Et on ne se moque de mes noms de fonctions et variables toujours franglaises /emoticons/biggrin@2x.png 2x" width="20" height="20">

[wydook66]

CURL, WGET, PERL, CYGWIN, DEGRIB, GRADS, PCgrib,

Faut installer quoi alors ? je suis perdu de chez perdu là

Je voudrais au moins les grandes ligne pour pouvoir réaliser des prévisions automatique pour ma ville Perpignan.

Avec quoi récupérer les grib automatiquement, comment les décoder, puis comment les envoyer sur ma page web automatiquement, le reste j'essayerais de me débrouiller.

Si quelqu'un peu m'aider en MP ca serait sympa, car la j'y comprend plus rien.

Merci beaucoup.

[snowman43]

Tu récupères les grib avec wget, tu fais une sortie txt avec wgrib et après tu te débrouilles avec ce fichier .txt (tu fais un programme PERL, C, Python ou PHP pour le lire).

[wydook66]

Tu récupères les grib avec wget, tu fais une sortie txt avec wgrib et après tu te débrouilles avec ce fichier .txt (tu fais un programme PERL, C, Python ou PHP pour le lire).

donc pas besoin de degrib ?

Wget faut un logiciel spécial pour l'ouvrir et pouvoir lui dire de récuperer les données sur tel ou tel site ?

désolé de vous ennuyer mais c'est que c'est compliqué tout ca pour quelqu'un qui n'a jamais fait ca à part php.

Merci

[snowman43]

http://www.gnu.org/software/wget/

Ce sont des logiciels ligne de commande, aucune interface pour tout ça.

[wydook66]

http://www.gnu.org/software/wget/

Ce sont des logiciels ligne de commande, aucune interface pour tout ça.

comment on fait pour lui donner à récupérer une url, si on peut rien ouvrir ?

[jeangab]

wget est un genre d'aspirateur a site, en ligne de commande (ou par un script), il faut lui indiquer l'url

la documentation de wget est la http://www.gnu.org/software/wget/manual/

par ex,

wget [url=http://nomad3.ncep.noaa.gov/pub/gfs_master/gfs20070130/gfs.t00z.master.grbf06]http://nomad3.ncep.noaa.gov/pub/gfs_master/gfs20070130/gfs.t00z.master.grbf06[/url]

copie le fichier gfs.t00z.master.grbf06 (complet !) sur ton ddpour cfr, on trouve l'inventaire des differentes variables gfs ici http://www.nco.ncep.noaa.gov/pmb/products/gfs/gfsinfo.shtml

heu, sinon, pour vflx tu est sur snowman? c'est le "momentum flux" en newton/m, qui a rapport avec la vague de gravité d'apres ce que j'ai pu comprendre.

j'ai bien trouver le duvb "uv-b downward solar flux" qui m'a l'air interressant, mais il est dans les fichiers gfs.tXXz.sfluxgrbfXX

[snowman43]

Je m'étais trompé JeanGab (j'avais donné ça de mémoire), j'utilise oPEVPR dans GFS 1° pour le flux solaire (je viens de vérifier).

[wydook66]

J'avais bien compris qu'il s'agit d'un aspirateur. Le problème est de savoir comment ouvrir, je suppose, un de ces fichier de wget pour lui dire de récupérer tel ou tel url. Je ne peut ouvrir aucun fichier, donc je ne peut pas lui indiquer de récupérer une url.

[snowman43]

wget est un logiciel ligne de commande ! tu n'ouvres pas le logiciel avec une petite icône.

En ligne de commande tu tapes wget http:// etc... et wget te récupère le fichier !

[wydook66]

wget est un logiciel ligne de commande ! tu n'ouvres pas le logiciel avec une petite icône.

En ligne de commande tu tapes wget http:// etc... et wget te récupère le fichier !

J'ai réussi à récupérer le fichier, il s'agit d'un fichier GRB je croit. C'est long à téléchargé mais je l'ai enfin 88514 ko le fichier. (je sait pas le temps qu'il doit falloir pour en récupérer plein )

Ensuite je fait quoi ? /emoticons/tongue@2x.png 2x" width="20" height="20">

j'ai télécharger wgrib que j'ai mis dans windows (je pense que c'est ca) et je fait quoi ensuite ?

[snowman43]

Tu tapes la commande que je t'ai fournit en page précédente.

wgrib gfs.t00z.pgrbf111 -lon 40.5 5.5 > out.txt

Si ca marche pas, tu fais wgrib -help pour connaître toutes les commandes dispos ou tu recherches sur le web de la documentation.

[wydook66]

Tu tapes la commande que je t'ai fournit en page précédente.

wgrib gfs.t00z.pgrbf111 -lon 40.5 5.5 > out.txt

Si ca marche pas, tu fais wgrib -help pour connaître toutes les commandes dispos ou tu recherches sur le web de la documentation.

J'ai :

C:/> wgrib -help

could not open file: -help

pourtant j'ai bien mis wgrib dans le dossier windows.

J'ai essayer wget -help et ca marche.