- 积分
- 17
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2017-9-10
- 最后登录
- 1970-1-1
|
发表于 2019-12-21 20:13:58
|
显示全部楼层
ts1 =f->TS
latt =f->lat
lonn =f->lon
lsdata = a->LSMASK
lsm = landsea_mask(lsdata,latt,lonn)
ts1 = mask(ts1,lsm.eq.1,False)
ts1 = mask(ts1,lsm.eq.3,False)
;ts1 =mask(ts1,lsm.eq.0, True)
;*********************************PDO**********************
sst =ts1(:,{20:60},{90:270})
siw =month_to_season(sst,"DJF")
lat =f->lat({20:60})
lon =f->lon({90:270})
X = sst(lat|:,lon|:,time|:)
neof=4
eof = eofcov( X, neof)
eof_ts = eofcov_ts( X, eof)
PDO =eof_ts(0,:)
PDO =runave(PDO,21,0)
printVarSummary(PDO)
PDO =dim_standardize_Wrap(PDO,0)
;print(PDO) |
|