Changeset 4744 for LMDZ6/trunk/libf/phylmd/dyn1d

Timestamp:

Nov 2, 2023, 10:09:59 AM (8 months ago)

Author:

jyg

Message:

Implementation of a two radii model for wake population dynamics.
It is activated with : iflag_wk_pop_dyn=3

Location:

LMDZ6/trunk/libf/phylmd/dyn1d

Files:

• old_lmdz1d.F90 (modified) (3 diffs)
• scm.F90 (modified) (3 diffs)

LMDZ6/trunk/libf/phylmd/dyn1d/old_lmdz1d.F90

@@ -15 +15 @@
        solsw, solswfdiff, t_ancien, q_ancien, u_ancien, v_ancien, rneb_ancien, &
        wake_delta_pbl_TKE, delta_tsurf, wake_fip, wake_pe, &
-       wake_deltaq, wake_deltat, wake_s, wake_dens, &
+       wake_deltaq, wake_deltat, wake_s, awake_s, wake_dens, &
        awake_dens, cv_gen, wake_cstar, &
        zgam, zmax0, zmea, zpic, zsig, &
@@ -905 +905 @@
         wake_pe = 0.
         wake_s = 0.
+        awake_s = 0.
         wake_dens = 0.
         awake_dens = 0.
@@ -938 +939 @@
 ! t_ancien,q_ancien,,frugs(:,is_oce),clwcon(:,1),rnebcon(:,1),ratqs(:,1)
 ! run_off_lic_0,pbl_tke(:,1:klev,nsrf), zmax0,f0,sig1,w01
-! wake_deltat,wake_deltaq,wake_s,wake_dens,awake_dens,cv_gen,wake_cstar,
+! wake_deltat,wake_deltaq,wake_s,awake_s,wake_dens,awake_dens,cv_gen,wake_cstar,
 ! wake_fip,wake_delta_pbl_tke(:,1:klev,nsrf)
 !

LMDZ6/trunk/libf/phylmd/dyn1d/scm.F90

@@ -11 +11 @@
        solsw, solswfdiff, t_ancien, q_ancien, u_ancien, v_ancien, rneb_ancien, &
        wake_delta_pbl_TKE, delta_tsurf, wake_fip, wake_pe, &
-       wake_deltaq, wake_deltat, wake_s, wake_dens, &
+       wake_deltaq, wake_deltat, wake_s, awake_s, wake_dens, &
        awake_dens, cv_gen, wake_cstar, &
        zgam, zmax0, zmea, zpic, zsig, &
@@ -663 +663 @@
         wake_pe = 0.
         wake_s = 0.
+        awake_s = 0.
         wake_dens = 0.
         awake_dens = 0.
@@ -696 +697 @@
 ! t_ancien,q_ancien,,frugs(:,is_oce),clwcon(:,1),rnebcon(:,1),ratqs(:,1)
 ! run_off_lic_0,pbl_tke(:,1:klev,nsrf), zmax0,f0,sig1,w01
-! wake_deltat,wake_deltaq,wake_s,wake_dens,awake_dens,cv_gen,wake_cstar,
+! wake_deltat,wake_deltaq,wake_s,awake_s,wake_dens,awake_dens,cv_gen,wake_cstar,
 ! wake_fip,wake_delta_pbl_tke(:,1:klev,nsrf)
 !