back out two answer changers

swensosc · Sep 11, 2023 · 3034da1 · 3034da1
1 parent 3f4b9a8
commit 3034da1
Show file tree

Hide file tree

Showing 2 changed files with 8 additions and 8 deletions.
diff --git a/src/biogeophys/SoilHydrologyMod.F90 b/src/biogeophys/SoilHydrologyMod.F90
@@ -2447,7 +2447,7 @@ subroutine SubsurfaceLateralFlow(bounds,  &
             do j = jwt(c)+1,1,-1
 
                ! ensure water is not added to frozen layers
-               if (zi(c,j) < frost_table(c)) then 
+!scs               if (zi(c,j) < frost_table(c)) then 
                   ! analytical expression for specific yield
                   s_y = watsat(c,j) &
                        * ( 1. - (1.+1.e3*zwt(c)/sucsat(c,j))**(-1./bsw(c,j)))
@@ -2466,7 +2466,7 @@ subroutine SubsurfaceLateralFlow(bounds,  &
                   else
                      zwt(c) = zi(c,j-1)
                   endif
-               endif
+!scs               endif
 
             enddo
 

diff --git a/src/biogeophys/SoilWaterMovementMod.F90 b/src/biogeophys/SoilWaterMovementMod.F90
@@ -1395,12 +1395,12 @@ subroutine soilwater_moisture_form(bounds, num_hydrologyc, &
          !  save number of adaptive substeps used during time step
          nsubsteps(c) = nsubstep
 
-         ! check for over-saturated layers                                      
-         do j = nlayers,2,-1
-            over_saturation   = max(h2osoi_liq(c,j)-(watsat(c,j)*m_to_mm*dz(c,j)),0._r8)
-            h2osoi_liq(c,j)   = min(watsat(c,j)*m_to_mm*dz(c,j), h2osoi_liq(c,j))
-            h2osoi_liq(c,j-1) = h2osoi_liq(c,j-1) + over_saturation
-         end do
+!!$         ! check for over-saturated layers                                      
+!!$         do j = nlayers,2,-1
+!!$            over_saturation   = max(h2osoi_liq(c,j)-(watsat(c,j)*m_to_mm*dz(c,j)),0._r8)
+!!$            h2osoi_liq(c,j)   = min(watsat(c,j)*m_to_mm*dz(c,j), h2osoi_liq(c,j))
+!!$            h2osoi_liq(c,j-1) = h2osoi_liq(c,j-1) + over_saturation
+!!$         end do
 
          ! check for negative moisture values
          do j = 2, nlayers