关于set gxout stat 问题

冬日情愫+

看到发贴人： 小金鱼  的帖子说 set gxout stat 命令可以查看数据的缺省值。看了grads的官网，上面说set gxout stat 是print statistical output to the terminal ，我用了这个命令（数据是2.5*2.5的nc文件）后出现Undef value = 32766 但是在q ctlinfo 命令查看数据描述文件时 undef -9.96921e+36 。到底应该是哪个呢？
ga-> q ctlinfo
dset e:/data/grads/air.1998.nc
title 4x daily NMC reanalysis (1998)
undef -9.96921e+36
xdef 144 linear 0 2.5
ydef 73 linear -90 2.5
zdef 17 levels 1000 925 850 700 600 500 400 300
250 200 150 100 70 50 30 20 10
tdef 1072 linear 00Z01JAN1998 360mn

ga-> set gxout stat
ga-> d air
Data Type = grid
Dimensions = 0 1
I Dimension = 1 to 145 Linear 0 2.5
J Dimension = 1 to 73 Linear -90 2.5
Sizes = 145 73 10585
Undef value = 32766
Undef count = 0  Valid count = 10585
Min, Max = 236 310
Cmin, cmax, cint = 240 310 5
Stats[sum,sumsqr,root(sumsqr),n]:     2.95504e+06 8.28165e+08 28777.9 10585
Stats[(sum,sumsqr,root(sumsqr))/n]:     279.172 78239.5 279.713
Stats[(sum,sumsqr,root(sumsqr))/(n-1)]: 279.199 78246.9 279.727
Stats[(sigma,var)(n)]:     17.3876 302.328
Stats[(sigma,var)(n-1)]:   17.3884 302.357
另外我在grads官网上没有找到关于Stats[sum,sumsqr,root(sumsqr),n]   和Stats[(sigma,var)(n)] 的详细介绍，不知道它们代表什么意思，其中的变量sum,sumsqr,root(sumsqr),sigma有是什么意思？

mofangbao

这个我也不清楚，不知道有没有朋友能解释，我分别用meteoinfo和grads打开nc文件，看到的undef值也不一样，确实很玄乎，也可能是我们没有理解其中的本质吧，继续探求中...

mofangbao

不知道是不是数据类型的长短造成的

冬日情愫+

我也不知道到底是什么原因。也在思索中。。。。

兰溪之水

@mofangbao我遇到过这个问题，用set gxout stat查看的那个是对的，至于q ctlinfo忘了哪里有简单的解释，反正我实践过，转换后应该设置为set gxout stat查看的那个缺省值~

mofangbao

兰溪之水 发表于 2011-11-19 20:22

                               
登录/注册后可看大图

@mofangbao我遇到过这个问题，用set gxout stat查看的那个是对的，至于q ctlinfo忘了哪里有简单的解释，反正 ...

下次涉及到的时候我也试试你说的

MeteoInfo

ctl文件是人为写的，写多少都可以
事实上nc文件中每个变量的missing value都可以设为不同，这是个自由国度的人弄的自由数据格式{:soso_e112:}

冬日情愫+

但是缺省值是不能随便定的吧，在数据计算的时候缺省值只是一个代号，系统是剔除不予计算的啊？

小金鱼

我用这个命令查过海温资料,结果缺测是32767,而且把资料转换成十进制进也可以看出是32767