基于opengrads脚本subplot.gs的一页多图的改进

普鲁斯特

一页多图不是一个新鲜的话题了，论坛里的帖子已经讲得很详细，也在此向前辈们致敬~

http://bbs.06climate.com/forum.php?mod=viewthread&tid=4586

http://bbs.06climate.com/forum.php?mod=viewthread&tid=5505

但我可能更喜欢适合自己的脚本，于是根据opengrads里的脚本subplot.gs修改了下。
附件里两个文件：subplot_x根据行来画，subplot_y按列来画。数据采用的是ncep reanalysis II 2010年的。色卡是用调色板调的
色标采用仿matlab格式：

http://bbs.06climate.com/forum.php?mod=viewthread&tid=9294

1. *
   
2. * Subplot(m,n,frame) - Like Matlab function subplot().
   
3. * Usage: subplot_x m n frame
   
4. *          m       :   the number of rows
   
5. *          n        :   the number of columns
   
6. *          frame  :   the figure to draw
   
7. * Examples:
   
8. *         'subplot_x 2 2 1'
   
9. *           |------| |------|
   
10. *           |plot1| |plot2|
    
11. *           |------| |------|
    
12. *           |plot3| |plot4|
    
13. *           |------| |------|      
    
14. *           'subplot_x 3 1 1'
    
15. * Notes:
    
16. *           plot 1st row first, then 2nd row, ...
    
17. * Modified by Nanshui. Sep 27, 2017.
    
18. 
    
19. function subplot(args)
    
20. 
    
21.    m     = subwrd(args,1)
    
22.    n     = subwrd(args,2)
    
23.    frame = subwrd(args,3)
    
24. 
    
25. *  find coordinates for this frame
    
26. *  -------------------------------
    
27.    ii = 0
    
28.     k = 0
    
29.    while ( ii < m )
    
30.       ii = ii + 1
    
31.       jj = 0
    
32.       while ( jj < n )
    
33.           jj = jj + 1
    
34.           k = k + 1
    
35.           if ( k = frame )
    
36.                i = ii
    
37.                j = jj
    
38.                break
    
39.           endif
    
40.       endwhile
    
41.    endwhile
    
42. 
    
43. * Determine page size
    
44. * -------------------
    
45.   'q gxinfo'
    
46.   line = sublin(result,2)
    
47.   xsize = subwrd(line,4)
    
48.   ysize = subwrd(line,6)
    
49. 
    
50.   say xsize ' ' ysize
    
51. 
    
52. * Determine xy in page coordinates
    
53. * --------------------------------
    
54.   dx = (xsize-2) / n
    
55.   dy = (ysize-1.5) / m
    
56. 
    
57.   xoff = 0.125 * dx
    
58.   yoff = 0.125 * dy
    
59. 
    
60.   dy = dy - yoff/4
    
61. 
    
62.   xa = (j-1) * dx + xoff + 0.5
    
63.   xb = xa + dx -  xoff
    
64.   ya = (m-i) * dy + yoff + 1
    
65.   yb = ya + dy - yoff
    
66. *我改动了dx dy xa ya 留出了画色标的地方
    
67. *  Note :
    
68. *        dx dy xa xb ya yb could be modified according to your own prefrence.
    
69. 
    
70.   say 'set parea ' xa ' ' xb ' ' ya ' ' yb
    
71.   'set parea ' xa ' ' xb ' ' ya ' ' yb
    
72. *加入了下面几行，设置xlab ylab在整个图最外侧标
    
73.   if (i < m)
    
74.     'set xlab off'
    
75.   else
    
76.     'set xlab on'
    
77.   endif
    
78.   
    
79.   if( j = 1 )
    
80.     'set ylab on'
    
81.   else
    
82.     'set ylab off'
    
83.   endif
    
84.   
    
85. 
    
86. return
    

复制代码

1. 'reinit'
   
2. 'open D:\Summer\F\NCEP-DOE-R2\vwnd\vwnd.year.ctl'
   
3. 'enable print D:\Summer\F\NCEP-DOE-R2\subplot_test2.gmf'
   
4. 
   
5. 
   
6. 'set grid off'
   
7. 'set lon 0 360'
   
8. 'set lat 0 90'
   
9. 'set z 10'
   
10. 'set map 1'
    
11. 'set mproj scaled'
    
12. 'set ylopts 1 5 0.15'
    
13. 'set xlopts 1 5 0.15'
    
14. 'set annot 1 5'
    
15. 
    
16. fig = abcdefghijklmn
    
17. 
    
18. i=1
    
19. while(i<=8)
    
20. 'set vpage 0 11 0 8.5'
    
21. 'set grads off'
    
22. 'subplot_y 4 2 'i''
    
23. 'set t 'i''
    
24. 'set ylint 20'
    
25. 'set xlint 60'
    
26. 
    
27. 'set gxout shaded'
    
28. 'color_16_BlDarkRe'
    
29. 'set rbcols 17 18 19 20 21 22 23 24 0 255 25 26 27 28 29 30 31 32'
    
30. 'set clevs -40 -35 -30 -25 -20 -15 -10 -5 0 5 10 15 20 25 30 35 40'
    
31. 'd vwnd'
    
32. 
    
33. * 'set gxout contour'
    
34. * 'set clevs -40 -35 -30 -25 -20 -15 -10 -5 5 10 15 20 25 30 35 40'
    
35. * 'd vwnd'
    
36. 
    
37. * 'draw_figstr2 TL 1.0 Fig.'i''
    
38. figi=substr(fig,i,1)
    
39. 'draw_figstr2 TL 1.0 ('figi')'
    
40. i=i+1
    
41. endwhile
    
42. 'cbar_matlab 2 1 0  5.20 0.5'
    
43. 'print'
    
44. 'disable print'
    
45. ;
    

复制代码

另外，在我的画图脚本里用到了@兰溪之水 的drawfigstr, 也根据自己的需要进行了修改。
http://bbs.06climate.com/forum.php?mod=viewthread&tid=15203
修改后可以在BL/BR/TR/TL标记长度不大于10的字符串，字符串外的框会随字符串长度调整。
但限于个人水平，修改之后使用比较受限，

兰溪原版使用如下：

*    1. drawfigstr BottomLeft (a)
*    2. drawfigstr TL (b)
*    3. drawfigstr BR (c) 1.5
修改后只能按照下面的格式
*    1. 'draw_figstr2 BL 1.0 (a)'
*    2. 'draw_figstr2 TL 1.2 (b)'
*    3. 'draw_figstr2 BR 1.5 phase 1'
BL/BR/TR/TL 不能再写全名，只能简写，放大因子必须保留一位小数，所标字符串需要与撇号紧挨，不然图上标记时字符后会有空白。
所用函数一并打包放附件。
附件.zip (1.31 MB, 下载次数: 91)

2017-9-26 14:48 上传

点击文件名下载附件

看云淡风清

楼主，这个很不错，有时间学习下！

river

好棒啊，图很漂亮。感谢分享

junyang517

学习学习～感谢分享

yqcxm

很不错，感谢分享

小其其格

厉害呀楼主，画出来的图也很漂亮！

天边的那边

grads图也很好看了，有些像ncl出来的

风叶ele

thank u very much!
have a try~

sandy175

楼主厉害了，学习学习

sandy175

楼主，cbar_matlab这个函数怎么获取？