In [92]: grid = np.meshgrid(
    ...:     np.linspace(-1, 1, 5),
    ...:     np.linspace(-1, 1, 4), indexing='ij')
    ...: 
In [93]: grid
Out[93]: 
[array([[-1. , -1. , -1. , -1. ],
        [-0.5, -0.5, -0.5, -0.5],
        [ 0. ,  0. ,  0. ,  0. ],
        [ 0.5,  0.5,  0.5,  0.5],
        [ 1. ,  1. ,  1. ,  1. ]]),
 array([[-1.        , -0.33333333,  0.33333333,  1.        ],
        [-1.        , -0.33333333,  0.33333333,  1.        ],
        [-1.        , -0.33333333,  0.33333333,  1.        ],
        [-1.        , -0.33333333,  0.33333333,  1.        ],
        [-1.        , -0.33333333,  0.33333333,  1.        ]])]

grid是一个包含两个数组的列表。第一个数组的数字来自第一个参数（有5个元素的那个）。第二个是第二个参数中的数字。
为什么np.shape(grid)是（5，4）？你想要什么样的布局？
np.shape(grid)实际上是np.array(grid).shape，这就是为什么有一个额外的维度。