- 积分
- 108
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2023-12-11
- 最后登录
- 1970-1-1
|
登录后查看更多精彩内容~
您需要 登录 才可以下载或查看,没有帐号?立即注册
x
我想用三个月滑动平均大于0.5,且维持5个月以上来判断厄尔尼诺年,但是我用这个条件算出来的厄尔尼诺年与有记录的厄尔尼诺年对不上。比如1976年,它是一个强厄尔尼诺年。但是我算出来的三个月滑动平均只有最后两个值是大于0.5的,其余值都小于0.5,都不能判定为一个厄尔尼诺年。想请问大家这是为什么啊,我该怎么改呢
#调取EINO3.4区的数据
sst34=f1['sst'].loc['1960-01-01':'2022-12-01'].loc[:,5:-5,190:240]
lon=f1['lon'].loc[190:240]
lat=f1['lat'].loc[5:-5]
#time=f1.variables['time']
#print(f1)
#print(sst34)#(744,5,26)
#将数据转换成四维数组(year,month,lat,lon)对年求平均场及距平场
sst34=np.array(sst34).reshape(63,12,5,26)
sst34_1=np.mean(np.mean(sst34,2),2)#平均场 (63,12)
sst34_2=np.mean(sst34_1,0) #(12)
sst34_jp=np.zeros([63,12])
for m in range(63):
for n in range(12):
sst34_jp[m,n]=sst34_1[m,n]-sst34_2[n]
#找出三个月滑动平均大于0.5,且维持5个月及以上的年份
sst34_move=np.zeros([63,12])
for m in range(63):
sst34_move[m,:]=gaussian_filter1d(sst34_jp[m,:], 3)
|
|