import numpy as np
import math
a=6.37e6
s=0.5
n=5
l=6
pi=3.1416
d=1e5
d=np.loadtxt('D:\实验数据\实验二\ddd.txt',dtype=int)
h=np.loadtxt('D:\实验数据\实验二\HHH.txt',dtype=int)
f=np.loadtxt('D:\实验数据\实验二\FF.txt',dtype=int)
f= np.concatenate(f, axis=0)
#风场订正
up=np.zeros((n,l),int)
vp=np.zeros((n,l),int)
for i in range(0,n):
for j in range(0,l):
up[i][j]=f*np.sin(d*pi/180-pi)
vp[i][j]=f*np.cos(d*pi/180-pi)
print('风速u分量\n',up)
print('风速v分量\n',vp)
f.shape = (30, )
np.sin(d * pi / 180 - pi).shape = (5, 6)
建议看一下是否忘了对d进行下标取值,如d[i][j]