Grads批处理nc文件

moerHuang

请各位大神指导下，我取了1979-2009年1月降水资料（.nc），想用grads批处理成二进制文件方便后续fortran的读取。我写了个gs文件，但grads一直报错，萌新实在不会改了，想求助下大家是哪里的问题。

*本程序提取（1979-2009年）1月nc格式的中国南方降水值
'reinit'
t5=1979

while(t5<=2009)
'set gxout fwrite'
'set fwrite F:\PT\Pdat\'%t5%P'.dat'
'sdfopen F:\PT\P\'%t5%P'.nc'

'set lon 110 125'
'set lat 25 35'
'set lev 0 0'

tt=1
while(tt<=62)
'set t 'tt''
'd tp'
tt=tt+1
endwhile

t5=t5+1

endwhile
'disable fwrite'
;




报错：All files closed; all defined objects released;
All GrADS attributes have been reinitialized
SET FWrite Error:  fwrite file is open
Use DISABLE FWRITE command to close file

谢谢大家

Wetter

用matlab的ncread系列函数搞吧~

oucdaqi

在't5=t5+1'前面加一句命令'close 1'

wet510

把disable fwrite命令与最后一个end while换下位置

moerHuang

谢谢楼上两位，是这样解决的

moerHuang

Wetter 发表于 2018-6-12 19:51
用matlab的ncread系列函数搞吧~

比较习惯用grads和fortran处理，已经解决啦，谢谢~