nx = 1;  flow = vertical;  plot = final
density = 1 g/cm3;  bulk_density = 2.3 g/cm3

rock sh
    A_perm = $AP;  B_perm = $BP log_darcy;  p_kxkz = 10
end_rock

X(sh) = 1;  reference = uncompacted;  heat_flow = 1 HFU

strat 'Basal unit'
    t_dep = $TIME0 m.y.
    thickness = 100 m

strat 'Shale deposited over run'
    t_dep = 0 yrs
    thickness = 10 km