def METHOD_NAME(data, region, spacing, expected_grid):
    """
    Run surface by passing in x, y, z numpy.ndarrays individually.
    """
    output = surface(
        x=data.x,
        y=data.y,
        z=data.z,
        spacing=spacing,
        region=region,
        verbose="e",  # Suppress warnings for IEEE 754 rounding
    )
    check_values(output, expected_grid)