Index: LMDZ6/branches/contrails/libf/phylmd/lmdz_cloudth.f90
===================================================================
--- LMDZ6/branches/contrails/libf/phylmd/lmdz_cloudth.f90	(revision 5418)
+++ LMDZ6/branches/contrails/libf/phylmd/lmdz_cloudth.f90	(revision 5489)
@@ -69,5 +69,4 @@
       REAL zpdf_a(ngrid),zpdf_b(ngrid),zpdf_e1(ngrid),zpdf_e2(ngrid) 
       REAL zqs(ngrid), qcloud(ngrid)
-      REAL erf
 
 
@@ -91,14 +90,14 @@
 ! Initialisation des variables r?elles
 !-------------------------------------------------------------------------------
-      sigma1(:,:)=0.
-      sigma2(:,:)=0.
-      qlth(:,:)=0.
-      qlenv(:,:)=0.  
-      qltot(:,:)=0.
-      rneb(:,:)=0.
+      sigma1(:,ind2)=0.
+      sigma2(:,ind2)=0.
+      qlth(:,ind2)=0.
+      qlenv(:,ind2)=0.  
+      qltot(:,ind2)=0.
+      rneb(:,ind2)=0.
       qcloud(:)=0.
-      cth(:,:)=0.
-      cenv(:,:)=0.
-      ctot(:,:)=0.
+      cth(:,ind2)=0.
+      cenv(:,ind2)=0.
+      ctot(:,ind2)=0.
       qsatmmussig1=0.
       qsatmmussig2=0.
@@ -317,19 +316,18 @@
       REAL zpdf_a(ngrid),zpdf_b(ngrid),zpdf_e1(ngrid),zpdf_e2(ngrid) 
       REAL zqs(ngrid), qcloud(ngrid)
-      REAL erf
 
 !------------------------------------------------------------------------------
 ! Initialisation des variables r?elles
 !------------------------------------------------------------------------------
-      sigma1(:,:)=0.
-      sigma2(:,:)=0.
-      qlth(:,:)=0.
-      qlenv(:,:)=0.  
-      qltot(:,:)=0.
-      rneb(:,:)=0.
+      sigma1(:,ind2)=0.
+      sigma2(:,ind2)=0.
+      qlth(:,ind2)=0.
+      qlenv(:,ind2)=0.  
+      qltot(:,ind2)=0.
+      rneb(:,ind2)=0.
       qcloud(:)=0.
-      cth(:,:)=0.
-      cenv(:,:)=0.
-      ctot(:,:)=0.
+      cth(:,ind2)=0.
+      cenv(:,ind2)=0.
+      ctot(:,ind2)=0.
       qsatmmussig1=0.
       qsatmmussig2=0.
@@ -644,5 +642,4 @@
       REAL zpdf_sig(ngrid),zpdf_k(ngrid),zpdf_delta(ngrid)
       REAL zpdf_a(ngrid),zpdf_b(ngrid),zpdf_e1(ngrid),zpdf_e2(ngrid) 
-      REAL erf
 
 
@@ -663,17 +660,17 @@
 ! Initialisation des variables r?elles
 !-------------------------------------------------------------------------------
-      sigma1(:,:)=0.
-      sigma2(:,:)=0.
-      qlth(:,:)=0.
-      qlenv(:,:)=0.  
-      qltot(:,:)=0.
-      rneb(:,:)=0.
+      sigma1(:,ind2)=0.
+      sigma2(:,ind2)=0.
+      qlth(:,ind2)=0.
+      qlenv(:,ind2)=0.  
+      qltot(:,ind2)=0.
+      rneb(:,ind2)=0.
       qcloud(:)=0.
-      cth(:,:)=0.
-      cenv(:,:)=0.
-      ctot(:,:)=0.
-      cth_vol(:,:)=0.
-      cenv_vol(:,:)=0.
-      ctot_vol(:,:)=0.
+      cth(:,ind2)=0.
+      cenv(:,ind2)=0.
+      ctot(:,ind2)=0.
+      cth_vol(:,ind2)=0.
+      cenv_vol(:,ind2)=0.
+      ctot_vol(:,ind2)=0.
       qsatmmussig1=0.
       qsatmmussig2=0.
@@ -878,5 +875,4 @@
       REAL zpdf_a(ngrid),zpdf_b(ngrid),zpdf_e1(ngrid),zpdf_e2(ngrid) 
       REAL zqs(ngrid), qcloud(ngrid)
-      REAL erf
 
       REAL rhodz(ngrid,klev)
@@ -895,17 +891,17 @@
 !------------------------------------------------------------------------------
 
-      sigma1(:,:)=0.
-      sigma2(:,:)=0.
-      qlth(:,:)=0.
-      qlenv(:,:)=0.  
-      qltot(:,:)=0.
-      rneb(:,:)=0.
+      sigma1(:,ind2)=0.
+      sigma2(:,ind2)=0.
+      qlth(:,ind2)=0.
+      qlenv(:,ind2)=0.  
+      qltot(:,ind2)=0.
+      rneb(:,ind2)=0.
       qcloud(:)=0.
-      cth(:,:)=0.
-      cenv(:,:)=0.
-      ctot(:,:)=0.
-      cth_vol(:,:)=0.
-      cenv_vol(:,:)=0.
-      ctot_vol(:,:)=0.
+      cth(:,ind2)=0.
+      cenv(:,ind2)=0.
+      ctot(:,ind2)=0.
+      cth_vol(:,ind2)=0.
+      cenv_vol(:,ind2)=0.
+      ctot_vol(:,ind2)=0.
       qsatmmussig1=0.
       qsatmmussig2=0.
@@ -1306,5 +1302,5 @@
       REAL qcloud(ngrid) !eau totale dans le nuage
         !Some arithmetic variables
-      REAL erf,pi,sqrt2,sqrt2pi
+      REAL  pi,sqrt2,sqrt2pi
         !Depth of the layer
       REAL dz(ngrid,klev)    !epaisseur de la couche en metre
@@ -1320,13 +1316,13 @@
 ! Initialization
 !------------------------------------------------------------------------------
-      qlth(:,:)=0.
-      qlenv(:,:)=0.  
-      qltot(:,:)=0.
-      cth_vol(:,:)=0.
-      cenv_vol(:,:)=0.
-      ctot_vol(:,:)=0.
-      cth_surf(:,:)=0.
-      cenv_surf(:,:)=0.
-      ctot_surf(:,:)=0.
+      qlth(:,ind2)=0.
+      qlenv(:,ind2)=0.  
+      qltot(:,ind2)=0.
+      cth_vol(:,ind2)=0.
+      cenv_vol(:,ind2)=0.
+      ctot_vol(:,ind2)=0.
+      cth_surf(:,ind2)=0.
+      cenv_surf(:,ind2)=0.
+      ctot_surf(:,ind2)=0.
       qcloud(:)=0.
       rdd=287.04
@@ -1579,5 +1575,4 @@
       REAL qlbef
       REAL dqsatenv(klon), dqsatth(klon)
-      REAL erf
       REAL zpdf_sig(klon),zpdf_k(klon),zpdf_delta(klon)
       REAL zpdf_a(klon),zpdf_b(klon),zpdf_e1(klon),zpdf_e2(klon) 
