Changeset 2911 for LMDZ5/trunk/libf/phylmd/cloudth_mod.F90

Timestamp:

Jun 13, 2017, 7:31:16 PM (7 years ago)

Author:

jbmadeleine

Message:

Changes in the vertical PDF used for shallow convective clouds.
Added vert_alpha_th which is the variance of the vertical PDF used for thermals.
If nothing is specified in the def files, it is set equal to vert_alpha.
Otherwise, it is now possible to have a different vertical PDF in the thermals
and in the environment.
Also added a fourth option iflag_cloudth_vert = 4 where
deltas=vert_alpha*sigmas instead of deltas=a*vert_alpha*sigmas.

File:

• LMDZ5/trunk/libf/phylmd/cloudth_mod.F90 (modified) (5 diffs)

LMDZ5/trunk/libf/phylmd/cloudth_mod.F90

@@ -860 +860 @@
       ! Change the width of the PDF used for vertical subgrid scale heterogeneity
       ! (J Jouhaud, JL Dufresne, JB Madeleine)
-      REAL,SAVE :: vert_alpha
-      !$OMP THREADPRIVATE(vert_alpha)
+      REAL,SAVE :: vert_alpha, vert_alpha_th
+      !$OMP THREADPRIVATE(vert_alpha, vert_alpha_th)
       LOGICAL, SAVE :: firstcall = .TRUE.
       !$OMP THREADPRIVATE(firstcall)
@@ -897 +897 @@
         CALL getin_p('cloudth_vert_alpha',vert_alpha)
         WRITE(*,*) 'cloudth_vert_alpha = ', vert_alpha
+        ! The factor used for the thermal is equal to that of the environment
+        !   if nothing is explicitly specified in the def file
+        vert_alpha_th=vert_alpha
+        CALL getin_p('cloudth_vert_alpha_th',vert_alpha_th)
+        WRITE(*,*) 'cloudth_vert_alpha_th = ', vert_alpha_th
         firstcall=.FALSE.
       ENDIF
@@ -1001 +1006 @@
       qltot(ind1,ind2)=fraca(ind1,ind2)*qlth(ind1,ind2)+(1.-1.*fraca(ind1,ind2))*qlenv(ind1,ind2)
 
-      ELSE IF (iflag_cloudth_vert == 3) THEN
+      ELSE IF (iflag_cloudth_vert >= 3) THEN
 
 !-------------------------------------------------------------------------------
@@ -1010 +1015 @@
 !      deltasenv=aenv*ratqs(ind1,ind2)*zqsatenv(ind1,ind2)
 !      deltasth=ath*ratqs(ind1,ind2)*zqsatth(ind1,ind2)
-      deltasenv=aenv*vert_alpha*sigma1s
-      deltasth=ath*vert_alpha*sigma2s
+      IF (iflag_cloudth_vert == 3) THEN
+        deltasenv=aenv*vert_alpha*sigma1s
+        deltasth=ath*vert_alpha_th*sigma2s
+      ELSE IF (iflag_cloudth_vert == 4) THEN
+        deltasenv=vert_alpha*sigma1s
+        deltasth=vert_alpha_th*sigma2s
+      ENDIF
       
       xenv1=-(senv+deltasenv)/(sqrt(2.)*sigma1s)
@@ -1053 +1063 @@
 
 
-      ENDIF ! of if (iflag_cloudth_vert==1 or 3)
+      ENDIF ! of if (iflag_cloudth_vert==1 or 3 or 4)
 
       if (cenv(ind1,ind2).lt.1.e-10.or.cth(ind1,ind2).lt.1.e-10) then