'reinit'
t1=1
'open swan.ctl'
'q file'
rec=sublin(result,5)
tt=subwrd(rec,12)


while(t1<=tt)
 'set vpage 2 5.5 4 7'
 'set grads off'
 'set mpdset hires'
 'set t 't1
 'q dims'
 line=sublin(result,5)
 time=subwrd(line,6)
 <!-- SET TIMESTEP HERE -->
 'set gxout shaded'
 'color 0 56 3 -kind darkblue->cyan->lime->yellow->orange->red->darkred->magenta->maroon->lavenderblush'
 'd const(htsgw*3.2808399,0,-u)'
 'run cbarn.gs'
 'set gxout vector'
 'set arrlab off'
 'set arrscl 0.2'
 'define rad = 4*atan2(1.0,1.0)/180'
 'define uwave = -1*sin(rad*wavedir)'
 'define vwave = -1*cos(rad*wavedir)'
 'set ccolor 1'
 'd skip(uwave,10,10);skip(vwave,10,10)'
 'set line 1 1 1'
 'draw shp marine_zones'
 'draw shp river_basins'
 'draw shp lakes'
 'set shpopts 15'
 'draw shp rivers'
 'draw shp zones'
 'draw title NWPS Significant Wave Height (ft) and Peak Wave Direction\Hour 't2' ('time')'

 'set vpage 2 5.5 .8 3.8'
 'set grads off'
 'set gxout shaded'
 'set mpdset hires'
 'color 0 135 10 -kind darkblue->dodgerblue->lime->yellow->orange->red->magenta->maroon->lavenderblush'
 'd const(wndspd,0,-u)'
 'run cbarn.gs'
 'set gxout barb'
 'set arrlab off'
 'set arrscl 0.2'
 'define rad = 4*atan2(1.0,1.0)/180'
 'define uwnd = -wndspd*sin(rad*wnddir)'
 'define vwnd = -wndspd*cos(rad*wnddir)'
 'set ccolor 1'
 'd skip(uwnd,10,10);skip(vwnd,10,10)'
 'set line 1 1 1'
 'draw shp marine_zones'
 'draw shp river_basins'
 'draw shp lakes'
 'set shpopts 15'
 'draw shp rivers'
 'draw shp zones'
 'draw title NWPS Wind (knots)\Hour 't2' ('time')'

 'set vpage 5.5 9 4 7'
 'set grads off'
 'set gxout shaded'
 'set mpdset hires'
 'color 0 25 2 -kind darkblue->dodgerblue->yellow->orange->red->magenta->maroon->lavenderblush'
 'd const(waveper,0,-u)'
 'run cbarn.gs'
 'set gxout vector'
 'set arrlab off'
 'set arrscl 0.2'
 'define rad = 4*atan2(1.0,1.0)/180'
 'define uwave = -1*sin(rad*wavedir)'
 'define vwave = -1*cos(rad*wavedir)'
 'set ccolor 1'
 'd skip(uwave,10,10);skip(vwave,10,10)'
 'set line 1 1 1'
 'draw shp marine_zones'
 'draw shp river_basins'
 'draw shp lakes'
 'set shpopts 15'
 'draw shp rivers'
 'draw shp zones'
 'draw title NWPS Peak Wave Period (s) and Direction\Hour 't2' ('time')'

 'set vpage 5.5 9 .8 3.8'
 'set grads off'
 'set gxout shaded'
 'set mpdset hires'
 'color 0 56 3 -kind darkblue->cyan->lime->yellow->orange->red->darkred->magenta->maroon->lavenderblush'
 'd const(swell*3.2808399,0,-u)'
 'run cbarn.gs'
 'set gxout vector'
 'set arrlab off'
 'set arrscl 0.2'
 'set line 1 1 1'
 'draw shp marine_zones'
 'draw shp river_basins'
 'draw shp lakes'
 'set shpopts 15'
 'draw shp rivers'
 'draw shp zones'
 'draw title NWPS Significant Swell Height (ft)\Hour 't2' ('time')'

 'set vpage off'
 'draw string 5.7 7.3 **EXPERIMENTAL**'

 if (t2 < 10)
  'printim swan_4panel_hr00't2'.png x1792 y1344 white'
  '!/usr/bin/convert swan_4panel_hr00't2'.png -draw "image over 270,180 0,0 noaa_logo.gif" swan_4panel_hr00't2'.png'
  '!/usr/bin/convert swan_4panel_hr00't2'.png -draw "image over 1430,180 0,0 nws_logo.gif" swan_4panel_hr00't2'.png'
 else
   if (t2 < 100)
    'printim swan_4panel_hr0't2'.png x1792 y1344 white'
    '!/usr/bin/convert swan_4panel_hr0't2'.png -draw "image over 270,180 0,0 noaa_logo.gif" swan_4panel_hr0't2'.png'
    '!/usr/bin/convert swan_4panel_hr0't2'.png -draw "image over 1430,180 0,0 nws_logo.gif" swan_4panel_hr0't2'.png' 
   else
    'printim swan_4panel_hr't2'.png x1792 y1344 white'
    '!/usr/bin/convert swan_4panel_hr't2'.png -draw "image over 270,180 0,0 noaa_logo.gif" swan_4panel_hr't2'.png'
    '!/usr/bin/convert swan_4panel_hr't2'.png -draw "image over 1430,180 0,0 nws_logo.gif" swan_4panel_hr't2'.png' 
   endif
 endif
 'c'
 t1=t1+1
endwhile

'quit'