Index: /LMDZ5/trunk/libf/phylmd/icefrac_lsc_mod.F90
===================================================================
--- /LMDZ5/trunk/libf/phylmd/icefrac_lsc_mod.F90	(revision 2960)
+++ /LMDZ5/trunk/libf/phylmd/icefrac_lsc_mod.F90	(revision 2961)
@@ -28,5 +28,5 @@
   REAL, DIMENSION(np), INTENT(OUT) :: icefrac
 
-  REAL :: sig0,www,tmin_tmp,icefrac_tmp
+  REAL :: sig0,www,tmin_tmp,liqfrac_tmp
   INTEGER :: ip
 
@@ -43,7 +43,11 @@
      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)
-     icefrac_tmp = MIN(MAX(icefrac_tmp,0.0),1.0)
-     icefrac(ip) = icefrac_tmp**exposant_glace
+     liqfrac_tmp=  (temp(ip)-tmin_tmp) / (t_glace_max-tmin_tmp)
+     liqfrac_tmp = MIN(MAX(liqfrac_tmp,0.0),1.0)
+     IF (iflag_t_glace.GE.3) THEN
+       icefrac(ip) = 1.0-liqfrac_tmp**exposant_glace
+     ELSE
+       icefrac(ip) = (1.0-liqfrac_tmp)**exposant_glace
+     ENDIF
   ENDDO