cam6_4_028: Dust emissions

.gitmodules

@@ -178,8 +178,8 @@ fxDONOTUSEurl = https://github.com/NCAR/ParallelIO
 
 [submodule "clm"]
 path = components/clm
-url = https://github.com/ESCOMP/CTSM
-fxtag = ctsm5.2.009
+url = https://github.com/ekluzek/CTSM.git
+fxtag = 8478b73
 fxrequired = ToplevelRequired
 fxDONOTUSEurl = https://github.com/ESCOMP/CTSM
 
@@ -189,4 +189,3 @@ url = https://github.com/ESCOMP/FMS_interface
 fxtag = fi_240516
 fxrequired = ToplevelRequired
 fxDONOTUSEurl = https://github.com/ESCOMP/FMS_interface
-

src/chemistry/bulk_aero/dust_model.F90

@@ -1,10 +1,13 @@
 !===============================================================================
 ! Dust for Bulk Aerosol Model
 !===============================================================================
-module dust_model 
+module dust_model
   use shr_kind_mod,    only: r8 => shr_kind_r8, cl => shr_kind_cl
   use spmd_utils,      only: masterproc
   use cam_abortutils,  only: endrun
+  use cam_logfile,     only: iulog
+  use shr_dust_emis_mod,only: is_dust_emis_zender, is_zender_soil_erod_from_atm
+  use cam_control_mod, only: aqua_planet
 
   implicit none
   private
@@ -36,6 +39,7 @@ module dust_model
 
   real(r8)          :: dust_emis_fact = -1.e36_r8        ! tuning parameter for dust emissions
   character(len=cl) :: soil_erod_file = 'soil_erod_file' ! full pathname for soil erodibility dataset
+
 contains
 
   !=============================================================================
@@ -44,8 +48,8 @@ module dust_model
   subroutine dust_readnl(nlfile)
 
     use namelist_utils,  only: find_group_name
-    use units,           only: getunit, freeunit
-    use mpishorthand
+    use spmd_utils,      only: mpicom, masterprocid, mpi_character, mpi_real8, mpi_success
+    use shr_dust_emis_mod, only: shr_dust_emis_readnl
 
     character(len=*), intent(in) :: nlfile  ! filepath for file containing namelist input
 
@@ -56,11 +60,11 @@ subroutine dust_readnl(nlfile)
     namelist /dust_nl/ dust_emis_fact, soil_erod_file
 
     !-----------------------------------------------------------------------------
+    if (aqua_planet) return
 
     ! Read namelist
     if (masterproc) then
-       unitn = getunit()
-       open( unitn, file=trim(nlfile), status='old' )
+       open( newunit=unitn, file=trim(nlfile), status='old' )
        call find_group_name(unitn, 'dust_nl', status=ierr)
        if (ierr == 0) then
           read(unitn, dust_nl, iostat=ierr)
@@ -69,14 +73,26 @@ subroutine dust_readnl(nlfile)
           end if
        end if
        close(unitn)
-       call freeunit(unitn)
     end if
 
-#ifdef SPMD
     ! Broadcast namelist variables
