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