- 积分
- 3711
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2018-6-12
- 最后登录
- 1970-1-1
|

楼主 |
发表于 2020-4-10 12:56:59
|
显示全部楼层
那个程序我后来没用了,下面这个是我用的,就是存储顺序跟你不太一样,我是把每一侯的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'
; |
|