WRF EMS - ems_install.pl - problème installation

[kevorage]

Avant d'utiliser le script de Kumar je te conseil vivement de lire la documentation sur grads dispo sur le net.

Essaye également de réaliser des cartes simples puis améliore les tranquillement en regardant la doc et en essayant les commandes disponibles.

Par exemple

'reinit''open lien_vers_ton_fichier_ctl''set background 0''c''set lat 46.09 47.02''set lon 5.91 7.16''set lev 1''set t 1''set display color white''set timelab off''d tmp2m -273.15''draw title Temperatures surface (2m)''printim tmp2m.png x800 y800'

A enregistrer dans un fichier avec l'extension .gs. N'oublies pas de modifier le lien vers le fichier .ctl et la latitude et longitude.

Ensuite tu lance grads en mode console et tu exécute le script. Si tout se passe bien tu aura une image nommée tmp2m.png enregistrée dans ton répertoire courant.

Merci, je vais essayer tout ceci. Moi qui pensais que tout était automatisé par l'autorun... Raté

EDIT : Tout ne se passe pas bien à priori ...

Pfiouuuuu...

Edit 2 : Je ne me suis pas encore penché sur la doc. Faut que j'y travaille, mais tout en Anglais, c'est titanesque.

[Tintin87]

Bonjour, ayant quelques connaissances sur GrADS :

set error, no files open yet

ça veut dire que dans ta première ligne de ton fichier .gs, GrADS n'a pas trouvé le fichier.ctl.Montre ton script gs pour voir d'où vient l'erreur.

Sinon, pour exécuter un script grads en mode console il faut faire :

run /home/.../script.gs

[ZHB]

Effectivement, le fichier ctl ne semble pas trouvé. Ensuite j'ai l'impression que tu essaye beaucoup de trucs sans trop savoir ce que tu fais !

Copies ceci et colle le dans ton fichier "test.gs"

'reinit''open /home/administrateur/wrfems/runs/MC/emsprd/grads/201108221000_arw_wrfout__d01.ctl''set background 0''set lat 45.47 48.50''set lon 0.85 4.60''set t 1''set display color white''set timelab off''d tmp2m -273.15''draw title Temperatures surface (2m)''printim tmp2m.png x800 y800'

Ensuite fais ceci pour lancer ton script grads :

# cd /home/administrateur/wrfems/runs/MC/

# grads -b

ga -> run test.gs

[ZHB]

Effectivement, le fichier ctl ne semble pas trouvé. Ensuite j'ai l'impression que tu essaye beaucoup de trucs sans trop savoir ce que tu fais !

Copies ceci et colle le dans ton fichier "test.gs"

'reinit''open /home/administrateur/wrfems/runs/MC/emsprd/grads/201108221000_arw_wrfout__d01.ctl''set background 0''set lat 45.47 48.50''set lon 0.85 4.60''set t 1''set display color white''set timelab off''d tmp2m -273.15''draw title Temperatures surface (2m)''printim tmp2m.png x800 y800'

Ensuite fais ceci pour lancer ton script grads :

# cd /home/administrateur/wrfems/runs/MC/

# grads -b

ga -> run test.gs

Et même si la doc est en anglais faut vraiment prendre la peine de lire la base.

Oups, je voulais faire modifier...

[kevorage]

Effectivement, le fichier ctl ne semble pas trouvé. Ensuite j'ai l'impression que tu essaye beaucoup de trucs sans trop savoir ce que tu fais !

Ce n'est pas une impression, j'y vais totalement à tâton depuis que j'ai réussi à créer mon domaine et à lancer un run. J'étais tellement content, et je pensais que j'avais fait le plus difficile

Donc effectivement, je n'ai pas tenté d'étudier la doc de GrADS aujourd'hui, car je pensais que l'autorun gérait la génération des cartes.

J'ai fait ce que tu m'as dis. Pour te prouver ma bonne foi, je comprend que ce que fais GrADS ici s'explique comme-suit :

Il examine le script en GS, qui le renvoie à un fichier CTL qui lui donne toutes les infos nécessaires à son fonctionnement.

Voici ce que j'obtiens en suivant ta manipulation, point par point et à la lettre :

