侯平均计算

软壳

onenuister 发表于 2020-4-9 11:08
tt=152
while(tt

出来的图不是一个值。
tt=151
while(tt<=240)
'define pre=ave(precip,t='tt',t='tt+4',1)'
'd pre'
tt=tt+5
endwhile
改成这样可以吗

软壳

onenuister 发表于 2020-4-9 11:08
tt=152
while(tt

楼主，想问你一个问题，你这样多年循环出来的数据，跟你用同样方法处理单独一年的数据，一样吗？我现在41年，年循环出来的图，只有第一年（1979年）的数据和单独处理1979年的图是一样的，其他年份的数据都不一样。

onenuister

软壳 发表于 2020-4-9 22:09
楼主，想问你一个问题，你这样多年循环出来的数据，跟你用同样方法处理单独一年的数据，一样吗？我现在41 ...

那个程序我后来没用了，下面这个是我用的，就是存储顺序跟你不太一样，我是把每一侯的40年存完再存下一侯，这样计算出来的数据跟同样方法处理单独一侯的数据是一样的
***1979-2018年7月第1侯-第6侯ESWJuwnd
'reinit'
'set fwrite f:\002\04\ESWJ_uwnd1979-2018.grd'
'set gxout fwrite'
it=183
while(it<=212)
year=1979
while(year<=2018)
'sdfopen f:\002\uwnd.'year'.nc'
'set lat -90 90'
'set lon 0 357.5'
'set lev 300'
***闰年
if(((math_fmod(year,4)=0)&(math_fmod(year,100)!=0))|(math_fmod(year,400)=0))
if(it=208)
'define u1=ave(uwnd,t='it',t='it+5')'
else
'define u1=ave(uwnd,t='it',t='it+4')'
endif
'define u2=aave(u1,lon=90,lon=110,lat=35,lat=40)'
'd u2'
else
***平年
if(it=208)
'define u1=ave(uwnd,t='it-1',t='it+4')'
else
'define u1=ave(uwnd,t='it-1',t='it+3')'
endif
'define u2=aave(u1,lon=90,lon=110,lat=35,lat=40)'
'd u2'
endif
year=year+1
'close 1'
endwhile
it=it+5
endwhile
'disable fwrite'
;

软壳

onenuister 发表于 2020-4-10 12:56
那个程序我后来没用了，下面这个是我用的，就是存储顺序跟你不太一样，我是把每一侯的40年存完再存下一侯 ...

我的问题也解决了，打扰了