HowTo : Calcul de la vraie pression à partir des niveaux hybrides sigma-pression de LMDZ : Différence entre versions
| Ligne 3 : | Ligne 3 : | ||
Les variables 4D (x,y,z,t) de LMDZ, comme ta, ua, va, etc sont sorties sur des niveaux hybrides sigma-pression. | Les variables 4D (x,y,z,t) de LMDZ, comme ta, ua, va, etc sont sorties sur des niveaux hybrides sigma-pression. | ||
| − | La | + | La variable 4D (x, y, z, t) "pres" contient les valeurs de pression de chaque maille pour chaque niveau vertical et pas de temps. |
| − | Lorsque "pres" n'est pas disponible, on peut utiliser la variable "geop" (qui correspond au géopotentiel de chaque maille pour chaque niveau vertical et pas de temps) et la pression à la surface pour calculer la pression sur la verticale. | + | Lorsque "pres" n'est pas disponible, on peut utiliser la variable "geop" (qui correspond au géopotentiel de chaque maille pour chaque niveau vertical et pas de temps) et la pression à la surface pour calculer la pression sur la verticale. |
Lorsque les variables "pres" ou "geop" ne sont pas disponibles, on peut reconstituer les vrais niveaux de pression du modele "true_presnivs", à l'aide des coefficients ap et bp et de la pression à la surface, appellée ici "preff". | Lorsque les variables "pres" ou "geop" ne sont pas disponibles, on peut reconstituer les vrais niveaux de pression du modele "true_presnivs", à l'aide des coefficients ap et bp et de la pression à la surface, appellée ici "preff". | ||
Version actuelle en date du 10 mars 2022 à 17:10
HowTo : Calcul de la vraie pression à partir des niveaux hybrides sigma-pression de LMDZ
Les variables 4D (x,y,z,t) de LMDZ, comme ta, ua, va, etc sont sorties sur des niveaux hybrides sigma-pression.
La variable 4D (x, y, z, t) "pres" contient les valeurs de pression de chaque maille pour chaque niveau vertical et pas de temps. Lorsque "pres" n'est pas disponible, on peut utiliser la variable "geop" (qui correspond au géopotentiel de chaque maille pour chaque niveau vertical et pas de temps) et la pression à la surface pour calculer la pression sur la verticale.
Lorsque les variables "pres" ou "geop" ne sont pas disponibles, on peut reconstituer les vrais niveaux de pression du modele "true_presnivs", à l'aide des coefficients ap et bp et de la pression à la surface, appellée ici "preff".
Dans la situation où les variables ap et bp sont données aux interfaces du modèle (ap et bp sont de dimension klevp1=80), on calcule "true_presnivs" par:
DO l = 1, klev
true_presnivs(l) = 0.5 *( ap(l)+bp(l)*preff + ap(l+1)+bp(l+1)*preff )
ENDDO
avec klev=79.
Dans le cas où ap et bp sont donnés aux milieux des couches (avec ap et bp de dimension klev), alors :
DO l = 1, klev
true_presnivs(l) = ap(l)+bp(l)*preff
ENDDO
La pression à la surface "preff" est également une variable 4D et est a priori incluse dans les fichiers 4D (comme ta*.nc, ua*.nc, va*.nc, etc).