Scanning description file: test.gs

Open Error: missing XDEF record

Open Error: missing YDEF record

Open Error: missing ZDEF record

Open Error: missing TDEF record

Open Error: missing UNDEF record

Open Error: missing DSET record

Open Error: missing VARS record

The data file was not opened.

[ZHB]

Je pense à une mauvaise configuration de ton domaine. Montre le contenu de ton fichier "201108221000_arw_wrfout__d01.ctl".

[kevorage]

Je pense à une mauvaise configuration de ton domaine. Montre le contenu de ton fichier "201108221000_arw_wrfout__d01.ctl".

Voici :

dset ^201108221200_arw_wrfout_d01.grb2index ^201108221200_arw_wrfout_d01.idxundef 9.999E+20title 201108221200_arw_wrfout_d01.grb2* produced by g2ctl v0.0.7* command line options: -verf -ts1hr 201108221200_arw_wrfout_d01.grb2 201108221200_arw_wrfout_d01.idx* griddef=1:0:(99 x 81):grid_template=0:winds(grid): lat-lon grid:(99 x 81) units 1e-06 input WE:SN output WE:SN res 56 lat 45.466000 to 48.588000 by 0.039000 lon 0.824000 to 4.650000 by 0.039000 #points=8019:winds(grid)dtype grib2ydef 81 linear 45.466000 0.039025xdef 99 linear 0.824000 0.0390408163265306tdef 25 linear 12Z22aug2011 1hr* PROFILE hPazdef 39 levels 101300 100000 97500 95000 92500 90000 87500 85000 82500 80000 77500 75000 72500 70000 67500 65000 62500 60000 57500 55000 52500 50000 47500 45000 42500 40000 37500 35000 32500 30000 27500 25000 22500 20000 15000 10000 5000 2500 1000options pascalsvars 142no4LFTX180_0mb  0,108,18000,0   0,7,193 ** 180-0 mb above ground Best (4 layer) Lifted Index [K]ACGRPLsfc   0,1,0   0,1,251,1 ** surface Accumulated grid scale Graupel [mm]ACHAILsfc   0,1,0   0,1,250,1 ** surface Accumulated grid scale Hail [mm]ACPCPsfc   0,1,0   0,1,10,1 ** surface Convective Precipitation [kg/m^2]APCPsfc   0,1,0   0,1,8,1 ** surface Total Precipitation [kg/m^2]CAPEsfc   0,1,0   0,7,6 ** surface Convective Available Potential Energy [J/kg]CAPE180_0mb  0,108,18000,0   0,7,6 ** 180-0 mb above ground Convective Available Potential Energy [J/kg]CAPE90_0mb  0,108,9000,0   0,7,6 ** 90-0 mb above ground Convective Available Potential Energy [J/kg]CAPE255_0mb  0,108,25500,0   0,7,6 ** 255-0 mb above ground Convective Available Potential Energy [J/kg]CFRZRsfc   0,1,0   0,1,193 ** surface Categorical Freezing Rain (yes=1; no=0) [non-dim]CICEPsfc   0,1,0   0,1,194 ** surface Categorical Ice Pellets (yes=1; no=0) [non-dim]CINsfc   0,1,0   0,7,7 ** surface Convective Inhibition [J/kg]CIN180_0mb  0,108,18000,0   0,7,7 ** 180-0 mb above ground Convective Inhibition [J/kg]CIN90_0mb  0,108,9000,0   0,7,7 ** 90-0 mb above ground Convective Inhibition [J/kg]CIN255_0mb  0,108,25500,0   0,7,7 ** 255-0 mb above ground Convective Inhibition [J/kg]CPRATsfc   0,1,0   0,1,196 ** surface Convective Precipitation Rate [kg/m^2/s]CRAINsfc   0,1,0   0,1,192 ** surface Categorical Rain (yes=1; no=0) [non-dim]CSNOWsfc   0,1,0   0,1,195 ** surface Categorical Snow (yes=1; no=0) [non-dim]DPT2m   0,103,2   0,0,6 ** 2 m above ground Dew Point Temperature [K]DPThy1   0,105,1   0,0,6 ** 1 hybrid level Dew Point Temperature [K]DPThy2   0,105,2   0,0,6 ** 2 hybrid level Dew Point Temperature [K]DPT30_0mb  0,108,3000,0   0,0,6 ** 30-0 mb above ground Dew Point Temperature [K]ELONsfc   0,1,0   0,191,193 ** surface East Longitude (0 - 360) [deg]GRPLMXsfc   0,1,0   0,1,245 ** surface Maximum column integrated graupel [kg/m^2]GUSTsfc   0,1,0   0,2,22 ** surface Wind Speed (Gust) [m/s]GUSTMX10m   0,103,10   0,2,248 ** 10 m above ground Maximum 10m wind gust between output time [m/s]HCDChcll   0,234,0   0,6,5 ** high cloud layer High Cloud Cover [%]HGTsfc   0,1,0   0,3,5 ** surface Geopotential Height [gpm]HGTprs    39,100  0,3,5 ** (1013 1000 975 950 925.. 150 100 50 25 10) Geopotential Height [gpm]HGThy1   0,105,1   0,3,5 ** 1 hybrid level Geopotential Height [gpm]HGThy2   0,105,2   0,3,5 ** 2 hybrid level Geopotential Height [gpm]HGTclb   0,2,0   0,3,5 ** cloud base Geopotential Height [gpm]HGTtop0C   0,204,0   0,3,5 ** highest tropospheric freezing level Geopotential Height [gpm]HGTceil   0,215,0   0,3,5 ** cloud ceiling Geopotential Height [gpm]HGTlwb0   0,245,0   0,3,5 ** lowest level of the wet bulb zero Geopotential Height [gpm]HGTclt   0,3,0   0,3,5 ** cloud top Geopotential Height [gpm]HGT0C   0,4,0   0,3,5 ** 0C isotherm Geopotential Height [gpm]HGTl5   0,5,0   0,3,5 ** level of adiabatic condensation from sfc Geopotential Height [gpm]HLCY3000_0m  0,103,3000,0   0,7,8 ** 3000-0 m above ground Storm Relative Helicity [m^2/s^2]HLCY1000_0m  0,103,1000,0   0,7,8 ** 1000-0 m above ground Storm Relative Helicity [m^2/s^2]HPBLsfc   0,1,0   0,3,196 ** surface Planetary Boundary Layer Height [m]ICECsfc   0,1,0   10,2,0 ** surface Ice Cover [Proportion]LANDsfc   0,1,0   2,0,0 ** surface Land Cover (1=land, 0=sea) [Proportion]LCDClcll   0,214,0   0,6,3 ** low cloud layer Low Cloud Cover [%]LFTXl100_100  0,100,50000,100000   0,7,192 ** 500-1000 mb Surface Lifted Index [K]MAXRH2m   0,103,2   0,1,253 ** 2 m above ground Maximum Relative Humidity [%]MCDCmcll   0,224,0   0,6,4 ** middle cloud layer Medium Cloud Cover [%]MINRH2m   0,103,2   0,1,252 ** 2 m above ground Minimum Relative Humidity [%]MSLETmsl   0,101,0   0,3,192 ** mean sea level MSLP (Eta model reduction) [Pa]MSTAV0_100cm  0,106,0,1   2,0,194 ** 0-1 m below ground Moisture Availability [%]MX10SP10m   0,103,10   0,2,245 ** 10 m above ground Maximum 10m wind speed between output times [m/s]MX10U10m   0,103,10   0,2,249 ** 10 m above ground U-Comp of Maximum 10m wind speed between output times [m/s]MX10V10m   0,103,10   0,2,250 ** 10 m above ground V-Comp of Maximum 10m wind speed between output times [m/s]MXCPRTsfc   0,1,0   0,1,248 ** surface Maximum cumulus total precipitation rate [mm/hr]MXGPRTsfc   0,1,0   0,1,249 ** surface Maximum grid scale total precipitation rate [mm/hr]MXSPRTsfc   0,1,0   0,1,246 ** surface Maximum grid scale snow/ice precipitation rate [Liq. Eqiv.] [mm/hr]MXTPRTsfc   0,1,0   0,1,247 ** surface Maximum total precipitation rate [mm/hr]NCPCPsfc   0,1,0   0,1,9,1 ** surface Large-Scale Precipitation (non-convective) [kg/m^2]NLATsfc   0,1,0   0,191,192 ** surface Latitude (-90 to +90) [deg]PLI30_0mb  0,108,3000,0   0,7,0 ** 30-0 mb above ground Parcel Lifted Index (to 500 mb) [K]PRATEsfc   0,1,0   0,1,7 ** surface Precipitation Rate [kg/m^2/s]PRESsfc   0,1,0   0,3,0 ** surface Pressure [Pa]PREShy1   0,105,1   0,3,0 ** 1 hybrid level Pressure [Pa]PREShy2   0,105,2   0,3,0 ** 2 hybrid level Pressure [Pa]PRESclb   0,2,0   0,3,0 ** cloud base Pressure [Pa]PRESclt   0,3,0   0,3,0 ** cloud top Pressure [Pa]PRESl5   0,5,0   0,3,0 ** level of adiabatic condensation from sfc Pressure [Pa]PRMSLmsl   0,101,0   0,3,1 ** mean sea level Pressure Reduced to MSL [Pa]PWAT30_0mb  0,108,3000,0   0,1,3 ** 30-0 mb above ground Precipitable Water [kg/m^2]PWATclm   0,200,0   0,1,3 ** entire atmosphere (considered as a single layer) Precipitable Water [kg/m^2]REFCclm   0,200,0   0,16,196 ** entire atmosphere (considered as a single layer) Composite reflectivity [dB]REFCMXsfc   0,1,0   0,16,246 ** surface Maximum simulated composite reflectivity between output times [dbZ]REFDprs    34,100  0,16,195 ** (1013 1000 975 950 925.. 300 275 250 225 200) Reflectivity [dB]REFD4000m   0,103,4000   0,16,195 ** 4000 m above ground Reflectivity [dB]REFD1000m   0,103,1000   0,16,195 ** 1000 m above ground Reflectivity [dB]REFDMX1000m   0,103,1000   0,16,245 ** 1000 m above ground Maximum simulated radar reflectivity between output times [dbZ]REFDMX3000m   0,103,3000   0,16,245 ** 3000 m above ground Maximum simulated radar reflectivity between output times [dbZ]RHprs    39,100  0,1,1 ** (1013 1000 975 950 925.. 150 100 50 25 10) Relative Humidity [%]RH2m   0,103,2   0,1,1 ** 2 m above ground Relative Humidity [%]RH30_0mb  0,108,3000,0   0,1,1 ** 30-0 mb above ground Relative Humidity [%]SNODsfc   0,1,0   0,1,11 ** surface Snow Depth [m]SNOMsfc   0,1,0   0,1,16 ** surface Snow Melt [kg/m^2]SNOWCsfc   0,1,0   0,1,201 ** surface Snow Cover [%]SOILL0_10cm  0,106,0,0.1   2,3,192 ** 0-0.1 m below ground Liquid Volumetric Soil Moisture (non Frozen) [Proportion]SOILM0_200cm  0,106,0,2   2,0,3 ** 0-2 m below ground Soil Moisture Content [kg/m^2]SOILW0_10cm  0,106,0,0.1   2,0,192 ** 0-0.1 m below ground Volumetric Soil Moisture Content [Fraction]SOILW10_40cm  0,106,0.1,0.4   2,0,192 ** 0.1-0.4 m below ground Volumetric Soil Moisture Content [Fraction]SOILW40_100cm  0,106,0.4,1   2,0,192 ** 0.4-1 m below ground Volumetric Soil Moisture Content [Fraction]SOILW100_200cm  0,106,1,2   2,0,192 ** 1-2 m below ground Volumetric Soil Moisture Content [Fraction]SOTYPsfc   0,1,0   2,3,0 ** surface Soil Type [-]SPFHprs    39,100  0,1,0 ** (1013 1000 975 950 925.. 150 100 50 25 10) Specific Humidity [kg/kg]SPFH2m   0,103,2   0,1,0 ** 2 m above ground Specific Humidity [kg/kg]SPFHhy1   0,105,1   0,1,0 ** 1 hybrid level Specific Humidity [kg/kg]SPFHhy2   0,105,2   0,1,0 ** 2 hybrid level Specific Humidity [kg/kg]SPFH30_0mb  0,108,3000,0   0,1,0 ** 30-0 mb above ground Specific Humidity [kg/kg]SRWEQsfc   0,1,0   0,1,12 ** surface Snowfall Rate Water Equivalent [kg/m^2/s]SSRUNsfc   0,1,0   1,0,193,1 ** surface Storm Surface Runoff [kg/m^2]TCDCprs    34,100  0,6,1 ** (1013 1000 975 950 925.. 300 275 250 225 200) Total Cloud Cover [%]TCDCclm   0,200,0   0,6,1 ** entire atmosphere (considered as a single layer) Total Cloud Cover [%]TCOLCclm   0,200,0   0,6,198 ** entire atmosphere (considered as a single layer) Total Column-Integrated Condensate [kg/m^2]TCOLIclm   0,200,0   0,6,197 ** entire atmosphere (considered as a single layer) Total Column-Integrated Cloud Ice [kg/m^2]TCOLRclm   0,200,0   0,1,204 ** entire atmosphere (considered as a single layer) Total Column Integrated Rain [kg/m^2]TCOLSclm   0,200,0   0,1,205 ** entire atmosphere (considered as a single layer) Total Column Integrated Snow [kg/m^2]TCOLWclm   0,200,0   0,6,196 ** entire atmosphere (considered as a single layer) Total Column-Integrated Cloud Water [kg/m^2]TMAX2m   0,103,2   0,0,4 ** 2 m above ground Maximum Temperature [K]TMIN2m   0,103,2   0,0,5 ** 2 m above ground Minimum Temperature [K]TMPsfc   0,1,0   0,0,0 ** surface Temperature [K]TMPprs    39,100  0,0,0 ** (1013 1000 975 950 925.. 150 100 50 25 10) Temperature [K]TMP2m   0,103,2   0,0,0 ** 2 m above ground Temperature [K]TMPhy1   0,105,1   0,0,0 ** 1 hybrid level Temperature [K]TMPhy2   0,105,2   0,0,0 ** 2 hybrid level Temperature [K]TMP30_0mb  0,108,3000,0   0,0,0 ** 30-0 mb above ground Temperature [K]TMPclt   0,3,0   0,0,0 ** cloud top Temperature [K]TSOIL0_10cm  0,106,0,0.1   2,0,2 ** 0-0.1 m below ground Soil Temperature [K]TSOIL10_40cm  0,106,0.1,0.4   2,0,2 ** 0.1-0.4 m below ground Soil Temperature [K]TSOIL40_100cm  0,106,0.4,1   2,0,2 ** 0.4-1 m below ground Soil Temperature [K]TSOIL100_200cm  0,106,1,2   2,0,2 ** 1-2 m below ground Soil Temperature [K]TSOIL300cm   0,106,3   2,0,2 ** 3 m underground Soil Temperature [K]UGRDprs    39,100  0,2,2 ** (1013 1000 975 950 925.. 150 100 50 25 10) U-Component of Wind [m/s]UGRD10m   0,103,10   0,2,2 ** 10 m above ground U-Component of Wind [m/s]UGRDhy1   0,105,1   0,2,2 ** 1 hybrid level U-Component of Wind [m/s]UGRDhy2   0,105,2   0,2,2 ** 2 hybrid level U-Component of Wind [m/s]UGRD30_0mb  0,108,3000,0   0,2,2 ** 30-0 mb above ground U-Component of Wind [m/s]UPHLsfc   0,1,0   0,7,197 ** surface Updraft Helicity [m^2/s^2]UPHLMXsfc   0,1,0   0,7,246 ** surface Maximum Updraft Helicity between output times [m^2/s^2]USTM6000_0m  0,103,6000,0   0,2,194 ** 6000-0 m above ground U-Component Storm Motion [m/s]VEGsfc   0,1,0   2,0,4 ** surface Vegetation [%]VGRDprs    39,100  0,2,3 ** (1013 1000 975 950 925.. 150 100 50 25 10) V-Component of Wind [m/s]VGRD10m   0,103,10   0,2,3 ** 10 m above ground V-Component of Wind [m/s]VGRDhy1   0,105,1   0,2,3 ** 1 hybrid level V-Component of Wind [m/s]VGRDhy2   0,105,2   0,2,3 ** 2 hybrid level V-Component of Wind [m/s]VGRD30_0mb  0,108,3000,0   0,2,3 ** 30-0 mb above ground V-Component of Wind [m/s]VGTYPsfc   0,1,0   2,0,198 ** surface Vegetation Type [integer(0-13)]VISsfc   0,1,0   0,19,0 ** surface Visibility [m]VISclb   0,2,0   0,19,0 ** cloud base Visibility [m]VSTM6000_0m  0,103,6000,0   0,2,195 ** 6000-0 m above ground V-Component Storm Motion [m/s]VVELprs    39,100  0,2,8 ** (1013 1000 975 950 925.. 150 100 50 25 10) Vertical Velocity (Pressure) [Pa/s]VWSH9800_0m  0,103,9800,0   0,2,192 ** 9800-0 m above ground Vertical speed sheer [1/s]WDNMXsfc   0,1,0   0,2,247 ** surface Maximum downward vertical motion between output times [m/s]WEASDsfc   0,1,0   0,1,13,1 ** surface Water Equivalent of Accumulated Snow Depth [kg/m^2]WTMPsfc   0,1,0   10,3,0 ** surface Water Temperature [K]WUPMXsfc   0,1,0   0,2,246 ** surface Maximum upward vertical motion between output times [m/s]ENDVARS

