Changeset 1338 for LMDZ4/branches/LMDZ4V5.0-dev/libf/phylmd/thermcell_plume.F90

Timestamp:

Apr 6, 2010, 2:49:00 PM (14 years ago)

Author:

idelkadi

Message:

Nouvelle version des thermiques

File:

• LMDZ4/branches/LMDZ4V5.0-dev/libf/phylmd/thermcell_plume.F90 (modified) (18 diffs)

LMDZ4/branches/LMDZ4V5.0-dev/libf/phylmd/thermcell_plume.F90

@@ -64 +64 @@
       REAL zqsatth(ngrid,klev)
       REAL zta_est(ngrid,klev)
+      REAL ztemp(ngrid),zqsat(ngrid)
       REAL zdw2
       REAL zw2modif
@@ -77 +78 @@
       real zdz,zfact,zbuoy,zalpha
       real zcor,zdelta,zcvm5,qlbef
-      real Tbef,qsatbef
-      real dqsat_dT,DT,num,denom
       REAL REPS,RLvCp,DDT0
       PARAMETER (DDT0=.01)
@@ -116 +115 @@
 !      write(lunout,*)'THERM 31H '
 
+      print*,'THERMCELL_PLUME OPTIMISE V0 '
 ! Initialisations des variables reeles
 if (1==0) then
@@ -239 +239 @@
 ! sans tenir compte du detrainement et de l'entrainement dans cette
 ! couche
+! C'est a dire qu'on suppose 
+! ztla(l)=ztla(l-1) et zqta(l)=zqta(l-1)
 ! Ici encore, on doit pouvoir ajouter entr_star (qui peut etre calculer
 ! avant) a l'alimentation pour avoir un calcul plus propre
 !---------------------------------------------------------------------------
 
-     call thermcell_condens(ngrid,active, &
-&          zpspsk(:,l),pplev(:,l),ztla(:,l-1),zqta(:,l-1),zqla_est(:,l))
+   ztemp(:)=zpspsk(:,l)*ztla(:,l-1)
+   call thermcell_qsat(ngrid,active,pplev(:,l),ztemp,zqta(:,l-1),zqsat)
+
 
     do ig=1,ngrid
         if(active(ig)) then
+        zqla_est(ig,l) = max(0.,zqta(ig,l-1)-zqsat(ig))
         ztva_est(ig,l) = ztla(ig,l-1)*zpspsk(ig,l)+RLvCp*zqla_est(ig,l)
         zta_est(ig,l)=ztva_est(ig,l)
@@ -366 +370 @@
     enddo
 
-   call thermcell_condens(ngrid,activetmp,zpspsk(:,l),pplev(:,l),ztla(:,l),zqta(:,l),zqla(:,l))
-
+   ztemp(:)=zpspsk(:,l)*ztla(:,l)
+   call thermcell_qsat(ngrid,activetmp,pplev(:,l),ztemp,zqta(:,l),zqsatth(:,l))
 
    do ig=1,ngrid
@@ -373 +377 @@
         if (prt_level.ge.20) print*,'coucou calcul detr 4512: ig, l', ig, l
 ! on ecrit de maniere conservative (sat ou non)
+           zqla(ig,l) = max(0.,zqta(ig,l)-zqsatth(ig,l))
 !          T = Tl +Lv/Cp ql
            ztva(ig,l) = ztla(ig,l)*zpspsk(ig,l)+RLvCp*zqla(ig,l)
@@ -382 +387 @@
 
 !on ecrit zqsat 
-           zqsatth(ig,l)=qsatbef  
 
 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
@@ -554 +558 @@
       REAL zqsatth(ngrid,klev)
       REAL zta_est(ngrid,klev)
+      REAL ztemp(ngrid),zqsat(ngrid)
       REAL zdw2
       REAL zw2modif
@@ -570 +575 @@
       real zcor,zdelta,zcvm5,qlbef,zdz2
       real betalpha,zbetalpha
