请选择 进入手机版 | 继续访问电脑版
爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

新浪微博登陆

只需一步, 快速开始

搜索
查看: 3601|回复: 5

[脚本编辑] Grads批处理nc文件

[复制链接]

新浪微博达人勋

发表于 2018-6-12 16:28:19 | 显示全部楼层 |阅读模式

登录后查看更多精彩内容~

您需要 登录 才可以下载或查看,没有帐号?立即注册 新浪微博登陆

x
本帖最后由 moerHuang 于 2018-6-12 16:28 编辑

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

*本程序提取(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

谢谢大家


gs.png
nc.png
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2018-6-12 19:51:01 | 显示全部楼层
用matlab的ncread系列函数搞吧~
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2018-6-13 08:30:14 | 显示全部楼层
在't5=t5+1'前面加一句命令'close 1'
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2018-6-13 14:20:06 | 显示全部楼层
把disable fwrite命令与最后一个end while换下位置
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2018-6-13 15:49:17 | 显示全部楼层
谢谢楼上两位,是这样解决的
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2018-6-13 15:50:43 | 显示全部楼层
Wetter 发表于 2018-6-12 19:51
用matlab的ncread系列函数搞吧~

比较习惯用grads和fortran处理,已经解决啦,谢谢~
密码修改失败请联系微信:mofangbao
您需要登录后才可以回帖 登录 | 立即注册 新浪微博登陆

本版积分规则

Copyright ©2011-2014 bbs.06climate.com All Rights Reserved.  Powered by Discuz! (京ICP-10201084)

本站信息均由会员发表,不代表气象家园立场,禁止在本站发表与国家法律相抵触言论

快速回复 返回顶部 返回列表