关于fwrite的一个问题，想将HadIsst.nc转化成二进制，就是不成功啊亲！！

happyhui · 发表于 2012-5-27 12:50:09

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

您需要登录才可以下载或查看，没有帐号？立即注册

x

最近在用HadIsst做项目，想用fwrite语句把nc格式转化成二进制然后用于fortran计算，但是就是不成功啊！！！以下是我的gs文件：
'reinit'
'sdfopen f:\sam\areacc\hadleycenter\HadISST.nc'
'set fwrite f:\sam\areacc\hadleycenter\sst.dat'
'set gxout fwrite'
'set lon -179.5 179.5'
'set lat -89.5 89.5'
'set z 1'
'set t 1 1706'
'd sst'
'reinit'

最神奇的地方是，这个gs我已经用过无数次了以前读NCEP的2.5精度文件完全没有问题，但这次转成了dat之后，用fortran读出来所有数据都是0。写了ctl之后用GRADS直接打开画图检验也都全是0值。抓狂了！！请大家帮助一下！项目已经让我因为这个事儿拖了一个礼拜进程了！非常感谢！~~~~

末端 · 发表于 2012-5-27 20:38:54

难道不是应该先set gxout fwrite
然后再set fwrite吗？

易小凯 · 发表于 2012-5-27 13:03:43

disable fwrite

mofangbao · 发表于 2012-5-27 13:28:10

不同资料的分辨率不同，不要使用set lon lat的方法，用set x 这种试试

happyhui · 发表于 2012-5-27 15:29:54

mofangbao 发表于 2012-5-27 13:28

登录/注册后可看大图

不同资料的分辨率不同，不要使用set lon lat的方法，用set x 这种试试

不行啊，刚才改成了set x的方法，还是不好用。真是太奇怪了。。。

happyhui · 发表于 2012-5-27 15:30:27

易小凯发表于 2012-5-27 13:03

登录/注册后可看大图

disable fwrite

disable过了。。没有效果还是没有读入值

做个霸气的木头 · 发表于 2012-5-27 15:49:22

disable必须有，再者把时间设成外层循环试试，再不出数据我就怀疑你文件路径或者原数据文件有问题，用原数据画画图看看

易小凯 · 发表于 2012-5-27 16:48:50

happyhui 发表于 2012-5-27 15:30

登录/注册后可看大图

disable过了。。没有效果还是没有读入值

用原数据画图看看

mofangbao · 发表于 2012-5-27 21:50:54

先不要提取，把输出的部分注释掉改为画图

冰枫 · 发表于 2012-5-28 08:27:01

没有提示什么错误嘛？可能此nc不是grads能读取的格式

		自动登录	找回密码
密码			立即注册

[分享资料] 关于fwrite的一个问题，想将HadIsst.nc转化成二进制，就是不成功啊亲！！

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

相关帖子