Discussion utilisateur:Ionela

HowTo : Calcul de la vraie pression à partir des niveaux hybrides sigma-pression de LMDZ

LMDZ utilise des niveaux hybrides sigma-pression sur la verticale.

Les variables 4D (x,y,z,t) de LMDZ, comme ta, ua, va, etc sont sorties sur ces niveaux hybrides sigma-pression.

La variable "pres" contient les valeurs de pression sur chaque niveau hybride.

Pour reconstituer les vrais niveaux de pression du modele "true_presnivs", il faut utiliser les coefficients ap et bp, ainsi que la pression a 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).