constructor Function

  pure function constructor(beta, melt_conversion, salinity_denom) result(this)
    real(r8), intent(in) :: beta
      !! The inverse stefan number, $$ \frac{c(T_a - T_m)}{L} $$
    real(r8), intent(in) :: melt_conversion
      !! The factor to convert between the scale for melt used by
      !! Dallaston et al. (2015) and that used in ISOFT, $$
      !! \frac{m_0x_0}{D_0U_0}, $$ where \(m_0\) is the melt scale
      !! used by Dalalston et al.
    real(r8), intent(in), optional :: salinity_denom
      !! The factor which, when used to divide the `melt_conversion`
      !! term, produces a conversion factor for the melt-terms in the
      !! salinity equation of Dallaston et al. (2015). It has the form
      !! $$ \frac{Q_{g0}}{D_0U_0}, $$ where \(Q_{g0}\) is the
      !! subglacial discharge across the grounding line.
    type(dallaston2015_melt) :: this
      !! The newly created object representing the melt relationship.
    this%coef = (beta + 1.0_r8)
    this%melt_conversion = melt_conversion
    if (present(salinity_denom)) this%salinity_denom = salinity_denom
  end function constructor