[Tintin87]

D'ailleurs, sur une région relativement petite, il met 9H pour générer 24H de prévis ...

Utilise le modèle NMM au lieu du ARW, tu devrais gagner 25% de temps.Autres paramètres très importants : les paramètres microphysiques : certains paramètres augmente le temps de 200% !

Scanning description file: test.gsOpen Error: missing XDEF recordOpen Error: missing YDEF recordOpen Error: missing ZDEF recordOpen Error: missing TDEF recordOpen Error: missing UNDEF recordOpen Error: missing DSET recordOpen Error: missing VARS recordThe data file was not opened.

c'est ce que tu obtiens quand tu tape en mode console ?

pour lancer GrADS

grads

puis une fois dans GrADS

'open /home/administrateur/wrfems/runs/MC/emsprd/grads/201108221000_arw_wrfout__d01.ctl'

[kevorage]

Utilise le modèle NMM au lieu du ARW, tu devrais gagner 25% de temps.

Autres paramètres très importants : les paramètres microphysiques : certains paramètres augmente le temps de 200% !

En réalité, je ne me suis pas encore penché sur la configuration plus poussée, encore une fois, j'ai tellement eu du mal à le mettre en place...

Si j'arrive enfin à générer une carte, j'attaque la configuration poussée, et je dump tout ça sur un disque externe

