diff --git a/MODFLOW-NWT/src/gwf2wel7_NWT.f b/MODFLOW-NWT/src/gwf2wel7_NWT.f
index 163c046..30570a1 100644
--- 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)