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

zhaoxiaotong · 发表于 2013-10-28 10:25:34

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

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

x

用哈德来中心的海温资料，求1970-2011年7-10月的全球海温平均，求出的应该为42年，空间为360*180
可是画出的图很奇怪啊， QQ截图20131028101908.jpg

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 · 发表于 2013-10-28 10:32:04

把缺省值改对了就行了

zhaoxiaotong · 发表于 2013-10-28 11:05:09

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

登录/注册后可看大图

把缺省值改对了就行了

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

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

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

river · 发表于 2013-10-28 12:09:03

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

ibeing · 发表于 2013-10-28 13:09:30

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

勇敢的悲伤 · 发表于 2013-10-28 15:43:22

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

jeffggg · 发表于 2013-10-28 18:41:13

可以肯定是缺测值的问题

zhaoxiaotong · 发表于 2013-10-28 21:34:56

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

登录/注册后可看大图

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

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

superangus · 发表于 2015-10-30 20:44:46

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

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

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

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