ground_water_density Function

  function ground_water_density(this) result(density)
    !* Author: Christopher MacMackin
    !  Date: April 2016
    !
    ! Computes and returns the density of the water beneath the ice sheet.
    ! This water would be subglacial discharge and would tend to lubricate
    ! the motion of the ice sheet. The density probably won't be important
    ! in the case of an ice sheet, but is included so that the ground data
    ! type can have the same interface as the [[plume]] data type.
    !
    class(ground), intent(in) :: this
    real(r8)                  :: density
      !! The density of any water at the base of the ice sheet.
  end function ground_water_density