Discussion utilisateur:Ionela : Différence entre versions

Version du 10 mars 2022 à 15:41

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 "pres" est une variable 4D (x, y, z, t) contennant 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 pa 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, variable 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).