关于Grads的fwrite

陈记祖

这是我的代码
'reinit'
'sdfopen E:\ncep\air.sig995.2009.nc'
i=1
while (i<=1460)
'set t 'i
'set lon 170 170'
'set lat 90 90'
'set gxout fwrite'
'set fwrite E:\air.dat'
i=i+1
endwhile
'disable fwrite'
结果出现
No hardcopy metafile open
All files closed;all defined objects released
all grads attributes have been reinitialized
但是运行后没有出现.dat文件，有没有人告诉我为什么，我是新手。

mofangbao

说明你最基础的手册都没有看过一遍，本板块有很多提取数据的代码，你自己找找原因吧，这个应该你自己解决

陈记祖

mofangbao 发表于 2012-7-8 09:45

                               
登录/注册后可看大图

说明你最基础的手册都没有看过一遍，本板块有很多提取数据的代码，你自己找找原因吧，这个应该你自己解决

我也是这几天采用，以前没学过，现在急用，您能帮我解决一下吗

言深深

我不懂grads，但是你是不是需要在endwhile前面，添加一句d air

陈记祖

言深深 发表于 2012-7-8 10:07

                               
登录/注册后可看大图

我不懂grads，但是你是不是需要在endwhile前面，添加一句d air

我刚加了d air，试了一下，出现了.dat文件，但是里面没有东西。{:eb303}

unice

感觉语句顺序有点别扭，不知道是不是这个原因，

陈记祖

unice 发表于 2012-7-8 10:37

                               
登录/注册后可看大图

感觉语句顺序有点别扭，不知道是不是这个原因，

我也不知道啊。是不是需要一个ctl说明文件，这个我真不懂，急死人了，

unice

如果我的话我习惯用这样的语序，楼主可以试一下，行不行告诉我哈
'reinit'
'sdfopen E:\ncep\air.sig995.2009.nc'
'set fwrite E:\air.dat'
'set gxout fwrite'
'set lon 170 170'
'set lat 90 90'
'set z 1'
i=1
while (i<=1460)
'set t 'i
'd air'
i=i+1
endwhile
'disable fwrite'
;
对了，好像还要设置一下层数，楼主试一下？

unice

nc文件是自描述文件，不需要ctl的，除非你这个是非标准的nc文件，好像是这样说的

做个霸气的木头

因为你要生成一个dat文件，那么i=i+1前面两行是放在外面的，另外如深哥所说，你必须有d air这行命令在循环里面才会有数据写到dat里面，gs语句最后一行应该加分号单独结尾，这是基础知识~