'set display color white'
* alters plot aspect ratio
'set vpage 0.5 8.5 0.5 8.5'
'clear'
'set grads off'
'set lev 1500'
'run rgbset.gs'
'set mproj off'

say 'Create gif images as well (1=yes ; 0=no)'
pull ans
frame = 1
'q file'
rec=sublin(result,5)
_endtime=subwrd(rec,12)
runscript = 1
* start at time 20
dis_t = 2
distm1 = 1

* movie loop
while(runscript)

'set t ' dis_t
'q dims'
rec=sublin(result,5)
_analysis=subwrd(rec,6)

say 'Time is ' _analysis

'clear'
'set grads off'
'set mproj off'
'set gxout shaded'
*'set xaxis 70 130 10'
*'set yaxis 0.2 2.2 0.4'
'set clevs   -2.56 -1.28 -.64 -.32 -.16 -.08 -.04 -.02 0 .02 .04 .08 .16 .32 .64 1.28 2.56'
'set ccols   49 48 47 46 45 44 43 42 41 0 61 62 63 64 65 66 67 68 69'
'd w*100'
'set strsiz .2'
'set string 1 l 4'
* hardwired 1800 sec plotting interval
'draw string .5 8.0 w (cm/s)  time = ' distm1*3 ' min'
*'cbarn 1.0 0 6.0'
'cbarn'
'set gxout contour'

if(ans)
'printim movie'frame'.gif gif '
frame=frame+1
endif
pull dummy
if ( dis_t=_endtime )
 runscript=0
endif 
dis_t = dis_t + 1
distm1 = distm1 + 1
endwhile