Changeset 2376

Timestamp:

Oct 20, 2015, 11:57:39 AM (9 years ago)

Author:

jyg

Message:

Add initialization of ice fraction below cloud
base in cv3_undilute2 in cv3_routines.F90.

File:

• LMDZ5/trunk/libf/phylmd/cv3_routines.F90 (modified) (3 diffs)

LMDZ5/trunk/libf/phylmd/cv3_routines.F90

@@ -1074 +1074 @@
   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 +1083 @@
   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 +1532 @@
           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