- 积分
- 4366
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2021-6-9
- 最后登录
- 1970-1-1
|
发表于 2021-10-2 23:38:12
|
显示全部楼层
*原始时间序列
'enable print d:\statistics\result2.gmf'
'set lat -5'
'set lon 190'
'set lev 1000'
'set t 1 30'
'define staave=aave(st,lon=190,lon=240,lat=-5,lat=5)'
'd tloop(staave)'
'draw title Nino3.4 SST time series'
'print'
'c'
*距平时间序列
'set lat -5'
'set lon 190'
'set lev 1000'
'set t 1 30'
'define staave1=ave(staave,t=1,t=30)'
'define staano=staave-staave1'
'd tloop(staano)'
'draw title Nino3.4 SST anomaly time series'
'print'
'c'
标准化时间序列
'set lat -5'
'set lon 190'
'set lev 1000'
'set t 1 30'
'define stamse=sqrt(ave(pow(staano,2),t=1,t=30))'
'define stz=staano/stamse'
'd tloop(stz)'
'draw title Nino3.4 SST standardization time series'
'print'
'disable print'
;
|
|