c'est ce que tu obtiens quand tu tape en mode console ?

pour lancer GrADS

grads

puis une fois dans GrADS

'open /home/administrateur/wrfems/runs/MC/emsprd/grads/201108221000_arw_wrfout__d01.ctl'

Si je réalise cette manipulation, avec donc le serveur X qui tourne, GrADS me dit qu'il a ouvert le fichier en tant que fichier 1, et il récupère latitude, longitude, le LEV ainsi que les valeurs de temps.

[Tintin87]

GrADS me dit qu'il a ouvert le fichier en tant que fichier 1

Donc il ouvre bien ton fichier ctlEt si ensuite tu tapes :

'q file'

Qu'est-ce que tu obtiens ?

[kevorage]

Donc il ouvre bien ton fichier ctl

Et si ensuite tu tapes :

'q file'

Qu'est-ce que tu obtiens ?

Il me "print" mon fichier ctl que j'ai posté plus au dessus.

[Tintin87]

Donc ton fichier a l'air bon.

Ressaye avec ton fichier test.gs en le simplifiant au maximum

'open /home/administrateur/wrfems/runs/MC/emsprd/grads/201108221000_arw_wrfout_d01.ctl''set background 0''set t 1''set gxout shaded''set csmooth on''d tmp2m -273.15''draw title Temperatures surface (2m)''printim tmp2m.png x800 y800'

