代码如下:
#计算风速
def fc(u,v):
n1=len(u[:,0])#行数,即纬度数
n2=len(u[0,:])#列数,即经度数
data=np.zeros((n1,n2))#创建与经纬度匹配的数组
#计算风速大小
for i in range(n1):
for j in range(n2):
data[i,j]=np.sqrt(u[i,j]**2+v[i,j]**2)#计算风速大小
return data
计算某一区域内各格点的风速大小。