Index: LMDZ6/trunk/libf/phylmd/lscp_ini_mod.F90
===================================================================
--- LMDZ6/trunk/libf/phylmd/lscp_ini_mod.F90	(revision 4562)
+++ LMDZ6/trunk/libf/phylmd/lscp_ini_mod.F90	(revision 4563)
@@ -17,5 +17,7 @@
   INTEGER,SAVE :: iflag_evap_prec=1             ! precipitation evaporation flag. 0: nothing, 1: "old way", 
                                                 ! 2: Max cloud fraction above to calculate the max of reevaporation
-                                                ! 4: LTP'method i.e. evaporation in the clear-sky fraction of the mesh only
+                                                ! >=4: LTP'method i.e. evaporation in the clear-sky fraction of the mesh only
+                                                ! pay attention that iflag_evap_prec=4 may lead to unrealistic and overestimated
+                                                ! evaporation. Use 5 instead
   !$OMP THREADPRIVATE(iflag_evap_prec)
 
Index: LMDZ6/trunk/libf/phylmd/lscp_mod.F90
===================================================================
--- LMDZ6/trunk/libf/phylmd/lscp_mod.F90	(revision 4562)
+++ LMDZ6/trunk/libf/phylmd/lscp_mod.F90	(revision 4563)
@@ -415,8 +415,8 @@
             IF (zrfl(i)+zifl(i).GT.0.) THEN
 
-            ! LudoTP: we only account for precipitation evaporation in the clear-sky (iflag_evap_prec=4).
+            ! LudoTP: we only account for precipitation evaporation in the clear-sky (iflag_evap_prec>=4).
             ! c_iso: likely important to distinguish cs from neb isotope precipitation
 
-                IF (iflag_evap_prec.EQ.4) THEN
+                IF (iflag_evap_prec.GE.4) THEN
                     zrfl(i) = zrflclr(i)
                     zifl(i) = ziflclr(i)
@@ -429,6 +429,11 @@
                 ENDIF
 
-                IF (iflag_evap_prec.EQ.4) THEN
+                IF (iflag_evap_prec.GT.4) THEN
+                ! Max evaporation not to saturate the clear sky precip fraction
+                ! i.e. the fraction where evaporation occurs
+                    zqev0 = MAX(0.0, (zqs(i)-zq(i))*znebprecipclr(i))
+                ELSEIF (iflag_evap_prec .EQ. 4) THEN
                 ! Max evaporation not to saturate the whole mesh
+                ! Pay attention -> lead to unrealistic and excessive evaporation
                     zqev0 = MAX(0.0, zqs(i)-zq(i))
                 ELSE
@@ -440,5 +445,5 @@
                 ! dP/dz=beta*(1-q/qsat)*sqrt(P)
                 ! formula from Sundquist 1988, Klemp & Wilhemson 1978
-                ! LTP: evaporation only in the clear sky part (iflag_evap_prec=4)
+                ! LTP: evaporation only in the clear sky part (iflag_evap_prec>=4)
 
                 IF (iflag_evap_prec.EQ.3) THEN
@@ -446,5 +451,5 @@
                     *SQRT(zrfl(i)/max(1.e-4,znebprecip(i)))        &
                     *(paprs(i,k)-paprs(i,k+1))/pplay(i,k)*zt(i)*RD/RG
-                ELSE IF (iflag_evap_prec.EQ.4) THEN
+                ELSE IF (iflag_evap_prec.GE.4) THEN
                      zqevt = znebprecipclr(i)*coef_eva*(1.0-zq(i)/qsl(i)) &
                     *SQRT(zrfl(i)/max(1.e-8,znebprecipclr(i))) &
@@ -463,5 +468,5 @@
                     *SQRT(zifl(i)/max(1.e-4,znebprecip(i))) &
                     *(paprs(i,k)-paprs(i,k+1))/pplay(i,k)*zt(i)*RD/RG
-                ELSE IF (iflag_evap_prec.EQ.4) THEN
+                ELSE IF (iflag_evap_prec.GE.4) THEN
                      zqevti = znebprecipclr(i)*coef_eva_i*(1.0-zq(i)/qsi(i)) &
                     *SQRT(zifl(i)/max(1.e-8,znebprecipclr(i))) &
@@ -518,5 +523,5 @@
                                            
 
-                IF (iflag_evap_prec.EQ.4) THEN
+                IF (iflag_evap_prec.GE.4) THEN
                     zrflclr(i) = zrfl(i)
                     ziflclr(i) = zifl(i)
@@ -536,5 +541,5 @@
 
                 ! update of rainfall and snowfall due to melting
-                IF (iflag_evap_prec.EQ.4) THEN
+                IF (iflag_evap_prec.GE.4) THEN
                     zrflclr(i)=zrflclr(i)+zmelt*ziflclr(i)
                     zrflcld(i)=zrflcld(i)+zmelt*ziflcld(i)
@@ -885,5 +890,5 @@
 
     ! LTP:
-    IF (iflag_evap_prec==4) THEN
+    IF (iflag_evap_prec .GE. 4) THEN
 
         !Partitionning between precipitation coming from clouds and that coming from CS
@@ -1077,5 +1082,5 @@
     DO i = 1, klon
             
-            IF (iflag_evap_prec.EQ.4) THEN
+            IF (iflag_evap_prec.GE.4) THEN
                 ziflprev(i)=ziflcld(i)
             ELSE
@@ -1119,5 +1124,5 @@
                     ! c_iso : mv here condensation of isotopes + redispatchage en precip 
 
-                IF (iflag_evap_prec.EQ.4) THEN
+                IF (iflag_evap_prec.GE.4) THEN
                     zrflcld(i) = zrflcld(i)+zqprecl(i) &
                     *(paprs(i,k)-paprs(i,k+1))/(RG*dtime)
@@ -1140,6 +1145,6 @@
 
     ! LTP: limit of surface cloud fraction covered by precipitation when the local intensity of the flux is below rain_int_min
-    ! if iflag_evap_pre=4
-    IF (iflag_evap_prec.EQ.4) THEN
+    ! if iflag_evap_prec>=4
+    IF (iflag_evap_prec.GE.4) THEN
 
         DO i=1,klon