Changeset 2201 for LMDZ5/trunk/libf/phylmd

Timestamp:

Feb 10, 2015, 3:23:19 PM (10 years ago)

Author:

crio

Message:

Expérience du tube de dentifrice: Ajout d'un terme de grande-échelle dans
la fermeture du schéma de convection. La fraction de la convergence
grande-échelle de masse au LFC ajoutée au flux de masse à la base calculé
par la fermeture en ALP est controlée par coef_clos_ls lu dans physiq.def
(valeur comprise entre 0 et 1, 0 par défaut).
Tube of toothpaste experiment: Additional large-scale term in the closure
formulation of the deep convection scheme. The fraction of the large-scale
convergence of mass at LFC additioned to the cloud-base mass-flux given by
the ALP closure is controlled by coef_clos_ls read in physiq.def (value
between 0 and 1, 0 by default).

Location:

LMDZ5/trunk/libf/phylmd

Files:

• YOMCST2.h (modified) (2 diffs)
• concvl.F90 (modified) (2 diffs)
• conf_phys_m.F90 (modified) (4 diffs)
• cv3a_compress.F90 (modified) (4 diffs)
• cv3p1_closure.F90 (modified) (5 diffs)
• cva_driver.F90 (modified) (6 diffs)

LMDZ5/trunk/libf/phylmd/YOMCST2.h

@@ -3 +3 @@
       REAL  gammas, alphas, betas, Fmax, qqa1, qqa2, qqa3, scut
       REAL  Qcoef1max,Qcoef2max,Supcrit1,Supcrit2
+      REAL coef_clos_ls
 !
       COMMON/YOMCST2/gammas,    alphas, betas, Fmax, scut,              &
@@ -8 +9 @@
      &               Qcoef1max,Qcoef2max,                               &
      &               Supcrit1, Supcrit2,                                &
-     &               choice,iflag_mix
+     &               choice,iflag_mix,coef_clos_ls
 !$OMP THREADPRIVATE(/YOMCST2/)
 !    --------------------------------------------------------------------

LMDZ5/trunk/libf/phylmd/concvl.F90

@@ -29 +29 @@
   USE dimphy
   USE infotrac, ONLY: nbtr
+  USE phys_local_var_mod, ONLY: omega
   IMPLICIT NONE
 ! ======================================================================
@@ -398 +399 @@
                     t, q, qs, t_wake, q_wake, qs_wake, s_wake, u, v, tra, &
                     em_p, em_ph, &
-                    Ale, Alp, &
+                    Ale, Alp, omega, &
                     em_sig1feed, em_sig2feed, em_wght, &
                     iflag, d_t, d_q, d_u, d_v, d_tra, rain, kbas, ktop, &

LMDZ5/trunk/libf/phylmd/conf_phys_m.F90

@@ -132 +132 @@
     integer,SAVE :: iflag_coupl_omp,iflag_clos_omp,iflag_wake_omp
     integer,SAVE :: iflag_cvl_sigd_omp
+    REAL, SAVE :: coef_clos_ls_omp
     REAL, SAVE :: supcrit1_omp, supcrit2_omp
     INTEGER, SAVE :: iflag_mix_omp
@@ -1343 +1344 @@
     iflag_clos_omp = 1
     call getin('iflag_clos',iflag_clos_omp)
+    !
+    !Config Key  = coef_clos_ls
+    !Config Desc =  
+    !Config Def  = 0
+    !Config Help = 
+    !
+    coef_clos_ls_omp = 0.
+    call getin('coef_clos_ls',coef_clos_ls_omp)
+
     !
     !Config Key  = iflag_cvl_sigd
@@ -1946 +1956 @@
     iflag_clos = iflag_clos_omp
     iflag_wake = iflag_wake_omp
+    coef_clos_ls = coef_clos_ls_omp
     alp_offset = alp_offset_omp
     iflag_cvl_sigd = iflag_cvl_sigd_omp
@@ -2135 +2146 @@
     write(lunout,*)' iflag_thermals_closure = ', iflag_thermals_closure
     write(lunout,*)' iflag_clos = ', iflag_clos
+    write(lunout,*)' coef_clos_ls = ', coef_clos_ls
     write(lunout,*)' type_run = ',type_run 
     write(lunout,*)' ok_cosp = ',ok_cosp

LMDZ5/trunk/libf/phylmd/cv3a_compress.F90

@@ -4 +4 @@
     th1_wake, tra1, h1, lv1, lf1, cpn1, p1, ph1, tv1, tp1, tvp1, clw1, &
     h1_wake, lv1_wake, lf1_wake, cpn1_wake, tv1_wake, sig1, w01, ptop21, &
-    ale1, alp1, iflag, nk, icb, icbs, plcl, tnk, qnk, gznk, hnk, unk, vnk, &
+    ale1, alp1, omega1, iflag, nk, icb, icbs, plcl, tnk, qnk, gznk, hnk, unk, vnk, &
     wghti, pbase, buoybase, t, q, qs, t_wake, q_wake, qs_wake, s_wake, u, v, &
     gz, th, th_wake, tra, h, lv, lf, cpn, p, ph, tv, tp, tvp, clw, h_wake, &
