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"
```

Thanks!

Joe