I'm seeking some advice on linking the residual of the limited heat equation to subglacial hydrology at the bed.

My model setup has a Dirichlet temperature condition at the bed (Pressure Melting Point) and, if I understand correctly, the heat sinks and sources imposed in satisfying this Dirichlet condition should be directly related to the melting/refreezing at the bed by the value of the latent heat of fusion.

From reading the documentation, I understand that

- Code: Select all
`Calculate Loads = Logical True`

can be set, which makes available the "Temperature Loads" variable.

My plan would be then to have the source term set such that:

- Code: Select all

$LHF = latentheatoffusion

Water Flux = Variable Temperature Loads

Real MATC "tx/$LHF"

Can anyone advise me if this will achieve the desired result or if I've missed something?

Thanks!

Joe