Changeset 5796 for LMDZ6/branches/contrails/libf/phylmd/phys_output_write_mod.F90

Timestamp:

Aug 4, 2025, 3:03:07 PM (11 days ago)

Author:

aborella

Message:

Additional diags for contrails + simplified coupling between deep conv and cirrus clouds + small modifsin RRTM for RF of contrails alone

File:

• LMDZ6/branches/contrails/libf/phylmd/phys_output_write_mod.F90 (modified) (7 diffs)

LMDZ6/branches/contrails/libf/phylmd/phys_output_write_mod.F90

@@ -245 +245 @@
          o_Tcritcont, o_qcritcont, o_potcontfraP, o_potcontfraNP, &
          o_AEI_contrails, o_AEI_surv_contrails, o_fsurv_contrails, o_section_contrails, &
+         o_nice_ygcont, o_iwc_ygcont, o_rvol_ygcont, o_tau_ygcont, &
+         o_nice_cont, o_iwc_cont, o_rvol_cont, o_tau_cont, &
          o_flight_dist, o_flight_fuel, o_qice_cont, &
          o_dcfcini, o_dqicini, o_dqtcini, o_dnicini, &
@@ -251 +253 @@
          o_dcfcsed, o_dqicsed, o_dqtcsed, o_dnicsed, &
          o_dcfcauto, o_dqicauto, o_dqtcauto, o_dnicauto, &
-         o_cldfra_nocont, o_cldtau_nocont, o_cldemi_nocont, o_cldh_nocont, &
-         o_contcov, o_conttau, o_contemi, o_iwp_nocont, o_iwc_nocont, o_ref_ice_nocont, &
+         o_cldfra_cont, o_cldtau_nocont, o_cldemi_nocont, o_cldh_nocont, &
+         o_contcov, o_conttau, o_contemi, o_iwp_cont, o_iwcon_cont, o_ref_ice_cont, &
          o_tops_nocont, o_topl_nocont, o_sols_nocont, o_soll_nocont, o_nettop_nocont, &
 !--interactive CO2
@@ -411 +413 @@
          Tcritcont, qcritcont, potcontfraP, potcontfraNP, &
          AEI_contrails, AEI_surv_contrails, fsurv_contrails, section_contrails, &
+         nice_ygcont, iwc_ygcont, rvol_ygcont, tau_ygcont, &
+         nice_cont, iwc_cont, rvol_cont, tau_cont, &
          flight_dist, flight_fuel, qice_cont, &
          dcfc_ini, dqic_ini, dqtc_ini, dnic_ini, &
@@ -417 +421 @@
          dcfc_sed, dqic_sed, dqtc_sed, dnic_sed, &
          dcfc_auto, dqic_auto, dqtc_auto, dnic_auto, &
-         cldfra_nocont, cldtau_nocont, cldemi_nocont, cldh_nocont, &
-         contcov, conttau, contemi, fiwp_nocont, fiwc_nocont, ref_ice_nocont, &
+         cldfra_cont, cldtau_nocont, cldemi_nocont, cldh_nocont, &
+         contcov, conttau, contemi, fiwp_cont, fiwc_cont, ref_ice_cont, &
          topsw_nocont, toplw_nocont, solsw_nocont, sollw_nocont, &
          alp_bl_det, alp_bl_fluct_m, alp_bl_conv, &
@@ -2463 +2467 @@
          CALL histwrite_phy(o_fsurv_contrails, fsurv_contrails)
          CALL histwrite_phy(o_section_contrails, section_contrails)
+         CALL histwrite_phy(o_nice_ygcont, nice_ygcont)
+         CALL histwrite_phy(o_iwc_ygcont, iwc_ygcont)
+         CALL histwrite_phy(o_rvol_ygcont, rvol_ygcont)
+         CALL histwrite_phy(o_tau_ygcont, tau_ygcont)
+         CALL histwrite_phy(o_nice_cont, nice_cont)
+         CALL histwrite_phy(o_iwc_cont, iwc_cont)
+         CALL histwrite_phy(o_rvol_cont, rvol_cont)
+         CALL histwrite_phy(o_tau_cont, tau_cont)
          CALL histwrite_phy(o_qice_cont, qice_cont)
          CALL histwrite_phy(o_dcfcini, dcfc_ini)
@@ -2485 +2497 @@
          CALL histwrite_phy(o_dqtcauto, dqtc_auto)
          CALL histwrite_phy(o_dnicauto, dnic_auto)
-         CALL histwrite_phy(o_cldfra_nocont, cldfra_nocont)
+         CALL histwrite_phy(o_cldfra_cont, cldfra_cont)
          CALL histwrite_phy(o_cldtau_nocont, cldtau_nocont)
          CALL histwrite_phy(o_cldemi_nocont, cldemi_nocont)
@@ -2492 +2504 @@
          CALL histwrite_phy(o_conttau, conttau)
          CALL histwrite_phy(o_contemi, contemi)
-         CALL histwrite_phy(o_iwp_nocont, fiwp_nocont)
-         CALL histwrite_phy(o_iwc_nocont, fiwc_nocont)
-         CALL histwrite_phy(o_ref_ice_nocont, ref_ice_nocont)
+         CALL histwrite_phy(o_iwp_cont, fiwp_cont)
+         CALL histwrite_phy(o_iwcon_cont, fiwc_cont)
+         CALL histwrite_phy(o_ref_ice_cont, ref_ice_cont)
          IF (ok_rad_contrail) THEN
            IF (vars_defined) zx_tmp_fi2d = topsw_nocont * swradcorr