run = steady;  start = 0 yrs
nx = 20;  delta_z = 200 m
temperature = 20 C;  fluid_density = 1 g/cm3;  bulk_density = 2.3 g/cm3

x_well(km)   0   50

strat 'Unit 1'
    t_dep = -1 m.y.;  X(rk1) = 100%
    thickness = 1 km;  water_depth = 0

strat 'Unit 2'
    t_dep = 0 yrs;  X(rk2) = 100%
    column  thickness(km)  water_depth(km)
    w(1)          2             -1
    w(2)          1              0

rock rk1
    bpor = 0;  A_perm = 0;  B_perm = $BP1 log_darcy;  p_kxkz = 10
rock rk2
    bpor = 0;  A_perm = 0;  B_perm = $BP2 log_darcy;  p_kxkz = 10
end_rock