Index: LMDZ5/trunk/libf/phylmd/cv3_routines.F90
===================================================================
--- LMDZ5/trunk/libf/phylmd/cv3_routines.F90	(revision 2374)
+++ LMDZ5/trunk/libf/phylmd/cv3_routines.F90	(revision 2376)
@@ -1074,4 +1074,5 @@
   REAL, DIMENSION (nloc, nd), INTENT (OUT)           :: ep, sigp, hp
   REAL, DIMENSION (nloc, nd), INTENT (OUT)           :: buoy
+  REAL, DIMENSION (nloc, nd), INTENT (OUT)           :: frac
 
 !local variables:
@@ -1082,5 +1083,4 @@
   REAL by, defrac, pden, tbis
   REAL ah0(nloc), cape(nloc), capem(nloc), byp(nloc)
-  REAL frac(nloc, nd)
   LOGICAL lcape(nloc)
   INTEGER iposit(nloc)
@@ -1532,4 +1532,12 @@
           hp(i, k) = hnk(i) + (lv(i,k)+(cpd-cpv)*t(i,k)+frac(i,k)*lf(i,k))* &
                               ep(i, k)*clw(i, k)
+        END IF
+      END DO
+    END DO
+! Below cloud base, set ice fraction to cloud base value
+    DO k = 1, nl
+      DO i = 1, ncum
+        IF (k<icb(i)) THEN
+          frac(i,k) = frac(i,icb(i))
         END IF
       END DO