Changeset 1496 for LMDZ5/trunk/libf/phylmd/thermcell_main.F90

Timestamp:

Mar 11, 2011, 10:44:05 AM (13 years ago)

Author:

Laurent Fairhead

Message:

Modified closure for the new physics package, new values for the iflag_coupl parameter
FH

---

Fermeture modifiée pour la nouvelle physique, nouvelles valeurs définies pour
le paramètre iflag_coupl
FH

File:

• LMDZ5/trunk/libf/phylmd/thermcell_main.F90 (modified) (8 diffs)

LMDZ5/trunk/libf/phylmd/thermcell_main.F90

@@ -8 +8 @@
      &                  ,fm0,entr0,detr0,zqta,zqla,lmax  &
      &                  ,ratqscth,ratqsdiff,zqsatth  &
-     &                  ,r_aspect,l_mix,tau_thermals,iflag_thermals_ed &
      &                  ,Ale_bl,Alp_bl,lalim_conv,wght_th &
      &                  ,zmax0, f0,zw2,fraca,ztv &
@@ -47 +46 @@
 #include "FCTTRE.h"
 #include "iniprint.h"
+#include "thermcell.h"
 
 !   arguments:
@@ -54 +54 @@
       INTEGER itap
 
-      INTEGER ngrid,nlay,w2di
-      real tau_thermals
-      integer iflag_thermals_ed
-      real ptimestep,l_mix,r_aspect
+      INTEGER ngrid,nlay
+      real ptimestep
       REAL pt(ngrid,nlay),pdtadj(ngrid,nlay)
       REAL pu(ngrid,nlay),pduadj(ngrid,nlay)
@@ -162 +160 @@
       real Ale_bl(klon)
       real Alp_bl(klon)
-      real alp_int(klon)
+      real alp_int(klon),dp_int(klon),zdp
       real ale_int(klon)
       integer n_int(klon)
@@ -474 +472 @@
 
 ! Appel avec la version seche
-      CALL thermcell_closure(ngrid,nlay,r_aspect,ptimestep,rho,  &
+      CALL thermcell_closure(ngrid,nlay,r_aspect_thermals,ptimestep,rho,  &
      &   zlev,lalim,alim_star_clos,f_star,zmax_sec,wmax_sec,f,lev_out)
 
@@ -480 +478 @@
 ! Appel avec les zmax et wmax tenant compte de la condensation
 ! Semble moins bien marcher
-!     CALL thermcell_closure(ngrid,nlay,r_aspect,ptimestep,rho,  &
+!     CALL thermcell_closure(ngrid,nlay,r_aspect_thermals,ptimestep,rho,  &
 !    &   zlev,lalim,alim_star,f_star,zmax,wmax,f,lev_out)
 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
@@ -702 +700 @@
       enddo
 
-!     print*,'AAAAAAA ',Alp_bl,Ale_bl,lmix
-
-
-! TEST. IL FAUT REECRIRE LES ALE et ALP
-!     Ale_bl(:)=0.5*wmax(:)*wmax(:)
-!     Alp_bl(:)=0.1*wmax(:)*wmax(:)*wmax(:)
-
 !test:calcul de la ponderation des couches pour KE
 !initialisations
-!      print*,'ponderation'
 
       fm_tot(:)=0.
@@ -743 +733 @@
       endif
       enddo
+
+!------------------------------------------------------------------------
+! Modif CR/FH 20110310 : Alp integree sur la verticale.
+! Integrale verticale de ALP.
+! wth3 etant aux niveaux inter-couches, on utilise d play comme masse des
+! couches
+!------------------------------------------------------------------------
+
+      alp_int(:)=0.
+      dp_int(:)=0.
+      do l=2,nlay
+        do ig=1,ngrid
+           if(l.LE.lmax(ig)) THEN
+           zdp=pplay(ig,l-1)-pplay(ig,l)
+           alp_int(ig)=alp_int(ig)+0.5*rhobarz(ig,l)*wth3(ig,l)*zdp
+           dp_int(ig)=dp_int(ig)+zdp
+           endif
+        enddo
+      enddo
+
+      if (iflag_coupl==3 .or. iflag_coupl==4) then
+      do ig=1,ngrid
+!valeur integree de alp_bl * 0.5:
+        if (dp_int(ig)>0.) then
+        Alp_bl(ig)=alp_int(ig)/dp_int(ig)
+        endif
+      enddo!
+      endif
+
+
+! Facteur multiplicatif sur Alp_bl
+      Alp_bl(:)=alp_bl_k*Alp_bl(:)
+
+!------------------------------------------------------------------------