• Solver Fortran File: SurfaceBoundaryEnthalpy.f90
• Solver Name: SurfBoundarySolver
• Required Output Variable(s): Surf Enth, Mass Balance, Densi, Firn, Melting, Refreeze, Accu, Rad_Fact, Rain, PotRad
• Required Input Variable(s): Depth and SurfGrad from the FlowDepth Solver
• Input Data: Daily air temperature timeserie
• Optional Input Data: Daily precipitation timeserie

SurfBoundarySolver is a pseudo-solver which compute surface mass balance and Dirichlet surface boundary condition for the Enthalpy solver. It takes into account firn heating processes by solving vertical melt-water percolation and refreezing.

The solver uses the provided air temperature (and precipitation) daily record to compute the associated mean surface characteristic of the glacier over the time period covered by the provided data time-serie. It outputs the following variables:

1. Mass Balance (m w.eq. yr^-1): Mean surface mass balance
2. Surf Enth (J kg^{-1}) : Enthalpy value bellow active layer. Can be use as a Dirichlet condition in the Enthalpy Solver
3. Densi (kg m^{-3}): Density field in 3D
4. Firn (m w. eq.) : Firn thickness
5. Melting (m w.eq. yr^{-1}) : Surface melting
6. Refreeze (m w.eq. yr^{-1}) : Amount of refreezing (superimposed ice)
7. Accu (m w.eq. yr^{-1}) : Snow accumulation
8. Rad_fact (m w.eq. (W m^{-2})^{-1}) : Melting factor for radiation
9. Rain (m w.eq. yr^{-1}) : Amount of rain
10. PotRad (W m-^{-2}) : Potential solar radiation