grads提取nc数据转化为dat格式的数据出现问题

云大小子

我的提取的gs文件是：
'reinit'
'sdfopen e:\nc\noaahaiyang\uwnd.mean.nc'
'set gxout fwrite'
'set fwrite e:\nc\noaahaiyang\uwnd.dat'

tt=169
while(tt<=180)
'set lon 55  167'
'set lat 5   55'
'set t 'tt''
'd uwnd'
say tt
tt=tt+1
endwhile
'disable fwrite'
'reinit'

然后用所得到的dat文件画图其中ctl文件是：
dset e:\nc\noaahaiyang\uwnd.dat
title ncep data
undef -9999
xdef 57 linear 55 2
ydef 26 linear 5 2
zdef 1  linear 1000   1
tdef   12 linear 00:00Z01JAN2005 1mo
vars 1
uwnd  1  99  data testing
endvars
用dat数据和nc数据画出的同一时刻 同一区域的图形不一样：
上面的是nc数据后面的是dat画出来的，好像只有0值的图形；

mofangbao

你 set gxout grid 填值看看转出来后的格点上的数据值如何

云大小子

mofangbao 发表于 2012-3-15 12:13

                               
登录/注册后可看大图

你 set gxout grid 填值看看转出来后的格点上的数据值如何

画出来是一片空白啊

云大小子

mofangbao 发表于 2012-3-15 12:13

                               
登录/注册后可看大图

你 set gxout grid 填值看看转出来后的格点上的数据值如何

我的gs文件写错了么，怎么提取出来的dat数据画图不合适啊，错在哪里啊

lotina

会不会是提取数据的那个GS错了，你在设置经度的后边在设置一下层次呢？'set lev 1000'
不知道你的NC资料是几层的!如果资料只是NC资料只是一层的话就应该没关系。
还有就是CTL中，层次不是非线性的，一般都是设置为  zdef 1 levels 1000，也不知道会不会是这问题。
一点小建议哈！

云大小子

lotina 发表于 2012-3-15 13:53

                               
登录/注册后可看大图

会不会是提取数据的那个GS错了，你在设置经度的后边在设置一下层次呢？'set lev 1000'
不知道你的NC资料是 ...

我按照你说的 尝试弄了一次 但是还是不行

lotina

云大小子 发表于 2012-3-15 14:22

                               
登录/注册后可看大图

我按照你说的 尝试弄了一次 但是还是不行

你的资料是一层？

godenflame135

sdf能打开的文件也是有规定的，你确定是规定的文件？
如果不是，还是fortran编程读取写成二进制网格点，再用grads画图吧！

godenflame135

sdf能打开的文件也是有规定的，你确定是规定的文件？
如果不是，还是fortran编程读取写成二进制网格点，再用grads画图吧！

云大小子

lotina 发表于 2012-3-15 14:27

                               
登录/注册后可看大图

你的资料是一层？

恩 一层数据