grads 批处理 输出文件问题

nickbsb

我想把nc资料转换成二进制grd的格式，我想提取hgt.1980.nc的数据直接转换输出成hgt.1980.grd的文件，然后提取hgt.1981.nc的数据直接转换输出成hgt.1981.grd的文件，以此类推，下文就是我写的，但是不对，请大神如何修改成我想要的。。

tt=1980
while(tt<=2013)
'set gxout fwrite'
'set fwrite I:\4_times_data\hgt_grd\hgt_'tt'.grd'
'sdfopen I:\4_times_data\hgt\hgt.'tt'.nc'
if(tt=1980|tt=1984|tt=1988|tt=1992|tt=1996|tt=2000|tt=2004|tt=2008|tt=2012)
tt0=1
while(tt0<=1464)
'set t 'tt0''
'set lat -90 90'
'set lon 0 357.5'
'd hgt'
tt0=tt0+1
endwhile
else
tt0=1
while(tt0<=1464)
'set t 'tt0''
'set lat -90 90'
'set lon 0 357.5'
'd hgt'
tt0=tt0+1
endwhile
endif
'close 1'
tt=tt+1
endwhile
'disable fwrite'
;

mofangbao

自己先检查  为什么觉得不对？

nickbsb

mofangbao 发表于 2014-2-26 12:03

                               
登录/注册后可看大图

自己先检查  为什么觉得不对？

请看我上传的图片

冰枫

你应该在循环中把fwrite文件也close。

mofangbao

至少有个错误是disable的放到循环里面去  

陌上花开123456

不知道你有木有解决，我这边有个gs是直接提取的，你可以参考看看

刷牙

我改了一下，可以运行，不知道对不对，你自己检查一下
'reinit'
tt=1982


while(tt<=2001)
'set fwrite D:\study\2012\compare\ncep\hgt_'tt'.grd'
'sdfopen D:\study\2012\compare\ncep\pres.tropp.'tt'.nc'
'set gxout fwrite'
'set lat 15 55'
'set lon 70 160'

if(tt=1984|tt=1988|tt=1992|tt=1996|tt=2000)
tt0=1
while(tt0<=1464)
'set t 'tt0''
'd pres'
tt0=tt0+1
endwhile

else
tt0=1
while(tt0<=1460)
'set t 'tt0''
'd pres'
tt0=tt0+1
endwhile
endif
'disable fwrite'
'reinit'
tt=tt+1
endwhile
;

水草

楼主后来这个问题怎么解决的呢？

Уēｓ、尜玉

请问楼主问题解决了没？？

凡不凡

我也碰到这样的错误