Index: LMDZ6/trunk/libf/phylmd/cloudth_mod.F90
===================================================================
--- LMDZ6/trunk/libf/phylmd/cloudth_mod.F90	(revision 4622)
+++ LMDZ6/trunk/libf/phylmd/cloudth_mod.F90	(revision 4623)
@@ -25,4 +25,5 @@
 
       INCLUDE "YOMCST.h"
+      INCLUDE "nuage.h" ! iflag_ratqs
       INCLUDE "YOETHF.h"
       INCLUDE "FCTTRE.h"
@@ -30,5 +31,4 @@
       INTEGER itap,ind1,ind2
       INTEGER ngrid,klev,klon,l,ig 
-!ym      INTEGER iflag_ratqs
       
       REAL ztv(ngrid,klev)
@@ -270,4 +270,5 @@
 
       INCLUDE "YOMCST.h"
+      INCLUDE "nuage.h" ! iflag_ratqs 
       INCLUDE "YOETHF.h"
       INCLUDE "FCTTRE.h"
@@ -275,5 +276,4 @@
       INTEGER itap,ind1,ind2
       INTEGER ngrid,klev,klon,l,ig 
-!ym      INTEGER iflag_ratqs
       
       REAL ztv(ngrid,klev)
@@ -610,4 +610,5 @@
 
       INCLUDE "YOMCST.h"
+      INCLUDE "nuage.h" ! iflag ratq
       INCLUDE "YOETHF.h"
       INCLUDE "FCTTRE.h"
@@ -615,5 +616,4 @@
       INTEGER itap,ind1,ind2
       INTEGER ngrid,klev,klon,l,ig 
-!ym      INTEGER iflag_ratqs
       
       REAL ztv(ngrid,klev)
@@ -833,4 +833,5 @@
 
       INCLUDE "YOMCST.h"
+      INCLUDE "nuage.h" ! iflag ratq
       INCLUDE "YOETHF.h"
       INCLUDE "FCTTRE.h"
@@ -838,5 +839,4 @@
       INTEGER itap,ind1,ind2
       INTEGER ngrid,klev,klon,l,ig 
-!ym      INTEGER iflag_ratqs
       
       REAL ztv(ngrid,klev)
@@ -1026,7 +1026,7 @@
       ENDIF
       sigma1s = sigma1s_fraca + sigma1s_ratqs
-!ym           IF (iflag_ratqs.eq.11) then
-!ym              sigma1s = ratqs(ind1,ind2)*po(ind1)*aenv
-!ym           ENDIF
+      IF (iflag_ratqs.eq.11) then
+         sigma1s = ratqs(ind1,ind2)*po(ind1)*aenv
+      ENDIF
       sigma2s=(sigma2s_factor*(((sth-senv)**2)**0.5)/((fraca(ind1,ind2)+0.02)**sigma2s_power))+0.002*zqta(ind1,ind2)
 !      tests
@@ -1566,4 +1566,5 @@
 
       INCLUDE "YOMCST.h"
+      INCLUDE "nuage.h" ! iflag_ratqs
       INCLUDE "YOETHF.h"
       INCLUDE "FCTTRE.h"
@@ -1902,7 +1903,7 @@
  
            sigma1s = sigma1s_fraca + sigma1s_ratqs
-!ym           IF (iflag_ratqs.eq.11) then
-!ym              sigma1s = ratqs(ind1,ind2)*po(ind1)*aenv
-!ym           ENDIF
+           IF (iflag_ratqs.eq.11) then
+              sigma1s = ratqs(ind1,ind2)*po(ind1)*aenv
+           ENDIF
            sigma2s=(sigma2s_factor*(((sth-senv)**2)**0.5)/((fraca(ind1,ind2)+0.02)**sigma2s_power))+0.002*zqta(ind1,ind2)
 
@@ -2147,10 +2148,10 @@
 
             sigma1s = sigma1s_fraca + sigma1s_ratqs
-!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
+            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
             deltasenv=aenvl*vert_alpha*sigma1s
             xenv1=-(senvl+deltasenv)/(sqrt(2.)*sigma1s)