Heat Flow and Bottom Boundary Condition

The bottom boundary condition is set by lbound and TDEEP.

Note

By default, the interfaces sets a fixed bottom temperature:

lbound = -1.//180.

Heat Flow

A geothermal heat flux can be set using lbound alone, with the heat flux given in mW/m2:

lbound = 35.

sets a heat flux of 35 mW/m2.

Insulation

 lbound = 0

insulates the bottom layer (no flux)

Set Bottom Temperature

The bottom temperature must be set first by TDEEP, then a flag is provided through lbound:

TDEEP = 180. #Bottom temperature set at 180.K
lbound = -1//TDEEP

Set the Initial Temperature of all the Layers

TDEEP = 180. #Bottom temperature set at 180.K
lbound = -2//TDEEP

All the layers are initially set at TDEEP, and the bottom layer is kept at TDEEP during the run.

Equilibrium

When selecting a lower boundary condition, make sure to verify that the model has been running long enough vs. depth enough vs. at adequate temperatures.