Skip to content

Commit

Permalink
Remove the ifs around different dust-emis methods for the Zender case…
Browse files Browse the repository at this point in the history
… as pointed out in the review
  • Loading branch information
ekluzek committed Jun 18, 2024
1 parent 8e69f12 commit e291534
Showing 1 changed file with 5 additions and 21 deletions.
26 changes: 5 additions & 21 deletions src/biogeochem/DustEmisZender2003.F90
Original file line number Diff line number Diff line change
Expand Up @@ -31,7 +31,6 @@ module DustEmisZender2003
use PatchType , only : patch
use ZenderSoilErodStreamType, only : soil_erod_stream_type
use DustEmisBase , only : dust_emis_base_type
use clm_varctl , only : dust_emis_method
!
! !PUBLIC TYPES
implicit none
Expand Down Expand Up @@ -160,6 +159,7 @@ end subroutine InitHistory

subroutine InitCold(this, bounds)
!
! Initialize values from a cold start
! !ARGUMENTS:
class (dust_emis_zender2003_type) :: this
type(bounds_type), intent(in) :: bounds
Expand All @@ -168,27 +168,11 @@ subroutine InitCold(this, bounds)
integer :: c,l
!-----------------------------------------------------------------------

! Set basin factor to 1 for now

if (dust_emis_method == 'Leung_2023') then
!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
call endrun( msg="Leung_2023 dust_emis_method is currently not available"//errMsg(sourcefile, __LINE__))
else if (dust_emis_method == 'Zender_2003') then
if ( this%soil_erod_stream%UseStreams() )then
call this%soil_erod_stream%CalcDustSource( bounds, &
this%mbl_bsn_fct_col(bounds%begc:bounds%endc) )
else
this%mbl_bsn_fct_col(:) = 1.0_r8
end if
if ( this%soil_erod_stream%UseStreams() )then
call this%soil_erod_stream%CalcDustSource( bounds, &
this%mbl_bsn_fct_col(bounds%begc:bounds%endc) )
else
write(iulog,*) 'dust_emis_method not recognized = ', trim(dust_emis_method)
call endrun( msg="dust_emis_method namelist item is not valid "//errMsg(sourcefile, __LINE__))
this%mbl_bsn_fct_col(:) = 1.0_r8
end if

end subroutine InitCold
Expand Down

0 comments on commit e291534

Please sign in to comment.