testing CESM3_SLH using cam_cesm2_2_rel_09_slh

RafaPedroFernandez · Oct 1, 2024 · 535c5a0 · 535c5a0
1 parent 3b1a221
commit 535c5a0
Show file tree

Hide file tree

Showing 15 changed files with 340 additions and 2,085 deletions.
diff --git a/bld/namelist_files/namelist_defaults_cam.xml b/bld/namelist_files/namelist_defaults_cam.xml
@@ -1684,7 +1684,9 @@
 <gas_wetdep_ice_uptake_list>'HNO3'</gas_wetdep_ice_uptake_list>
 
 <!-- dry deposition method for gas-phase species -->
-<drydep_method>xactive_lnd</drydep_method>
+<!-- rpf_CESM3_SLH testing cam_cesm2_2_rel_09_slh -->
+<!-- <drydep_method>xactive_lnd</drydep_method> -->
+<!-- rpf_CESM3_SLH testing cam_cesm2_2_rel_09_slh -->
 
 
 <!-- SLH scaling factors for different compsets -->

diff --git a/bld/namelist_files/namelist_definition.xml b/bld/namelist_files/namelist_definition.xml
@@ -6680,7 +6680,8 @@ List of species that undergo dry deposition.
 Default: set by build-namelist.
 </entry>
 
-<entry id="drydep_method" type="char*16" category="dry_deposition"
+<!-- rpf_CESM3_SLH testing cam_cesm2_2_rel_09_slh -->
+<!--<entry id="drydep_method" type="char*16" category="dry_deposition"
        group="drydep_inparm" valid_values="xactive_atm,xactive_lnd,table" >
 Dry deposition method used.  This specifies the method used to calculate dry
 deposition velocities of gas-phase chemical species.  The available methods
@@ -6689,7 +6690,17 @@ are:
  'xactive_atm' - interactive method in CAM
  'xactive_lnd' - interactive method in CLM
 Default: set by build-namelist
+</entry>-->
+<!-- rpf_CESM3_SLH testing cam_cesm2_2_rel_09_slh -->
+
+<!-- rpf_CESM3_SLH testing cam_cesm2_2_rel_09_slh -->
+<entry id="dep_data_file" type="char*500" input_pathname="abs" category="dry_deposition"
+       group="drydep_inparm" valid_values="" >
+Full pathname of file containing gas phase deposition data including effective
+Henry's law coefficients.
+Default: set by build-namelist.
 </entry>
+<!-- rpf_CESM3_SLH testing cam_cesm2_2_rel_09_slh -->
 
 <entry id="rxn_rate_sums" type="char*256(200)" category="cam_chem"
        group="rxn_rate_diags_nl" valid_values="" >

diff --git a/bld/namelist_files/use_cases/2000_cam4_trop_chem.xml b/bld/namelist_files/use_cases/2000_cam4_trop_chem.xml
@@ -33,7 +33,9 @@
 <prescribed_ozone_cycle_yr> 2000                                  </prescribed_ozone_cycle_yr>
 
 <chem_rad_passive>.true.</chem_rad_passive>
-<drydep_method>'xactive_lnd'</drydep_method>
+<!-- rpf_CESM3_SLH testing cam_cesm2_2_rel_09_slh -->
+<!-- <drydep_method>'xactive_lnd'</drydep_method> -->
+<!-- rpf_CESM3_SLH testing cam_cesm2_2_rel_09_slh -->
 
 <!-- sim_year used for CLM datasets -->
 <sim_year>2000</sim_year>

diff --git a/bld/namelist_files/use_cases/sd_waccm_sulfur.xml b/bld/namelist_files/use_cases/sd_waccm_sulfur.xml
@@ -23,7 +23,9 @@
 <ncdata dyn="fv"  hgrid="1.9x2.5" >atm/waccm/ic/f40.2008-2011.2deg.wcm.sd.carma.sulf.002.cam2.i.2010-01-01-00000.nc</ncdata>
 
 <!-- Deposition Velocities -->
-<drydep_method>'xactive_atm'</drydep_method>
+<!-- rpf_CESM3_SLH testing cam_cesm2_2_rel_09_slh -->
+<!-- <drydep_method>'xactive_atm'</drydep_method> -->
+<!-- rpf_CESM3_SLH testing cam_cesm2_2_rel_09_slh -->
 
 
 <!-- Emissions DEFAULTS OKAY for: airpl_emis_file,ch2o_emis_file,co_emis_file,nox_emis_file-->

diff --git a/chem_proc b/chem_proc
diff --git a/src/atmos_phys b/src/atmos_phys
diff --git a/src/chemistry/mozart/mo_chemini.F90 b/src/chemistry/mozart/mo_chemini.F90
@@ -54,8 +54,10 @@ subroutine chemini &
     use mo_sulf,           only : sulf_inti
     use mo_photo,          only : photo_inti
     use mo_lightning,      only : lightning_inti
+!rpf_CESM3_SLH testing cam_cesm2_2_rel_09_slh
     use mo_drydep,         only : drydep_inti
-    use seq_drydep_mod,    only : DD_XLND, DD_XATM, drydep_method
+!     use seq_drydep_mod,    only : DD_XLND, DD_XATM, drydep_method
+!rpf_CESM3_SLH testing cam_cesm2_2_rel_09_slh
     use mo_imp_sol,        only : imp_slv_inti
     use mo_exp_sol,        only : exp_sol_inti
     use spmd_utils,        only : iam
@@ -193,11 +195,17 @@ subroutine chemini &
     !-----------------------------------------------------------------------
     !	... initialize the dry deposition module
     !-----------------------------------------------------------------------
-    if ( drydep_method == DD_XATM .or. drydep_method == DD_XLND ) then
-       call drydep_inti(depvel_lnd_file, clim_soilw_file, season_wes_file )
-    else
-       call drydep_inti( depvel_file )
-    endif
+!rpf_CESM3_SLH testing cam_cesm2_2_rel_09_slh
+!     if ( drydep_method == DD_XATM .or. drydep_method == DD_XLND ) then
+!        call drydep_inti(depvel_lnd_file, clim_soilw_file, season_wes_file )
+!     else
+!        call drydep_inti( depvel_file )
+!     endif
+!rpf_CESM3_SLH testing cam_cesm2_2_rel_09_slh
+
+!rpf_CESM3_SLH testing cam_cesm2_2_rel_09_slh
+    call drydep_inti(depvel_lnd_file)
+!rpf_CESM3_SLH testing cam_cesm2_2_rel_09_slh
 
     if (masterproc) write(iulog,*) 'chemini: after drydep_inti on node ',iam