anchor.diff:    end if
   
     return
  -end subroutine enforce_moderate_free_surface
  +end subroutine enforce_depth_thickness_consistency
   
   ! ---------------------------------------------------------------------------
   !> Ensure that layer heights do not fall below the prescribed minimum
  
positive.code: ! ---------------------------------------------------------------------------
  !> Check that the free surface anomaly and layer thicknesses are consistent
  
  subroutine enforce_moderate_free_surface(h, eta, depth, &
      freesurfFac, thickness_error, nx, ny, layers)
    implicit none
  
  
negative.code: 
      ! Check that the supplied free surface anomaly and layer
      ! thicknesses are consistent
      call enforce_depth_thickness_consistency(h, eta, depth, &
          freesurfFac, thickness_error, nx, ny, layers)
  
    end if
  