Index: /LMDZ5/trunk/libf/phylmd/fisrtilp.F90
===================================================================
--- /LMDZ5/trunk/libf/phylmd/fisrtilp.F90	(revision 2506)
+++ /LMDZ5/trunk/libf/phylmd/fisrtilp.F90	(revision 2507)
@@ -623,5 +623,5 @@
                     if (iflag_t_glace.eq.0) then
                     zdelta = MAX(0.,SIGN(1.,t_glace_min_old-Tbef(i)))
-                    else if (iflag_t_glace.eq.1) then
+                    else if (iflag_t_glace.ge.1) then
                     zdelta = MAX(0.,SIGN(1.,t_glace_min-Tbef(i)))
                     endif
@@ -673,5 +673,5 @@
                  else
                  ! Iteration pour convergence avec qsat(T)
-                 if (iflag_t_glace.eq.1) then
+                 if (iflag_t_glace.ge.1) then
                  CALL icefrac_lsc(klon,zt(:),pplay(:,k)/paprs(:,1),zfice(:))
                  endif
@@ -687,5 +687,5 @@
                  endif
                  
-                 if (iflag_t_glace.eq.1) then
+                 if (iflag_t_glace.ge.1) then
                  dzfice(i)= exposant_glace * zfice(i)**(exposant_glace-1) / (t_glace_min - t_glace_max)
                  endif
@@ -811,5 +811,5 @@
         endif
      ELSE
-         if (iflag_t_glace.eq.1) then
+         if (iflag_t_glace.ge.1) then
             CALL icefrac_lsc(klon,zt(:),pplay(:,k)/paprs(:,1),zfice(:))
          endif
Index: /LMDZ5/trunk/libf/phylmd/icefrac_lsc_mod.F90
===================================================================
--- /LMDZ5/trunk/libf/phylmd/icefrac_lsc_mod.F90	(revision 2506)
+++ /LMDZ5/trunk/libf/phylmd/icefrac_lsc_mod.F90	(revision 2507)
@@ -14,6 +14,7 @@
   ! Doutriaux-Boucher & Quaas 2004, section 2.2.)
   !
-  ! (JBM 3/14 8/14)
+  ! (JBM 3/14 8/14 5/16)
   
+  USE print_control_mod, ONLY: lunout, prt_level
   INCLUDE "nuage.h"
 
@@ -33,5 +34,12 @@
 
   DO ip=1,np
-     www=(max(sig(ip)-sig0,0.))/(1.-sig0) ! w=1 at the surface and 0 for sig < sig0
+     IF (iflag_t_glace.EQ.1) THEN
+       ! Transition to ice close to surface for T<Tmax
+       ! w=1 at the surface and 0 for sig < sig0
+       www=(max(sig(ip)-sig0,0.))/(1.-sig0) 
+     ELSEIF (iflag_t_glace.GE.2) THEN
+       ! No convertion to ice close to surface
+       www = 0.
+     ENDIF
      tmin_tmp=www*t_glace_max+(1.-www)*t_glace_min
      icefrac_tmp= 1.0 - (temp(ip)-tmin_tmp) / (t_glace_max-tmin_tmp)