问一下nc文件中undef一般设置为多少？

artmunich

我下载了1980年每天四次的再分析uwnd资料，然后写了它的ctl文件
dset f:/data/1980/uwnd.1980.nc
undef -9.99E33
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  1464 linear 0Z01JAN1980 4hr
vars 1
uwnd 17 99 4xDaily U-wind
endvars
我想试一下能不能用，就set gxout contour，画了uwnd，结果图像是这样的：


还说了：
contouring：-1.2e+38 to 1.2e+38 interval 3e+37,数值明显不对啊
我猜测可能是undef的问题？
做过的麻烦指点一下，不胜感激啊。。。

MeteoInfo

GrADS中netCDF数据可以通过三种方式打开：（1）sdfopen，只限于符合COARDS约定（http://ferret.wrc.noaa.gov/noaa_coop/coop_cdf_profile.html）的nc数据。netCDF格式（还有HDF格式）非常灵活，用于气象数据时需要有一些约定，最常用的约定是CF（COARDS可以看作CF约定的子集），对于维、变量、属性有详细的规定，这样以来软件才能通过约定对数据进行正确的判读。不知道为什么GrADS不支持CF约定（COARDS约定过于简单，无法描述复杂的数据），因此sdfopen能打开的数据十分有限。（2）xdfopen，对于非COARDS约定的nc数据可以尝试xdfopen来打开，需要一个简单的control文件，“传说中的谁”给出了一个例子，具体可以看官方说明http://www.iges.org/grads/gadoc/gradcomdxdfopen.html。（3）open，如同打开GrADS二进制文件，需要一个完整的control文件，详见http://www.iges.org/grads/gadoc/SDFdescriptorfile.html。

第2和3种打开nc文件的方法都需要control文件，是为了软件能够正确读取维和变量等信息。MeteoInfo打开nc文件的方式类似sdfopen，无需control文件，不过MeteoInfo支持CF约定（包含了COARDS约定）、WRF Out约定（WRF模式输出）和IOAPI约定（用于CMAQ模式）。因此可以在不需要control文件的情况下打开更多的nc数据文件。

传说中的谁

artmunich 发表于 2011-7-10 00:54

                               
登录/注册后可看大图

哦，对的，改过之后还是错误的，不过ls@随缘 说不需要ctl，可能解释

你可以试试另一个办法，就是用sdfopen之后，用fwrite将数据写成二进制，然后再编写ctl画图。

随缘

nc资料自带ctl说明文件的啊，干什么还要写呢？
直接sdfopen命令打开资料画图就可以了啊

artmunich

随缘 发表于 2011-7-9 23:59

                               
登录/注册后可看大图

nc资料自带ctl说明文件的啊，干什么还要写呢？
直接sdfopen命令打开资料画图就可以了啊

因为下载的资料较多，有些是一天四次的，一些是一天一次的，每次查都挺费事，就写ctl文件方便记忆
nc文件哪里带了ctl说明？我没见过呀

随缘

artmunich 发表于 2011-7-10 00:02

                               
登录/注册后可看大图

因为下载的资料较多，有些是一天四次的，一些是一天一次的，每次查都挺费事，就写ctl文件方便记忆
nc文件 ...

nc资料都是自带ctl说明文件的，在每个字段的头部都有文件信息的

你sdfopen命令打开nc资料后，用q ctlinfo命令查看就可以了

传说中的谁

参照FNL的？undef 9.999E+20

artmunich

随缘 发表于 2011-7-10 00:05

                               
登录/注册后可看大图

nc资料都是自带ctl说明文件的，在每个字段的头部都有文件信息的

你sdfopen命令打开nc资料后，用q ctli ...

谢谢您的帮助，我把里面的数值都改成了头文件中的信息，基本上改后的ctl文件完全就是nc里的头文件了
这样应该是可以的，但是我一方面用sdfopen打开nc文件，画一张图，如下

直接打开nc文件

再用open打开nc文件对应的ctl文件，同样的设置画一张图，如下

ctl作图

难道出现了bug？？？？

artmunich

传说中的谁 发表于 2011-7-10 00:11

                               
登录/注册后可看大图

参照FNL的？undef 9.999E+20

看了头文件里，undef是-9.96921e+36，不过改了也没用，画图照样出现上面的错误啊:Q

传说中的谁

artmunich 发表于 2011-7-10 00:27

                               
登录/注册后可看大图

看了头文件里，undef是-9.96921e+36，不过改了也没用，画图照样出现上面的错误啊

没试过写CTL来画NC的资料，我看看吧。

传说中的谁

artmunich 发表于 2011-7-10 00:27

                               
登录/注册后可看大图

看了头文件里，undef是-9.96921e+36，不过改了也没用，画图照样出现上面的错误啊

小小提醒一下你的CTL里面TDEF有点小错误，应该是6hr才对吧？

artmunich

传说中的谁 发表于 2011-7-10 00:34

                               
登录/注册后可看大图

小小提醒一下你的CTL里面TDEF有点小错误，应该是6hr才对吧？

哦，对的，改过之后还是错误的，不过ls@随缘 说不需要ctl，可能解释