-    call mpibcast(dust_emis_fact, 1,                   mpir8,   0, mpicom)
-    call mpibcast(soil_erod_file, len(soil_erod_file), mpichar, 0, mpicom)
-#endif
+    call mpi_bcast(soil_erod_file, len(soil_erod_file), mpi_character, masterprocid, mpicom, ierr)
+    if (ierr/=mpi_success) then
+       call endrun(subname//' MPI_BCAST ERROR: soil_erod_file')
+    end if
+    call mpi_bcast(dust_emis_fact, 1, mpi_real8, masterprocid, mpicom, ierr)
+    if (ierr/=mpi_success) then
+       call endrun(subname//' MPI_BCAST ERROR: dust_emis_fact')
+    end if
+
+    call shr_dust_emis_readnl(mpicom, 'drv_flds_in')
+
+    if (masterproc) then
+       write(iulog,*) subname,': soil_erod_file = ',trim(soil_erod_file)
+       write(iulog,*) subname,': dust_emis_fact = ',dust_emis_fact
+       write(iulog,*) subname,': is_dust_emis_zender : ',is_dust_emis_zender()
+       write(iulog,*) subname,': is_zender_soil_erod_from_atm : ',is_zender_soil_erod_from_atm()
+    end if
 
   end subroutine dust_readnl
 
@@ -89,13 +105,17 @@ subroutine dust_init()
 
     integer :: n
 
+    if (aqua_planet) return
+
     do n = 1, dust_nbin
        call cnst_get_ind(dust_names(n), dust_indices(n),abort=.false.)
     end do
     dust_active = any(dust_indices(:) > 0)
     if (.not.dust_active) return
 
-    call  soil_erod_init( dust_emis_fact, soil_erod_file )
+    if (is_zender_soil_erod_from_atm()) then
+       call  soil_erod_init( dust_emis_fact, soil_erod_file )
+    endif
 
     call dust_set_params( dust_nbin, dust_dmt_grd, dust_dmt_vwr, dust_stk_crc )
 
@@ -115,25 +135,36 @@ subroutine dust_emis( ncol, lchnk, dust_flux_in, cflx, soil_erod )
 
    ! local vars
     integer :: i, m, idst
+    real(r8) :: erodfctr(ncol)
     real(r8), parameter :: dust_emis_sclfctr(dust_nbin) &
          = (/ 0.011_r8/0.032456_r8, 0.087_r8/0.174216_r8, 0.277_r8/0.4085517_r8, 0.625_r8/0.384811_r8 /)
 
+    if (.not.dust_active) return
+
     ! set dust emissions
 
-    col_loop: do i =1,ncol
+    if (is_zender_soil_erod_from_atm()) then
+       do i =1,ncol
+          ! adjust emissions based on soil erosion
+          soil_erod(i) = soil_erodibility( i, lchnk )
+          erodfctr(i) = soil_erod(i)/soil_erod_fact*1.15_r8
+       end do
+    else
+       erodfctr(:) = 1._r8
+    end if
+
+    ! rebin dust emissons
 
-       soil_erod(i) = soil_erodibility( i, lchnk )
+    do i =1,ncol
 
-       ! adjust emissions based on soil erosion 
        do m = 1,dust_nbin
 
           idst = dust_indices(m)
-          cflx(i,idst) = -dust_flux_in(i,m) &
-               * dust_emis_sclfctr(m)*soil_erod(i)/soil_erod_fact*1.15_r8
+          cflx(i,idst) = -dust_flux_in(i,m) * dust_emis_sclfctr(m) * erodfctr(i)
 
        enddo
 
-    end do col_loop
+    end do
 
   end subroutine dust_emis
 

src/chemistry/modal_aero/dust_model.F90

@@ -6,6 +6,9 @@ module dust_model
   use spmd_utils,       only: masterproc
   use cam_abortutils,   only: endrun
   use modal_aero_data,  only: ntot_amode, ndst=>nDust
+  use cam_logfile,      only: iulog
+  use shr_dust_emis_mod,only: is_dust_emis_zender, is_zender_soil_erod_from_atm
+  use cam_control_mod, only: aqua_planet
 
   implicit none
   private
@@ -30,8 +33,8 @@ module dust_model
   real(r8), allocatable :: dust_dmt_vwr(:)
   real(r8), allocatable :: dust_stk_crc(:)
 
-  real(r8)          :: dust_emis_fact = -1.e36_r8        ! tuning parameter for dust emissions
-  character(len=cl) :: soil_erod_file = 'soil_erod_file' ! full pathname for soil erodibility dataset
+  real(r8)          :: dust_emis_fact = -huge(1._r8)       ! tuning parameter for dust emissions
+  character(len=cl) :: soil_erod_file = 'soil_erod_file'   ! full pathname for soil erodibility dataset
 
   logical :: dust_active = .false.
 
@@ -43,8 +46,8 @@ module dust_model
   subroutine dust_readnl(nlfile)
 
     use namelist_utils,  only: find_group_name
-    use units,           only: getunit, freeunit
-    use mpishorthand
+    use spmd_utils,      only: mpicom, masterprocid, mpi_character, mpi_real8, mpi_success
+    use shr_dust_emis_mod, only: shr_dust_emis_readnl
 
     character(len=*), intent(in) :: nlfile  ! filepath for file containing namelist input
 
@@ -55,11 +58,11 @@ subroutine dust_readnl(nlfile)
     namelist /dust_nl/ dust_emis_fact, soil_erod_file
 
     !-----------------------------------------------------------------------------
+    if (aqua_planet) return
 
     ! Read namelist
     if (masterproc) then
-       unitn = getunit()
-       open( unitn, file=trim(nlfile), status='old' )
+       open( newunit=unitn, file=trim(nlfile), status='old' )
        call find_group_name(unitn, 'dust_nl', status=ierr)
        if (ierr == 0) then
           read(unitn, dust_nl, iostat=ierr)
@@ -68,14 +71,26 @@ subroutine dust_readnl(nlfile)
           end if
        end if
        close(unitn)
-       call freeunit(unitn)
     end if
 
-#ifdef SPMD
     ! Broadcast namelist variables
-    call mpibcast(dust_emis_fact, 1,                   mpir8,   0, mpicom)
-    call mpibcast(soil_erod_file, len(soil_erod_file), mpichar, 0, mpicom)
-#endif
+    call mpi_bcast(soil_erod_file, len(soil_erod_file), mpi_character, masterprocid, mpicom, ierr)
+    if (ierr/=mpi_success) then
+       call endrun(subname//' MPI_BCAST ERROR: soil_erod_file')
+    end if
+    call mpi_bcast(dust_emis_fact, 1, mpi_real8, masterprocid, mpicom, ierr)
+    if (ierr/=mpi_success) then
+       call endrun(subname//' MPI_BCAST ERROR: dust_emis_fact')
+    end if
+
+    call shr_dust_emis_readnl(mpicom, 'drv_flds_in')
+
+    if (masterproc) then
+       write(iulog,*) subname,': soil_erod_file = ',trim(soil_erod_file)
+       write(iulog,*) subname,': dust_emis_fact = ',dust_emis_fact
+       write(iulog,*) subname,': is_dust_emis_zender : ',is_dust_emis_zender()
+       write(iulog,*) subname,': is_zender_soil_erod_from_atm : ',is_zender_soil_erod_from_atm()
+    end if
 
   end subroutine dust_readnl
 
@@ -91,6 +106,8 @@ subroutine dust_init()
     character(len=32) :: spec_name
     integer, parameter :: mymodes(7) = (/ 2, 1, 3, 4, 5, 6, 7 /) ! tricky order ...
 
+    if (aqua_planet) return
+
     dust_nbin = ndst
     dust_nnum = ndst
 
@@ -131,7 +148,9 @@ subroutine dust_init()
     dust_active = any(dust_indices(:) > 0)
     if (.not.dust_active) return
 
-    call  soil_erod_init( dust_emis_fact, soil_erod_file )
+    if (is_zender_soil_erod_from_atm()) then
+       call  soil_erod_init( dust_emis_fact, soil_erod_file )
+    end if
 
     call dust_set_params( dust_nbin, dust_dmt_grd, dust_dmt_vwr, dust_stk_crc )
 
@@ -155,32 +174,31 @@ subroutine dust_emis( ncol, lchnk, dust_flux_in, cflx, soil_erod )
     integer :: i, m, idst, inum
     real(r8) :: x_mton
     real(r8),parameter :: soil_erod_threshold = 0.1_r8
+    real(r8) :: erodfctr(ncol)
 
     ! set dust emissions
 
-    col_loop: do i =1,ncol
-
-       soil_erod(i) = soil_erodibility( i, lchnk )
+    if (is_zender_soil_erod_from_atm()) then
+       col_loop: do i =1,ncol
+          soil_erod(i) = soil_erodibility( i, lchnk )
+          if( soil_erod(i) .lt. soil_erod_threshold ) soil_erod(i) = 0._r8
+          erodfctr(i) = soil_erod(i)/soil_erod_fact*1.15_r8
+       end do col_loop
+    else
+       erodfctr(:) = 1._r8
+    end if
 
-       if( soil_erod(i) .lt. soil_erod_threshold ) soil_erod(i) = 0._r8
+    ! rebin dust emissons
 
-       ! rebin and adjust dust emissons..
+    do i = 1,ncol
        do m = 1,dust_nbin
-
           idst = dust_indices(m)
-
-          cflx(i,idst) = sum( -dust_flux_in(i,:) ) &
-               * dust_emis_sclfctr(m)*soil_erod(i)/soil_erod_fact*1.15_r8
-
+          cflx(i,idst) = sum( -dust_flux_in(i,:) ) * dust_emis_sclfctr(m) * erodfctr(i) ! mass mixing ratio
           x_mton = 6._r8 / (pi * dust_density * (dust_dmt_vwr(m)**3._r8))
-
           inum = dust_indices(m+dust_nbin)
-
-          cflx(i,inum) = cflx(i,idst)*x_mton
-
-       enddo
-
-    end do col_loop
+          cflx(i,inum) = cflx(i,idst)*x_mton ! number mixing ratio
+       end do
+    end do
 
   end subroutine dust_emis