Changeset 5491 for LMDZ6/trunk/libf/phylmd/cv3p1_closure.f90

Timestamp:

Jan 19, 2025, 6:48:10 PM (5 hours ago)

Author:

jyg

Message:

New outputs :

+ coef_clos = [conv mass flux given by Alp closure]/[conv mass flux given by Emanuel scheme closure]
+ coef_clos_eff = effective coefficient used in the convective scheme.

File:

• LMDZ6/trunk/libf/phylmd/cv3p1_closure.f90 (modified) (4 diffs)

LMDZ6/trunk/libf/phylmd/cv3p1_closure.f90

@@ -4 +4 @@
 SUBROUTINE cv3p1_closure(nloc, ncum, nd, icb, inb, pbase, plcl, p, ph, tv, &
     tvp, buoy, supmax, ok_inhib, ale, alp, omega,sig, w0, ptop2, cape, cin, m, &
-    iflag, coef, plim1, plim2, asupmax, supmax0, asupmaxmin, cbmf, plfc, &
-    wbeff)
+    iflag, coef, coeftrue, plim1, plim2, asupmax, supmax0, asupmaxmin, &
+    cbmf, plfc, wbeff)
 
 
@@ -48 +48 @@
   REAL, DIMENSION (nloc), INTENT (OUT)               :: cape, cin
   REAL, DIMENSION (nloc, nd), INTENT (OUT)           :: m
+  REAL, DIMENSION (nloc), INTENT (OUT)               :: coef, coeftrue
   REAL, DIMENSION (nloc), INTENT (OUT)               :: plim1, plim2
   REAL, DIMENSION (nloc, nd), INTENT (OUT)           :: asupmax
@@ -74 +75 @@
   REAL cbmflim(nloc), cbmf1(nloc), cbmfmax(nloc)
   REAL cbmflast(nloc)
-  REAL coef(nloc)
   REAL xp(nloc), xq(nloc), xr(nloc), discr(nloc), b3(nloc), b4(nloc)
   REAL theta(nloc), bb(nloc)
@@ -598 +598 @@
   DO il = 1, ncum
     coef(il) = (cbmf(il)+1.E-10)/(cbmflim(il)+1.E-10)
+    coeftrue(il) = coef(il)
   END DO
   IF (prt_level>=20) PRINT *, 'cv3p1_param apres coef_plantePLUS'