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