• Solver Fortran File: GlaDSCoupledSolver.F90 and GlaDSchannelSolver.F90
• Solver Name: GlaDSCoupled, GlaDSsheetThickDummy and GlaDSchannelOut
• Required Output Variable(s): Hydraulic Potential, Sheet Thickness and Channel Area
• Required Input Variable(s): None
• Optional Output Variable(s): Vclose, Wopen, Water Pressure, Effective Pressure, Sheet Discharge, Sheet Storage
• Optional Input Variable(s): Zb

The complete description of the equations solved by the GlaDS solver can be found in Werder et al. (2013). The implementation follows exactly these equations, except that optionally the hydraulic potential can be computed at the top of the water sheet instead than at the bed (keyword: Neglect Sheet Thickness in Potential).

The GlaDS solver solves for the hydraulic potential, the water sheet thickness and the cross-sectional area of the channels. Whereas the two first variables are nodal variable and define continuous fields, the Channel area is a discrete field only defined on the edge of the elements.

The GlaDS model is composed of three solvers:

The required keywords in the SIF file for the Porous Solver are:

! Define some useful parameters using MATC
$yearinsec = 365.25*24*60*60

An example using the GlaDS Solver can be found in [ELMER_TRUNK]/elmerice/examples/GlaDS.

The description of the GlaDS model is in:
Werder M.A., I.J. Hewitt, C.G. Schoof and G.E. Flowers, 2013. Modeling channelized and distributed subglacial drainage in two dimensions. Journal of Geophysical Research: Earth Surface, 118(4), 2140-2158.