run = steady; start = 0 yrs
nx = 30;  delta_z = 100 m
temperature = full; surface_temp = 20 C

rock ss
    A_perm = 0;  B_perm = -2 log_darcy;  p_kxkz = 10
    phi0 = 30%; phi1 = 0; bpor = 0
end_rock

X(ss) = 1
x_well(km)  0 10 20

strat 'Basal unit'
    t_dep = -1 m.y.
    thickness = 2 km

strat 'Topography'
    t_dep = 0 yrs
    column    water_depth(m)  thickness(m)  heat_flow
    w(1)        $TOPO           $THICK           4
    w(2)            0                0           2
    w(3)            0                0           1