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

'reinit'
'open real_z.ctl'

** need input data on height levels

'set mpdset hires'
'set display color white'
'run rgbset.gs'

say 'Create gif images as well (1=yes ; 0=no)'
pull ans
frame = 1

'q file'
rec=sublin(result,5)
_endtime=subwrd(rec,12)
_endlev=subwrd(rec,9)

runscript = 1
dis_t = 1

while(runscript)
dis_z = 1
'set t ' dis_t
'q dims'
rec=sublin(result,5)
_analysis=subwrd(rec,6)
say 'Time is ' _analysis

runlevel = 1
while(runlevel)
'set z ' dis_z
'q dims'
rec=sublin(result,4)
newz=subwrd(rec,6)
say 'Level is ' newz

** Looking for 2 km data
if (newz=2)
say 'Plotting'
'c'
'set grads off'
'set gxout shaded'
'set clevs 10 20 30 40 50 60 70 80 90 100'        
'set ccols 0 0 0  0 33 34 35 36 37 38  39'
'd rh'
'set gxout contour'
'set ccolor 4'
'set cthick 6'
'set cint 2'
'd p'
'set ccolor 2'
'set cthick 6'
'set cint 2'
'd tc'
'set strsiz .2'
'set string 1 l 6'
'draw string 2.7 8.35  Height (m,blue), T(C,red), RH(color)'
'set strsiz .15'
'set string 1 l 3'
'draw string 7.5 8.0 ' _analysis
'draw string 1.5 8.0 ' newz ' km'
'run cbar.gs'
if(ans)
'printim realz'frame'.gif gif '
frame=frame+1
endif
pull dummy
endif

** Looking for 5 km data
if (newz=5)
say 'Plotting'
'c'
'set ccolor 4'
'set cthick 6'
'set cint 2'
'd p'
'set ccolor 2'
'set cthick 6'
'set cint 2'
'd tc'
'set strsiz .2'
'set string 1 l 6'
'draw string 3.7 8.35  Height (m,blue), T(C,red)'
'set strsiz .15'
'set string 1 l 3'
'draw string 7.5 8.0 ' _analysis
'draw string 1.5 8.0 ' newz ' km'
if(ans)
'printim realz'frame'.gif gif '
frame=frame+1
endif
pull dummy
endif

** Looking for 10 km data
if (newz=10)
say 'Plotting'
'c'
'set ccolor 4'
'set cthick 6'
'set cint 2'
'd p'
'set ccolor 2'
'set cthick 6'
'set cint 2'
'd tc'
'set strsiz .2'
'set string 1 l 6'
'draw string 3.7 8.35  Height (m,blue), T(C,red)'
'set strsiz .15'
'set string 1 l 3'
'draw string 7.5 8.0 ' _analysis
'draw string 1.5 8.0 ' newz ' km'
if(ans)
'printim realz'frame'.gif gif '
frame=frame+1
endif
pull dummy
endif

** Looking for 16 km data
if (newz=16)
say 'Plotting'
'c'
'set ccolor 4'
'set cthick 6'
'set cint 2'
'd p'
'set ccolor 2'
'set cthick 6'
'set cint 2'
'd tc'
'set strsiz .2'
'set string 1 l 6'
'draw string 3.7 8.35  Height (m,blue), T(C,red)'
'set strsiz .15'
'set string 1 l 3'
'draw string 7.5 8.0 ' _analysis
'draw string 1.5 8.0 ' newz ' km'
if(ans)
'printim realz'frame'.gif gif '
frame=frame+1
endif
pull dummy
endif

if ( dis_z=_endlev )
 runlevel=0
endif
dis_z = dis_z + 1
endwhile

if ( dis_t=_endtime )
 runscript=0
endif 
dis_t = dis_t + 1
endwhile