-    lv_wake, lf_wake, cpn_wake, tv_wake, sig, w0, ptop2, ale, alp)
+    lv_wake, lf_wake, cpn_wake, tv_wake, sig, w0, ptop2, ale, alp, omega)
   ! **************************************************************
   ! *
@@ -40 +40 @@
   REAL sig1(len, nd), w01(len, nd), ptop21(len)
   REAL ale1(len), alp1(len)
+  REAL omega1(len,nd)
 
   ! outputs:
@@ -60 +61 @@
   REAL sig(len, nd), w0(len, nd), ptop2(len)
   REAL ale(len), alp(len)
+  REAL omega(len,nd)
 
   ! local variables:
@@ -102 +104 @@
         sig(nn, k) = sig1(i, k)
         w0(nn, k) = w01(i, k)
+        omega(nn, k) = omega1(i, k)
       END IF
     END DO

LMDZ5/trunk/libf/phylmd/cv3p1_closure.F90

@@ -3 +3 @@
 
 SUBROUTINE cv3p1_closure(nloc, ncum, nd, icb, inb, pbase, plcl, p, ph, tv, &
-    tvp, buoy, supmax, ok_inhib, ale, alp, sig, w0, ptop2, cape, cin, m, &
+    tvp, buoy, supmax, ok_inhib, ale, alp, omega,sig, w0, ptop2, cape, cin, m, &
     iflag, coef, plim1, plim2, asupmax, supmax0, asupmaxmin, cbmf, plfc, &
     wbeff)
@@ -37 +37 @@
   LOGICAL ok_inhib ! enable convection inhibition by dryness
   REAL ale(nloc), alp(nloc)
+  REAL omega(nloc,nd)
 
   ! input/output:
@@ -52 +53 @@
 
   ! local variables:
-  INTEGER il, i, j, k, icbmax, i0(nloc)
+  INTEGER il, i, j, k, icbmax, i0(nloc), klfc
   REAL deltap, fac, w, amu
   REAL rhodp
@@ -523 +524 @@
   END DO
 
+!CR:Compute k at plfc
+  DO k=1,nl
+     DO il=1,ncum
+        if ((plfc(il).lt.ph(il,k)).and.(plfc(il).ge.ph(il,k+1))) then
+           klfc=k
+        endif
+     ENDDO
+  ENDDO
+!RC
 
   DO il = 1, ncum
@@ -528 +538 @@
     ! c       cbmf1(il) = alp2(il)/(0.5*wb*wb-Cin(il))
     cbmf1(il) = alp2(il)/(2.*wbeff(il)*wbeff(il)-cin(il))
+!CR: Add large-scale component to the mass-flux
+!encore connu sous le nom "Experience du tube de dentifrice"
+    if (coef_clos_ls.gt.0.) then 
+       cbmf1(il) = cbmf1(il) - coef_clos_ls*min(0.,1./RG*omega(il,klfc))
+    endif
+!RC
     IF (cbmf1(il)==0 .AND. alp2(il)/=0.) THEN
       WRITE (lunout, *) 'cv3p1_closure cbmf1=0 and alp NE 0 il alp2 alp cin ' &

LMDZ5/trunk/libf/phylmd/cva_driver.F90

@@ -7 +7 @@
                       u1, v1, tra1, &
                       p1, ph1, &
-                      Ale1, Alp1, &
+                      Ale1, Alp1, omega1, &
                       sig1feed1, sig2feed1, wght1, &
                       iflag1, ft1, fq1, fu1, fv1, ftra1, &
@@ -178 +178 @@
   REAL Ale1(len)
   REAL Alp1(len)
+  REAL omega1(len,nd)
   REAL sig1feed1 ! pressure at lower bound of feeding layer
   REAL sig2feed1 ! pressure at upper bound of feeding layer
@@ -455 +456 @@
   REAL supmax(nloc, klev)
   REAL Ale(nloc), Alp(nloc), coef_clos(nloc)
+  REAL omega(nloc,klev)
   REAL sigd(nloc)
 ! real mp(nloc,klev), qp(nloc,klev), up(nloc,klev), vp(nloc,klev)
@@ -776 +778 @@
                          h1_wake, lv1_wake, lf1_wake, cpn1_wake, tv1_wake, &
                          sig1, w01, ptop21, &
-                         Ale1, Alp1, &
+                         Ale1, Alp1, omega1, &
                          iflag, nk, icb, icbs, &
                          plcl, tnk, qnk, gznk, hnk, unk, vnk, &
@@ -786 +788 @@
                          h_wake, lv_wake, lf_wake, cpn_wake, tv_wake, &
                          sig, w0, ptop2, &
-                         Ale, Alp)
+                         Ale, Alp, omega)
 
 ! print*,'tv ',tv
@@ -877 +879 @@
         CALL cv3p1_closure(nloc, ncum, nd, icb, inb, &         ! na->nd
                            pbase, plcl, p, ph, tv, tvp, buoy, &
-                           supmax, ok_inhib, Ale, Alp, &
+                           supmax, ok_inhib, Ale, Alp, omega, &
                            sig, w0, ptop2, cape, cin, m, iflag, coef_clos, &
                            Plim1, plim2, asupmax, supmax0, &