分享一个grads批量出图的脚本

kongfeng0824

1. 如题是我批量出图用的。是将处理好的数据拿来出图。没有任何对数据的处理。唯一值得提的是对批量数据的读取。文件的名称是随时间变得。所以按照循环写的。用了通配符把数据批量出图出来的。  'c'
   
2. 'reinit'
   
3. 'open f:\1.ctl'
   
4. year=2008
   
5. **设定开始的年
   
6. month=1
   
7. day=1
   
8. hour=0
   
9. t=1
   
10. while(year<2009)
    
11. while(month<13)
    
12. while(day<32)
    
13. while(hour<24)
    
14. if(month<10)
    
15.     if(day<10)
    
16.         if(hour<10)
    
17. ***  设定输出的字符，输出的格式是4位年+2位月
    
18. aaa=year%0%month%0%day%0%hour
    
19. ***如果月小于10,则在月前加0，保证月是两位
    
20.         else
    
21.           aaa=year%0%month%0%day%hour
    
22.         endif
    
23.     else
    
24.      aaa=year%0%month%day%hour
    
25.     endif
    
26. else
    
27. aaa=year%month%day%hour
    
28. endif
    
29. **如果月大于等于10，则原样输出，因为它已经是2位了
    
30. 
    
31. *'set lon 78 99'
    
32. *'set lat 26.5 36.5'
    
33. 'set t 't   
    
34. *这里就是设置时间，每次循环都会改变t，也就是输出不同时刻的图，一定要注意第一个t后面有一个空格
    
35. *'q dim'
    
36. *it=sublin(result,5)
    
37. *it=subwrd(it,6)
    
38. 
    
39. 
    
40. 'set gxout shaded'
    
41. 'set mpdset cnworld'
    
42. *'set parea 0.5 9.5 0.1 8.0'
    
43. 'set parea 0.7 10.2 0.3 8.5'
    
44. *'set mpdset hires cnworld'
    
45. *which result in overlay of different boudary layer
    
46. * cnworld and another file should be copied from pcgrads to the dat directory in program file
    
47. 'set grads off'
    
48. *30-36:  red (rr%: 0, 600%)
    
49. 'set rgb    30 255 215 225'
    
50. 'set rgb    31 255 190 190'
    
51. 'set rgb    32 255 160 150'
    
52. 'set rgb    33 255 130 115'
    
53. 'set rgb     34 255 100 85'
    
54. 'set rgb      35 255 70 50'
    
55. 'set rgb       36 255 40 0'
    
56. 'set rgb      37 161 1   1'
    
57. 'set rgb      38 121  1  1'
    
58. 'set rgb      39 80   1  1'
    
59. *40-46:green (rr%: 0, 600%)
    
60. 'set rgb 40 225 255 225'
    
61. 'set rgb 41 190 255 190'
    
62. 'set rgb 42 150 255 150'
    
63. 'set rgb 43 115 255 115'
    
64. 'set rgb 44 75  255 75 '
    
65. 'set rgb 45 45  255 45 '
    
66. 'set rgb 46 0   255   0'
    
67. 
    
68. *50-56:  blue (rr%: 0,600%)
    
69. 'set rgb   50  200 254 254'
    
70. 'set rgb   51  178 251 254'
    
71. 'set rgb   52  155 242 254'
    
72. 'set rgb   53  130 230 254'
    
73. 'set rgb   54  108 211 254'
    
74. 'set rgb    55  90 191 254'
    
75. 'set rgb   56   70 166 254'
    
76. 'set rgb   57   42 136 254'
    
77. 'set rgb   58   24 105 254'
    
78. 'set rgb     59   0 64 254'
    
79. 
    
80. *40-46: yello (rr%: 0,600%)
    
81. 'set rgb   60  254 254 230'
    
82. 'set rgb   61  254 254 205'
    
83. 'set rgb   62  254 254 180'
    
84. 'set rgb   63  254 253 155'
    
85. 'set rgb   64  254 253 130'
    
86. 'set rgb   65  254 250 105'
    
87. 'set rgb   66  254 240  80'
    
88. 'set rgb   67  254 238  50'
    
89. 
    
90. 
    
