Water table

Fillable pore volume, retained storage, and water-table rise.

Water-table rise approximations based on mass balance.

cumulative_simpson_midpoint

cumulative_simpson_midpoint(y: ndarray, y_midpoint: ndarray, x: ndarray) -> np.ndarray

Integrate values using Simpson's rule and explicit interval midpoints.

fillable_volume_curve(z: ndarray, theta: ndarray, *, theta_s: float) -> np.ndarray

Compute fillable pore volume from each depth to the column base.

V_f(z) = integral_z^z_w0 [theta_s - theta(xi)] dxi

water_table_rise_from_volume(available_volume: float, z: ndarray, fillable_volume: ndarray, *, initial_water_table_depth: float) -> float

Estimate water-table rise by inverting fillable pore volume.