ground_solve Subroutine

  subroutine ground_solve(this, ice_thickness, ice_density, &
                         ice_temperature, time, success)
    !* Author: Chris MacMackin
    !  Date: April 2017
    !
    ! Solves the state of the ground for the specified ice properties,
    ! at the specified time.
    !
    class(ground), intent(inout)    :: this
    class(scalar_field), intent(in) :: ice_thickness
      !! Thickness of the ice above the basal surface
    real(r8), intent(in)            :: ice_density
      !! The density of the ice above the basal surface, assumed uniform
    real(r8), intent(in)            :: ice_temperature
      !! The temperature of the ice above the basal surface, assumed uniform
    real(r8), intent(in)            :: time
      !! The time to which the basal surface should be solved
    logical, intent(out)            :: success
      !! True if the solver is successful, false otherwise
  end subroutine ground_solve