I would like to suggest that the code for area calculation be changed from:
R = 6.37e6
# we know already that the spacing of the points is one degree latitude
dϕ = np.deg2rad(1.)
dλ = np.deg2rad(1.)
dA = R**2 * dϕ * dλ * np.cos(np.deg2rad(ds.lat))
dA.plot()
I would like to suggest that the code for area calculation be changed from:
to
as the data seems to suggest that it 2 degree grid