CALL PlumeSolver(Zi, Xi, Ta, Sa, Q0, ZOutput, MROutput)

SUBROUTINE PlumeSolver (Depth, Front, TempA, SalA, Discharge, DepthOutput, MeltOutput)

REAL(KIND=dp), ALLOCATABLE :: Depth(:), Front(:), TempA(:), SalA(:)
REAL(KIND=dp) :: Discharge
REAL(KIND=dp), POINTER :: DepthOutput(:), MeltOutput(:)

INTERFACE
  SUBROUTINE PlumeSolver
    USE TYPES
    REAL(KIND=dp), ALLOCATABLE :: Depth(:), Front(:), TempA(:), SalA(:)
    REAL(KIND=dp) :: Discharge
    REAL(KIND=dp), POINTER :: DepthOutput(:), MeltOutput(:)
  END SUBROUTINE PlumeSolver
END INTERFACE