利用GrADS将nc文件转化成为grd文件出现的错误

Daisy°L

这学期刚刚开始学习GrADS 进行上机实习过程中出现了问题 麻烦大家帮忙看下
所用资料是起始时间为1854年1月开始的海温资料
要用GrADS转化其中的纬度在-20-60ºN，经度在90ºE-200Eº，时间在1971年1月至2000年12月的资料
编写的gs文件如下
'renint'
'sdfopen f:\data\sst.mnmean.nc'
'set gxout fwrite'
'set fwrite f:\data\sst.mnmean.grd'
'set lon 90 200'
'set lat -20 60'
'set lev 0'
it=1405
while(it<=1764)
'set t 'it''
'd sst'
it=it+1
endwhile
'disable fwrite'
然后编写的ctl文件如下
dset f:\data\sst.mnmean.grd
undef -9.96921e+36
xdef 56 linear 90 2
ydef 41 linear -20 2
zdef 1 levels 0
tdef 360 linear JAN1971 1mo
vars 1
sst 0 -999 Monthly Means of Sea Surface Temperature
endvars
感觉没有什么问题，但是画出来的图和直接用nc画出来的图不一样，麻烦大家看一下啦

river

怎么不一样了，得有个对比吧。还有你那ctl里的缺测值怎么来的？明显不太正常

伊伊

图怎么不一样啊

LOSER

看下缺测值对不对

Daisy°L

river 发表于 2013-12-5 22:19

                               
登录/注册后可看大图

怎么不一样了，得有个对比吧。还有你那ctl里的缺测值怎么来的？明显不太正常

是这样的，一开始我写的是-9999，但是图画不出来 问老师 老师说缺测值不对 让我在GrADS里面打开nc文件 再输入q  ctlinfo 然后出来缺测值是那个 我就改成那个了

Daisy°L

伊伊 发表于 2013-12-5 22:32

                               
登录/注册后可看大图

图怎么不一样啊

那个 我晚上回去看能不能把图贴出来 准备考研呢 最近有些忙

mofangbao

输出之后缺测值都是 -9.99e+08
还有不要写
'set lon 90 200'
'set lat -20 60'
请改为对应的x y  否则容易出现偏移  除非你非常确定这些点刚好落在网格点上

LOSER

mofangbao 发表于 2013-12-6 09:03

                               
登录/注册后可看大图

输出之后缺测值都是 -9.99e+08
还有不要写
'set lon 90 200'

清风大神
sst.mnmean.nc 的 undef  -9.99e+33
hgt.mon.mean.nc 的 undef  -9.96921e+36
为什么您说输出之后就是-9.99e+08？有点不懂啊~

mofangbao

LOSER 发表于 2013-12-6 09:40

                               
登录/注册后可看大图

清风大神
sst.mnmean.nc 的 undef  -9.99e+33
hgt.mon.mean.nc 的 undef  -9.96921e+36

你可以set gxout fwrite　输出后用fortran打开看看的 一般都会变为 -9.99e+08  

mofangbao

Daisy°L 发表于 2013-12-6 08:59

                               
登录/注册后可看大图

那个 我晚上回去看能不能把图贴出来 准备考研呢 最近有些忙

额。。。别人是在帮助你呢  别人及时回答你了你却说可能还不能及时交流。。。太忙的话还是以考研为重  有空了再来发帖提问比较好的~