puis en ligne de commande, en ayant quitté GrADS :

grads -blc 'run /home/chemin de ton fichiertest/test.gs'

Normalement ça devrait marcherPS : quelle version de GrADS tu as ?

EDIT : regarde bien que dans ton fichier test que le nom 201108221000_arw_wrfout__

d01.ctl s'écrive bien comme ça (présence de 2 tirets)

[kevorage]

Donc ton fichier a l'air bon.

Ressaye avec ton fichier test.gs en le simplifiant au maximum

'open /home/administrateur/wrfems/runs/MC/emsprd/grads/201108221000_arw_wrfout_d01.ctl''set background 0''set t 1''set gxout shaded''set csmooth on''d tmp2m -273.15''draw title Temperatures surface (2m)''printim tmp2m.png x800 y800'

puis en ligne de commande, en ayant quitté GrADS :

grads -blc 'run /home/chemin de ton fichiertest/test.gs'

Normalement ça devrait marcherPS : quelle version de GrADS tu as ?

Grid Analysis and Display System (GrADS) Version 2.0.a9

Copyright © 1988-2010 by Brian Doty and the

Institute for Global Environment and Society (IGES)

GrADS comes with ABSOLUTELY NO WARRANTY

See file COPYRIGHT for more information

Config: v2.0.a9 little-endian readline printim grib2 netcdf hdf4-sds hdf5 opendap-grids,stn geotiff shapefile

