Move the dust emission source function soil erodibility for the Zender scheme from CAM to CTSM

bld/CLMBuildNamelist.pm

@@ -1560,6 +1560,7 @@ sub process_namelist_inline_logic {
   setup_logic_snowpack($opts,  $nl_flags, $definition, $defaults, $nl);
   setup_logic_fates($opts,  $nl_flags, $definition, $defaults, $nl);
   setup_logic_misc($opts, $nl_flags, $definition, $defaults, $nl);
+  setup_logic_zender_soilerod($opts, $nl_flags, $definition, $defaults, $nl); # dmleung added 11 Mar 2023
 
   #########################################
   # namelist group: atm2lnd_inparm
@@ -4187,6 +4188,15 @@ sub setup_logic_misc {
 
 #-------------------------------------------------------------------------------
 
+sub setup_logic_zender_soilerod {
+  # dmleung added on 31 Dec 2022 for reading Prigent's roughness stream file
+  my ($opts, $nl_flags, $definition, $defaults, $nl) = @_;
+  add_default($opts, $nl_flags->{'inputdata_rootdir'}, $definition, $defaults, $nl, 'stream_fldfilename_zendersoilerod' );
+  add_default($opts, $nl_flags->{'inputdata_rootdir'}, $definition, $defaults, $nl, 'stream_meshfile_zendersoilerod' );
+}
+
+#-------------------------------------------------------------------------------
+
 sub write_output_files {
   my ($opts, $nl_flags, $defaults, $nl) = @_;
 
@@ -4235,6 +4245,7 @@ sub write_output_files {
   push @groups, "ch4finundated";
   push @groups, "soilbgc_decomp";
   push @groups, "clm_canopy_inparm";
+  push @groups, "zendersoilerod"; # dmleung 11 Mar 2023
   if (remove_leading_and_trailing_quotes($nl->get_value('snow_cover_fraction_method')) eq 'SwensonLawrence2012') {
      push @groups, "scf_swenson_lawrence_2012_inparm";
   }

bld/namelist_files/namelist_defaults_ctsm.xml

@@ -2537,6 +2537,18 @@ lnd/clm2/surfdata_map/release-clm5.0.30/surfdata_ne0np4.CONUS.ne30x8_hist_78pfts
 <stream_meshfile_ch4finundated finundation_method="ZWT_inversion"
 >lnd/clm2/paramdata/finundated_inversiondata_0.9x1_ESMFmesh_cdf5_130621.nc</stream_meshfile_ch4finundated>
 
+<!-- ===========================================  -->
+<!-- Defaults for dust namelist items             -->
+<!-- Defaults for Zender's soil erodibility data  -->
+<!-- dmleung                                      -->
+<!-- 11 Mar 2023                                  -->
+<!-- ===========================================  -->
+
+<stream_fldfilename_zendersoilerod
+>/gpfs/fs1/work/dleung/Zender_dust_source/dst_source2x2tunedcam6-2x2-forCLM_c230311.nc</stream_fldfilename_zendersoilerod>
+<stream_meshfile_zendersoilerod
+>/gpfs/fs1/work/dleung/Zender_dust_source/lnd_mesh.nc</stream_meshfile_zendersoilerod>
+
 <!-- =========================================  -->
 <!-- Defaults for different BGC/decomp modes    -->
 <!-- =========================================  -->

bld/namelist_files/namelist_definition_ctsm.xml

@@ -1693,6 +1693,21 @@ Mapping method from Nitrogen deposition input file to the model resolution
     copy     = copy using the same indices
 </entry>
 
+<!-- ========================================================================================  -->
+<!-- Dust namelist and Zender soil erodibility stream -->
+<!-- dmleung added 11 Mar 2023 -->
+<!-- ========================================================================================  -->
+
+<entry id="stream_fldfilename_zendersoilerod" type="char*256" category="datasets"
+       input_pathname="abs" group="zendersoilerod" valid_values="" >
+Filename of input stream data for Zender's soil erodibility source function
+</entry>
+
+<entry id="stream_meshfile_zendersoilerod" type="char*256" category="datasets"
+       input_pathname="abs" group="zendersoilerod" valid_values="" >
+mesh filename of input stream data for Zender's soil erodibility source function
+</entry>
+
 <!-- ========================================================================================  -->
 <!-- finundated stream -->
 <!-- ========================================================================================  -->

src/biogeochem/DUSTMod.F90

@@ -30,6 +30,7 @@ module DUSTMod
   use LandunitType         , only : lun
   use ColumnType           , only : col
   use PatchType            , only : patch
+  use ZenderSoilErodStreamType,  only : zendersoilerodstream_type   ! dmleung added 11 Mar 2023
   !  
   ! !PUBLIC TYPES
   implicit none
@@ -59,7 +60,9 @@ module DUSTMod
      real(r8), pointer, private :: vlc_trb_2_patch           (:)   ! turbulent deposition velocity 2(m/s)
      real(r8), pointer, private :: vlc_trb_3_patch           (:)   ! turbulent deposition velocity 3(m/s)
      real(r8), pointer, private :: vlc_trb_4_patch           (:)   ! turbulent deposition velocity 4(m/s)
-     real(r8), pointer, private :: mbl_bsn_fct_col           (:)   ! basin factor
+     !########### added by dmleung 11 Mar 2023 ########################################################################
+     type(zendersoilerodstream_type), private :: zendersoilerodstream      ! Zender soil erodibility stream data
+     real(r8), pointer, private :: mbl_bsn_fct_col           (:)   ! [dimensionless] basin factor, or soil erodibility, time-constant
 
    contains
 
@@ -78,13 +81,17 @@ module DUSTMod
 contains
 
   !------------------------------------------------------------------------
-  subroutine Init(this, bounds)
+  !subroutine Init(this, bounds)
+  !##### dmleung edited for initializing stream files 11 Mar 2023  ########
+  subroutine Init(this, bounds, NLFilename)
 
     class(dust_type) :: this
     type(bounds_type), intent(in) :: bounds  
+    character(len=*),  intent(in) :: NLFilename   ! dmleung added 11 Mar 2023
 
     call this%InitAllocate (bounds)
     call this%InitHistory  (bounds)
+    call this%zendersoilerodstream%Init( bounds, NLFilename )  ! dmleung added 11 Mar 2023
     call this%InitCold     (bounds)
     call this%InitDustVars (bounds)
 
@@ -129,9 +136,11 @@ subroutine InitHistory(this, bounds)
     !
     ! !LOCAL VARIABLES:
     integer :: begp,endp
+    integer :: begc,endc  ! dmleung added 11 Mar 2023
     !------------------------------------------------------------------------
 
     begp = bounds%begp; endp = bounds%endp
+    begc = bounds%begc; endc = bounds%endp
 
     this%flx_mss_vrt_dst_tot_patch(begp:endp) = spval
     call hist_addfld1d (fname='DSTFLXT', units='kg/m2/s',  &
@@ -158,6 +167,13 @@ subroutine InitHistory(this, bounds)
          avgflag='A', long_name='turbulent deposition velocity 4', &
          ptr_patch=this%vlc_trb_4_patch, default='inactive')
 
+    !#####added by dmleung 11 Mar 2023 ########################################
+    this%mbl_bsn_fct_col(begc:endc) = spval
+    call hist_addfld1d (fname='LND_MBL', units='fraction',  &
+         avgflag='A', long_name='Soil erodibility factor', &
+         ptr_col=this%mbl_bsn_fct_col)
+    !##########################################################################
+
   end subroutine InitHistory
 
   !-----------------------------------------------------------------------
@@ -173,13 +189,18 @@ subroutine InitCold(this, bounds)
 
     ! Set basin factor to 1 for now
 
-    do c = bounds%begc, bounds%endc
-       l = col%landunit(c)
+    !if (dust_emission_scheme == 'Leung2023')
+    !   do c = bounds%begc, bounds%endc
+    !      l = col%landunit(c)
 
-       if (.not.lun%lakpoi(l)) then
-          this%mbl_bsn_fct_col(c) = 1.0_r8
-       end if
-    end do
+    !      if (.not.lun%lakpoi(l)) then
+    !         this%mbl_bsn_fct_col(c) = 1.0_r8
+    !      end if
+    !   end do
+    !else if (dust_emission_scheme == 'Zender2003')
+       call this%zendersoilerodstream%CalcDustSource( bounds, &
+                                this%mbl_bsn_fct_col(bounds%begc:bounds%endc) )
+    !end if
 
   end subroutine InitCold
 
@@ -262,7 +283,7 @@ subroutine DustEmission (bounds, &
          fv                  => frictionvel_inst%fv_patch            , & ! Input:  [real(r8) (:)   ]  friction velocity (m/s) (for dust model)          
          u10                 => frictionvel_inst%u10_patch           , & ! Input:  [real(r8) (:)   ]  10-m wind (m/s) (created for dust model)          
 
-         mbl_bsn_fct         => dust_inst%mbl_bsn_fct_col            , & ! Input:  [real(r8) (:)   ]  basin factor                                      
+         mbl_bsn_fct         => dust_inst%mbl_bsn_fct_col            , & ! Input:  [real(r8) (:)   ]  basin factor  ; dmleung changed the value in InitCold 11 Mar 2023       
          flx_mss_vrt_dst     => dust_inst%flx_mss_vrt_dst_patch      , & ! Output: [real(r8) (:,:) ]  surface dust emission (kg/m**2/s)               
          flx_mss_vrt_dst_tot => dust_inst%flx_mss_vrt_dst_tot_patch    & ! Output: [real(r8) (:)   ]  total dust flux back to atmosphere (pft)
          )