ctl文件出现的问题

逸儿~~

运行出sdaa.ctl文件配套的sdaa.grd的gs文件如下：
'reinit'
'sdfopen d:\plan4\netcdf00.nc'
'sdfopen d:\plan4\netcdf06.nc'
'sdfopen d:\plan4\netcdf12.nc'
'sdfopen d:\plan4\netcdf18.nc'

'set fwrite d:\plan4\sdaa.grd'
'set gxout fwrite'
'set lon 70 105'
'set lat 25 40'
'set lev 0 0'
i=1
while(i<=16436)
'set t 'i''
'define sda=sd+sd.2(time+6hr)+sd.3(time+12hr)+sd.4(time+18hr)'
'sdaa=sda/4'
'd sdaa'
i=i+1
endwhile
'disable fwrite'
;

gs文件中4个nc资料分别为每日00,06,12,18时刻的雪深，通过文件里的循环求得日平均，也就是sdaa.grd~nc资料的详细信息如下：


sdaa.grd的配套sdaa.ctl文件如下：
dset d:\plan4\sdaa.grd
title snow depth of everyday
undef -9.96921e+36
xdef 15 linear 70 2.5
ydef 7  linear 47 2.5
zdef 1  linear 0  0
tdef 16436  linear 00Z01SEP1957 1440mn
vars 1
sd 0 -999 Snow depth
endvars
在grads里运行的错误提示如下：



垂直方向nc资料是一层lev o o,进行数据处理时也是set lev 0 0,ctl里zdef是 1linear 0 0,最后为什么会出现那样的提示呢？？实在找不到问题，还请各位大神指教！！

river

逸儿~~ 发表于 2014-11-17 13:43
这个是积雪深度资料，截图1中是原始4个不同时刻nc资料的描述，上面zdef 1 linear 0 0，我以为这是表示地 ...

zdef 1 linear 0 0 这是NC自描述文件的一部分，和平时咱们用的那个ctl 虽然看上去是一样的，但其实还是隐含了很多其他的东西。因为NC资料的格式有其复杂性，具体的可以看一下王老师的帖子http://bbs.06climate.com/forum.php?mod=viewthread&tid=1267。

又说你说了半天，那么改了以后是对了吗？

river

你不是有一层么，怎么还是zdef 1  linear 0  0
改成zdef 1  linear 1  1试试吧还有缺测值也不对，经过grads提取，缺测值一般是-9.99E+8

逸儿~~

river 发表于 2014-11-17 12:09
你不是有一层么，怎么还是zdef 1  linear 0  0
改成zdef 1  linear 1  1试试吧还有缺测值也不对，经过grad ...

这个是积雪深度资料，截图1中是原始4个不同时刻nc资料的描述，上面zdef 1 linear 0 0，我以为这是表示地表呢，不知道地表垂直方向算不算一层，不然可能是我叙述错误了。。有问题的ctl文件描述的grd文件只是对这4个不同时刻的nc资料取了个平均，x,y,z方向没有变化，所以ctl中zdef还是按原始资料取了zdef 1 linear 0 0，没想到grads中提示错误。。。

逸儿~~

river 发表于 2014-11-17 14:53
zdef 1 linear 0 0 这是NC自描述文件的一部分，和平时咱们用的那个ctl 虽然看上去是一样的，但其实还是隐 ...

对了对了对了！！非常感谢，另外也要谢谢你宝贵的时间~~

逸儿~~

river 发表于 2014-11-17 14:53
zdef 1 linear 0 0 这是NC自描述文件的一部分，和平时咱们用的那个ctl 虽然看上去是一样的，但其实还是隐 ...

怎么可以酱紫，理由写了一堆，却没显示全。。。后面的大体意思是说一看见你回复的最后一句话就想笑，感觉好像犬夜叉的哥哥杀生丸明明是为别人做了好事却让别人感觉自己像坏人，当然不完全是这样喽，杀生丸也很好玩儿的~~多耽误了你一点时间，把问题弄得更清楚啦~总之，多谢咯O(∩_∩)O~~

river

逸儿~~ 发表于 2014-11-17 15:23
怎么可以酱紫，理由写了一堆，却没显示全。。。后面的大体意思是说一看见你回复的最后一句话就 ...

你太客气了。我的贡献够用，你自己留着下载东西吧。想多给几个，但是上限就是2······

逸儿~~

river 发表于 2014-11-17 18:24
你太客气了。我的贡献够用，你自己留着下载东西吧。想多给几个，但是上限就是2······

已经很好喽~