求大神指教！GRADS转哈德来中心海温资料，画出图来很奇怪

zhaoxiaotong

用哈德来中心的海温资料，求1970-2011年7-10月的全球海温平均，求出的应该为42年，空间为360*180
可是画出的图很奇怪啊，
GRADS显示范围 。明显不对啊。
gs如下：
'reinit'
'sdfopen E:\HadISST_sst.nc'
'set fwrite E:\sst_7-10_4monave_globle.grd'
'set gxout fwrite'
'set x 1 360'
'set y 1 180'
'set z 1'

tt=1207
while(tt<=1704)
tt1=tt+3
'set t ' tt
'define aa=ave(sst,t='tt',t='tt1')'
'd aa'
tt=tt+12
endwhile

  'disable fwrite'
  'reinit'   
    ;

ctl文件：
dset E:\sst_7-10_4monave_globle.grd
undef -1e+30
xdef 360 linear -179.5 1
ydef 180 linear -89.5 1
zdef 1 levels 1000
tdef 42 linear JAN1970 1yr
vars 1
sst 0 99 ***
endvars


求大神帮忙啊~

lqouc

把缺省值改对了就行了

zhaoxiaotong

lqouc 发表于 2013-10-28 10:32

                               
登录/注册后可看大图

把缺省值改对了就行了

原本的缺省就是-1e+30.我就算进行了平均，缺省应该不会变的吧

lqouc

显然变了，否则怎么会出错，你改成-9.99e8试试。

river

缺测值问题，楼上说的如果不行，再试试-9.99e+33

ibeing

我以前也碰到过同样的情况，应该就是缺测值的问题

勇敢的悲伤

绝对是缺测值的问题。你先把数据变成TXT格式，看一看，把CTL文件的缺省值改了就可以了

jeffggg

可以肯定是缺测值的问题

zhaoxiaotong

lqouc 发表于 2013-10-28 11:22

                               
登录/注册后可看大图

显然变了，否则怎么会出错，你改成-9.99e8试试。

我换了都不行~最后发现是grads版本的问题。一开始用的整合版，后来换成grads1.9算出来就对的。至今不知道为什么。。。

superangus

楼主，请问你的哈德来海温资料是从哪里下载的呢？怎么会是nc格式,我下载是txt格式的，不知道怎么读