我的求水汽通量的gs文件运行时出了这个问题，请大家帮帮忙

469482462

我的求水汽通量的gs文件运行时出了这个问题，请大家帮帮忙，看看哪里有错吧
'reinit'
iyear=1981
while(iyear<=2010)
'sdfopen D:\lw2011\huanliu\uwnd.'iyear'.nc'
'sdfopen D:\lw2011\huanliu\vwnd.'iyear'.nc'
'sdfopen D:\lw2011\huanliu\rhum.'iyear'.nc'
'sdfopen D:\lw2011\huanliu\air.'iyear'.nc'
'set gxout fwrite'
'set fwrite D:\tmax\huanliu_daily\sj\q850-'iyear'.grd'
'set x 1 144'
'set y 1 73'
'set lev 850'
k=1
while(k<=365)
'set t ' k
'define prs=lev'
'define g=9.8'
'define es=(6.112*exp((17.67*(air.4-273.15))/(air.4-29.65)))'
'define qs=(0.62197*es/(prs-0.378*es))'
'define q=qs*rhum.3/100'
'define aa=mag(q*uwnd/g,q*vwnd.2/g)*1000'
k=k+1
endwhile

'define zz=ave(aa,time=30jun'iyear',time=14aug'iyear')'
'd zz '
'disable fwrite'
'close 1'
iyear=iyear+1
endwhile
'reinit'

才仔

看提示就是文件关闭错误，close 1那里，就是不知道怎么改

wet510

将close 1 去掉试试，或者加上close 2 close 3 close 4，这个命令关掉了打开的第一个文件，也就是纬向风的nc资料，如果要关闭资料的话应该都关闭。还有你已经disable fwrite 了，输出已停止。

sun_shine_Xia

你在gs文件里面一共打开了4个数据文件，顺序号分别是1,2,3,4！你在close的时候只能从4开始close！然后到最后close 1 ~

river

你这个就是把close的问题解决了也出来不了正确的东西。楼主把自己想要得到什么东西说清楚，因为这个涉及到set fwrite命令的位置

469482462

river 发表于 2013-8-27 19:38

                               
登录/注册后可看大图

你这个就是把close的问题解决了也出来不了正确的东西。楼主把自己想要得到什么东西说清楚，因为这个涉及到s ...

我想把资料处理成每年6月30号到8月14号的平均水汽输送通量，求指导

river

469482462 发表于 2013-9-2 15:47

                               
登录/注册后可看大图

我想把资料处理成每年6月30号到8月14号的平均水汽输送通量，求指导

在'define zz=ave(aa,time=30jun'iyear',time=14aug'iyear')'前面加上'set t 1'；在'close 1'前面加上'close 3' 'close 2'