Index: LMDZ5/trunk/libf/phylmd/cloudth_mod.F90
===================================================================
--- LMDZ5/trunk/libf/phylmd/cloudth_mod.F90	(revision 2959)
+++ LMDZ5/trunk/libf/phylmd/cloudth_mod.F90	(revision 2960)
@@ -876,5 +876,6 @@
       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,4 +930,7 @@
         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,5 +992,9 @@
      &                (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)