added check for psiramp > 1 to make it absolute rather than fraction …

…of cell thickness
rniswon · Jun 27, 2024 · bb1751c · bb1751c
1 parent 1dc941c
commit bb1751c
Showing 1 changed file with 1 addition and 0 deletions.
diff --git a/MODFLOW-NWT/src/gwf2wel7_NWT.f b/MODFLOW-NWT/src/gwf2wel7_NWT.f
@@ -807,6 +807,7 @@ DOUBLE PRECISION FUNCTION smooth3(H,T,B,dQ)
       DOUBLE PRECISION, INTENT(OUT) :: dQ
       smooth3 = 0.0D0
       s = PSIRAMP
+      IF ( PSIRAMP > 1. ) s = PSIRAMP/(T-B)  !for Farag
       s = s*(T-B)   ! puming rate begins to be ramped down.
       x = (H-B)
       aa = -6.0d0/(s**3.0d0)