一个GrADS的bug: 一行内level数值过多图上会出现细线

StephenTien

分享一点今天捉虫（debug）的经历，给大家一个参考。

今天在使用GrADS2.1版绘制CMIP5 near-surface wind speed 场时，得到的结果有数条细线，如图1。

我深感困惑，遂用matlab绘制的相同的数据，结果未发现任何异常，排除数据问题。检查gs文件后，也基本排除脚本问题。

后来对ctl做了简单的修改，解决了问题，如图2。过程如下：

因为此模式数据水平网格不是线性的，所以用levels 来表示，结构示意如下：

dset ...
title ...
undef ...
xdef 192 levels  0.9375 2.8125 4.6875 6.5625 8.4375 10.3125 12.1875 14.0625 15.9375 17.8125 19.6875 25.3125 ...        
ydef 144 levels  -89.375 -88.125  -86.875 -85.625 -84.375 -83.125 -81.875 -80.625 -79.375 -78.125 -76.875 -75.625 ...
zdef 1 levels 0 1
tdef 1 linear 00Z01JAN1850 1mo
vars 1
sfcwind 1 99
endvars

修改后示意如下：

dset ...
title ...
undef ...
xdef 192 levels  0.9375 2.8125 4.6875 6.5625 8.4375 10.3125
12.1875 14.0625 15.9375 17.8125 19.6875 21.5625 23.4375        25.3125
27.1875        29.0625        30.9375        32.8125        34.6875        36.5625        38.4375        40.3125        
42.1875        44.0625        45.9375        47.8125        49.6875        51.5625        53.4375        55.3125        
57.1875        59.0625        60.9375        62.8125        64.6875        66.5625        68.4375        70.3125 ...
ydef 144 levels  -89.375 -88.125        -86.875 -85.625        -84.375 -83.125
-81.875 -80.625        -79.375 -78.125        -76.875 -75.625        -74.375 -73.125        
-71.875 -70.625        -69.375 -68.125        -66.875 -65.625        -64.375 -63.125        
-61.875 -60.625        -59.375 -58.125        -56.875 -55.625        -54.375 -53.125        
-51.875 -50.625        -49.375 -48.125        -46.875 -45.625        -44.375 -43.125        ...
zdef 1 levels 0 1
tdef 1 linear 00Z01JAN1850 1mo
vars 1
sfcwind 1 99
endvars

差别在于开始我把经度192个level，纬度144个level都写在了一行上，后来我把经纬度的level都分成了很多行。这应该是GrADS的一个bug，在此提醒大家。


        

qixiangren34

一份耕耘一分收获

StephenTien

qixiangren34 发表于 2017-11-28 23:26
一份耕耘一分收获

谢谢您的鼓励

香草拿铁

赞！！！！
有收获

guzhong

你这数据xy是线性的啊，怎么不是了。。
xdef  192 linear 0.9375 1.875
ydef  144 linear  -89.375 1.25
不能写成这样的么。。
不过若是必须穷举的话，楼主的经验可能会有用吧。。

兰北

有收获  谢谢

随缘

真是这样么？去试试去，好多图在生成eps格式的时候有细线，其它格式输出还没发现

自导自演

这倒是没想到，但一般ctl不用这么列举吧

rabin_xu

我倒是发现过有横线的问题，但多是把图放到文档里后才会出现

StephenTien

guzhong 发表于 2017-11-29 11:55
你这数据xy是线性的啊，怎么不是了。。
xdef  192 linear 0.9375 1.875
ydef  144 linear  -89.375 1.25
...

仔细看了一下还真是，惭愧惭愧，多谢指教