np meshgrid

1>>> nx, ny = (3, 2)
2>>> x = np.linspace(0, 1, nx)
3>>> y = np.linspace(0, 1, ny)
4>>> xv, yv = np.meshgrid(x, y)
5>>> xv
6array([[0. , 0.5, 1. ],
7       [0. , 0.5, 1. ]])
8>>> yv
9array([[0.,  0.,  0.],
10       [1.,  1.,  1.]])
11>>> xv, yv = np.meshgrid(x, y, sparse=True)  # make sparse output arrays
12>>> xv
13array([[0. ,  0.5,  1. ]])
14>>> yv
15array([[0.],
16       [1.]])
17