def METHOD_NAME(tmp_path):
    """Test square neighbourhooding with a mask"""
    kgo_dir = acc.kgo_root() / "nbhood/mask"
    kgo_path = kgo_dir / "kgo_masked.nc"
    input_path = kgo_dir / "input_masked.nc"
    output_path = tmp_path / "output.nc"
    args = [
        input_path,
        "--neighbourhood-output",
        "probabilities",
        "--neighbourhood-shape",
        "square",
        "--radii",
        "20000",
        "--output",
        output_path,
    ]
    run_cli(args)
    acc.compare(output_path, kgo_path)