Changeset 3500

Timestamp:

May 13, 2019, 7:50:31 PM (6 years ago)

Author:

jyg

Message:

Change in the computation of the ice fraction in unsaturated dondrafts, for backward compatibility

File:

• LMDZ6/trunk/libf/phylmd/cv3_routines.F90 (modified) (1 diff)

LMDZ6/trunk/libf/phylmd/cv3_routines.F90

@@ -2962 +2962 @@
 
     IF (cvflag_ice) THEN                                                                                !!jygprl
-      DO il = 1, ncum                                                                                   !!jygprl
-        IF (i<=inb(il) .AND. lwork(il)) THEN                                                            !!jygprl
-          frac(il, i) = (frac_a(il,i)*wdtrainA(il,i)+frac_s(il,i)*(wdtrainS(il,i)+wdtrainM(il,i))) / &  !!jygprl
-                        max(wdtrainA(il,i)+wdtrainS(il,i)+wdtrainM(il,i),smallestreal)                  !!jygprl
-          fraci(il, i) = frac(il, i)                                                                    !!jygprl
-        END IF                                                                                          !!jygprl
-      END DO                                                                                            !!jygprl
+      IF (cvflag_prec_eject) THEN
+        DO il = 1, ncum                                                                                   !!jygprl
+          IF (i<=inb(il) .AND. lwork(il)) THEN                                                            !!jygprl
+            frac(il, i) = (frac_a(il,i)*wdtrainA(il,i)+frac_s(il,i)*(wdtrainS(il,i)+wdtrainM(il,i))) / &  !!jygprl
+                          max(wdtrainA(il,i)+wdtrainS(il,i)+wdtrainM(il,i),smallestreal)                  !!jygprl
+            fraci(il, i) = frac(il, i)                                                                    !!jygprl
+          END IF                                                                                          !!jygprl
+        END DO                                                                                            !!jygprl
+      ELSE  ! (cvflag_prec_eject)
+        DO il = 1, ncum                                                                                   !!jygprl
+          IF (i<=inb(il) .AND. lwork(il)) THEN                                                            !!jygprl
+          frac(il, inb(il)) = 1. - (t(il,inb(il))-243.15)/(263.15-243.15)
+          frac(il, inb(il)) = min(max(frac(il,inb(il)),0.), 1.)
+!       Ice fraction computed again here as a function of the temperature seen by unsaturated downdraughts
+!       (i.e. the cold pool temperature) for compatibility with earlier versions.
+          fraci(il, i) = 1. - (t(il,i)-243.15)/(263.15-243.15)
+          fraci(il, i) = min(max(fraci(il,i),0.0), 1.0)
+          END IF                                                                                          !!jygprl
+        END DO                                                                                            !!jygprl
+      ENDIF  ! (cvflag_prec_eject)
     END IF                                                                                              !!jygprl