Index: LMDZ6/trunk/libf/phylmd/cloudth_mod.F90
===================================================================
--- LMDZ6/trunk/libf/phylmd/cloudth_mod.F90	(revision 4613)
+++ LMDZ6/trunk/libf/phylmd/cloudth_mod.F90	(revision 4617)
@@ -30,5 +30,5 @@
       INTEGER itap,ind1,ind2
       INTEGER ngrid,klev,klon,l,ig 
-      INTEGER iflag_ratqs
+!ym      INTEGER iflag_ratqs
       
       REAL ztv(ngrid,klev)
@@ -275,5 +275,5 @@
       INTEGER itap,ind1,ind2
       INTEGER ngrid,klev,klon,l,ig 
-      INTEGER iflag_ratqs
+!ym      INTEGER iflag_ratqs
       
       REAL ztv(ngrid,klev)
@@ -615,5 +615,5 @@
       INTEGER itap,ind1,ind2
       INTEGER ngrid,klev,klon,l,ig 
-      INTEGER iflag_ratqs
+!ym      INTEGER iflag_ratqs
       
       REAL ztv(ngrid,klev)
@@ -838,5 +838,5 @@
       INTEGER itap,ind1,ind2
       INTEGER ngrid,klev,klon,l,ig 
-      INTEGER iflag_ratqs
+!ym      INTEGER iflag_ratqs
       
       REAL ztv(ngrid,klev)
@@ -1026,7 +1026,7 @@
       ENDIF
       sigma1s = sigma1s_fraca + sigma1s_ratqs
-           IF (iflag_ratqs.eq.11) then
-              sigma1s = ratqs(ind1,ind2)*po(ind1)*aenv
-           ENDIF
+!ym           IF (iflag_ratqs.eq.11) then
+!ym              sigma1s = ratqs(ind1,ind2)*po(ind1)*aenv
+!ym           ENDIF
       sigma2s=(sigma2s_factor*(((sth-senv)**2)**0.5)/((fraca(ind1,ind2)+0.02)**sigma2s_power))+0.002*zqta(ind1,ind2)
 !      tests
@@ -1607,5 +1607,5 @@
 
       INTEGER itap,ind1,l,ig,iter,k 
-      INTEGER iflag_ratqs
+!      INTEGER iflag_ratqs
       INTEGER iflag_topthermals, niter
       LOGICAL falseklon(klon)
@@ -1902,7 +1902,7 @@
  
            sigma1s = sigma1s_fraca + sigma1s_ratqs
-           IF (iflag_ratqs.eq.11) then
-              sigma1s = ratqs(ind1,ind2)*po(ind1)*aenv
-           ENDIF
+!ym           IF (iflag_ratqs.eq.11) then
+!ym              sigma1s = ratqs(ind1,ind2)*po(ind1)*aenv
+!ym           ENDIF
            sigma2s=(sigma2s_factor*(((sth-senv)**2)**0.5)/((fraca(ind1,ind2)+0.02)**sigma2s_power))+0.002*zqta(ind1,ind2)
 
@@ -2147,10 +2147,10 @@
 
             sigma1s = sigma1s_fraca + sigma1s_ratqs
-           IF (iflag_ratqs.eq.11) then
-              sigma1s = ratqs(ind1,ind2)*po(ind1)*aenv
-           ENDIF
-           IF (iflag_ratqs.eq.11) then
-              sigma1s = ratqs(ind1,ind2)*po(ind1)*aenvl
-           ENDIF
+!ym           IF (iflag_ratqs.eq.11) then
+!ym              sigma1s = ratqs(ind1,ind2)*po(ind1)*aenv
+!ym           ENDIF
+!ym           IF (iflag_ratqs.eq.11) then
+!ym             sigma1s = ratqs(ind1,ind2)*po(ind1)*aenvl
+!ym           ENDIF
             deltasenv=aenvl*vert_alpha*sigma1s
             xenv1=-(senvl+deltasenv)/(sqrt(2.)*sigma1s)