-      real Tbef,qsatbef,b1,eps, afact
-      real dqsat_dT,DT,num,denom,m
+      real eps, afact
       REAL REPS,RLvCp,DDT0
       PARAMETER (DDT0=.01)
@@ -589 +593 @@
 
 !      print*,'THERM 31B'
+      print*,'THERMCELL_PLUME OPTIMISE V1 CCC '
 
 ! Initialisations des variables reeles
@@ -624 +629 @@
       linter(:)=1.
 !     linter(:)=1.
-      b1=2.
 ! Initialisation des variables entieres
       lmix(:)=1
@@ -631 +635 @@
       lalim(:)=1
 
-!      print*,'THERMCELL PLUME ARNAUD DEDANS'
+   print*,'THERMCELL PLUME QSAT2 NDDDDN'
 
 !-------------------------------------------------------------------------
@@ -708 +712 @@
 ! sans tenir compte du detrainement et de l'entrainement dans cette
 ! couche
+! C'est a dire qu'on suppose 
+! ztla(l)=ztla(l-1) et zqta(l)=zqta(l-1)
 ! Ici encore, on doit pouvoir ajouter entr_star (qui peut etre calculer
 ! avant) a l'alimentation pour avoir un calcul plus propre
 !---------------------------------------------------------------------------
 
-     call thermcell_condens(ngrid,active, &
-&          zpspsk(:,l),pplev(:,l),ztla(:,l-1),zqta(:,l-1),zqla_est(:,l))
-
+   ztemp(:)=zpspsk(:,l)*ztla(:,l-1)
+   call thermcell_qsat(ngrid,active,pplev(:,l),ztemp,zqta(:,l-1),zqsat(:))
 
     do ig=1,ngrid 
 !       print*,'active',active(ig),ig,l
         if(active(ig)) then 
+        zqla_est(ig,l)=max(0.,zqta(ig,l-1)-zqsat(ig))
         ztva_est(ig,l) = ztla(ig,l-1)*zpspsk(ig,l)+RLvCp*zqla_est(ig,l)
         zta_est(ig,l)=ztva_est(ig,l)
@@ -747 +753 @@
 !-------------------------------------------------
 
+     print*,'THERM V1 SANS DQ'
      do ig=1,ngrid
         if (active(ig)) then
@@ -758 +765 @@
           zalpha=f0(ig)*f_star(ig,l)/sqrt(w_est(ig,l+1))/rhobarz(ig,l)
           zdqt(ig,l)=max(zqta(ig,l-1)-po(ig,l),0.)/po(ig,l)
+          zdqt(ig,l)=0.
 
           
@@ -800 +808 @@
     enddo
 
-   call thermcell_condens(ngrid,activetmp,zpspsk(:,l),pplev(:,l),ztla(:,l),zqta(:,l),zqla(:,l))
-
+   ztemp(:)=zpspsk(:,l)*ztla(:,l)
+   call thermcell_qsat(ngrid,activetmp,pplev(:,l),ztemp,zqta(:,l),zqsatth(:,l))
 
    do ig=1,ngrid
@@ -808 +816 @@
 ! on ecrit de maniere conservative (sat ou non)
 !          T = Tl +Lv/Cp ql
+           zqla(ig,l)=max(0.,zqta(ig,l)-zqsatth(ig,l))
            ztva(ig,l) = ztla(ig,l)*zpspsk(ig,l)+RLvCp*zqla(ig,l)
            ztva(ig,l) = ztva(ig,l)/zpspsk(ig,l)
@@ -814 +823 @@
            ztva(ig,l) = ztva(ig,l)*(1.+RETV*(zqta(ig,l)  &
      &              -zqla(ig,l))-zqla(ig,l))
-
-!on ecrit zqsat 
-           zqsatth(ig,l)=qsatbef  
-
            zbuoy(ig,l)=RG*(ztva(ig,l)-ztv(ig,l))/ztv(ig,l)
            zdz=zlev(ig,l+1)-zlev(ig,l)