Changeset 6112 for LMDZ6/trunk/libf

Timestamp:

Mar 17, 2026, 6:38:15 PM (2 weeks ago)

Author:

asima

Message:

For Simple Plume aerosols (flag_aerosol=7) : possibility to output dNovrN factor (multiplying the preind. could droplet number cdnc_pi to mimic the Twomey aerosol indirect effect)

Location:

LMDZ6/trunk/libf/phylmd

Files:

• phys_output_ctrlout_mod.F90 (modified) (1 diff)
• phys_output_write_mod.F90 (modified) (3 diffs)

LMDZ6/trunk/libf/phylmd/phys_output_ctrlout_mod.F90

@@ -1758 +1758 @@
   TYPE(ctrl_out), SAVE :: o_od550SPaer = ctrl_out((/ 2, 6, 10, 10, 10, 10, 11, 11, 11, 11/), &
     'od550SPaer', 'Simple Plume aerosol optical depth', '-', (/ ('', i=1, 10) /))
+  TYPE(ctrl_out), SAVE :: o_dNovrN = ctrl_out((/ 6, 6, 10, 10, 10, 10, 11, 11, 11, 11/), &
+    'od550SPaer', 'Simple Plume multiplication factor of preind cloud droplet number', '-', (/ ('', i=1, 10) /))
   TYPE(ctrl_out), SAVE :: o_sconcso4 = ctrl_out((/ 2, 6, 10, 10, 10, 10, 11, 11, 11, 11/), &
     'sconcso4', 'Surface Concentration of Sulfate ', 'kg/m3', (/ ('', i=1, 10) /))

LMDZ6/trunk/libf/phylmd/phys_output_write_mod.F90

@@ -124 +124 @@
          o_solldown, o_dtsvdfo, o_dtsvdft, &
          o_dtsvdfg, o_dtsvdfi, o_z0m, o_z0h,  o_od443aer, o_od550aer, &
-         o_dryod550aer, o_od865aer, o_abs550aer, o_od550lt1aer, o_od550SPaer,  &
+         o_dryod550aer, o_od865aer, o_abs550aer, o_od550lt1aer, o_od550SPaer, o_dNovrN,  &
          o_sconcso4, o_sconcno3, o_sconcoa, o_sconcbc, &
          o_sconcss, o_sconcdust, o_concso4, o_concno3, &
@@ -402 +402 @@
          pmflxr, pmflxs, prfl, psfl,bsfl, re, fl, rh2m, &
          qsat2m, tpote, tpot, d_ts, od443aer, od550aer, dryod550aer, &
-         od865aer, abs550aer, od550lt1aer, od550SPaer, sconcso4, sconcno3, &
+         od865aer, abs550aer, od550lt1aer, od550SPaer, dNovrN,  sconcso4, sconcno3, &
          sconcoa, sconcbc, sconcss, sconcdust, concso4, concno3, &
          concoa, concbc, concss, concdust, loadso4, &
@@ -2095 +2095 @@
              CALL histwrite_phy(o_dryod550aer, dryod550aer)
              IF (flag_aerosol.EQ.7) THEN
-                ! (ASima) od550 for Simple Plume aerosols
+                ! (ASima) od550 for Simple Plume aerosols, and dNovrN
                 CALL histwrite_phy(o_od550SPaer, od550SPaer)
+                CALL histwrite_phy(o_dNovrN, dNovrN)
              ENDIF
              DO naero = 1, naero_tot-1