import numpy as np
R = 6371e3
dy = R * np.deg2rad(dlat)
for l in range(len(selected_sst.lon)):
sst_gradient[:-1, l] = (selected_sst[1:, l].values - selected_sst[:-1, l].values) / float(dy)
sst_gradient[-1, :] = (selected_sst[-1, :].values - selected_sst[-2, :].values) / float(-dy)
v=np.array(selected_sst.values, dtype=np.float)
sst_gradient, _ = np.gradient(v) #error in this line code here
我想找到海面温度梯度。