Changeset 2960 for LMDZ5/trunk/libf

Timestamp:

Jul 19, 2017, 12:19:26 PM (7 years ago)

Author:

fhourdin

Message:

Un lievre leve, sans doute.
Un tout petit seuil sur les largeurs de distributions sous maille
du schéma de nuage, qui produit un gros effet sur le contraste
strato-cumulus / cumulus ...

File:

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

LMDZ5/trunk/libf/phylmd/cloudth_mod.F90

@@ -876 +876 @@
       REAL,SAVE :: sigma1s_factor=1.1
       REAL,SAVE :: sigma1s_power=0.6
-      !$OMP THREADPRIVATE(sigma1s_factor,sigma1s_power)
+      REAL,SAVE :: cloudth_ratqsmin=-1.
+      !$OMP THREADPRIVATE(sigma1s_factor,sigma1s_power,cloudth_ratqsmin)
       INTEGER, SAVE :: iflag_cloudth_vert_noratqs=0
       !$OMP THREADPRIVATE(iflag_cloudth_vert_noratqs)
@@ -929 +930 @@
         CALL getin_p('cloudth_sigma1s_power',sigma1s_power)
         WRITE(*,*) 'cloudth_sigma1s_power = ', sigma1s_power
+        ! Minimum value for the environmental air subgrid water distrib
+        CALL getin_p('cloudth_ratqsmin',cloudth_ratqsmin)
+        WRITE(*,*) 'cloudth_ratqsmin = ', cloudth_ratqsmin
         ! Remove the dependency to ratqs from the variance of the vertical PDF
         CALL getin_p('iflag_cloudth_vert_noratqs',iflag_cloudth_vert_noratqs)
@@ -988 +992 @@
      &                (1-fraca(ind1,ind2))*((sth-senv)**2)**0.5
 !     sigma1s_fraca = (1.1**0.5)*(fraca(ind1,ind2)**0.6)/(1-fraca(ind1,ind2))*((sth-senv)**2)**0.5
-      sigma1s_ratqs = ratqs(ind1,ind2)*po(ind1)
+      IF (cloudth_ratqsmin>0.) THEN
+         sigma1s_ratqs = cloudth_ratqsmin*po(ind1)
+      ELSE
+         sigma1s_ratqs = ratqs(ind1,ind2)*po(ind1)
+      ENDIF
       sigma1s = sigma1s_fraca + sigma1s_ratqs
       sigma2s=(0.09*(((sth-senv)**2)**0.5)/((fraca(ind1,ind2)+0.02)**0.5))+0.002*zqta(ind1,ind2)