A more straight forward fix

emorway-usgs · Jul 24, 2024 · 9d5c093 · 9d5c093
1 parent b9e7a00
commit 9d5c093
Showing 1 changed file with 2 additions and 5 deletions.
diff --git a/src/Model/ModelUtilities/UzfCellGroup.f90 b/src/Model/ModelUtilities/UzfCellGroup.f90
@@ -859,9 +859,7 @@ subroutine simgwet(this, igwetflag, icell, hgwf, trhs, thcof, det)
     real(DP), intent(inout) :: det
     ! -- local
     real(DP) :: s, x, c, b, et
-    real(DP) :: signflip
     !
-    signflip = DONE
     this%gwet(icell) = DZERO
     trhs = DZERO
     thcof = DZERO
@@ -877,12 +875,11 @@ subroutine simgwet(this, igwetflag, icell, hgwf, trhs, thcof, det)
                       this%celtop(icell), this%celbot(icell))
     else if (igwetflag == 2) then
       et = etfunc_nlin(s, x, c, det, trhs, thcof, hgwf)
-      signflip = -DONE
     end if
     ! this%gwet(icell) = et * this%uzfarea(icell)
     trhs = trhs * this%uzfarea(icell)
     thcof = thcof * this%uzfarea(icell)
-    this%gwet(icell) = (trhs - (thcof * hgwf)) * signflip
+    this%gwet(icell) = trhs - (thcof * hgwf)
     ! write(99,*)'in group', icell, this%gwet(icell)
     !
     ! -- Return
@@ -971,7 +968,7 @@ function etfunc_nlin(s, x, c, det, trhs, thcof, hgwf)
     range = DEM3 * x
     call sCubic(depth, range, det, scale)
     etgw = etgw * scale
-    trhs = -etgw
+    trhs = etgw
     det = -det * etgw
     etfunc_nlin = etgw
     !