求时间序列

龙井

我用的资料是NECP上月平均温度

'reinit'
'sdfopen  D:\Thesis\download\month\air.nc'
'enable print D:\Thesis\download\month\tempe.gmf'
'set grid off'
'set grads off'
'set x 1'
'set y 1'
'set z 6'
'set t 373 816'
'tt=tloop(aave(air,lon=75,lon=100,lat=25,lat=40))'
'ttt=tloop(aave(air,lon=75,lon=120,lat=25,lat=40))'
'd tt-ttt'
'print'
'c'
'disable print'
'reinit'

balfulosa

在定义tt和ttt之前先求年平均，然后对年平均变量计算tt和ttt试试呢？

龙井

balfulosa 发表于 2016-5-1 23:11
在定义tt和ttt之前先求年平均，然后对年平均变量计算tt和ttt试试呢？

file:///C:/Users/Administrator/Desktop/QQ截图20160502082450
改为这样，但还是画不出来图，有错误，显示如图
i=1
while(i<=432)
     'qw=ave(air,t='i+372',t='i+384')'
     'tt=tloop(aave(qw,lon=75,lon=100,lat=25,lat=40))'
     i=i+12
endwhile

'set lon 75 120'
'set lat 25 40'
'set z 6'
*'set t 373 816'
i=1
while(i<=432)
     'wd=ave(air,t='i+372',t='i+384')'
     'ttt=tloop(aave(wd,lon=75,lon=120,lat=25,lat=40))'
     i=i+12
endwhile
     'd tt-ttt'

balfulosa

龙井 发表于 2016-5-2 08:27
改为这样，但还是画不出来图，有错误，显示如图
i=1
while(i

这样求是不行的，因为每一个循环tt只有一个值并不是时间序列。想了一下，还真没用grads算年平均，倒是可以但应该很复杂。建议先用FORTRAN处理一下然后在计算。

balfulosa

龙井 发表于 2016-5-2 08:27
改为这样，但还是画不出来图，有错误，显示如图
i=1
while(i

这样求是不行的，因为每一个循环tt只有一个值并不是时间序列。想了一下，还真没用grads算年平均，倒是可以但应该很复杂。建议先用FORTRAN处理一下然后在计算。