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

Timestamp:

Aug 4, 2025, 3:03:07 PM (10 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_ctrlout_mod.F90 (modified) (2 diffs)

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

@@ -2714 +2714 @@
   TYPE(ctrl_out), SAVE :: o_section_contrails = ctrl_out((/ 11, 11, 11, 11, 11, 11, 11, 11, 11, 11/), &
     'section_cont', 'Cross section of newly formed contrails', 'm2', (/ ('', i=1, 10) /))
-  TYPE(ctrl_out), SAVE :: o_cldfra_nocont = ctrl_out((/ 11, 11, 11, 11, 11, 11, 11, 11, 11, 11/), &
-    'cldfra_nocont', 'Cloud fraction w/o contrails', '-', (/ ('', i=1, 10) /))
+  TYPE(ctrl_out), SAVE :: o_nice_ygcont = ctrl_out((/ 11, 11, 11, 11, 11, 11, 11, 11, 11, 11/), &
+    'nice_ygcont', 'Ice particle number concentration in young contrails', '#/cm3', (/ ('', i=1, 10) /))
+  TYPE(ctrl_out), SAVE :: o_iwc_ygcont = ctrl_out((/ 11, 11, 11, 11, 11, 11, 11, 11, 11, 11/), &
+    'iwc_ygcont', 'Ice water content in young contrails', 'g/m3', (/ ('', i=1, 10) /))
+  TYPE(ctrl_out), SAVE :: o_rvol_ygcont = ctrl_out((/ 11, 11, 11, 11, 11, 11, 11, 11, 11, 11/), &
+    'rvol_ygcont', 'Ice crystals volumic radius in young contrails', 'microns', (/ ('', i=1, 10) /))
+  TYPE(ctrl_out), SAVE :: o_tau_ygcont = ctrl_out((/ 11, 11, 11, 11, 11, 11, 11, 11, 11, 11/), &
+    'tau_ygcont', 'Young contrails optical depth', '-', (/ ('', i=1, 10) /))
+  TYPE(ctrl_out), SAVE :: o_nice_cont = ctrl_out((/ 11, 11, 11, 11, 11, 11, 11, 11, 11, 11/), &
+    'nice_cont', 'Ice particle number concentration in contrails', '#/cm3', (/ ('', i=1, 10) /))
+  TYPE(ctrl_out), SAVE :: o_iwc_cont = ctrl_out((/ 11, 11, 11, 11, 11, 11, 11, 11, 11, 11/), &
+    'iwc_cont', 'Ice water content in contrails', 'g/m3', (/ ('', i=1, 10) /))
+  TYPE(ctrl_out), SAVE :: o_rvol_cont = ctrl_out((/ 11, 11, 11, 11, 11, 11, 11, 11, 11, 11/), &
+    'rvol_cont', 'Ice crystals volumic radius in contrails', 'microns', (/ ('', i=1, 10) /))
+  TYPE(ctrl_out), SAVE :: o_tau_cont = ctrl_out((/ 11, 11, 11, 11, 11, 11, 11, 11, 11, 11/), &
+    'tau_cont', 'Contrails optical depth', '-', (/ ('', i=1, 10) /))
+  TYPE(ctrl_out), SAVE :: o_cldfra_cont = ctrl_out((/ 11, 11, 11, 11, 11, 11, 11, 11, 11, 11/), &
+    'cldfra_cont', 'Cloud fraction of contrails', '-', (/ ('', i=1, 10) /))
   TYPE(ctrl_out), SAVE :: o_cldtau_nocont = ctrl_out((/ 11, 11, 11, 11, 11, 11, 11, 11, 11, 11/), &
     'cldtau_nocont', 'Cloud optical thickness w/o contrails', '1', (/ ('', i=1, 10) /))
@@ -2728 +2744 @@
   TYPE(ctrl_out), SAVE :: o_contemi = ctrl_out((/ 11, 11, 11, 11, 11, 11, 11, 11, 11, 11/), &
     'contemi', 'Contrails optical emissivity', '1', (/ ('', i=1, 10) /))
-  TYPE(ctrl_out), SAVE :: o_iwp_nocont = ctrl_out((/ 11, 11, 11, 11, 11, 11, 11, 11, 11, 11/), &
-    'iwp_nocont', 'Cloud ice water path w/o contrails', 'kg/m2', (/ ('', i=1, 10) /))
-  TYPE(ctrl_out), SAVE :: o_iwc_nocont = ctrl_out((/ 11, 11, 11, 11, 11, 11, 11, 11, 11, 11/), &
-    'iwc_nocont', 'Cloud ice water content seen by radiation w/o contrails', 'kg/kg', (/ ('', i=1, 10) /))
-  TYPE(ctrl_out), SAVE :: o_ref_ice_nocont = ctrl_out((/ 4, 10, 10, 10, 10, 10, 11, 11, 11, 11/), &
-    'ref_ice_nocont', 'Effective radius of ice crystals w/o contrails', 'microns', (/ ('', i=1, 10) /))
+  TYPE(ctrl_out), SAVE :: o_iwp_cont = ctrl_out((/ 11, 11, 11, 11, 11, 11, 11, 11, 11, 11/), &
+    'iwp_cont', 'Cloud ice water path of contrails', 'kg/m2', (/ ('', i=1, 10) /))
+  TYPE(ctrl_out), SAVE :: o_iwcon_cont = ctrl_out((/ 11, 11, 11, 11, 11, 11, 11, 11, 11, 11/), &
+    'iwcon_cont', 'Cloud ice water content seen by radiation of contrails', 'kg/kg', (/ ('', i=1, 10) /))
+  TYPE(ctrl_out), SAVE :: o_ref_ice_cont = ctrl_out((/ 4, 10, 10, 10, 10, 10, 11, 11, 11, 11/), &
+    'ref_ice_cont', 'Effective radius of ice crystals of contrails', 'microns', (/ ('', i=1, 10) /))
   TYPE(ctrl_out), SAVE :: o_tops_nocont = ctrl_out((/ 11, 11, 11, 11, 11, 11, 11, 11, 11, 11/), &
     'tops_nocont', 'Solar rad. at TOA w/o contrails', 'W/m2', (/ ('', i=1, 10) /))