- Fichiers de commandes à essayer
- Variables de givrage par niveau de pression pour l'aviation
- Variables de turbulence par niveau de pression pour l'aviation
- Probabilité de givrage, pression, vents et température aux niveaux de vol standards
- Paramètres aéronautiques divers
- Température, humidité relative, hauteurs, et vents horizontaux et verticaux
- Glace et eau nuageuses, et condensat total
- Couvert nuageux, eau précipitable, humidité relative, et hydrométéores précipitants et nuages totaux dans la colonne
- Variables près de la surface, y compris visibilité, vents, humidité relative, paramètres de stabilité, courants ascendants et descendants, et autres variables en surface
- Variables de la couche limite planétaire (CLP), y compris qualité de l'air et indices de danger d'incendie, et température, humidité relative, pression et vents moyens dans la CLP et par couche de 30 hPa
- Scripts GrADS
Afin de vous aider à prévoir le temps en utilisant le modèle WRF-EMS, cette annexe contient des renseignements et des instructions qui vous serviront à extraire les paramètres que le modèle calcule pour l'aéronautique. Accédez au modèle WRF-EMS et aux instructions vous permettant de l'utiliser localement à l'adresse http://strc.comet.ucar.edu/software/newrems/
Pour utiliser les scripts figurant dans ce document, assurez-vous d'avoir bien installé et configuré le modèle EMS et que celui-ci génère des fichiers GRIB2. Pour de plus amples renseignements sur ces sujets, consultez la documentation à l'adresse http://strc.comet.ucar.edu/software/newrems/
Révisez les chapitres 2 et 9.
Pour représenter graphiquement les données des fichiers GRIB2, utilisez le logiciel GrADS, qui est installé parallèlement au modèle EMS. Pour que GrADS ait accès aux données et pour qu'il les traite, vous devez configurer les fichiers index et les fichiers de commandes, puis les lier avec les données.
Pour de plus amples renseignements sur la façon de créer des fichiers de commandes qui lisent les données GRIB1 et GRIB2, consultez la documentation à l'adresse http://grads.iges.org/grads/gadoc/grib.html
Fichiers de commandes à essayer
Cette section contient les lignes de code à inclure dans les fichiers de commandes, et ce, pour divers paramètres aéronautiques. Vous pouvez les adapter à vos propres besoins d'affichage.
Exemple de code que contient un des fichiers:
Ligne 1: dset [nom et emplacement du fichier de données]
Ligne 2: index [nom et emplacement du fichier index]
Ligne 3: undef 9.999E+20 [valeur assignée aux variables non définies]
Ligne 4: title [titre pertinent, p. ex. Paramètres aéronautiques aux niveaux de vol, début de la prévision <heure>UTC <date><mois><année>]
Ligne 5: dtype grib2 [indique le format des données de prévision, dans ce cas-ci, GRIB2]
Ligne 6: pdef 700 500 lccr -15.236000 19.971000 1 1 -2.000000 -2.000000 39.000000 6000.000000 6000.000000 [pdef est la projection utilisée pour les données du modèle; lccr est la projection de Lambert. Cette ligne de code indique au logiciel GrADS de quelle façon interpoler les données du modèle aux latitudes/longitudes qu'il utilise]
Ligne 7: xdef 674 linear 19.971631 0.0565324828044194 [la ligne xdef définit la longitude du début du tracé d'ouest en est, le nombre de points de l'ouest à l'est et la distance entre les points en x]
Ligne 8: ydef 474 linear -15.342380 0.0565324828044194 [la ligne ydef définit la latitude du début du tracé du sud au nord, le nombre de points du sud au nord et la distance entre les points en y]
Ligne 9: tdef 31 linear 12Z11apr2005 1hr [tdef définit l'heure de début (12 UTC, le 11 avril 2005, dans ce cas), le nombre de prévisions pour la période (31) et l'intervalle temporel entre les prévisions (ici, 1 heure)]
Ligne 10: zdef 15 levels 30 50 80 100 305 457 610 914 1524 1829 2134 2743 3658 4572 6000 [zdef indique le nombre de niveaux dans la verticale et les hauteurs à afficher (ici, les niveaux de vol sont en mètre, à partir du niveau inférieur)]
Ligne 11: vars 5 [vars est le nombre de variables à lire à partir du fichier de données. Ici, il y a cinq variables. Leur emplacement GRIB et leur description suivent]
Ligne 12: ICPRBm15,102 0,19,233 ** (6000 4572 3658 ... 80 50 30) Probabilité de givrage [non dim.]
[première des cinq variables du fichier de données, disponible sur les 15 niveaux et autres paramètres décrivant la variable. Ces paramètres varieront selon que la variable est une valeur instantanée, une valeur moyenne ou une valeur statistique comme une probabilité, etc.]
Ligne 13: PRESm15,102 0,3,0 ** (6000 4572 3658 ... 80 50 30) Pression [Pa] [comme ci-dessus, mais pour la pression à chaque niveau de vol. Notez que les unités apparaissent entre crochets, ici, en Pascal]
Ligne 14: TMPm15,102 0,0,0 ** (6000 4572 3658 ... 80 50 30) Température [K] [comme ci-dessus, mais pour la température à chaque niveau de vol, en Kelvin]
Ligne 15: UGRDm15,102 0,2,2 ** (6000 4572 3658 ... 80 50 30) Composante U du vent [m/s] [comme ci-dessus, mais pour le vent zonal, en mètre par seconde]
Ligne 16: VGRDm15,102 0,2,3 ** (6000 4572 3658 ... 80 50 30) Composante V du vent [m/s] [comme ci-dessus, mais pour le vent méridional, en mètre par seconde]
Ligne 17: ENDVARS [obligatoire à la fin de chaque fichier de commandes]
Notez que la ligne «dset» devrait pointer vers les données à lire par le fichier de commandes et que la ligne «index» devrait pointer vers le fichier qui relie le fichier de commande aux données. Ces «pointeurs» doivent être insérés dans le fichier de commandes où se trouve du [texte en italique entre crochets].
Les autres éléments du fichier de commandes dépendent des dimensions du domaine du modèle EMS (latitude, longitude, hauteur et heure, ainsi que de la projection servant à tracer les données. Notez que les latitudes, les longitudes et l'heure de début dans chacun de ces fichiers de commandes sont les mêmes, car elles proviennent toutes de la même simulation du modèle. Les lignes de code du premier fichier de commandes sont annotées pour permettre la compréhension de la structure du fichier.
Variables de givrage par niveau de pression pour l'aviation
Nom: [aaaammjjhhmm]_plevs_aviation_icing_d01.ctl
Les lignes commençant par un astérisque (*) sont des commentaires explicatifs et l'ordinateur ne les exécute pas.
dset [nom et emplacement du fichier de données]
index [nom et emplacement du fichier index]
undef 9.999E+20 [valeur assignée aux variables non définies]
title [titre pertinent, p. ex. Givrage aux niveaux de pression, début de la prévision <heure>UTC <date><mois><année>]
dtype grib2
*
* «pdef» est la projection utilisée pour les données du modèle; «lccr» est la projection de Lambert, mais d'autres options
* pourraient mieux convenir à votre domaine.
*
pdef 700 500 lccr -15.236000 19.971000 1 1 -2.000000 -2.000000 39.000000 6000.000000 6000.000000
*
* Les lignes «xdef» et «ydef» définissent l'emplacement du début du tracé de l'ouest à l'est et du sud au nord, et la résolution en x et en y.
*
xdef 674 linear 19.971631 0.0565324828044194
ydef 474 linear -15.342380 0.0565324828044194
*
* «tdef» définit l'heure de début (12 UTC, le 11 avril 2005, dans ce cas), le nombre de prévisions pour la période (31)
* et l'intervalle temporel entre les prévisions (ici, 1 heure).
*
tdef 31 linear 12Z11apr2005 1hr
*
* «zdef» indique le nombre de niveaux dans la verticale et les niveaux à afficher (ici, la pression est en Pascal).
* PROFIL hPa
*
zdef 45 levels 85000 84000 83000 82500 82000 81000 80000 79000 78000 77500 77000 76000 75000 74000 73000 72500 72000 71000 70000 69000 68000 67500 67000 66000 65000 64000 63000 62500 62000 61000 60000 59000 58000 57500 57000 56000 55000 54000 53000 52500 52000 51000 50000 47500 45000
*
* L'option «pascals» indique à GrADS que les niveaux sont donnés en Pascal.
*
options pascals
*
* Le fichier de commandes contient deux variables.
*
vars 2
*
* Nom des variables, nombre de niveaux dans la verticale (45), adresses dans le fichier GRIB (0,19,233 et 0,1,206).
*
ICPRBprs45,100 0,19,233 ** (850 840 830 825 820.. 520 510 500 475 450) Probabilité de givrage [non dim.]
TIPDprs45,100 0,1,206 ** (850 840 830 825 820.. 520 510 500 475 450) Diagnostic de givrage potentiel total [non dim.]
ENDVARS [obligatoire à la fin du fichier de commandes]
Variables de turbulence par niveau de pression pour l'aviation
Nom: [aaaammjjhhmm]_plevs_aviation_turb_d01.ctl
dset [nom et emplacement du fichier de données]
index [nom et emplacement du fichier index]
undef 9.999E+20
title [titre pertinent, p. ex. Givrage aux niveaux de pression, début de la prévision <heure>UTC <date><mois><année>]
* produit par g2ctl v0.0.7
* griddef=1:0:(700 x 500):grid_template=30:vents(grille): Lambert: (700 x 500) entrée OE:SN sortie OE:SN res 56 Lat1 -15.236000 Lon1 19.971000 LoV 39.000000 LatD -2.000000 Latin1 -2.000000 Latin2 -2.000000 LatSP -90.000000 LonSP 0.000000 Sud
dtype grib2
*
* (les derniers nombres représentent les résolutions en x et en y, en mètre).
*
pdef 700 500 lccr -15.236000 19.971000 1 1 -2.000000 -2.000000 39.000000 6000.000000 6000.000000
xdef 674 linear 19.971631 0.0565324828044194
ydef 474 linear -15.342380 0.0565324828044194
tdef 31 linear 12Z11apr2005 1hr
* PROFIL hPa
zdef 45 levels 85000 84000 83000 82500 82000 81000 80000 79000 78000 77500 77000 76000 75000 74000 73000 72500 72000 71000 70000 69000 68000 67500 67000 66000 65000 64000 63000 62500 62000 61000 60000 59000 58000 57500 57000 56000 55000 54000 53000 52500 52000 51000 50000 47500 45000
options pascals
vars 2
TPFIprs45,100 0,19,219 ** (850 840 830 825 820.. 520 510 500 475 450) Indice de turbulence potentielle prévue [-]
TURBINDXprs45,100 0,19,254 ** (850 840 830 825 820.. 520 510 500 475 450) Valeurs de turbulence utilisées pour déterminer la catégorie de l'indice de turbulence potentielle [-]
ENDVARS
Probabilité de givrage, pression, vents et température aux niveaux de vol standards
Nom: [aaaammjjhhmm]_aviation_fdlevs_d01.ctl
dset [nom et emplacement du fichier de données]
index [nom et emplacement du fichier index]
undef 9.999E+20
title [titre pertinent, p. ex. Givrage aux niveaux de pression, début de la prévision <heure>UTC <date><mois><année>]
* produit par g2ctl v0.0.7
dtype grib2
pdef 700 500 lccr -15.236000 19.971000 1 1 -2.000000 -2.000000 39.000000 6000.000000 6000.000000
xdef 674 linear 19.971631 0.0565324828044194
ydef 474 linear -15.342380 0.0565324828044194
tdef 31 linear 12Z11apr2005 1hr
* PROFIL m au-dessus du NMM
zdef 15 levels 30 50 80 100 305 457 610 914 1524 1829 2134 2743 3658 4572 6000
vars 5
ICPRBm15,102 0,19,233 ** (6000 4572 3658 2743 2134.. 305 100 80 50 30) Probabilité de givrage [non dim.]
PRESm15,102 0,3,0 ** (6000 4572 3658 2743 2134.. 305 100 80 50 30) Pression [Pa]
TMPm15,102 0,0,0 ** (6000 4572 3658 2743 2134.. 305 100 80 50 30) Température [K]
UGRDm15,102 0,2,2 ** (6000 4572 3658 2743 2134.. 305 100 80 50 30) Composante U du vent [m/s]
VGRDm15,102 0,2,3 ** (6000 4572 3658 2743 2134.. 305 100 80 50 30) Composante V du vent [m/s]
ENDVARS
Paramètres aéronautiques divers
Nom: [aaaammjjhhmm]_aviation_other_d01.ctl
dset [nom et emplacement du fichier de données]
index [nom et emplacement du fichier index]
undef 9.999E+20
title [titre pertinent, p. ex. Givrage aux niveaux de pression, début de la prévision <heure>UTC <date><mois><année>]
* produit par g2ctl v0.0.7
* griddef=1:0:(700 x 500):grid_template=30:vents(grille): Lambert: (700 x 500) entrée OE:SN sortie OE:SN res 56 Lat1 -15.236000 Lon1 19.971000 LoV 39.000000 LatD -2.000000 Latin1 -2.000000 Latin2 -2.000000 LatSP -90.000000 LonSP 0.000000 Sud
dtype grib2
pdef 700 500 lccr -15.236000 19.971000 1 1 -2.000000 -2.000000 39.000000 6000.000000 6000.000000
xdef 674 linear 19.971631 0.0565324828044194
ydef 474 linear -15.342380 0.0565324828044194
tdef 31 linear 12Z11apr2005 1hr
zdef 1 linear 1 1
vars 15
CEILsfc 0,1,0 0,6,13 ** Plafond près du sol [m]
FLGHTsfc 0,1,0 0,19,205 ** Catégorie de vol en surface [non dimensionnelle]
HGTtop0C 0,204,0 0,3,5 ** Hauteur géopotentielle du plus haut niveau de congélation troposphérique [gpm]
HGT0C 0,4,0 0,3,5 ** Hauteur géopotentielle de l'isotherme 0 degré C [gpm]
HGTmwl 0,6,0 0,3,5 ** Hauteur géopotentielle du vent maximal [gpm]
ICAHTmwl 0,6,0 0,3,3 ** Hauteur de référence (atmosphère standard de l'OACI) du vent maximal [m]
ICAHTtrop 0,7,0 0,3,3 ** Hauteur de référence (atmosphère standard de l'OACI) de la tropopause [m]
PREStop0C 0,204,0 0,3,0 ** Pression au niveau de congélation troposphérique le plus haut [Pa]
PRESmwl 0,6,0 0,3,0 ** Niveau de pression du vent maximal [Pa]
RHtop0C 0,204,0 0,1,1 ** Humidité relative au plus haut niveau troposphérique de congélation [%]
RH0C 0,4,0 0,1,1 ** Humidité relative à l'isotherme 0 degré C [%]
TMPmwl 0,6,0 0,0,0 ** Température au niveau du vent maximal [K]
UGRDmwl 0,6,0 0,2,2 ** Composante U du vent maximal [m/s]
VGRDmwl 0,6,0 0,2,3 ** Composante V du vent maximal [m/s]
VWSH610m 0,103,610 0,2,192 ** Cisaillement vertical de vitesse du vent, 610 m au-dessus du sol [1/s]
ENDVARS
Température, humidité relative, hauteurs, et vents horizontaux et verticaux
Nom: [aaaammjjhhmm]_plevs_depvars_d01.ctl
dset [nom et emplacement du fichier de données]
index [nom et emplacement du fichier index]
undef 9.999E+20
title [titre pertinent, p. ex. Givrage aux niveaux de pression, début de la prévision <heure>UTC <date><mois><année>]
* produit par g2ctl v0.0.7
* griddef=1:0:(700 x 500):grid_template=30:vents(grille): Lambert: (700 x 500) entrée OE:SN sortie OE:SN res 56 Lat1 -15.236000 Lon1 19.971000 LoV 39.000000 LatD -2.000000 Latin1 -2.000000 Latin2 -2.000000 LatSP -90.000000 LonSP 0.000000 Sud
dtype grib2
pdef 700 500 lccr -15.236000 19.971000 1 1 -2.000000 -2.000000 39.000000 6000.000000 6000.000000
xdef 674 linear 19.971631 0.0565324828044194
ydef 474 linear -15.342380 0.0565324828044194
tdef 31 linear 12Z11apr2005 1hr
* PROFIL hPa
zdef 35 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 25000 20000 15000 10000 5000
options pascals
vars 6
HGTprs35,100 0,3,5 ** (1013 1000 975 950 925.. 250 200 150 100 50) Hauteur géopotentielle [gpm]
RHprs35,100 0,1,1 ** (1013 1000 975 950 925.. 250 200 150 100 50) Humidité relative [%]
TMPprs35,100 0,0,0 ** (1013 1000 975 950 925.. 250 200 150 100 50) Température [K]
UGRDprs35,100 0,2,2 ** (1013 1000 975 950 925.. 250 200 150 100 50) Composante U du vent [m/s]
VGRDprs35,100 0,2,3 ** (1013 1000 975 950 925.. 250 200 150 100 50) Composante V du vent [m/s]
VVELprs35,100 0,2,8 ** (1013 1000 975 950 925.. 250 200 150 100 50) Vitesse du vent vertical (Pression) [Pa/s]
ENDVARS
Glace et eau nuageuses, et condensat total
Nom: [aaaammjjhhmm]_plevs_moistvars_d01.ctl
dset [nom et emplacement du fichier de données]
index [nom et emplacement du fichier index]
undef 9.999E+20
title [titre pertinent, p. ex. Givrage aux niveaux de pression, début de la prévision <heure>UTC <date><mois><année>]
* produit par g2ctl v0.0.7
* griddef=1:0:(700 x 500):grid_template=30:vents(grille): Lambert: (700 x 500) entrée OE:SN sortie OE:SN res 56 Lat1 -15.236000 Lon1 19.971000 LoV 39.000000 LatD -2.000000 Latin1 -2.000000 Latin2 -2.000000 LatSP -90.000000 LonSP 0.000000 Sud
dtype grib2
pdef 700 500 lccr -15.236000 19.971000 1 1 -2.000000 -2.000000 39.000000 6000.000000 6000.000000
xdef 674 linear 19.971631 0.0565324828044194
ydef 474 linear -15.342380 0.0565324828044194
tdef 31 linear 12Z11apr2005 1hr
* PROFIL hPa
zdef 35 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 25000 20000 15000 10000 5000
options pascals
vars 3
CICEprs35,100 0,6,0 ** (1013 1000 975 950 925.. 250 200 150 100 50) Glace nuageuse [kg/m^2]
CLWMRprs35,100 0,1,22 ** (1013 1000 975 950 925.. 250 200 150 100 50) Rapport de mélange des nuages [kg/kg]
TCONDprs35,100 0,6,195 ** (1013 1000 975 950 925.. 250 200 150 100 50) Condensat total [kg/kg]
ENDVARS
Couvert nuageux, eau précipitable, humidité relative, et hydrométéores précipitants et nuages totaux dans la colonne
Nom: [aaaammjjhhmm]_clouds_column_d01.ctl
dset [nom et emplacement du fichier de données]
index [nom et emplacement du fichier index]
undef 9.999E+20
title [titre pertinent, p. ex. Givrage aux niveaux de pression, début de la prévision <heure>UTC <date><mois><année>]
* produit par g2ctl v0.0.7
* griddef=1:0:(700 x 500):grid_template=30:vents(grille): Lambert: (700 x 500) entrée OE:SN sortie OE:SN res 56 Lat1 -15.236000 Lon1 19.971000 LoV 39.000000 LatD -2.000000 Latin1 -2.000000 Latin2 -2.000000 LatSP -90.000000 LonSP 0.000000 Sud
dtype grib2
pdef 700 500 lccr -15.236000 19.971000 1 1 -2.000000 -2.000000 39.000000 6000.000000 6000.000000
xdef 674 linear 19.971631 0.0565324828044194
ydef 474 linear -15.342380 0.0565324828044194
tdef 31 linear 12Z11apr2005 1hr
* PROFIL hPa
zdef 65 levels 101300 101000 100000 99000 98000 97500 97000 96000 95000 94000 93000 92500 92000 91000 90000 89000 88000 87500 87000 86000 85000 84000 83000 82500 82000 81000 80000 79000 78000 77500 77000 76000 75000 74000 73000 72500 72000 71000 70000 69000 68000 67500 67000 66000 65000 64000 63000 62500 62000 61000 60000 59000 58000 57500 57000 56000 55000 54000 53000 52500 52000 51000 50000 47500 45000
options pascals
vars 20
HCDChcll 0,234,0 0,6,5 ** Étendue de la couche nuageuse supérieure [%]
LCDClcll 0,214,0 0,6,3 ** Étendue de la couche nuageuse inférieure [%]
MCDCmcll 0,224,0 0,6,4 ** Étendue de la couche nuageuse moyenne [%]
MCONVclm 0,200,0 0,1,197 ** Divergence horizontale de l'humidité dans toute l'atmosphère (considérée comme une seule couche) [kg/kg/s]
PWATclm 0,200,0 0,1,3 ** Eau précipitable dans toute l'atmosphère (considérée comme une seule couche) [kg/m^2]
PWATMAXclm 0,200,0 0,1,240 ** Eau précipitable maximale dans toute l'atmosphère (considérée comme une seule couche) [kg/m^2]
RHclm 0,200,0 0,1,1 ** Humidité relative dans toute l'atmosphère (considérée comme une seule couche) [%]
TCDCprs65,100 0,6,1 ** (1013 1010 1000 990 980.. 520 510 500 475 450) Couvert nuageux total [%]
TCDCclm 0,200,0 0,6,1 ** Couvert nuageux total dans toute l'atmosphère (considérée comme une seule couche) [%]
TCLSWclm 0,200,0 0,1,209 ** Eau liquide surfondue totale sur toute la colonne [kg/m^2]
TCOLCclm 0,200,0 0,6,198 ** Condensat total sur toute la colonne [kg/m^2]
TCOLGclm 0,200,0 0,1,244 ** Neige roulée totale sur toute la colonne [kg/m^2]
TCOLGMXclm 0,200,0 0,1,243 ** Neige roulée totale maximale sur toute la colonne [kg/m^2]
TCOLHclm 0,200,0 0,1,242 ** Grêle totale sur toute la colonne [kg/m^2]
TCOLHMXclm 0,200,0 0,1,241 ** Grêle totale maximale sur toute la colonne [kg/m^2]
TCOLIclm 0,200,0 0,6,197 ** Glace nuageuse totale sur toute la colonne [kg/m^2]
TCOLMclm 0,200,0 0,1,210 ** Glace fondante totale sur toute la colonne [kg/m^2]
TCOLRclm 0,200,0 0,1,204 ** Pluie totale sur toute la colonne [kg/m^2]
TCOLSclm 0,200,0 0,1,205 ** Neige totale sur toute la colonne [kg/m^2]
TCOLWclm 0,200,0 0,6,196 ** Eau nuageuse totale sur toute la colonne [kg/m^2]
ENDVARS
Variables près de la surface, y compris visibilité, vents, humidité relative, paramètres de stabilité, courants ascendants et descendants, et autres variables en surface
Nom: [aaaammjjhhmm]_fields_near_surface_d01.ctl
dset [nom et emplacement du fichier de données]
index [nom et emplacement du fichier index]
undef 9.999E+20
title [titre pertinent, p. ex. Givrage aux niveaux de pression, début de la prévision <heure>UTC <date><mois><année>]
* produit par g2ctl v0.0.4i
dtype grib2
pdef 700 500 lccr -15.236000 19.971000 1 1 -2.000000 -2.000000 39.000000 6000.000000 6000.000000
xdef 673 linear 19.971631 0.056532482804954
ydef 491 linear -15.342380 0.0545454545454545
tdef 31 linear 12Z11apr2005 1hr
zdef 1 linear 1 1
vars 47
no4LFTX180_0mb 0,108,18000,0 0,7,193 ** Meilleur indice de soulèvement de 180 à 0 hPa au-dessus du sol (4 couches) [K]
no5WAVA1000_0m 0,103,1000,0 0,3,197 ** Anomalie de hauteur géopotentielle de 1000 à 0 hPa au-dessus du sol, 5 ondes [gpm]
CAPEsfc 0,1,0 0,7,6 ** Énergie potentielle de convection disponible en surface [J/kg]
CINsfc 0,1,0 0,7,7 ** Énergie d'inhibition de la convection en surface [J/kg]
DPT2m 0,103,2 0,0,6 ** Température du point de rosée 2 m au-dessus du sol [K]
ELONsfc 0,1,0 0,191,193 ** Longitude est en surface (0 à 360) [°]
GUST10m 0,103,10 0,2,22 ** Vitesse du vent (rafale) à 10 m au-dessus du sol [m/s]
HGTsfc 0,1,0 0,3,5 ** Hauteur géopotentielle en surface [gpm]
HGTeql 0,247,0 0,3,5 ** Hauteur géopotentielle du niveau d'équilibre [gpm]
IMGD2m 0,103,2 255,255,255 ** Données d'image à 2 m au-dessus du sol [-]
LANDsfc 0,1,0 2,0,0 ** Couverture du sol en surface (0=mer, 1=terre) [proportion]
LFTXl100_100 0,100,50000,100000 0,7,192 ** Indice de soulèvement en surface (500 à 1000 hPa) [K]
LRGHR6000_0m 0,103,6000,0 0,0,195 ** Taux de réchauffement du condensat à grande échelle, de 6000 à 0 m au-dessus du sol [K/s]
LTNGsfc 0,1,0 0,17,192 ** Foudre en surface [non dim.]
MAXDVVl100_100 0,100,100000,40000 0,2,221,2 ** Maximum horaire de la vitesse du vent vertical descendant entre 1000 et 400 hPa [m/s]
MAXRH2m 0,103,2 0,1,27 ** Humidité relative maximale à 2 m au-dessus du sol [%]
MAXUVVl100_100 0,100,100000,40000 0,2,220,2 ** Maximum horaire de la vitesse du vent vertical ascendant entre 1000 et 400 hPa [m/s]
MINRH2m 0,103,2 0,1,198 ** Humidité relative minimale à 2 m au-dessus du sol [%]
MIXR2m 0,103,2 0,1,2 ** Rapport de mélange de l'humidité à 2 m au-dessus du sol [kg/kg]
MSLETmsl 0,101,0 0,3,192 ** Pression au niveau moyen de la mer (réduction du modèle Eta) [Pa]
MXUPHL5000_2000m 0,103,5000,2000 0,7,199,2 ** Maximum horaire de l'hélicité du vent ascendant situé entre 5000 et 2000 m au-dessus du sol [m^2/s^2]
NLATsfc 0,1,0 0,191,192 ** Latitude en surface (-90 à 90) [°]
PLI30_0mb 0,108,3000,0 0,7,0 ** Indice de soulèvement de la particule de 30 à 0 hPa au-dessus du sol (jusqu'à 500 hPa) [K]
PRESsfc 0,1,0 0,3,0 ** Pression en surface [Pa]
PRES2m 0,103,2 0,3,0 ** Pression à 2 m au-dessus du sol [Pa]
PRMSLmsl 0,101,0 0,3,1 ** Pression au niveau moyen de la mer réduite au niveau moyen de la mer [Pa]
RH2m 0,103,2 0,1,1 ** Humidité relative à 2 m au-dessus du sol [%]
SPFH2m 0,103,2 0,1,0 ** Humidité spécifique à 2 m au-dessus du sol [kg/kg]
TMPsfc 0,1,0 0,0,0 ** Température en surface [K]
TMP2m 0,103,2 0,0,0 ** Température à 2 m au-dessus du sol [K]
UGRD10m 0,103,10 0,2,2 ** Composante U du vent à 10 m au-dessus du sol [m/s]
UPHL5000_2000m 0,103,5000,2000 0,7,197 ** Hélicité du vent ascendant entre 5000 et 2000 m au-dessus du sol [m^2/s^2]
VGRD10m 0,103,10 0,2,3 ** Composante V du vent à 10 m au-dessus du sol [m/s]
VISsfc 0,1,0 0,19,0 ** Visibilité en surface [m]
WTMPsfc 0,1,0 10,3,0 ** Température superficielle de l'eau [K]
t2mean2m 0,103,2 0,0,254,0 ** Température moyenne à 2 m au-dessus du sol [K]
rhmean2m 0,103,2 0,1,254,0 ** Humidité relative moyenne à 2 m au-dessus du sol [%]
vvmeanl100_100 0,100,100000,40000 0,2,246,0 ** Vitesse moyenne du vent vertical de 1000 à 400 hPa [m/s]
v10max10m 0,103,10 0,2,247 ** Composante V du vent maximal à 10 m au-dessus du sol [m/s]
u10max10m 0,103,10 0,2,248 ** Composante U du vent maximal à 10 m au-dessus du sol [m/s]
s10max10m 0,103,10 0,2,249 ** Magnitude du vent maximal à 10 m au-dessus du sol [m/s]
g10max10m 0,103,10 0,2,250 ** Rafale maximale à 10 m au-dessus du sol [m/s]
v10mean10m 0,103,10 0,2,252,0 ** Composante V du vent moyen à 10 m au-dessus du sol [m/s]
u10mean10m 0,103,10 0,2,253,0 ** Composante U du vent moyen à 10 m au-dessus du sol [m/s]
s10mean10m 0,103,10 0,2,254,0 ** Magnitude du vent moyen à 10 m au-dessus du sol [m/s]
ENDVARS
Variables de la couche limite planétaire (CLP), y compris qualité de l'air et indices de danger d'incendie, et température, humidité relative, pression et vents moyens dans la CLP et par couche de 30 hPa
Nom: [aaaammjjhhmm]_fields_pbl_d01.ctl
dset [nom et emplacement du fichier de données]
index [nom et emplacement du fichier index]
undef 9.999E+20
title [titre pertinent, p. ex. Givrage aux niveaux de pression, début de la prévision <heure>UTC <date><mois><année>]
* produit par g2ctl v0.0.7
* griddef=1:0:(700 x 500):grid_template=30:vents(grille): Lambert: (700 x 500) entrée OE:SN sortie OE:SN res 56 Lat1 -15.236000 Lon1 19.971000 LoV 39.000000 LatD -2.000000 Latin1 -2.000000 Latin2 -2.000000 LatSP -90.000000 LonSP 0.000000 Sud
dtype grib2
pdef 700 500 lccr -15.236000 19.971000 1 1 -2.000000 -2.000000 39.000000 6000.000000 6000.000000
xdef 674 linear 19.971631 0.0565324828044194
ydef 474 linear -15.342380 0.0565324828044194
tdef 31 linear 12Z11apr2005 1hr
zdef 1 linear 1 1
vars 39
FIREDLAsfc 0,1,0 2,4,2 ** Région potentielle d'incendie dû à la foudre sèche en surface [%]
HGTpbl 0,220,0 0,3,5 ** Hauteur géopotentielle de la couche limite planétaire [gpm]
HGTl5 0,5,0 0,3,5 ** Hauteur géopotentielle du niveau de condensation par soulèvement adiabatique à partir de la surface [gpm]
HPBLsfc 0,1,0 0,3,196 ** Hauteur de la couche limite planétaire à partir de la surface [m]
PBLREGsfc 0,1,0 0,19,12 ** Régime de la couche limite planétaire [-]
PRES30_0mb 0,108,3000,0 0,3,0 ** Pression entre la surface et 30 hPa au-dessus du sol [Pa]
PRES60_30mb 0,108,6000,3000 0,3,0 ** Pression entre 60 et 30 hPa au-dessus du sol [Pa]
PRES90_60mb 0,108,9000,6000 0,3,0 ** Pression entre 90 et 60 hPa au-dessus du sol [Pa]
PRES120_90mb 0,108,12000,9000 0,3,0 ** Pression entre 120 et 90 hPa au-dessus du sol [Pa]
PRES150_120mb 0,108,15000,12000 0,3,0 ** Pression entre 150 et 120 hPa au-dessus du sol [Pa]
PRES180_150mb 0,108,18000,15000 0,3,0 ** Pression entre 180 et 150 hPa au-dessus du sol [Pa]
PRESl5 0,5,0 0,3,0 ** Pression au niveau de condensation par soulèvement adiabatique à partir de la surface [Pa]
RH30_0mb 0,108,3000,0 0,1,1 ** Humidité relative entre la surface et 30 hPa au-dessus du sol [%]
RH60_30mb 0,108,6000,3000 0,1,1 ** Humidité relative entre 60 et 30 hPa au-dessus du sol [%]
RH90_60mb 0,108,9000,6000 0,1,1 ** Humidité relative entre 90 et 60 hPa au-dessus du sol [%]
RH120_90mb 0,108,12000,9000 0,1,1 ** Humidité relative entre 120 et 90 hPa au-dessus du sol [%]
RH150_120mb 0,108,15000,12000 0,1,1 ** Humidité relative entre 150 et 120 hPa au-dessus du sol [%]
RH180_150mb 0,108,18000,15000 0,1,1 ** Humidité relative entre 180 et 150 hPa au-dessus du sol [%]
TMP30_0mb 0,108,3000,0 0,0,0 ** Température entre la surface et 30 hPa au-dessus du sol [K]
TMP60_30mb 0,108,6000,3000 0,0,0 ** Température entre 60 et 30 hPa au-dessus du sol [K]
TMP90_60mb 0,108,9000,6000 0,0,0 ** Température entre 90 et 60 hPa au-dessus du sol [K]
TMP120_90mb 0,108,12000,9000 0,0,0 ** Température entre 120 et 90 hPa au-dessus du sol [K]
TMP150_120mb 0,108,15000,12000 0,0,0 ** Température entre 150 et 120 hPa au-dessus du sol [K]
TMP180_150mb 0,108,18000,15000 0,0,0 ** Température entre 180 et 150 hPa au-dessus du sol [K]
UGRD30_0mb 0,108,3000,0 0,2,2 ** Composante U du vent de la surface à 30 hPa au-dessus du sol [m/s]
UGRD60_30mb 0,108,6000,3000 0,2,2 ** Composante U du vent de 60 à 30 hPa au-dessus du sol [m/s]
UGRD90_60mb 0,108,9000,6000 0,2,2 ** Composante U du vent de 90 à 60 hPa au-dessus du sol [m/s]
UGRD120_90mb 0,108,12000,9000 0,2,2 ** Composante U du vent de 120 à 90 hPa au-dessus du sol [m/s]
UGRD150_120mb 0,108,15000,12000 0,2,2 ** Composante U du vent de 150 à 120 hPa au-dessus du sol [m/s]
UGRD180_150mb 0,108,18000,15000 0,2,2 ** Composante U du vent de 180 à 150 hPa au-dessus du sol [m/s]
UGRDpbl 0,220,0 0,2,2 ** Composante U du vent dans la couche limite planétaire [m/s]
VGRD30_0mb 0,108,3000,0 0,2,3 ** Composante V du vent de la surface à 30 hPa au-dessus du sol [m/s]
VGRD60_30mb 0,108,6000,3000 0,2,3 ** Composante V du vent de 60 à 30 hPa au-dessus du sol [m/s]
VGRD90_60mb 0,108,9000,6000 0,2,3 ** Composante V du vent de 90 à 60 hPa au-dessus du sol [m/s]
VGRD120_90mb 0,108,12000,9000 0,2,3 ** Composante V du vent de 120 à 90 hPa au-dessus du sol [m/s]
VGRD150_120mb 0,108,15000,12000 0,2,3 ** Composante V du vent de 150 à 120 hPa au-dessus du sol [m/s]
VGRD180_150mb 0,108,18000,15000 0,2,3 ** Composante V du vent de 180 à 150 hPa au-dessus du sol [m/s]
VGRDpbl 0,220,0 0,2,3 ** Composante V du vent dans la couche limite planétaire [m/s]
VRATEpbl 0,220,0 0,2,224 ** Taux de ventilation dans la couche limite planétaire [m^2/s]
ENDVARS
Scripts GrADS
Quand les résultats du modèle EMS sont prêts, vous pouvez utiliser les scripts GrADS pour les représenter graphiquement. Cette action peut s'effectuer sur l'ordinateur à partir de la ligne de commande de la fenêtre GrADS, après que les fichiers GRIB2 ont été écrits et reliés aux fichiers de commandes. Dans ce cas, copiez les scripts sur l'ordinateur, de préférence dans un répertoire réservé aux scripts GrADS. Ces scripts utilisent les fichiers de commandes figurant ci-dessus et les variables doivent porter le même nom que dans les fichiers de commandes.
Vous trouverez ci-dessous une liste de scripts utiles. Le nom du script apparaît entre parenthèses dans le titre. Vous verrez aussi un exemple de graphique créé par ce script. Les lignes de texte contenant un astérisque (*) expliquent à l'utilisateur la fonction qui les suit.
Parallèlement à l'examen des scripts, vous pouvez aussi consulter la documentation GrADS.
Série temporelle de paramètres à une station (pointts.gs)
*
* Ouvrir le fichier GrADS contenant le paramètre voulu. Ici, le plafond
* nuageux
*
'open [aaaammjjhhmm]_aviation_other_d01.ctl'
*
* Régler les coordonnées de longitude et de latitude. Les longitudes est et les latitudes nord sont positives.
*
'set lon 36.92 * Longitude de Nairobi, Kenya
'set lat -1.32 * Latitude de Nairobi, Kenya
*
* Entrer la durée de la série temporelle voulue. Ici, la durée est réglée pour afficher
* la prévision en entier, du début à la fin de la période de 30 heures.
*
'set t 1 31'
*
* Régler la longueur de l'axe vertical (y) pour l'affichage du plafond. La valeur maximale est de 5000 pieds, puisqu'au-dessus de celle-ci
* les règles de vol à vue (VFR) s'appliquent.
*
'set vrange 0 5000'
*
* Régler les seuils des catégories de vol, ici, ceux des plafonds. Les valeurs vont de
* 200 pieds (VLIFR) à 500 pieds (LIFR) à 1000 pieds (IFR) à 3000 pieds (MVFR), avec un maximum
* à 5000 pieds (VFR). Ces valeurs serviront à tracer des lignes horizontales sur le graphique.
*
'set vlifr=200'
'set lifr=500'
'set ifr=1000'
'set mvfr=3000'
*
* Définir les couleurs pour chaque ligne horizontale marquant les limites du plafond des catégories de vol
* et tracer chaque ligne. Régler «cmark» à 0 informe GrADS de ne pas tracer de marques pour chaque heure de validité de la prévision.
*
'set cthick 1' * Un trait gris très fin
'set ccolor 15'
'set cmark 0'
'd vlifr'
'set cthick 7'
'set cmark 0'
'd lifr' * Un trait gris épais
'set cthick 5'
'set ccolor 2'
'd ifr' * Un trait rouge épais
'set cthick 5'
'set ccolor 7'
'd mvfr' * Un trait jaune épais
*
* Définir la couleur et tracer le plafond prévu par l'EMS pour Nairobi, Kenya.
*
'set ccolor 4'
'set cmark 0'
'set cthick 3'
'd ceilsfc'
*
* Écrire le titre et le nom des axes. D'abord, obtenir les valeurs qui serviront au titre.
* Obtenir l'heure de lancement de la prévision en réglant l'heure à 1 (heure de début), en demandant cette heure,
* en «extrayant» la ligne contenant l'heure et en choisissant le «mot» représentant l'heure.
*
'set t 1'
'q time'
fwrd=subwrd(result,3)
fhr=substr(fwrd,1,2)
fdate=substr(fwrd,4,2)
fmo=substr(fwrd,6,3)
fyear=substr(fwrd,9,4)
*
* Maintenant, écrire le titre à partir des variables fhr, fdate, fmo et fyear, ci-dessus.
* Le saut de ligne est représenté par «\».
*
'draw title Jomo Kenyatta International Airport Cloud Ceiling (nearest grid point)\
EMS Forecast From 'fhr' UTC 'fdate' 'fmo' 'fyear
*
* Écrire le nom des axes x et y.
*
'draw xlab Forecast Valid Time'
'draw ylab Ceiling Height (feet)'
*
* Créer un fichier image en format .gif, avec un fond blanc.
*
'printim 'fhr'z'fdate''fmo''fyear'_JomoKenyatta_ceiling.gif gif white'
Cisaillement de vent à basse altitude (vwsh610mplot.gs)
*
* Ouvrir le fichier GrADS contenant divers paramètres aéronautiques.
*
'open [aaaammjjhhmm]_aviation_other_d01.ctl'
*
* Démarrer la boucle de traçage. i est un index servant à lire tous les fichiers de prévision. Commencer avec i = 1.
*
i = 1
*
* Entrer la valeur maximale de «i», donnée par le nombre d'heures de validité + 1. Ici, il existe 30 fichiers horaires commençant à
* T + 00 (heure de départ) et finissant à T + 30 (prévision de 30 heures, à intervalle horaire
* entre les prévisions).
*
while (i < 32)
*
* L'option «hires» de «mpdset» affiche les frontières politiques et les caractéristiques géographiques les plus détaillées.
*
'set mpdset hires'
'set map 1 1 7'
*
* Si vous désirez ne tracer qu'une sous-section du domaine du modèle EMS, vous devrez
* entrer une longitude et une latitude minimales et maximales. Voir ci-dessous.
*
* minlon et maxlon sont les longitudes minimale et maximale du domaine à afficher.
* minlat et maxlat sont les latitudes minimale et maximale du domaine à afficher.
*
'set lon <minlon> <maxlon>'
'set lat <minlon> <maxlon>''
'set t 'i
'q time'
*
* Utiliser le langage du script pour extraire du contenu servant à créer le titre et le nom du fichier. <cyc> est le cycle de prévision,
* <rundate> est le jour du mois, <runmo> est le mois et <runyear> est l'année de l'exécution du modèle.
* Les définir au début de l'exécution.
*
if (i = 1)
runwrd=subwrd(result,3)
cyc=substr(runwrd,1,2)
say 'cyc = 'cyc
rundate=substr(runwrd,4,2)
runmo=substr(runwrd,6,3)
runyear=substr(runwrd,9,4)
'q files'
runlin=sublin(result,1)
say 'runlin = 'runlin
runwrd=subwrd(runlin,4)
runfile=substr(runwrd,1,10)
endif
*
* Maintenant, obtenir les variables pour le titre du graphique. <fdate> est le jour du mois, <fmo> est le mois et
* <fyear> est l'année composant la date de validité.
*
'q time'
fwrd=subwrd(result,3)
fhr=substr(fwrd,1,2)
fdate=substr(fwrd,4,2)
fmo=substr(fwrd,6,3)
fyear=substr(fwrd,9,4)
'set gxout grfill'
'set grads off'
'set clevs 15 20 25 30 35 40 45 50'
'set ccols 0 9 4 5 3 10 7 8 2 6'
'd vwsh610m'
'cbarn'
*
* <ghr> est le délai de prévision en heure. Les fichiers pour ce script se succèdent à intervalles d'une heure,
* il s'agit donc de soustraire un pour que le premier graphique soit à T + 00, l'heure de début de la prévision.
*
ghr = i - 1
*
'draw title Lowest 2000 Ft. Wind Shear, EMS 'cyc' UTC 'rundate' 'runmo' 'runyear' Forecast\Valid Time 'fhr' UTC 'fdate' 'fmo' 'fyear
*
* Nom du fichier à tracer. Notez le zéro non significatif pour les délais de prévision de moins de 10. Vous pouvez
* modifier le nom pour qu'il soit pertinent. Ici, le nom de la variable à afficher est utilisé.
*
if (ghr < 10)
'printim vwsh601m_'runfile'f0'ghr'.gif gif white'
else
'printim vwsh610m_'runfile'f'ghr'.gif gif white'
endif
*
* Suppression du contenu et préparation du graphique suivant.
*
'clear'
i = i + 1
endwhile
Visibilité (visplot.gs)
'open [aaaammjjhhmm]_aviation_other_d01.ctl'
i = 1
while (i < 32)
say 'i = 'i
'set mpdset hires'
'set map 1 1 7'
'set t 'i
'q time'
if (i = 1)
runwrd=subwrd(result,3)
cyc=substr(runwrd,1,2)
say 'cyc = 'cyc
rundate=substr(runwrd,4,2)
runmo=substr(runwrd,6,3)
runyear=substr(runwrd,9,4)
'q files'
runlin=sublin(result,1)
say 'runlin = 'runlin
runwrd=subwrd(runlin,4)
runfile=substr(runwrd,1,10)
endif
'q time'
fwrd=subwrd(result,3)
fhr=substr(fwrd,1,2)
fdate=substr(fwrd,4,2)
fmo=substr(fwrd,6,3)
fyear=substr(fwrd,9,4)
'set gxout grfill'
'set grads off'
'set clevs 1 3 5'
'set ccols 15 2 7 0'
'd vissfc/1609.34'
'cbarn'
ghr = i - 1
'draw title Visibility (statute miles), EMS 'cyc' UTC 'rundate' 'runmo' 'runyear' Forecast\Valid Time 'fhr' UTC 'fdate' 'fmo' 'fyear
if (ghr < 10)
'printim vissfc_'runfile'f0'ghr'zoom.gif gif white'
else
'printim vissfc_'runfile'f'ghr'zoom.gif gif white'
endif
'clear'
i = i + 1
endwhile
Plafond nuageux (cldbase.gs)
'open [aaaammjjhhmm]_aviation_other_d01.ctl'
i = 1
while (i < 32)
'set mpdset hires'
'set map 1 1 7'
'set t 'i
'q time'
if (i = 1)
runwrd=subwrd(result,3)
cyc=substr(runwrd,1,2)
say 'cyc = 'cyc
rundate=substr(runwrd,4,2)
runmo=substr(runwrd,6,3)
runyear=substr(runwrd,9,4)
'q files'
runlin=sublin(result,1)
say 'runlin = 'runlin
runwrd=subwrd(runlin,4)
runfile=substr(runwrd,1,10)
endif
'q time'
fwrd=subwrd(result,3)
fhr=substr(fwrd,1,2)
fdate=substr(fwrd,4,2)
fmo=substr(fwrd,6,3)
fyear=substr(fwrd,9,4)
'set gxout grfill'
'set grads off'
'set clevs 500 1000 3000'
'set ccols 15 2 7 0'
'd ceilsfc/3.28084'
'cbarn'
ghr = i - 1
'draw title Cloud Ceiling (feet), WRF-EMS 'cyc' UTC 'rundate' 'runmo' 'runyear' Forecast\Valid Time 'fhr' UTC 'fdate' 'fmo' 'fyear
if (ghr < 10)
'printim ceilsfc_'runfile'f0'ghr'.gif gif white'
else
'printim ceilsfc_'runfile'f'ghr'.gif gif white'
endif
'clear'
i = i + 1
endwhile
Catégorie de vol (fltsfc.gs)
'open [aaaammjjhhmm]_aviation_other_d01.ctl'
i = 1
while (i < 32)
say 'i = 'i
'set mpdset hires'
'set map 1 1 7'
'set t 'i
'q time'
if (i = 1)
runwrd=subwrd(result,3)
cyc=substr(runwrd,1,2)
say 'cyc = 'cyc
rundate=substr(runwrd,4,2)
runmo=substr(runwrd,6,3)
runyear=substr(runwrd,9,4)
'q files'
runlin=sublin(result,1)
say 'runlin = 'runlin
runwrd=subwrd(runlin,4)
runfile=substr(runwrd,1,10)
endif
'q time'
fwrd=subwrd(result,3)
fhr=substr(fwrd,1,2)
fdate=substr(fwrd,4,2)
fmo=substr(fwrd,6,3)
fyear=substr(fwrd,9,4)
'set gxout grfill'
'set grads off'
'set clevs 1 2 3'
'set ccols 6 8 7 0'
'd flghtsfc'
'cbarn'
ghr = i - 1
'draw title Flight Condition, WRF-EMS 'cyc' UTC 'rundate' 'runmo' 'runyear' Forecast\Valid Time 'fhr' UTC 'fdate' 'fmo' 'fyear
if (ghr < 10)
'printim flghtsfc_'runfile'f0'ghr'.gif gif white'
else
'printim flghtsfc_'runfile'f'ghr'.gif gif white'
endif
'clear'
i = i + 1
endwhile
Rapport de mélange de l'eau nuageuse et température (clwmrplot.gs)
'open [aaaammjjhhmm]_plevs_moistvars_d01.ctl'
'open [aaaammjjhhmm]_plevs_depvars_d01.ctl'
i = 1
say 'what level do you want?' * Quel niveau voulez-vous?
pull level
'set lev 'level
while (i < 32)
say 'i = 'i
'set mpdset hires'
'set map 1 1 7'
'set t 'i
'q time'
if (i = 1)
runwrd=subwrd(result,3)
cyc=substr(runwrd,1,2)
say 'cyc = 'cyc
rundate=substr(runwrd,4,2)
runmo=substr(runwrd,6,3)
runyear=substr(runwrd,9,4)
'q files'
runlin=sublin(result,1)
say 'runlin = 'runlin
runwrd=subwrd(runlin,4)
runfile=substr(runwrd,1,10)
endif
'q time'
fwrd=subwrd(result,3)
fhr=substr(fwrd,1,2)
fdate=substr(fwrd,4,2)
fmo=substr(fwrd,6,3)
fyear=substr(fwrd,9,4)
'set gxout grfill'
'set grads off'
'set clevs 0.01 0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4'
'set ccols 0 9 14 4 10 3 7 12 8 2 6'
'd clwmrprs*1000.'
'cbarn'
'set gxout contour'
'set ccolor 1'
'set cthick 7'
'set cint 3'
'd tmpprs.2-273.15'
ghr = i - 1
'draw title 'level'-hPa Cld Wat (g/Kg) and Tmp, WRF-EMS 'cyc' UTC 'rundate' 'runmo' 'runyear' Fcst.\Valid Time 'fhr' UTC 'fdate' 'fmo' 'fyear
if (ghr < 10)
'printim WRF_'level'Tclwmr_'runfile'f0'ghr'.gif gif white'
else
'printim WRF_'level'Tclwmr_'runfile'f'ghr'.gif gif white'
endif
'clear'
i = i + 1
endwhile
Probabilité de givrage aux niveaux de vol standards (icfltevplot.gs)
'open [aaaammjjhhmm]_aviation_fdlevs_d01.ctl'
'q file'
say 'what level do you want?' * Quel niveau voulez-vous? (6000, 4753, etc.)'
pull level
i = 1
while (i < 32)
#while (i < 13)
say 'i = 'i
'set mpdset hires'
'set map 1 1 7'
'set t 'i
'set lev 'level
'q time'
if (i = 1)
runwrd=subwrd(result,3)
cyc=substr(runwrd,1,2)
say 'cyc = 'cyc
rundate=substr(runwrd,4,2)
runmo=substr(runwrd,6,3)
runyear=substr(runwrd,9,4)
'q files'
runlin=sublin(result,1)
say 'runlin = 'runlin
runwrd=subwrd(runlin,4)
runfile=substr(runwrd,1,10)
endif
'q time'
fwrd=subwrd(result,3)
fhr=substr(fwrd,1,2)
fdate=substr(fwrd,4,2)
fmo=substr(fwrd,6,3)
fyear=substr(fwrd,9,4)
'set gxout grfill'
'set grads off'
'set clevs 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9'
'set ccols 0 9 4 5 3 10 7 12 8 2'
'd icprbm'
'cbarn'
ghr = i - 1
'draw title Icing Prob. (fraction) at 'level' meters, WRF-EMS 'cyc' UTC 'rundate' 'runmo' 'runyear'\Forecast Valid Time 'fhr' UTC 'fdate' 'fmo' 'fyear
if (ghr < 10)
'printim icbm'level'_'runfile'f0'ghr'.gif gif white'
else
'printim icbm'level'_'runfile'f'ghr'.gif gif white'
endif
say 'continue?'
pull dummy
'clear'
i = i + 1
endwhile
Meilleur indice de soulèvement (LI, LIsfcplot.gs)
'open [aaaammjjhhmm]_fields_near_surface_d01.ctl'
i = 1
while (i < 32)
say 'i = 'i
'set mpdset hires'
'set map 1 1 7'
'set t 'i
'q time'
if (i = 1)
runwrd=subwrd(result,3)
cyc=substr(runwrd,1,2)
say 'cyc = 'cyc
rundate=substr(runwrd,4,2)
runmo=substr(runwrd,6,3)
runyear=substr(runwrd,9,4)
'q files'
runlin=sublin(result,1)
say 'runlin = 'runlin
runwrd=subwrd(runlin,4)
runfile=substr(runwrd,1,10)
endif
'q time'
fwrd=subwrd(result,3)
fhr=substr(fwrd,1,2)
fdate=substr(fwrd,4,2)
fmo=substr(fwrd,6,3)
fyear=substr(fwrd,9,4)
'set gxout grfill'
'set grads off'
'set clevs -8 -6 -4 -2 0 2 4 6'
'set ccols 9 4 5 3 10 0 7 8 2 6'
'd no4LFTX180_0mb'
'cbarn'
ghr = i - 1
'draw title Best Lifted Index (0-180hPa above sfc), WRF-EMS 'cyc' UTC 'rundate' 'runmo' 'runyear'\Forecast Valid Time 'fhr' UTC 'fdate' 'fmo' 'fyear
if (ghr < 10)
'printim bestLI_'runfile'f0'ghr'.gif gif white'
else
'printim bestLI_'runfile'f'ghr'.gif gif white'
endif
'clear'
i = i + 1
endwhile
Énergie potentielle de convection disponible (EPCD) et énergie d'inhibition de la convection (EIC) (capesfcplot.gs)
'open [aaaammjjhhmm]_fields_near_surface_d01.ctl'
i = 1
while (i < 32)
say 'i = 'i
'set mpdset hires'
'set map 1 1 7'
'set t 'i
'q time'
if (i = 1)
runwrd=subwrd(result,3)
cyc=substr(runwrd,1,2)
say 'cyc = 'cyc
rundate=substr(runwrd,4,2)
runmo=substr(runwrd,6,3)
runyear=substr(runwrd,9,4)
'q files'
runlin=sublin(result,1)
say 'runlin = 'runlin
runwrd=subwrd(runlin,4)
runfile=substr(runwrd,1,10)
endif
'q time'
fwrd=subwrd(result,3)
fhr=substr(fwrd,1,2)
fdate=substr(fwrd,4,2)
fmo=substr(fwrd,6,3)
fyear=substr(fwrd,9,4)
'set gxout shaded'
'set grads off'
'set clevs 0 500 1000 1500 2000 2500 3000 3500 4000 5000 6000'
'set ccols 0 9 4 5 3 10 7 12 8 2 6 15'
'd capesfc'
'cbarn'
'set gxout contour'
'set clevs -500 -400 -300 -200 -100 -50 0 50 100 200 300 400 500'
'set ccolor 1'
'set cthick 7'
'set clopts -1 7 0.1'
'd cinsfc'
ghr = i - 1
'draw title Surface CAPE (shaded) and CIN (contour),\WRF-EMS 'cyc' UTC 'rundate' 'runmo' 'runyear'\Forecast Valid Time 'fhr' UTC 'fdate' 'fmo' 'fyear
if (ghr < 10)
'printim capecin_'runfile'f0'ghr'.gif gif white'
else
'printim capecin_'runfile'f'ghr'.gif gif white'
endif
'clear'
i = i + 1
endwhile
Température de luminance infrarouge (IR) (brtmptoaplot.gs)
'open 200610080000_fields_precip_d01.ctl'
i = 1
while (i < 32)
say 'i = 'i
'set mpdset hires'
'set map 1 1 7'
'set t 'i
'q time'
if (i = 1)
runwrd=subwrd(result,3)
cyc=substr(runwrd,1,2)
say 'cyc = 'cyc
rundate=substr(runwrd,4,2)
runmo=substr(runwrd,6,3)
runyear=substr(runwrd,9,4)
'q files'
runlin=sublin(result,1)
say 'runlin = 'runlin
runwrd=subwrd(runlin,4)
runfile=substr(runwrd,1,10)
endif
'q time'
fwrd=subwrd(result,3)
fhr=substr(fwrd,1,2)
fdate=substr(fwrd,4,2)
fmo=substr(fwrd,6,3)
fyear=substr(fwrd,9,4)
'set gxout grfill'
'set grads off'
*
* Réglage des couleurs de l'illustration pour correspondance avec les couleurs EUMETSAT du canal infrarouge.
*
'set clevs -80 -79 -78 -77 -76 -75 -74 -73 -72 -71 -70 -69 -68 -67 -66 -65 -64 -63 -62 -61 -60 -59 -58 -57 -56 -55 -54 -53 -52 -51 -50 -49 -48 -47 -46 -45 -44 -43 -42 -41 -40 -39 -38 -37 -34 -31 -27 -23 -19 -15 -11 -7 -3 1 5 9 13 17 21 25 29 33 37 40'
'set ccols 153 151 149 148 147 146 145 144 143 142 141 140 139 138 136 134 133 132 129 128 127 126 125 124 123 122 120 119 118 117 116 115 114 113 112 111 110 109 108 107 106 105 104 103 102 101 100 99 98 97 96 95 94 93 92 91 90 89 88 87 86 85 84 83 82 81'
'd brtmptoa'
'cbarn'
ghr = i - 1
'draw title TOA Brightness Temperature (K), WRF-EMS 'cyc' UTC 'rundate' 'runmo' 'runyear' Forecast\Valid Time 'fhr' UTC 'fdate' 'fmo' 'fyear
if (ghr < 10)
'printim brtmptoa_'runfile'f0'ghr'.gif gif white'
else
'printim brtmptoa_'runfile'f'ghr'.gif gif white'
endif
'clear'
i = i + 1
endwhile
Réflectivité radar composite (plotcompradarrefl.gs)
'open 200908311200_fields_precip_d01.ctl'
'open 200908311200_fields_near_surface_d01.ctl'
'set lon -7 7'
'set lat 7 17'
'set mpdset hires'
'set map 1 1 5'
'q files'
line=sublin(result,1)
word=subwrd(line,4)
runmo=substr(word,5,2)
runyr=substr(word,1,4)
rundt=substr(word,7,2)
runhr=substr(word,9,2)
say 'run month = 'runmo
say 'run year = 'runyr
say 'run date = 'rundt
say 'run hour = 'runhr
'set gxout grfill'
'run radarcols.gs' * ce script règle les couleurs
t = 1
'set t 't
while (t < 37)
'set t 't
'set clevs 5 10 15 20 25 30 35 40 45 50 55 60 65 70 75'
'set ccols 0 21 22 23 24 25 26 27 28 29 30 31 32 33 34 0'
'set grads off'
'd refcclm'
'cbarn'
fhr=(t-1)
'set gxout barb'
'set digsiz 0.07'
'set ccolor 1'
'set cthick 7'
'd skip(ugrd10m.2,9)*1.94;vgrd10m.2*1.94'
'draw title 'fhr'-hr Forecast Sim Composite Radar Refl/10m Wind\Run from 'runhr' UTC 'runmo'/'rundt'/'runyr
'printim simcomprefl_'word'f'fhr'zoom.gif gif white'
t=t+1
'clear'
endwhile
Échelle de couleur des réflectivités (radarcols.gs)
*
* Réglage des couleurs de réflectivité radar pour correspondance avec les couleurs des images radar américaines.
*
'set rgb 21 0 236 236'
'set rgb 22 1 160 246'
'set rgb 23 0 0 246'
'set rgb 24 0 255 0'
'set rgb 25 0 200 0'
'set rgb 26 0 144 0'
'set rgb 27 255 255 0'
'set rgb 28 231 192 0'
'set rgb 29 255 144 0'
'set rgb 30 255 0 0'
'set rgb 31 214 0 0'
'set rgb 32 192 0 0'
'set rgb 33 255 0 255'
'set rgb 34 153 85 201'
Exemple d'illustration issue du script ci-dessus, pour le domaine complet couvrant l'Afrique occidentale, centré près du Burkina Faso.
Sommet des nuages par radar (retopclmplot.gs)
'open 200908311200_fields_precip_d01.ctl'
'set lon -7 7'
'set lat 7 17'
i = 1
while (i < 32)
#while (i < 18)
say 'i = 'i
'set mpdset hires'
'set map 1 1 7'
'set t 'i
'q time'
if (i = 1)
runwrd=subwrd(result,3)
cyc=substr(runwrd,1,2)
say 'cyc = 'cyc
rundate=substr(runwrd,4,2)
runmo=substr(runwrd,6,3)
runyear=substr(runwrd,9,4)
'q files'
runlin=sublin(result,1)
say 'runlin = 'runlin
runwrd=subwrd(runlin,4)
runfile=substr(runwrd,1,10)
endif
'q time'
fwrd=subwrd(result,3)
fhr=substr(fwrd,1,2)
fdate=substr(fwrd,4,2)
fmo=substr(fwrd,6,3)
fyear=substr(fwrd,9,4)
'set gxout grfill'
'set grads off'
'set clevs 0 25000 30000 35000 40000 45000 50000 55000 60000'
'set ccols 0 9 4 5 3 7 8 2 6'
'd retopclm*3.28'
'cbarn'
ghr = i - 1
'draw title Simulated Echo Cloud Tops, WRF-EMS 'cyc' UTC 'rundate' 'runmo' 'runyear'\Forecast Valid Time 'fhr' UTC 'fdate' 'fmo' 'fyear
if (ghr < 10)
'printim retopclm_'runfile'f0'ghr'zoom.gif gif white'
else
'printim retopclm_'runfile'f'ghr'zoom.gif gif white'
endif
say 'hit enter to continue' * Touche de retour pour continuer
pull dummy
'clear'
i = i + 1
endwhile
Exemple d'illustration issue du script ci-dessus, pour un sous-domaine couvrant l'Afrique occidentale, centré près de Ouagadougou, Burkina Faso.
Indice de turbulence (niveaux de pression, turbprsplot.gs)
'open 200908311200_plevs_aviation_turb_d01.ctl'
'q file'
say 'what pressure level do you want? * Quel niveau de pression voulez-vous? (450, 475, 500, 510, 520, 530, etc.)'
pull level
i = 1
while (i < 32)
#while (i < 13)
say 'i = 'i
'set mpdset hires'
'set map 1 1 7'
'set t 'i
'set lev 'level
'q time'
if (i = 1)
runwrd=subwrd(result,3)
cyc=substr(runwrd,1,2)
say 'cyc = 'cyc
rundate=substr(runwrd,4,2)
runmo=substr(runwrd,6,3)
runyear=substr(runwrd,9,4)
'q files'
runlin=sublin(result,1)
say 'runlin = 'runlin
runwrd=subwrd(runlin,4)
runfile=substr(runwrd,1,10)
endif
'q time'
fwrd=subwrd(result,3)
fhr=substr(fwrd,1,2)
fdate=substr(fwrd,4,2)
fmo=substr(fwrd,6,3)
fyear=substr(fwrd,9,4)
'set gxout grfill'
'set grads off'
'set clevs 1 2 3'
'set ccols 0 7 8 2'
'd tpfiprs'
'cbarn'
ghr = i - 1
'draw title Turb. Index at 'level'-hPa pressure, WRF-EMS 'cyc' UTC 'rundate' 'runmo' 'runyear'\Forecast Valid Time 'fhr' UTC 'fdate' 'fmo' 'fyear
if (ghr < 10)
'printim tpfi'level'_'runfile'f0'ghr'.gif gif white'
else
'printim tpfi'level'_'runfile'f'ghr'.gif gif white'
endif
say 'continue?'
pull dummy
'clear'
i = i + 1
endwhile
Vitesse du vent vertical (niveaux de pression, vvelplot.gs)
'open 200610080000_plevs_depvars_d01.ctl'
i = 1
say 'what level do you want?' * Quel niveau voulez-vous?
pull level
'set lev 'level
while (i < 32)
say 'i = 'i
'set mpdset hires'
'set map 1 1 7'
'set t 'i
'q time'
if (i = 1)
runwrd=subwrd(result,3)
cyc=substr(runwrd,1,2)
say 'cyc = 'cyc
rundate=substr(runwrd,4,2)
runmo=substr(runwrd,6,3)
runyear=substr(runwrd,9,4)
'q files'
runlin=sublin(result,1)
say 'runlin = 'runlin
runwrd=subwrd(runlin,4)
runfile=substr(runwrd,1,10)
endif
'q time'
fwrd=subwrd(result,3)
fhr=substr(fwrd,1,2)
fdate=substr(fwrd,4,2)
fmo=substr(fwrd,6,3)
fyear=substr(fwrd,9,4)
'set gxout grfill'
'set grads off'
'set clevs -2.5 -2 -1.5 -1 -0.5 0.5 1 1.5 2 2.5'
'set ccols 6 2 8 12 7 0 3 10 4 14 9'
'd vvelprs'
'cbarn'
ghr = i - 1
'draw title 'level'-hPa Vert Vel (Pa/sec), WRF-EMS 'cyc' UTC 'rundate' 'runmo' 'runyear' Forecast\Valid Time 'fhr' UTC 'fdate' 'fmo' 'fyear
if (ghr < 10)
'printim WRF_'level'vvel_'runfile'f0'ghr'.gif gif white'
else
'printim WRF_'level'vvel_'runfile'f'ghr'.gif gif white'
endif
'clear'
i = i + 1
endwhile
Exemple d'illustration issue du script ci-dessus, pour le domaine complet couvrant l'Afrique australe, centré près du sud du Mozambique.