Issue 'q config' command for more detailed configuration information

GX Package Initialization: Size = 11 8.5

Running in Batch mode

Open error on tmp2m.png

PRINTIM error: open error

EDIT : regarde bien que dans ton fichier test que le nom 201108221000_arw_wrfout__

d01.ctl s'écrive bien comme ça (présence de 2 tirets)

J'ai vérifié, et modifié le contenu du script en conséquence, il n'y a qu'un tiret.

[Tintin87]

On avance, on avance, car apparemment il ouvre bien le fichier ctl.

Open error on tmp2m.png

PRINTIM error: open error

A mon avis tu lui a indiqué un dossier qui n'existe pas.

Creer un dossier /home/.../cartes/tmp2m/

dans ton fichier test, à la fin tape :

'printim /home/.../cartes/tmp2m/tmp2m.png x800 y800'

J'ai vérifié, et modifié le contenu du script en conséquence, il n'y a qu'un tiret.

C'est pour ça qu'il ne pouvais pas ouvrir le fichier ctl...

[kevorage]

On avance, on avance, car apparemment il ouvre bien le fichier ctl.

A mon avis tu lui a indiqué un dossier qui n'existe pas.

Creer un dossier /home/.../cartes/tmp2m/

dans ton fichier test, à la fin tape :

