Merge branch 'mingxuanwupnnl/atm/MAM5_tuning_forv3HR' (PR #7795)

Add tunable aerosol diameter thresholds for renaming to strat coarse mode in MAM5

Add tunable aerosol diameter thresholds for renaming to strat coarse mode in MAM5

When using CMIP7 volcanic SO2 emission, the model will produce larger stratospheric
AOD and sulfate burden due to increased emission amount as well as changing injection
height. The stratospheric sulfate burden during Pinatubo eruption when using CMIP7
volcanic SO2 emission is larger than HIRS observations.

For v3.HR simulations, we identified two tunable aerosol mode diameter thresholds
for renaming accumulation mode (mode 1) sulfate to stratospheric coarse mode (mode 5).
By changing them from 0.4 um to 1 um, changing the dgnumhi in mode 1 physprop file,
and changing dgnumlo in mode 5 physprop file, the stratospheric AOD and sulfate
burden during Pinatubo are reduced, the RESTOM and SW flux at TOA are also reduced
to the level when using default CMIP6 volcanic emission.

In this PR, we only introduce two namelist variables for the two tunable thresholds
and keep values the same as v3.LR simulations.

[BFB]
[NML] Two new namelist variables introduced affecting all tests using EAMv3

* origin/mingxuanwupnnl/atm/MAM5_tuning_forv3HR:
  Fix add_defaults for the new parameters only when defaults are available
  revise and add description for the two namelist variable according to review comments
  revise the description for the two threshold diameter
  Add missing args for phys_getopts
  add tunable aerosol mode diameter threholds for renaming to strat coarse mode in MAM5

Author: wlin7

components/eam/bld/build-namelist

@@ -3736,6 +3736,11 @@ add_default($nl, 'is_output_interactive_volc');
 add_default($nl, 'mam_amicphys_optaa'); #for enabling amicphys code in eam
 add_default($nl, 'n_so4_monolayers_pcage'); #number of so4(+nh4) monolayers needed to "age" a carbon particle
 
+#cutoff diameter for larger accum mode particles renamed to the strat coarse mode 
+add_default($nl, 'dp_cut_accum_rename') if (get_default_value('dp_cut_accum_rename')); 
+#threshold diameter for all accum mode particles transferred to the strat coarse mode
+add_default($nl, 'dp_xferall_thresh_accum_rename') if (get_default_value('dp_xferall_thresh_accum_rename')); 
+
 # Microphysics scheme
 add_default($nl, 'use_subcol_microp');
 add_default($nl, 'microp_scheme');

components/eam/bld/namelist_files/namelist_defaults_eam.xml

@@ -1917,6 +1917,8 @@ with se_tstep, dt_remap_factor, dt_tracer_factor set to -1
 <p3_wbf_coeff phys="default" clubb_sgs="1" microphys="p3"> 1.0      </p3_wbf_coeff>
 <p3_nc_autocon_expon  phys="default" clubb_sgs="1" microphys="p3">  -1.10D0   </p3_nc_autocon_expon>
 <n_so4_monolayers_pcage  phys="default" clubb_sgs="1" microphys="p3">  3.0D0   </n_so4_monolayers_pcage>
+<dp_cut_accum_rename  phys="default">  4.0D-7   </dp_cut_accum_rename>
+<dp_xferall_thresh_accum_rename phys="default">  4.0D-7  </dp_xferall_thresh_accum_rename>
 <clubb_C1 phys="default"> 2.4        </clubb_C1>
 <clubb_C11 phys="default"> 0.70       </clubb_C11>
 <clubb_C11b phys="default"> 0.20       </clubb_C11b>

components/eam/bld/namelist_files/namelist_definition.xml

@@ -4201,6 +4201,18 @@ number of so4(+nh4) monolayers needed to "age" a carbon particle
 Default: Set by namelist
 </entry>
 
+<entry id="dp_cut_accum_rename" type="real"  category="scavenging"
+       group="phys_ctl_nl" valid_values="" >
+cutoff diameter for larger accum mode particles renamed to the strat coarse mode
+Default: Set by namelist
+</entry>
+
+<entry id="dp_xferall_thresh_accum_rename" type="real"  category="scavenging"
+       group="phys_ctl_nl" valid_values="" >
+threshold diameter for all accum mode particles transferred to the strat coarse mode
+Default: Set by namelist
+</entry>
+
 
 <!-- Diagnostics -->
 

components/eam/src/chemistry/modal_aero/modal_aero_amicphys.F90

@@ -52,6 +52,14 @@ module modal_aero_amicphys
 ! aging criterion is approximate so do not try to distinguish
 !    sulfuric acid, bisulfate, ammonium sulfate
 
+  real(r8), public, protected :: dp_cut_accum_rename = huge(1.0_r8)
+! cutoff diameter for larger accum mode particles renamed to the strat coarse mode
+
+  real(r8), public, protected :: dp_xferall_thresh_accum_rename = huge(1.0_r8)
+! threshold diameter for all accum mode particles transferred to the strat coarse mode
+
+  real(r8), public, protected :: micro_mg_accre_enhan_fac = huge(1.0_r8)
+
 #if ( defined( CAMBOX_ACTIVATE_THIS ) )
   integer, public :: cldy_rh_sameas_clear = 0
 ! this is only used for some specific box model tests
@@ -4550,11 +4558,11 @@ subroutine mam_rename_1subarea_strat(                               &
          shrinkaa(n) = .false. 
          shrinkbb = .false. ! set initial value shrink switch: if calculate shrink fraction
          if ( mtoo == ncrsf ) then ! stratosphere renaming
-              dp_cut(mfrm) = 4.0e-7_r8
+              dp_cut(mfrm) = dp_cut_accum_rename
               lndp_cut(mfrm) = log( dp_cut(mfrm) )
               dp_belowcut(mfrm) = 0.99*dp_cut(mfrm)
               dp_xfernone_thresh(mfrm) = 1.0e-7_r8
-              dp_xferall_thresh(mfrm)    = 4.0e-7_r8
+              dp_xferall_thresh(mfrm)    = dp_xferall_thresh_accum_rename
          end if
          if (1>2) then ! turn off the shrink path for now
          if ((mfrm == ncrsf .and. ncrsf > 0)) then
@@ -5994,7 +6002,8 @@ end subroutine gas_aer_uptkrates_1box1gas
 
 !--------------------------------------------------------------------------------
 !--------------------------------------------------------------------------------
-      subroutine modal_aero_amicphys_init( imozart, species_class,n_so4_monolayers_pcage_in)
+      subroutine modal_aero_amicphys_init( imozart, species_class,n_so4_monolayers_pcage_in, &
+                                           dp_cut_accum_rename_in, dp_xferall_thresh_accum_rename_in)
 
 !-----------------------------------------------------------------------
 !
@@ -6034,6 +6043,8 @@ subroutine modal_aero_amicphys_init( imozart, species_class,n_so4_monolayers_pca
    integer, intent(in)  :: imozart
    integer, intent(in)  :: species_class(:)
    real(r8), intent(in) :: n_so4_monolayers_pcage_in
+   real(r8), intent(in) :: dp_cut_accum_rename_in
+   real(r8), intent(in) :: dp_xferall_thresh_accum_rename_in
 
 !-----------------------------------------------------------------------
 ! local
@@ -6059,6 +6070,8 @@ subroutine modal_aero_amicphys_init( imozart, species_class,n_so4_monolayers_pca
 n_so4_monolayers_pcage  = n_so4_monolayers_pcage_in
 dr_so4_monolayers_pcage = n_so4_monolayers_pcage * 4.76e-10
 
+dp_cut_accum_rename = dp_cut_accum_rename_in
+dp_xferall_thresh_accum_rename = dp_xferall_thresh_accum_rename_in
 
 #if ( defined( CAMBOX_ACTIVATE_THIS ) )
       ldiag82  = .true.  ; lun82  = 82

components/eam/src/chemistry/modal_aero/modal_aero_initialize_data.F90

@@ -425,6 +425,8 @@ subroutine modal_aero_initialize(pbuf2d, imozart, species_class)
        integer  :: numaerosols     ! number of bulk aerosols in climate list
        character(len=20) :: bulkname
        real(r8) :: pi,n_so4_monolayers_pcage_in
+       real(r8) :: dp_cut_accum_rename_in
+       real(r8) :: dp_xferall_thresh_accum_rename_in
        complex(r8), pointer  :: refindex_aer_sw(:), &
             refindex_aer_lw(:)
        real(r8), pointer :: qqcw(:,:)
@@ -441,7 +443,9 @@ subroutine modal_aero_initialize(pbuf2d, imozart, species_class)
        call phys_getopts(convproc_do_gas_out = convproc_do_gas, &
             convproc_do_aer_out = convproc_do_aer, &
             mam_amicphys_optaa_out = mam_amicphys_optaa, &
-            n_so4_monolayers_pcage_out = n_so4_monolayers_pcage_in) 
+            n_so4_monolayers_pcage_out = n_so4_monolayers_pcage_in, &
+            dp_cut_accum_rename_out = dp_cut_accum_rename_in, &
+            dp_xferall_thresh_accum_rename_out = dp_xferall_thresh_accum_rename_in)
 
 
        ! Mode specific properties.
@@ -628,7 +632,9 @@ subroutine modal_aero_initialize(pbuf2d, imozart, species_class)
        if ( mam_amicphys_optaa > 0 ) then
           call modal_aero_calcsize_init( pbuf2d, species_class )
           call modal_aero_newnuc_init( mam_amicphys_optaa )
-          call modal_aero_amicphys_init( imozart, species_class,n_so4_monolayers_pcage_in )
+          call modal_aero_amicphys_init( imozart, species_class,n_so4_monolayers_pcage_in, &
+                                         dp_cut_accum_rename_in, &
+                                         dp_xferall_thresh_accum_rename_in)
        else
           call modal_aero_rename_init
           !   calcsize call must follow rename call

components/eam/src/physics/cam/phys_control.F90

@@ -125,6 +125,8 @@ module phys_control
                                                        !                                    newnuc, and coag routines) 
                                                        !  > 0 -- use new microphysics code (single call to amicphys routine)
 real(r8)          :: n_so4_monolayers_pcage = huge(1.0_r8) ! number of so4(+nh4) monolayers needed to "age" a carbon particle
+real(r8)          :: dp_cut_accum_rename = huge(1.0_r8) ! cutoff diameter for larger accum mode particles renamed to the strat coarse mode
+real(r8)          :: dp_xferall_thresh_accum_rename = huge(1.0_r8) ! threshold diameter for all accum mode particles transferred to the strat coarse mode 
 real(r8)          :: micro_mg_accre_enhan_fac = huge(1.0_r8) !!Accretion enhancement factor
 logical           :: liqcf_fix            = .false.    ! liq cld fraction fix calc.                     
 logical           :: regen_fix            = .false.    ! aerosol regeneration bug fix for ndrop.F90 
@@ -261,7 +263,8 @@ subroutine phys_ctl_readnl(nlfile)
       cld_macmic_num_steps, micro_do_icesupersat, &
       fix_g1_err_ndrop, ssalt_tuning, resus_fix, convproc_do_aer, &
       convproc_do_gas, convproc_method_activate, liqcf_fix, regen_fix, demott_ice_nuc, pergro_mods, pergro_test_active, &
-      mam_amicphys_optaa, n_so4_monolayers_pcage,micro_mg_accre_enhan_fac, &
+      mam_amicphys_optaa, n_so4_monolayers_pcage, dp_cut_accum_rename, &
+      dp_xferall_thresh_accum_rename, micro_mg_accre_enhan_fac, &
       cflx_cpl_opt, &
       l_tracer_aero, l_vdiff, l_rayleigh, l_gw_drag, l_ac_energy_chk, &
       l_bc_energy_fix, l_dry_adj, l_st_mac, l_st_mic, l_rad, prc_coef1,prc_exp,prc_exp1,cld_sed,mg_prc_coeff_fix, &
@@ -389,6 +392,8 @@ subroutine phys_ctl_readnl(nlfile)
    call mpibcast(convproc_method_activate,        1 , mpiint,  0, mpicom)
    call mpibcast(mam_amicphys_optaa,              1 , mpiint,  0, mpicom)
    call mpibcast(n_so4_monolayers_pcage,          1 , mpir8,   0, mpicom)
+   call mpibcast(dp_cut_accum_rename,             1 , mpir8,   0, mpicom)
+   call mpibcast(dp_xferall_thresh_accum_rename,  1 , mpir8,   0, mpicom)
    call mpibcast(micro_mg_accre_enhan_fac,        1 , mpir8,   0, mpicom)
    call mpibcast(liqcf_fix,                       1 , mpilog,  0, mpicom)
    call mpibcast(regen_fix,                       1 , mpilog,  0, mpicom)
@@ -620,6 +625,7 @@ subroutine phys_getopts(deep_scheme_out, shallow_scheme_out, eddy_scheme_out, &
                         cld_macmic_num_steps_out, micro_do_icesupersat_out, &
                         fix_g1_err_ndrop_out, ssalt_tuning_out,resus_fix_out,convproc_do_aer_out,  &
                         convproc_do_gas_out, convproc_method_activate_out, mam_amicphys_optaa_out, n_so4_monolayers_pcage_out, &
+                        dp_cut_accum_rename_out, dp_xferall_thresh_accum_rename_out, &
                         micro_mg_accre_enhan_fac_out, liqcf_fix_out, regen_fix_out,demott_ice_nuc_out, pergro_mods_out, pergro_test_active_out &
                        ,cflx_cpl_opt_out &
                        ,l_tracer_aero_out, l_vdiff_out, l_rayleigh_out, l_gw_drag_out, l_ac_energy_chk_out  &
@@ -712,6 +718,8 @@ subroutine phys_getopts(deep_scheme_out, shallow_scheme_out, eddy_scheme_out, &
    integer,           intent(out), optional :: convproc_method_activate_out 
    integer,           intent(out), optional :: mam_amicphys_optaa_out
    real(r8),          intent(out), optional :: n_so4_monolayers_pcage_out
+   real(r8),          intent(out), optional :: dp_cut_accum_rename_out
+   real(r8),          intent(out), optional :: dp_xferall_thresh_accum_rename_out
    real(r8),          intent(out), optional :: micro_mg_accre_enhan_fac_out
    logical,           intent(out), optional :: liqcf_fix_out       
    logical,           intent(out), optional :: regen_fix_out       
@@ -818,6 +826,8 @@ subroutine phys_getopts(deep_scheme_out, shallow_scheme_out, eddy_scheme_out, &
    if ( present(convproc_method_activate_out ) ) convproc_method_activate_out = convproc_method_activate
    if ( present(mam_amicphys_optaa_out  ) ) mam_amicphys_optaa_out  = mam_amicphys_optaa
    if ( present(n_so4_monolayers_pcage_out  ) ) n_so4_monolayers_pcage_out = n_so4_monolayers_pcage
+   if ( present(dp_cut_accum_rename_out ) ) dp_cut_accum_rename_out  = dp_cut_accum_rename
+   if ( present(dp_xferall_thresh_accum_rename_out ) ) dp_xferall_thresh_accum_rename_out = dp_xferall_thresh_accum_rename
    if ( present(micro_mg_accre_enhan_fac_out)) micro_mg_accre_enhan_fac_out = micro_mg_accre_enhan_fac
    if ( present(liqcf_fix_out           ) ) liqcf_fix_out            = liqcf_fix      
    if ( present(regen_fix_out           ) ) regen_fix_out            = regen_fix