Skip to content

Commit

Permalink
remove write statement
Browse files Browse the repository at this point in the history
  • Loading branch information
emorway-usgs committed Nov 8, 2024
1 parent 0ae1d10 commit 60363e5
Showing 1 changed file with 2 additions and 2 deletions.
4 changes: 2 additions & 2 deletions autotest/TestUzfEtUtil.f90
Original file line number Diff line number Diff line change
Expand Up @@ -124,6 +124,7 @@ subroutine test_etfunc_nlin(error)
real(DP) :: celtop !< elevation of the top of the cell

! water table exactly in the middle of the extinction depth, should return pET
atol = 1d-12
deriv_et = DZERO
extdp = DONE
pET = DEM1
Expand All @@ -134,6 +135,7 @@ subroutine test_etfunc_nlin(error)
deriv_et = DZERO
trhs = DZERO
thcof = DZERO

rate1 = etfunc_nlin(celtop, extdp, pET, deriv_et, trhs, thcof, hgwf)
call check(error, is_close(rate1, pET))
if (allocated(error)) return
Expand All @@ -152,7 +154,6 @@ subroutine test_etfunc_nlin(error)
range = DEM3 * extdp ! an intermediate calc
hgwf = celtop - extdp + (range * DHALF)
rate1 = etfunc_nlin(celtop, extdp, pET, deriv_et, trhs, thcof, hgwf)
atol = 1d-12
call check(error, is_close(rate1, pET * DHALF, atol=atol))
! however, if water table is above (or in this case, at) the calculated "range",
! there should be no scaling
Expand All @@ -164,7 +165,6 @@ subroutine test_etfunc_nlin(error)
! in no gwet
hgwf = celtop - extdp
rate1 = etfunc_nlin(celtop, extdp, pET, deriv_et, trhs, thcof, hgwf)
write (*, *) 'error: ', rate1 - DZERO
call check(error, is_close(rate1, DZERO, atol=atol))

end subroutine test_etfunc_nlin
Expand Down

0 comments on commit 60363e5

Please sign in to comment.