Changeset 4803 for LMDZ6/trunk/libf/phylmdiso

Timestamp:

Feb 5, 2024, 10:16:07 PM (10 months ago)

Author:

evignon

Message:

implementation sous flag des premiers changements
concernant le traitement des precipitations grande echelle
dans le cadre de l'atelier nuages
Audran, Lea, Niels, Gwendal et Etienne

Location:

LMDZ6/trunk/libf/phylmdiso

Files:

• lmdz_lscp_poprecip.F90 (added)
• phys_local_var_mod.F90 (modified) (3 diffs)
• phys_output_ctrlout_mod.F90 (modified) (1 diff)
• physiq_mod.F90 (modified) (3 diffs)

LMDZ6/trunk/libf/phylmdiso/phys_local_var_mod.F90

@@ -639 +639 @@
       !$OMP THREADPRIVATE(fcontrP)
 
+!--POPRECIP variables
+      REAL, SAVE, ALLOCATABLE :: dqreva(:,:)
+      !$OMP THREADPRIVATE(dqreva)
+      REAL, SAVE, ALLOCATABLE :: dqrauto(:,:)
+      !$OMP THREADPRIVATE(dqrauto)
+      REAL, SAVE, ALLOCATABLE :: dqrcol(:,:)
+      !$OMP THREADPRIVATE(dqrcol)
+      REAL, SAVE, ALLOCATABLE :: dqrmelt(:,:)
+      !$OMP THREADPRIVATE(dqrmelt)
+       REAL, SAVE, ALLOCATABLE :: dqrfreez(:,:)
+      !$OMP THREADPRIVATE(dqrfreez)
+      REAL, SAVE, ALLOCATABLE :: dqssub(:,:)
+      !$OMP THREADPRIVATE(dqssub)
+      REAL, SAVE, ALLOCATABLE :: dqsauto(:,:)
+      !$OMP THREADPRIVATE(dqsauto)
+      REAL, SAVE, ALLOCATABLE :: dqsagg(:,:)
+      !$OMP THREADPRIVATE(dqsagg)
+      REAL, SAVE, ALLOCATABLE :: dqsrim(:,:)
+      !$OMP THREADPRIVATE(dqsrim)
+      REAL, SAVE, ALLOCATABLE :: dqsmelt(:,:)
+      !$OMP THREADPRIVATE(dqsmelt)
+      REAL, SAVE, ALLOCATABLE :: dqsfreez(:,:)
+      !$OMP THREADPRIVATE(dqsfreez)
+
+
+
+
+
 #ifdef CPP_StratAer
 !
@@ -1085 +1113 @@
       ALLOCATE(zqsatl(klon,klev), zqsats(klon,klev))
       ALLOCATE(Tcontr(klon,klev), qcontr(klon,klev), qcontr2(klon,klev), fcontrN(klon,klev), fcontrP(klon,klev))
+
+!--POPRECIP variables
+      ALLOCATE(dqreva(klon,klev),dqssub(klon,klev))
+      ALLOCATE(dqrauto(klon,klev), dqrcol(klon,klev), dqrmelt(klon,klev), dqrfreez(klon,klev))
+      ALLOCATE(dqsauto(klon,klev), dqsagg(klon,klev), dqsrim(klon,klev), dqsmelt(klon,klev), dqsfreez(klon,klev))
 
 #ifdef CPP_StratAer
@@ -1444 +1477 @@
       DEALLOCATE(Tcontr, qcontr, qcontr2, fcontrN, fcontrP)
 
+!--POPRECIP variables
+      DEALLOCATE(dqreva,dqssub)
+      DEALLOCATE(dqrauto,dqrcol,dqrmelt,dqrfreez)
+      DEALLOCATE(dqsauto,dqsagg,dqsrim,dqsmelt,dqsfreez)
+
 #ifdef CPP_StratAer
 ! variables for strat. aerosol CK

LMDZ6/trunk/libf/phylmdiso/phys_output_ctrlout_mod.F90

@@ -1477 +1477 @@
   TYPE(ctrl_out), SAVE :: o_pfracld = ctrl_out((/ 11, 11, 11, 11, 11, 11, 11, 11, 11, 11/), &
     'pfracld', 'LS precipitation fraction cloudy part', '-', (/ ('', i=1, 10) /))
