关于nc文件的问题 求指导~~纠结好久啦但是没有解决的问题

dzxconan

大家好~我在处理提取nc资料时候，如果直接画图看上去挺正确（就是用以下gs——*语句部分）
'reinit'
'sdfopen d:\paper\160\sst1854-2005.mon.mean.nc'
'set gxout fwrite'
'set fwrite d:\paper\160\sst.grd'
*'enable print d:\paper\160\sst.gmf'
'set t 1 1814'
n=1164
while(n<=1814)
'set t 'n''
'set lon 0 360'
'set lat -88 88'
'd sst'
*'print'
*'c'
n=n+12
endwhile
'disable fwrite'
*'disable print'
;
但是如果用上述gs先写出到grd中 然后配如下ctl
dset d:\paper\160\sst.grdtitle sstundef -9.96921e+36xdef 180 linear 0 2ydef 89 linear -88 2zdef 1 levels 0tdef 55 linear dec1950 1yrvars 1sst 0 99 sstendvars却得到下面这样的图啦

一直弄不好~~求助la~~~

river

你能不能仔细看看nc资料的缺测值啊，是你ctl里面的那个值吗？

dzxconan

river 发表于 2013-4-25 15:31

                               
登录/注册后可看大图

你能不能仔细看看nc资料的缺测值啊，是你ctl里面的那个值吗？

就是这个呢~~~

dzxconan

如果同时打开两个 还会出现这个 不懂啊~

river

dzxconan 发表于 2013-4-25 15:50

                               
登录/注册后可看大图

就是这个呢~~~

首先你提取资料循环就是时间循环，那在前面设置'set t 1 1814'是没有必要的，去掉。剩下就是ctl也改改试试
dset d:\paper\160\sst.grd
title sst
undef -9.96921e+36
xdef 180 linear 0 2
ydef 89 linear -88 2
zdef 1 levels 0     改成和原来一样的zdef 1 linear 0 0 试试
tdef 55 linear dec1950 1yr
vars 1
sst 0 99 ss
tendvars

dzxconan

river 发表于 2013-4-25 16:23

                               
登录/注册后可看大图

首先你提取资料循环就是时间循环，那在前面设置'set t 1 1814'是没有必要的，去掉。剩下就是ctl也改改试试 ...

表示ctl改成 0 0提示错误 我就改成zdef 1 linear 0 1 可是还是不行 图片一样呢~~

river

dzxconan 发表于 2013-4-25 16:49

                               
登录/注册后可看大图

表示ctl改成 0 0提示错误 我就改成zdef 1 linear 0 1 可是还是不行 图片一样呢~~

你还是看看缺测值吧，可以输出到txt里看看是不是那个缺测值。图画成那样一般都是缺测值的缘故

/猪头哼哼

缺测值的问题

对雪吟歌

应该是缺测值。用 ！ncdump -h 路径加文件 试试，可以看到真实的nc文件信息

一水天下

缺测值不对啊，用fortran把生成的grd转成txt看看缺测值多少，应该和原来的nc缺测不一样