Changeset 5491 for LMDZ6/trunk/libf/phylmd/cva_driver.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/cva_driver.f90 (modified) (9 diffs)

LMDZ6/trunk/libf/phylmd/cva_driver.f90

@@ -21 +21 @@
                       ftd1, fqd1, &
                       Plim11, Plim21, asupmax1, supmax01, asupmaxmin1, &
+                      coef_clos1, coef_clos_eff1, &
                       lalim_conv1, & 
 !!                      da1,phi1,mp1,phi21,d1a1,dam1,sigij1,clw1, &        ! RomP
@@ -249 +250 @@
   REAL, DIMENSION (len), INTENT (OUT)                :: supmax01
   REAL, DIMENSION (len), INTENT (OUT)                :: asupmaxmin1
+  REAL, DIMENSION (len), INTENT (OUT)                :: coef_clos1, coef_clos_eff1
   REAL, DIMENSION (len, nd), INTENT (OUT)            :: qtc1    ! in cloud water content (intensive)   ! cld
   REAL, DIMENSION (len, nd), INTENT (OUT)            :: sigt1   ! fract. cloud area (intensive)        ! cld
@@ -495 +497 @@
   REAL elij(nloc, nd, nd)
   REAL supmax(nloc, nd)
-  REAL Ale(nloc), Alp(nloc), coef_clos(nloc)
+  REAL Ale(nloc), Alp(nloc), coef_clos(nloc), coef_clos_eff(nloc)
   REAL omega(nloc,nd)
   REAL sigd(nloc)
@@ -679 +681 @@
   DO il = 1, nloc
     coef_clos(il) = 1.
+    coef_clos_eff(il) = 1.
   END DO
 
@@ -1003 +1006 @@
                            pbase, plcl, p, ph, tv, tvp, buoy, &
                            supmax, ok_inhib, Ale, Alp, omega, &
-                           sig, w0, ptop2, cape, cin, m, iflag, coef_clos, &
+                           sig, w0, ptop2, cape, cin, m, iflag, &
+                           coef_clos_eff, coef_clos, &
                            Plim1, plim2, asupmax, supmax0, &
                            asupmaxmin, cbmf, plfc, wbeff)
@@ -1016 +1020 @@
                            pbase, plcl, p, ph, tv, tvp, buoy, &
                            supmax, ok_inhib, Ale, Alp, omega, &
-                           sig, w0, ptop2, cape, cin, m, iflag, coef_clos, &
+                           sig, w0, ptop2, cape, cin, m, iflag, coef_clos_eff, &
                            Plim1, plim2, asupmax, supmax0, &
                            asupmaxmin, cbmf, plfc, wbeff)
@@ -1091 +1095 @@
                      th_wake, tv_wake, lv_wake, lf_wake, cpn_wake, &
                      ep, sigp, clw, frac_s, qpreca, frac_a, qta, &                    !!jygprl
-                     m, ment, elij, delt, plcl, coef_clos, &
+                     m, ment, elij, delt, plcl, coef_clos_eff, &
                      mp, qp, up, vp, trap, wt, water, evap, fondue, ice, &
                      faci, b, sigd, &
@@ -1218 +1222 @@
                            Plim1, plim2, asupmax, supmax0, &
                            asupmaxmin, &
+                           coef_clos, coef_clos_eff, &
                            da, phi, mp, phi2, d1a, dam, sigij, &         ! RomP
                            qta, clw, elij, evap, ep, epmlmMm, eplaMm, &  ! RomP
@@ -1231 +1236 @@
                            Plim11, plim21, asupmax1, supmax01, &
                            asupmaxmin1, &
+                           coef_clos1, coef_clos_eff1, &
                            da1, phi1, mp1, phi21, d1a1, dam1, sigij1,  &       ! RomP
                            qta1, clw1, elij1, evap1, ep1, epmlmMm1, eplaMm1, & ! RomP