grads折线图y轴左右两边设置不同刻度

曾记否

grads的神级坐标大家要多学习清风大侠的神帖http://bbs.06climate.com/forum.php?mod=viewthread&tid=10110
我这里只是随便弄个简单的供大家参考参考
比如有2组数：（ -446.8047      -289.0048      -21.10474       360.4953      -71.70477     -1.304749      -21.80475       75.59528       26.19525      -144.2048     157.5953       296.4953      -177.0048      -63.80475      -267.9048     -19.00476       318.5953      -197.2048       415.5953       20.49524    49.79523）

（77.44   78.74   80.35   82.09   83.51   84.86  86.07   87.52   88.64   89.84   99.08  106.48   107.07   117.51   120.52  123.18
126.05   128.71   131.25   133.65   136.20）


这两组数据用同一个坐标肯定看起来不舒服，因此就想第一组数用左边的坐标，第二组用右边的坐标。把这两个数放入同一个grd，写成ctl
dset D:\test\popul_preci.grd
undef -9999.0
xdef 1 linear 1 1
ydef 1 linear 1 1
zdef 1 linear 1 1
tdef 21 linear 01JAN1990 1yr
vars 2
m 0 99 liiang
n 0 99 popolation
endvars


gs程序如下：
'reinit'
'open D:\test\zxt.ctl'
'set x 1'
'set y 1'
'set t 1 21'
'set grads off'
'set grid off'
'enable print D:\test\popul_preci1.gmf'

'set gxout line'
'set cthick 10'
'set cstyle 1'
'set vrange -500 500'
'set ylint 100'
'set ccolor 15'
'set cmark 4'
'd m'
'set ylpos 0.0 r'
'set vrange 70 140'
'set ylint 10'
'set cmark 3'
'd n'
'print'
'disable print'
'reinit'
;

注意红色部分！'set ylpos 0.0 r' 0.0表示离开坐标轴的尺寸，0.0表示不离开，这样看起来好一点，r表示在右边坐标轴标值。'set vrange 70 140'我试了一下可以，但set yaxis命令就不行了。图中黄色的线就是看右边坐标了。x轴上方改变就设置'set xlpos 0.0 u'

njzqxt

顶楼主了，真心漂亮啊。
敢问楼主，能不能把  “把这两个数放入同一个grd，”  的程序也贴上来呢？

lqouc

顶一下~

排山倒海

不错的教程

mofangbao

嗯 这个也不错，GrADS自带函数的挖掘潜力很大

qxtlyf

太强了，真心漂亮啊。

北五味

唔，Grads功能还真强啊，就是指令太多了，看着头疼，不下功夫不行啊

曾记否

njzqxt 发表于 2013-4-27 14:37

                               
登录/注册后可看大图

顶楼主了，真心漂亮啊。
敢问楼主，能不能把  “把这两个数放入同一个grd，”  的程序也贴上来呢？{:eb317 ...

这个用fortran转一下，设一个二维数组读取，应该不困难。你可以自己试一下

曾记否

北五味 发表于 2013-4-27 15:59

                               
登录/注册后可看大图

唔，Grads功能还真强啊，就是指令太多了，看着头疼，不下功夫不行啊

嗯，只要有点印象就好，到时候可以再翻手册

故乡的云