'printim /home/.../cartes/tmp2m/tmp2m.png x800 y800'

C'est pour ça qu'il ne pouvais pas ouvrir le fichier ctl...

Gagné, ça fonctionne. Bon j'ai une carte vierge avec les données, mais il génère bien.

Etape suivante, modifier le script Kumar pour que ça fonctionne ?

Ou d'abord m'occuper de la carte de fond ?

Merci

[ZHB]

Désolé pour les erreurs que j'ai glissé dans le script

Prochaine étape comprendre comment générer les cartes que tu veux en utilisant grads. Je te conseil de laisser le script Kumar de côté et de faire par toi même pour comprendre un minimum la génération de cartes avec Grads. De plus ça te permettra d'avoir des cartes personnalisées selon tes envies.

[kevorage]

Désolé pour les erreurs que j'ai glissé dans le script

Prochaine étape comprendre comment générer les cartes que tu veux en utilisant grads. Je te conseil de laisser le script Kumar de côté et de faire par toi même pour comprendre un minimum la génération de cartes avec Grads. De plus ça te permettra d'avoir des cartes personnalisées selon tes envies.

Pas de problèmes, faut dire que j'ai du en faire de bonnes aussi jusqu'ici

D'accord.

Je vais tâcher de commenter le script :

'open /home/administrateur/wrfems/runs/MC/emsprd/grads/201108221200_arw_wrfout_d01.ctl'