+  TYPE(ctrl_out), SAVE :: o_dqreva = ctrl_out((/ 11, 11, 11, 11, 11, 11, 11, 11, 11, 11/), &
+    'dqreva', 'LS rain tendency due to evaporation', 'kg/m2/s', (/ ('', i=1, 10) /))
+   TYPE(ctrl_out), SAVE :: o_dqrauto = ctrl_out((/ 11, 11, 11, 11, 11, 11, 11, 11, 11, 11/), &
+    'dqrauto', 'LS rain tendency due to autoconversion', 'kg/m2/s', (/ ('', i=1, 10) /))
+  TYPE(ctrl_out), SAVE :: o_dqrcol = ctrl_out((/ 11, 11, 11, 11, 11, 11, 11, 11, 11, 11/), &
+    'dqrcol', 'LS rain tendency due to collection', 'kg/m2/s', (/ ('', i=1, 10) /))
+  TYPE(ctrl_out), SAVE :: o_dqrmelt = ctrl_out((/ 11, 11, 11, 11, 11, 11, 11, 11, 11, 11/), &
+    'dqrmelt', 'LS rain tendency due to melting', 'kg/m2/s', (/ ('', i=1, 10) /))
+  TYPE(ctrl_out), SAVE :: o_dqrfreez = ctrl_out((/ 11, 11, 11, 11, 11, 11, 11, 11, 11, 11/), &
+    'dqrfreez', 'LS rain tendency due to freezing', 'kg/m2/s', (/ ('', i=1, 10) /))
+  TYPE(ctrl_out), SAVE :: o_dqssub = ctrl_out((/ 11, 11, 11, 11, 11, 11, 11, 11, 11, 11/), &
+    'dqssub', 'LS snow tendency due to sublimation', 'kg/m2/s', (/ ('', i=1, 10) /))
+  TYPE(ctrl_out), SAVE :: o_dqsauto = ctrl_out((/ 11, 11, 11, 11, 11, 11, 11, 11, 11, 11/), &
+    'dqsauto', 'LS snow tendency due to autoconversion', 'kg/m2/s', (/ ('', i=1, 10) /))
+  TYPE(ctrl_out), SAVE :: o_dqsagg = ctrl_out((/ 11, 11, 11, 11, 11, 11, 11, 11, 11, 11/), &
+    'dqsagg', 'LS snow tendency due to aggragation', 'kg/m2/s', (/ ('', i=1, 10) /))
+  TYPE(ctrl_out), SAVE :: o_dqsrim = ctrl_out((/ 11, 11, 11, 11, 11, 11, 11, 11, 11, 11/), &
+    'dqsrim', 'LS snow tendency due to riming', 'kg/m2/s', (/ ('', i=1, 10) /))
+  TYPE(ctrl_out), SAVE :: o_dqsmelt = ctrl_out((/ 11, 11, 11, 11, 11, 11, 11, 11, 11, 11/), &
+    'dqsmelt', 'LS snow tendency due to melting', 'kg/m2/s', (/ ('', i=1, 10) /))
+  TYPE(ctrl_out), SAVE :: o_dqsfreez = ctrl_out((/ 11, 11, 11, 11, 11, 11, 11, 11, 11, 11/), &
+    'dqsfreez', 'LS snow tendency due to freezing', 'kg/m2/s', (/ ('', i=1, 10) /))
   TYPE(ctrl_out), SAVE :: o_rhum = ctrl_out((/ 2, 5, 10, 10, 10, 10, 11, 11, 11, 11/), &
     'rhum', 'Relative humidity', '-', (/ ('', i=1, 10) /))

LMDZ6/trunk/libf/phylmdiso/physiq_mod.F90

@@ -229 +229 @@
        east_gwstress,west_gwstress, &
        d_q_ch4, &
+       ! proprecip
+       dqreva, dqssub, &
+       dqrauto,dqrcol,dqrmelt,dqrfreez, &
+       dqsauto,dqsagg,dqsrim,dqsmelt,dqsfreez, &
        !  Special RRTM
        ZLWFT0_i,ZSWFT0_i,ZFLDN0,  &
@@ -4876 +4880 @@
     CALL lscp(klon,klev,phys_tstep,missing_val,paprs,pplay, &
          t_seri, q_seri,ptconv,ratqs, &
-         d_t_lsc, d_q_lsc, d_ql_lsc, d_qi_lsc, rneb, rneblsvol, rneb_seri, &
+         d_t_lsc, d_q_lsc, d_ql_lsc, d_qi_lsc, rneb, rneblsvol, rneb_seri, & 
          pfraclr,pfracld, &
          radocond, picefra, rain_lsc, snow_lsc, &
@@ -4885 +4889 @@
          qclr, qcld, qss, qvc, rnebclr, rnebss, gamma_ss, &
          Tcontr, qcontr, qcontr2, fcontrN, fcontrP , &
-         cloudth_sth,cloudth_senv,cloudth_sigmath,cloudth_sigmaenv)
+         cloudth_sth,cloudth_senv,cloudth_sigmath,cloudth_sigmaenv, &
+         dqreva,dqssub,dqrauto,dqrcol,dqrmelt,dqrfreez,dqsauto,dqsagg,dqsrim,dqsmelt,dqsfreez)