source: trunk/MESOSCALE_DEV/SRC/ARWpost/scripts/zlevels.gs @ 897

Last change on this file since 897 was 207, checked in by aslmd, 14 years ago

MESOSCALE: A GENERAL CLEAN-UP FOLLOWING UPDATING THE USER MANUAL. EVERYTHING ESSENTIAL IS IN MESOSCALE (much lighter than before). EVERYTHING FOR DEVELOPPERS OR EXPERTS IS IN MESOSCALE_DEV.

  • Property svn:executable set to *
File size: 2.5 KB
Line 
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
10say 'Create gif images as well (1=yes ; 0=no)'
11pull ans
12frame = 1
13
14'q file'
15rec=sublin(result,5)
16_endtime=subwrd(rec,12)
17_endlev=subwrd(rec,9)
18
19runscript = 1
20dis_t = 1
21
22while(runscript)
23dis_z = 1
24'set t ' dis_t
25'q dims'
26rec=sublin(result,5)
27_analysis=subwrd(rec,6)
28say 'Time is ' _analysis
29
30runlevel = 1
31while(runlevel)
32'set z ' dis_z
33'q dims'
34rec=sublin(result,4)
35newz=subwrd(rec,6)
36say 'Level is ' newz
37
38** Looking for 2 km data
39if (newz=2)
40say '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'
64if(ans)
65'printim realz'frame'.gif gif '
66frame=frame+1
67endif
68pull dummy
69endif
70
71** Looking for 5 km data
72if (newz=5)
73say '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'
90if(ans)
91'printim realz'frame'.gif gif '
92frame=frame+1
93endif
94pull dummy
95endif
96
97** Looking for 10 km data
98if (newz=10)
99say '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'
116if(ans)
117'printim realz'frame'.gif gif '
118frame=frame+1
119endif
120pull dummy
121endif
122
123** Looking for 16 km data
124if (newz=16)
125say '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'
142if(ans)
143'printim realz'frame'.gif gif '
144frame=frame+1
145endif
146pull dummy
147endif
148
149if ( dis_z=_endlev )
150 runlevel=0
151endif
152dis_z = dis_z + 1
153endwhile
154
155if ( dis_t=_endtime )
156 runscript=0
157endif
158dis_t = dis_t + 1
159endwhile
160
Note: See TracBrowser for help on using the repository browser.