Skip to content

Commit

Permalink
Address bug #36, effective LAI, total SAI and effective SAI were bein…
Browse files Browse the repository at this point in the history
…g zero'd strangely in ed_clm_link.

It appears that this part of the code is responsible for zero'ing the bare-ground portion of the code.  There is also the
assumption that patches that are not within the perview of ED should be zero'd here too. This would be crop patches, as well as the
bear-ground patches (grrr) just mentioned.
  • Loading branch information
rgknox committed May 13, 2016
1 parent f11a43f commit e23aeb7
Showing 1 changed file with 7 additions and 6 deletions.
13 changes: 7 additions & 6 deletions main/EDCLMLinkMod.F90
Original file line number Diff line number Diff line change
Expand Up @@ -1040,12 +1040,13 @@ subroutine ed_clm_link( this, bounds, ed_allsites_inst, ed_phenology_inst, &

clmpatch%is_veg(begp:endp) = .false.
clmpatch%is_bareground(begp:endp) = .false.

tlai(begp:endp) = 0.0_r8
elai(firstsoilpatch(g)) = 0.0_r8
tsai(firstsoilpatch(g)) = 0.0_r8
esai(firstsoilpatch(g)) = 0.0_r8
htop(begp:endp) = 0.0_r8
hbot(begp:endp) = 0.0_r8
elai(begp:endp) = 0.0_r8
tsai(begp:endp) = 0.0_r8
esai(begp:endp) = 0.0_r8

do g = begg,endg

Expand All @@ -1056,9 +1057,9 @@ subroutine ed_clm_link( this, bounds, ed_allsites_inst, ed_phenology_inst, &
! Zero the bare ground tile BGC variables.
! ============================================================================

tlai(firstsoilpatch(g)) = 0.0_r8
htop(firstsoilpatch(g)) = 0.0_r8
hbot(firstsoilpatch(g)) = 0.0_r8
!tlai(firstsoilpatch(g)) = 0.0_r8 ! This is redundant (test removal, RGK)
!htop(firstsoilpatch(g)) = 0.0_r8 ! This is redundant (test removal, RGK)
!hbot(firstsoilpatch(g)) = 0.0_r8 ! This is redundant (test removal, RGK)

patchn = 0
total_bare_ground = 0.0_r8
Expand Down

0 comments on commit e23aeb7

Please sign in to comment.