Ouverture du run 12Z de WRF ARW.

'set background 0'

Aucun fond de carte (d'ailleurs c'est un souci )

'set t 1'

Il s'agit de l'heure étudiée (ici certainement 13H)

'set gxout shaded'

Aucune idée...

'set csmooth on'

Bis...

'd tmp2m -273.15'

Encore une fois

'draw title Temperatures surface (2m)''printim /home/administrateur/wrfems/runs/MC/cartes/tmp2m/tmp2m.png x800 y800'

Afficher un titre particulier sur la carte, et la générer à tel endroit, avec comme dimension carrée 800pxlJ'ai bon ?

[Tintin87]

Bon j'ai une carte vierge avec les données

Quésaco ?

Je te conseil de laisser le script Kumar de côté

Oui tu peux très bien faire sans, c'est mon cas.

'set gxout shaded'

C'est pour ton type de graphique : à lire

'set csmooth on'

C'est pour "arrondir" les traits en gros.

'd tmp2m -273.15'

la commande d sert à "dessiner" tes donnéestmp2m : variable température en Kelvin

Je te conseille de suivre ça

[kevorage]

Quésaco ?

Oui tu peux très bien faire sans, c'est mon cas.

J'ai une carte de mon domaine, mais sans aucun fond (normal, le paramètre 'set background' est à 0).

Mais j'ai une carte

[ZHB]

'set background' => http://www.iges.org/grads/gadoc/commandsatt.html => http://www.iges.org/grads/gadoc/gradcomdsetbackground.html => set background color to color or color index => défini la couleur de fond.

Tu vois, la doc c'est utile

Petit bonus pour avoir un fond de carte, mais après faut vraiment que tu prennes la peine de cherche par toi même. Après si t'arrive vraiment pas y a pas de soucis pour poser des questions.

'set mproj nps''set mpdset hires'

Cependant si tu veux plus de détail sur ta carte il va falloir passer par les shapefiles, mais c'est une autre histoire.

[Tintin87]

Ah oui;

alors remplace

'set background 0'

par

'set mpdset hires'

et tout devrait rentrer dans l'ordre.EDIT : ZHB m'a devancé

[kevorage]

Ah oui;

alors remplace

'set background 0'

par

'set mpdset hires'

et tout devrait rentrer dans l'ordre.EDIT : ZHB m'a devancé

Plop, dernière chose et après, j'attaque le tuto

Ca n'a rien changé sur ma génération

[Moumoune]

Je crois que ma question s'est faite ensevelir par tous les derniers messages donc je me permet de la remettre:

En ces temps orageux je regarde pour tester des cartes de CAPE.

J'ai vu ce sujet: http://forums.infocl...-mucape-sbcape/

Et sur keraunos je vois qu'ils parlent souvent de MUCAPE. Lundi, leur carte de MUCAPE de WRF 8km était bien différente de la carte de CAPE de GFS 0.5 sur meteociel et plus représentative du risque orageux sur ma zone du nord IDF.

Dans mes fichiers ctl que j'arrive à sortir, en faisant un q file j'ai la liste des paramètres dispos et en cape j'ai ceci:

capesfc

cape180_0mb

cape90_0mb

cape255_0mb

Si j'ai compris le capesfc correspond au SBCAPE mais comment faire pour sortir une carte de MUCAPE?

merci

[Tintin87]

Bonjour, pour la MUCAPE : cape180_0mb

[Moumoune]

si c'est aussi simple que cela

Et une autre question, il y a des paramètres comme pwat30_0mb qui est l'eau précipitable sur les 30 premiers hPa de l'atmosphère si je comprends bien? (donc en gros 30x8 = 240m)

Ce n'est pas doublon avec pwatclm pour toute l'épaisseur de l'atmosphère? Que vaut il mieux utiliser en prévision des temps instables?

J'ai vue que les cartes de theta 'w sont biens pour la prévision de certains phénomènes comme les tempêtes où les gros fronts orageux. Ce paramètre n'est pas dispo dans la liste, faut il le créer à partir d'autres présents avec des formules par exemple?