grads提取nc数据，gs运行出错

Rolle0524

我想用grads把nc数据的某个经纬度的范围内的数据（所要的数据是wvap，表示蒸发量。。。。。。这个nc数据只有四个参数，经度、纬度、时间和蒸发量）提取出来，输入为txt，在论坛上看了别人的帖子，自己改了下，编了一个gs文件。如下：'reinit'
'sdfopen E:\RDATA\wvap-morton_Amon_CNRM-CM5_historical_r1i1p1_native.nc'
i=1
while(i<=1872)
'set t' i
'set lon 89.30 92.11'
'set lat 29.42 32.22'
'set gxout fwrite'
'set fwrite E:\RDATA\wvap-morton_Amon_CNRM-CM5_historical_r1i1p1_native.txt'
'd wvap'
i=i+1
endwhile
'disable fwrite'

我把这个gs文件放在桌面上，打开grads2.0，输入run，然后把桌面上的gs拖到grads里面，结果运行界面都是红色字体，程序结束后是如图所示：
虽然在指定的路径下生成了txt文档，但是打开都都是乱码的。
对grads不算太熟，希望大家指点下。某不甚感激。

river

建议楼主去看一下grads实用手册

Rolle0524

手册浏览了一遍,不知道是哪出的问题,
真心急着用这个结果,
知道的话麻烦告知下吧
谢了

Rolle0524

river 发表于 2014-9-15 19:23
建议楼主去看一下grads实用手册

手册浏览了一遍,不知道是哪出的问题,
真心急着用这个结果,
知道的话麻烦告知下吧
谢了

一水天下

set gxout fwrite出来的是二进制文件，不能当文本打开的

Rolle0524

一水天下 发表于 2014-9-15 21:24
set gxout fwrite出来的是二进制文件，不能当文本打开的

那红丝字体提示的set fwrite error：fwrite is open是啥意思呢?
这个脚本文件写的对吗？

一水天下

应该写成'set t 'i，还得把set fwrite写在循环外。
具体的多看看手册吧，应该有例子

传说中的谁

槽点太多。。。。。。。首先fwrite是不能输出你所谓的txt文件的，其次你set fwrite 放在循环里面，disable fwrite放在循环外面

Rolle0524

传说中的谁 发表于 2014-9-15 21:52
槽点太多。。。。。。。首先fwrite是不能输出你所谓的txt文件的，其次你set fwrite 放在循环里面，disable  ...

我搞错了，生成的文件应该是dat格式的。
不过我感觉我的gs文件都是模仿你写的这个http://bbs.06climate.com/forum.php?mod=viewthread&tid=944来的，
除了我的脚本中没有高度值的设置之外，其他感觉都一样呢。
怎么回事呢

小淼在梦游№

'reinit'
'sdfopen E:\RDATA\wvap-morton_Amon_CNRM-CM5_historical_r1i1p1_native.nc'
'set gxout fwrite'
'set fwrite E:\RDATA\wvap-morton_Amon_CNRM-CM5_historical_r1i1p1_native.dat'
'set lon 89.30 92.11'
'set lat 29.42 32.22'
i=1
while(i<=1872)
'set t' i
'd wvap'
i=i+1
endwhile
'disable fwrite‘
；
fwrite应该写在循环外面吧