set c_dst

swensosc · Nov 28, 2023 · dbf6845 · dbf6845
1 parent 5491aae
commit dbf6845
Showing 1 changed file with 5 additions and 2 deletions.
diff --git a/src/biogeophys/SoilHydrologyMod.F90 b/src/biogeophys/SoilHydrologyMod.F90
@@ -2297,13 +2297,16 @@ subroutine SubsurfaceLateralFlow(bounds,  &
             ! should a warning be used instead?
             head_gradient = min(max(head_gradient,-2._r8),2._r8)
 
-            ! Calculate transmissivity of source column
+            ! Determine source and destination columns
             if (head_gradient >= 0._r8) then
                c_src = c
+               c_dst = col%cold(c)
             else
                c_src = col%cold(c)
+               c_dst = c
             endif
-
+
+            ! Calculate transmissivity of source column
             transmis = 0._r8
             if(c_src /= ispval) then 
                ! transmissivity non-zero only when saturated conditions exist