Function threshold_gradient 

fn threshold_gradient(
    field: &[f64],
    grid: &Grid,
    axis: usize,
    threshold: f64,
) -> Vec<f64>

1.0 where |∂_axis field| < threshold, else 0.0.