1 | 'reinit' |
---|
2 | 'open real_z.ctl' |
---|
3 | |
---|
4 | ** need input data on height levels |
---|
5 | |
---|
6 | 'set mpdset hires' |
---|
7 | 'set display color white' |
---|
8 | 'run rgbset.gs' |
---|
9 | |
---|
10 | say 'Create gif images as well (1=yes ; 0=no)' |
---|
11 | pull ans |
---|
12 | frame = 1 |
---|
13 | |
---|
14 | 'q file' |
---|
15 | rec=sublin(result,5) |
---|
16 | _endtime=subwrd(rec,12) |
---|
17 | _endlev=subwrd(rec,9) |
---|
18 | |
---|
19 | runscript = 1 |
---|
20 | dis_t = 1 |
---|
21 | |
---|
22 | while(runscript) |
---|
23 | dis_z = 1 |
---|
24 | 'set t ' dis_t |
---|
25 | 'q dims' |
---|
26 | rec=sublin(result,5) |
---|
27 | _analysis=subwrd(rec,6) |
---|
28 | say 'Time is ' _analysis |
---|
29 | |
---|
30 | runlevel = 1 |
---|
31 | while(runlevel) |
---|
32 | 'set z ' dis_z |
---|
33 | 'q dims' |
---|
34 | rec=sublin(result,4) |
---|
35 | newz=subwrd(rec,6) |
---|
36 | say 'Level is ' newz |
---|
37 | |
---|
38 | ** Looking for 2 km data |
---|
39 | if (newz=2) |
---|
40 | say 'Plotting' |
---|
41 | 'c' |
---|
42 | 'set grads off' |
---|
43 | 'set gxout shaded' |
---|
44 | 'set clevs 10 20 30 40 50 60 70 80 90 100' |
---|
45 | 'set ccols 0 0 0 0 33 34 35 36 37 38 39' |
---|
46 | 'd rh' |
---|
47 | 'set gxout contour' |
---|
48 | 'set ccolor 4' |
---|
49 | 'set cthick 6' |
---|
50 | 'set cint 2' |
---|
51 | 'd p' |
---|
52 | 'set ccolor 2' |
---|
53 | 'set cthick 6' |
---|
54 | 'set cint 2' |
---|
55 | 'd tc' |
---|
56 | 'set strsiz .2' |
---|
57 | 'set string 1 l 6' |
---|
58 | 'draw string 2.7 8.35 Height (m,blue), T(C,red), RH(color)' |
---|
59 | 'set strsiz .15' |
---|
60 | 'set string 1 l 3' |
---|
61 | 'draw string 7.5 8.0 ' _analysis |
---|
62 | 'draw string 1.5 8.0 ' newz ' km' |
---|
63 | 'run cbar.gs' |
---|
64 | if(ans) |
---|
65 | 'printim realz'frame'.gif gif ' |
---|
66 | frame=frame+1 |
---|
67 | endif |
---|
68 | pull dummy |
---|
69 | endif |
---|
70 | |
---|
71 | ** Looking for 5 km data |
---|
72 | if (newz=5) |
---|
73 | say 'Plotting' |
---|
74 | 'c' |
---|
75 | 'set ccolor 4' |
---|
76 | 'set cthick 6' |
---|
77 | 'set cint 2' |
---|
78 | 'd p' |
---|
79 | 'set ccolor 2' |
---|
80 | 'set cthick 6' |
---|
81 | 'set cint 2' |
---|
82 | 'd tc' |
---|
83 | 'set strsiz .2' |
---|
84 | 'set string 1 l 6' |
---|
85 | 'draw string 3.7 8.35 Height (m,blue), T(C,red)' |
---|
86 | 'set strsiz .15' |
---|
87 | 'set string 1 l 3' |
---|
88 | 'draw string 7.5 8.0 ' _analysis |
---|
89 | 'draw string 1.5 8.0 ' newz ' km' |
---|
90 | if(ans) |
---|
91 | 'printim realz'frame'.gif gif ' |
---|
92 | frame=frame+1 |
---|
93 | endif |
---|
94 | pull dummy |
---|
95 | endif |
---|
96 | |
---|
97 | ** Looking for 10 km data |
---|
98 | if (newz=10) |
---|
99 | say 'Plotting' |
---|
100 | 'c' |
---|
101 | 'set ccolor 4' |
---|
102 | 'set cthick 6' |
---|
103 | 'set cint 2' |
---|
104 | 'd p' |
---|
105 | 'set ccolor 2' |
---|
106 | 'set cthick 6' |
---|
107 | 'set cint 2' |
---|
108 | 'd tc' |
---|
109 | 'set strsiz .2' |
---|
110 | 'set string 1 l 6' |
---|
111 | 'draw string 3.7 8.35 Height (m,blue), T(C,red)' |
---|
112 | 'set strsiz .15' |
---|
113 | 'set string 1 l 3' |
---|
114 | 'draw string 7.5 8.0 ' _analysis |
---|
115 | 'draw string 1.5 8.0 ' newz ' km' |
---|
116 | if(ans) |
---|
117 | 'printim realz'frame'.gif gif ' |
---|
118 | frame=frame+1 |
---|
119 | endif |
---|
120 | pull dummy |
---|
121 | endif |
---|
122 | |
---|
123 | ** Looking for 16 km data |
---|
124 | if (newz=16) |
---|
125 | say 'Plotting' |
---|
126 | 'c' |
---|
127 | 'set ccolor 4' |
---|
128 | 'set cthick 6' |
---|
129 | 'set cint 2' |
---|
130 | 'd p' |
---|
131 | 'set ccolor 2' |
---|
132 | 'set cthick 6' |
---|
133 | 'set cint 2' |
---|
134 | 'd tc' |
---|
135 | 'set strsiz .2' |
---|
136 | 'set string 1 l 6' |
---|
137 | 'draw string 3.7 8.35 Height (m,blue), T(C,red)' |
---|
138 | 'set strsiz .15' |
---|
139 | 'set string 1 l 3' |
---|
140 | 'draw string 7.5 8.0 ' _analysis |
---|
141 | 'draw string 1.5 8.0 ' newz ' km' |
---|
142 | if(ans) |
---|
143 | 'printim realz'frame'.gif gif ' |
---|
144 | frame=frame+1 |
---|
145 | endif |
---|
146 | pull dummy |
---|
147 | endif |
---|
148 | |
---|
149 | if ( dis_z=_endlev ) |
---|
150 | runlevel=0 |
---|
151 | endif |
---|
152 | dis_z = dis_z + 1 |
---|
153 | endwhile |
---|
154 | |
---|
155 | if ( dis_t=_endtime ) |
---|
156 | runscript=0 |
---|
157 | endif |
---|
158 | dis_t = dis_t + 1 |
---|
159 | endwhile |
---|
160 | |
---|