(已解决)nc资料转换成dat资料的时候出问题，求助

sosbbxc · 发表于 2012-11-15 20:09:05

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

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

x

本帖最后由 sosbbxc 于 2012-11-21 21:01 编辑

本人是个grads小菜鸟，才学习grads没多久，最近想把nc资料转换成dat资料，结果转出来的dat资料画出来的图和原先用nc资料画出来的图不一样。。研究了很长时间都没什么进展，所以在此求助。。。我用的nc资料是温度场，选取了500hPa的高度，转换的gs文件是这样的：
'reinit'
'sdfopen F:/air.mon.mean.nc'
'set gxout fwrite'
'set fwrite F:/air.dat'
it=7
while(it<=775)
'set t 'it''
'set lev 500'
'd air'
it=it+12
endwhile
'disable fwrite'
我想提取的是每年（1948年-2012年）7月份500hpa的温度数据。然后用nc资料画出来的图是这样的：

nc资料出的图

然后用提取的dat资料画出来的图是这样的：
dat出的图.jpg

两张图完全不一样。dat资料的ctl文件我是这样写的：
DSET f:\air.dat
title tezhengxiangliang
undef -9.99e+33
xdef 144 linear 0 2.5
ydef 73 linear -90 2.5
zdef 1 linear 500 1
tdef 65 linear 00Z01Jul1948 1yr
vars 1
air 0 99 air mon
endvars
很郁闷。不知道是哪里出问题了。。求高手指点。。感激不尽！！

Aires · 发表于 2012-11-15 20:26:56

我也不大懂你想做500hpa的 ctl中的变量air 后面应该是1吧

wc619lf · 发表于 2012-11-16 12:33:08

grads可以转成dat？试试grd格式吧

三江魂 · 发表于 2012-11-16 15:30:53

你资料没给全目前看应该是DAT的经纬度设置有问题

river · 发表于 2012-11-16 18:42:45

把set x 1 73和set y 1 144加上，还有你的ctl里面时间应该是间隔12个月，要写成12mo，要不然时间对不上了

做个霸气的木头 · 发表于 2012-11-16 23:01:43

看gs没什么问题，是不是你两张图的时间对应不上，或者你格点设置有问题

sosbbxc · 发表于 2012-11-17 21:58:34

river 发表于 2012-11-16 18:42

登录/注册后可看大图

把set x 1 73和set y 1 144加上，还有你的ctl里面时间应该是间隔12个月，要写成12mo，要不然时间对不上了

太谢谢你了。。按照你的方法问题就解决了。。没想到是这些微小的问题导致了我的错误。谢谢！

river · 发表于 2012-11-17 22:47:13

本帖最后由 river 于 2012-11-17 22:53 编辑

sosbbxc 发表于 2012-11-17 21:58

登录/注册后可看大图

太谢谢你了。。按照你的方法问题就解决了。。没想到是这些微小的问题导致了我的错误。谢谢！

解决就好

。使用set fwrite这个命令提取资料的时候最好是使用set x /set y这样的命令表示范围，使用set lat /set lon这样的命令很容易出现偏移。哦，还有，编辑一下帖子标题添加【已解决】吧，这样别人就不会来回答了。

shirly · 发表于 2012-11-22 18:29:08

学到新东西了，谢谢！！

wygohappy · 发表于 2013-4-9 10:19:15

受益匪浅谢谢~~！

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

[分享资料] (已解决)nc资料转换成dat资料的时候出问题，求助

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

相关帖子

点评