Function iterate_noise_column_with_aquifer 

pub(crate) fn iterate_noise_column_with_aquifer<N: DimensionNoises>(
    cache: &mut N::ColumnCache,
    noises: &N,
    aquifer: &mut Aquifer<N>,
    block_x: i32,
    block_z: i32,
    ocean_floor: bool,
) -> i32

Matches vanilla’s iterateNoiseColumn: iterates by Y cells, evaluating inner density functions at 8 cell corners, trilinearly interpolating each channel independently, then applying outer operations (squeeze, min, etc.) per-block via combine_interpolated.

Returns getBaseHeight (= getFirstFreeHeight = first Y above surface).