x = np.linspace(0, 2*np.pi, nx, endpoint=True) # X-points y = np.linspace(0, 2*np.pi, ny, endpoint=True) # Y-points dzdx = np.cos(X)*np.sin(Y) # partial derivative (x ...