Returns the density of the ice in the shelf, which is assumed to be uniform across its domain.
| Type | Intent | Optional | Attributes | Name | ||
|---|---|---|---|---|---|---|
| class(ice_shelf), | intent(in) | :: | this | 
The ice density.
  pure function shelf_temperature(this) result(temperature)
    !* Author: Christopher MacMackin
    !  Date: April 2016
    !
    ! Returns the density of the ice in the shelf, which is assumed to be
    ! uniform across its domain.
    !
    class(ice_shelf), intent(in) :: this
    real(r8)                     :: temperature !! The ice density.
    temperature = -15.0_r8 !TODO: Will probably want to change this at some point.
#ifdef DEBUG
    call logger%debug('ice_shelf%temperature','Ice shelf has temperature '// &
                      trim(str(temperature)))
#endif
  end function shelf_temperature