Skip to content

Commit

Permalink
Skip the methane balance check if dynamics lakes are on and it's the …
Browse files Browse the repository at this point in the history
…beginning of the year see ESCOMP#1356
  • Loading branch information
ekluzek committed Aug 29, 2021
1 parent 55ad5ab commit a9eef2c
Showing 1 changed file with 11 additions and 0 deletions.
11 changes: 11 additions & 0 deletions src/biogeochem/ch4Mod.F90
Original file line number Diff line number Diff line change
Expand Up @@ -1684,6 +1684,8 @@ subroutine ch4 (bounds, num_soilc, filter_soilc, num_lakec, filter_lakec, &
use ch4varcon , only : replenishlakec, allowlakeprod, ch4offline
use clm_varcon , only : secspday
use ch4varcon , only : finundation_mtd, finundation_mtd_h2osfc
use clm_time_manager, only : is_beg_curr_year
use dynSubgridControlMod, only : get_do_transient_lakes
!
! !ARGUMENTS:
type(bounds_type) , intent(in) :: bounds
Expand Down Expand Up @@ -2318,7 +2320,16 @@ subroutine ch4 (bounds, num_soilc, filter_soilc, num_lakec, filter_lakec, &
ch4_inst%totcolch4_grc(begg:endg), &
c2l_scale_type= 'unity', l2g_scale_type='unity' )

!
! Gricell level balance
!

! Skip the check if it's the beginning of a new year and dynamic lakes are on
! See (https://github.com/ESCOMP/CTSM/issues/1356#issuecomment-905963583)
!
if ( is_beg_curr_year() .and. get_do_transient_lakes() )then
ch4_first_time_grc(begg:endg) = .true.
end if

do g = begg, endg
if (.not. ch4_first_time_grc(g)) then
Expand Down

0 comments on commit a9eef2c

Please sign in to comment.