grads画出的等值线相交

yyyy

求助：各位大侠，有谁知道在什么情况下grads画出的等值线会相交的？

我将19720101.atmos_hgt_slp.nc用grads转成dat后
写了gs：
'reinit'
'open F:\month\one\epmonth\1972.ctl'
'enable print F:\month\one\epmonth\1972.gmf'
'set x 1 144'
'set y 1 90'
'set z 1'
'set t 1'
'set csmooth on'
'd ucomp'
'print'

'reinit'

结果就成这样了，有谁知道可能是什么原因吗？

river

ctl呢，ctl才是最重要的最基本的东西。只有ctl编写正确，可以真正描述实际的资料的样子，grads才能画出正确的图

yyyy

river 发表于 2013-5-12 22:49

                               
登录/注册后可看大图

ctl呢，ctl才是最重要的最基本的东西。只有ctl编写正确，可以真正描述实际的资料的样子，grads才能画出正确 ...

这是我的ctl，谢谢指教：
dset F:\month\one\epmonth\19720101.dat
title U
undef -9.9900000E+08
xdef 144 linear 1.25 2.5
ydef 90 levels -89.4944 -87.9775 -85.9551 -83.9326 -81.9101 -79.8876 -77.8652 -75.8427 -73.8202 -71.7978 -69.7753 -67.7528 -65.7303 -63.7079 -61.6854 -59.6629 -57.6404 -55.618 -53.5955 -51.573 -49.5506 -47.5281 -45.5056 -43.4831 -41.4607 -39.4382 -37.4157 -35.3933 -33.3708 -31.3483 -29.3258 -27.3034 -25.2809 -23.2584 -21.236 -19.2135 -17.191 -15.1685 -13.1461 -11.1236 -9.10112 -7.07865 -5.05618 -3.03371 -1.01124 1.01124 3.03371 5.05618 7.07865 9.10112 11.1236 13.1461 15.1685 17.191 19.2135 21.236 23.2584 25.2809 27.3034 29.3258 31.3483 33.3708 35.3933 37.4157 39.4382 41.4607 43.4831 45.5056 47.5281 49.5506 51.573 53.5955 55.618 57.6404 59.6629 61.6854 63.7079 65.7303 67.7528 69.7753 71.7978 73.8202 75.8427 77.8652 79.8876 81.9101 83.9326 85.9551 87.9775 89.4944
zdef 17 levels 1000 925 850 700 600 500 400 300 250 200 150 100 70 50 30 20 10
tdef 1 linear 12Z16DEC1971 1mn
vars 1
ucomp   17  999 *t,z,y,x  zonal wind
endvars

river

yyyy 发表于 2013-5-12 23:39

                               
登录/注册后可看大图

这是我的ctl，谢谢指教：
dset F:\month\one\epmonth\19720101.dat
title U

缺测值好像不太对，你把转换出来的数据用fortran在屏幕上输出一下，看看缺测值是不是-9.9900000E+08

yyyy

river 发表于 2013-5-13 07:25

                               
登录/注册后可看大图

缺测值好像不太对，你把转换出来的数据用fortran在屏幕上输出一下，看看缺测值是不是-9.9900000E+08

输出过了，这就是用fortran输出后改的，原来不是这个。。。。

yyyy

river 发表于 2013-5-13 07:25

                               
登录/注册后可看大图

缺测值好像不太对，你把转换出来的数据用fortran在屏幕上输出一下，看看缺测值是不是-9.9900000E+08

没改以前，画出来的线都挤在一起。。。。。。

river

yyyy 发表于 2013-5-13 08:31

                               
登录/注册后可看大图

没改以前，画出来的线都挤在一起。。。。。。

那提取资料的gs也贴上来看看吧

yyyy

river 发表于 2013-5-13 12:37

                               
登录/注册后可看大图

那提取资料的gs也贴上来看看吧

谢谢你了，我好像知道了，一个研究生学姐说，是q ctlinfo时，从原来nc文件里读出的纬度信息是错的，得改一下。。。不知为什么，grads读的ctl信息好多都是错的。。。谢谢你啦

丹徵

yyyy 发表于 2013-5-13 14:38

                               
登录/注册后可看大图

谢谢你了，我好像知道了，一个研究生学姐说，是q ctlinfo时，从原来nc文件里读出的纬度信息是错的，得改一 ...

那你后来正确读了没？我也读过，而且我配的ctl是直接复制q ctlinfo里的结果，但是画的图比你还过分的，相交，高斯格点那个貌似纬度排列要求很高啊。。。