def METHOD_NAME(left, right, resolution):
    """
    Create a grid from the coordinates of two opposite corners, and also resolution.
    The grid coordinate arrays are also returned.
    """
    # coords are for cell centers
    x, y, z = [np.arange(left[dim] + resolution[dim] / 2., right[dim], resolution[dim])
               for dim in range(3)]
    return x, y, z, np.stack(np.meshgrid(x, y, z), axis=-1)