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

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).