Changeset 5609 for LMDZ6/branches/contrails/libf/phylmd/dyn1d

Timestamp:

Apr 13, 2025, 7:10:19 PM (3 months ago)

Author:

aborella

Message:

• changed treatment of prognostic variables for prognostic clouds
• adapted sedimentation and autoconversion for prognostic cirrus clouds
• cloud mixing, ice sedimentation and ISSR diagnosis are now consistent with the water vapor PDF
• simplified assumptions for ice crystals deposition / sublimation
• first version of the coupling between prognostic cirrus clouds and deep convection
• added persistent contrail cirrus clouds in radiative diagnostics

Location:

LMDZ6/branches/contrails/libf/phylmd/dyn1d

Files:

• old_lmdz1d.f90 (modified) (3 diffs)
• scm.f90 (modified) (6 diffs)

LMDZ6/branches/contrails/libf/phylmd/dyn1d/old_lmdz1d.f90

@@ -10 +10 @@
       USE ioipsl, only: ju2ymds, ymds2ju, ioconf_calendar,getin
    USE phys_state_var_mod, ONLY : phys_state_var_init, phys_state_var_end, &
-       clwcon, detr_therm, &
+       cwcon, clwcon, detr_therm, &
        qsol, fevap, z0m, z0h, agesno, &
        du_gwd_rando, du_gwd_front, entr_therm, f0, fm_therm, &
@@ -23 +23 @@
        zstd, zthe, zval, ale_bl, ale_bl_trig, alp_bl, &
        ql_ancien, qs_ancien, qbs_ancien, &
-       cf_ancien, qvc_ancien, cfa_ancien, qta_ancien, &
+       cf_ancien, qvc_ancien, cfa_ancien, pcf_ancien, qva_ancien, qia_ancien, &
        prlw_ancien, prsw_ancien, prbsw_ancien, prw_ancien, &
        u10m,v10m,ale_wake,ale_bl_stat
@@ -874 +874 @@
           cf_ancien = 0.
           qvc_ancien = 0.
+          cwcon = 0.
         ENDIF
         IF ( ok_plane_contrail ) THEN
           cfa_ancien = 0.
-          qta_ancien = 0.
+          pcf_ancien = 0.
+          qva_ancien = 0.
+          qia_ancien = 0.
         ENDIF
 !jyg<

LMDZ6/branches/contrails/libf/phylmd/dyn1d/scm.f90

@@ -6 +6 @@
    USE ioipsl, only: ju2ymds, ymds2ju, ioconf_calendar,getin
    USE phys_state_var_mod, ONLY : phys_state_var_init, phys_state_var_end, &
-       clwcon, detr_therm, &
+       cwcon, clwcon, detr_therm, &
        qsol, fevap, z0m, z0h, agesno, &
        du_gwd_rando, du_gwd_front, entr_therm, f0, fm_therm, &
@@ -19 +19 @@
        zstd, zthe, zval, ale_bl, ale_bl_trig, alp_bl, &
        ql_ancien, qs_ancien, qbs_ancien, &
-       cf_ancien, qvc_ancien, cfa_ancien, qta_ancien, &
+       cf_ancien, qvc_ancien, cfa_ancien, pcf_ancien, qva_ancien, qia_ancien, &
        prlw_ancien, prsw_ancien, prbsw_ancien, prw_ancien, &
        u10m,v10m,ale_wake,ale_bl_stat, ratqs_inter_
@@ -616 +616 @@
           cf_ancien = 0.
           qvc_ancien = 0.
+          cwcon = 0.
         ENDIF
         IF ( ok_plane_contrail ) THEN
           cfa_ancien = 0.
-          qta_ancien = 0.
+          pcf_ancien = 0.
+          qva_ancien = 0.
+          qia_ancien = 0.
         ENDIF
         rain_fall=0.
@@ -820 +823 @@
       ENDIF
       
+      ! calculation of potential temperature for the advection
+      teta=temp*(pzero/play)**rkappa
+
       ! vertical tendencies computed as d X / d t = -W  d X / d z
 
@@ -829 +835 @@
         d_v_vert_adv(l)=-w_adv(l)*(v(l+1)-v(l-1))/(z_adv(l+1)-z_adv(l-1))
         ! d theta / dt = -W d theta / d z, transformed into d temp / d t dividing by (pzero/play(l))**rkappa
-        d_t_vert_adv(l)=-w_adv(l)*(temp(l+1)-temp(l-1))/(z_adv(l+1)-z_adv(l-1)) 
+        d_t_vert_adv(l)=-w_adv(l)*(teta(l+1)-teta(l-1))/(z_adv(l+1)-z_adv(l-1)) / (pzero/play(l))**rkappa
         d_q_vert_adv(l,:)=-w_adv(l)*(q(l+1,:)-q(l-1,:))/(z_adv(l+1)-z_adv(l-1))
       enddo
@@ -842 +848 @@
           d_u_vert_adv(l)=-w_adv(l)*(u(l)-u(l-1))/(z_adv(l)-z_adv(l-1))
           d_v_vert_adv(l)=-w_adv(l)*(v(l)-v(l-1))/(z_adv(l)-z_adv(l-1))
-          d_t_vert_adv(l)=-w_adv(l)*(temp(l)-temp(l-1))/(z_adv(l)-z_adv(l-1))
+          ! d theta / dt = -W d theta / d z, transformed into d temp / d t dividing by (pzero/play(l))**rkappa
+          d_t_vert_adv(l)=-w_adv(l)*(teta(l)-teta(l-1))/(z_adv(l)-z_adv(l-1)) / (pzero/play(l))**rkappa
           d_q_vert_adv(l,:)=-w_adv(l)*(q(l,:)-q(l-1,:))/(z_adv(l)-z_adv(l-1))
         ELSE
           d_u_vert_adv(l)=-w_adv(l)*(u(l+1)-u(l))/(z_adv(l+1)-z_adv(l))
           d_v_vert_adv(l)=-w_adv(l)*(v(l+1)-v(l))/(z_adv(l+1)-z_adv(l))
-          d_t_vert_adv(l)=-w_adv(l)*(temp(l+1)-temp(l))/(z_adv(l+1)-z_adv(l))
+          ! d theta / dt = -W d theta / d z, transformed into d temp / d t dividing by (pzero/play(l))**rkappa
+          d_t_vert_adv(l)=-w_adv(l)*(teta(l+1)-teta(l))/(z_adv(l+1)-z_adv(l)) / (pzero/play(l))**rkappa
           d_q_vert_adv(l,:)=-w_adv(l)*(q(l+1,:)-q(l,:))/(z_adv(l+1)-z_adv(l))
         ENDIF