- Timestamp:
- Jul 22, 2024, 9:29:09 PM (2 months ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ6/branches/Amaury_dev/libf/dyn3d_common/exner_milieu_m.F90
r5082 r5099 6 6 7 7 SUBROUTINE exner_milieu ( ngrid, ps, p, pks, pk, pkf ) 8 ! 8 9 9 ! Auteurs : F. Forget , Y. Wanherdrick 10 10 ! P.Le Van , Fr. Hourdin . 11 11 ! .......... 12 ! 12 13 13 ! .... ngrid, ps,p sont des argum.d'entree au sous-prog ... 14 14 ! .... pks,pk,pkf sont des argum.de sortie au sous-prog ... 15 ! 15 16 16 ! ************************************************************************ 17 17 ! Calcule la fonction d'Exner pk = Cp * (p/preff) ** kappa , aux milieux des … … 21 21 ! .. N.B : Au sommet de l'atmosphere, p(llm+1) = 0. , et ps et pks sont 22 22 ! la pression et la fonction d'Exner au sol . 23 ! 23 24 24 ! WARNING : CECI est une version speciale de exner_hyb originale 25 25 ! Utilise dans la version martienne pour pouvoir … … 28 28 ! energie totale/ interne / potentielle (F.Forget 2001) 29 29 ! ( voir note de Fr.Hourdin ) , 30 ! 31 ! 30 31 32 32 USE comconst_mod, ONLY: jmp1, cpp, kappa, r 33 33 USE comvert_mod, ONLY: preff … … 100 100 ! .... Calcul de pk pour la couche l 101 101 ! -------------------------------------------- 102 ! 102 103 103 dum1 = cpp * (2*preff)**(-kappa) 104 104 DO l = 1, llm-1
Note: See TracChangeset
for help on using the changeset viewer.