python经向温度梯度

975474549

求教大佬python经向温度梯度的程序怎么写
我自己写了一个，但是算出来的结果不对
#梯度计算(94,192)
r0=6371
r=np.zeros((94))
for i in range(94):
    r[i]=r0*(np.sin((90-lat[i])*(np.pi/180)))

graddsat=np.zeros((40,94,192))
for y in range(40):
    for lati in range(94):
        for loni in range(191):
            graddsat[y,lati,loni]=(dsat[y,lati,loni+1]-dsat[y,lati,loni])/((2*np.pi*r[lati])/192)
        graddsat[y,lati,191]=(dsat[y,lati,0]-dsat[y,lati,191])/((2*np.pi*r[lati])/192)