91. *50-56:purple (rr%: 0, 600%)  
    
92. 'set rgb 70 225 225 255'
    
93. 'set rgb 71 190 190 255'
    
94. 'set rgb 72 150 150 255'
    
95. 'set rgb 73 115 115 255'
    
96. 'set rgb 74 75 75 255'
    
97. 'set rgb 75 45 45 255'
    
98. 'set rgb 76 0 0 255'
    
99. 
    
100. 'set clevs 0.001 0.2 0.4 0.6 0.8 1.0 1.2 1.4 1.6 1.8 2.0'
     
101. *'set clevs 0.01 0.05 0.2 0.4 0.6 0.8 1.0 1.2 1.4 1.6 1.8'
     
102. *'set ccols 0 14 4 11 5 13 3 10 7 12 8 2 6'
     
103. 'set ccols 0 71 51 54 55 42 45 62 64 66 33 35'
     
104. 'set csmooth on'
     
105. 'd cmorph'
     
106. *'set gxout contour'
     
107. *'d aod'
     
108. 'draw title CMORPH_'%aaa%''
     
109. 'cbarn.gs 1 0 5.5 0.3 '
     
110. *'cbarn.gs 1 1 10 4'
     
111. *'cbarn.gs'
     
112. 'enable print f:\1\CMORPH_'%aaa%'.gmf'
     
113. *这里的'%aaa%'是把数字转化为字符。得到一个完整的输出文件名
     
114. 'printim  f:\1\CMORPH_'%aaa%'.png x800 y600 white'
     
115. *gif x800 y600 white'
     
116. 'print'
     
117. 'disable print'
     
118. 'c'
     
119. 
     
120. 
     
121. t=t+1
     
122. 
     
123. hour=hour+1
     
124. endwhile
     
125. 
     
126. hour=0
     
127. day=day+1
     
128. endwhile
     
129. 
     
130. day=1
     
131. month=month+1
     
132. endwhile
     
133. 
     
134. month=1
     
135. year=year+1
     
136. endwhile
     
137. 
     
138. 'reinit'

复制代码

mofangbao

如果是楼主自己写的话可以说明一下是原创的,然后介绍一下脚本的作用以及使用的注意点.
如果是转载的话建议楼主能给出原始链接,算是对作者的尊重,谢谢分享!

kongfeng0824

mofangbao 发表于 2013-2-21 08:48

                               
登录/注册后可看大图

如果是楼主自己写的话可以说明一下是原创的,然后介绍一下脚本的作用以及使用的注意点.
如果是转载的话建议 ...

这个是我自己写的。如题是我批量出图用的。是将处理好的数据拿来出图。没有任何对数据的处理。唯一值得提的是对批量数据的读取。文件的名称是随时间变得。所以按照循环写的。用了通配符把数据批量出图出来的。

kongfeng0824

mofangbao 发表于 2013-2-21 08:48

                               
登录/注册后可看大图

如果是楼主自己写的话可以说明一下是原创的,然后介绍一下脚本的作用以及使用的注意点.
如果是转载的话建议 ...

请问清风版主如何设置成原创呢？

mofangbao

kongfeng0824 发表于 2013-2-21 09:30

                               
登录/注册后可看大图

请问清风版主如何设置成原创呢？

原创是我们版主设置的啦 你帖子里提一下就好~

红河岸

1.ctl里是什么也一起分享下吧。

kongfeng0824

红河岸 发表于 2013-2-21 12:24

                               
登录/注册后可看大图

1.ctl里是什么也一起分享下吧。

对于ctl是自己写的。可以根据现有的ctl改写，比如你的数据，你就可以改写。

kongfeng0824

mofangbao 发表于 2013-2-21 11:50

                               
登录/注册后可看大图

原创是我们版主设置的啦 你帖子里提一下就好~

多谢清风版主！！！！！

红河岸

kongfeng0824 发表于 2013-2-22 00:03

                               
登录/注册后可看大图

对于ctl是自己写的。可以根据现有的ctl改写，比如你的数据，你就可以改写。

不是，我是想知道你的1.ctl里面是什么？本身就是画一张固定的图还是里面也是批量出图的那种= =

river

不错，很有用的批处理方法