### Linking subglacial hydrology to Heat Equation

Posted:

**20 May 2013, 14:28**Hi all,

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

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

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

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

Thanks!

Joe

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