function ensemble(params)
'reinit'
******************************************************
*
* USAGE: run wwiii.gs _tind timestamp product map
*
* map		    map index number
* prod          product number
*               
*******************************************************

*Script location:::

_model = "ENSEMBLE"

_scriptData = "/home/amec/gradscripts/ENSEMBLE/white/"
say "Script directory: "_scriptData

_WRFDir = "/home/amec/models/WRF/GEAST"
_GEMREGDir = "/home/amec/models/GEMREG"
_XGEMGLBDir = "/home/amec/models/XGEMGLB"
_XGFSDir = "/home/amec/models/XGFS"
_NOGAPSDir = "/home/amec/models/NOGAPS"

*Get the arguments from the function call:::

_tind = subwrd(params,1);
print "tind: "_tind
_timeStamp = subwrd(params,2);
print "timeStamp: "_timeStamp
_product = subwrd(params,3);
print "product: "_product
_map = subwrd(params, 4);
print "map: "_map

********************
* The variable names are defined and modified for units here. cmap 
* will convert the model names into these standard names.
********************
* for winds in knots (marine forecasts), use following conversions
gfssurfacewindu = "ugrd10m.3*1.944"
gfssurfacewindv = "vgrd10m.3*1.944"
nogapssfcwindu = "ugrd10m.5*1.944"
nogapssfcwindv = "vgrd10m.5*1.944"
regsurfacewindu = "ugrd10m.1*1.944"
regsurfacewindv = "vgrd10m.1*1.944"
glbsurfacewindu = "ugrd10m.4*1.944"
glbsurfacewindv = "vgrd10m.4*1.944"
wrfsurfacewindu = "ugrd10m.2*1.944"
wrfsurfacewindv = "vgrd10m.2*1.944"
********************
* Set the surface level altitude.
********************

sfclevel = 0

********************
* Select gif files wanted
if (_product = 2); SURFACEWIND = 1; else SURFACEWIND = 0; endif

********************
* Open the file.
********************
GEMREG_file = 'GEMREG.ctl'
WRF_file = 'WRF.ctl'
XGFS_file = 'XGFS.ctl'
XGEMGLB_file = 'XGEMGLB.ctl'
NOGAPS_file = 'NOGAPS1.ctl'

full_file = _GEMREGDir%'/'%_timeStamp%'/'%GEMREG_file
say 'OPENING 'full_file
'open ' full_file

full_file = _WRFDir%'/'%_timeStamp%'/'%WRF_file
say 'OPENING 'full_file
'open ' full_file

full_file = _XGFSDir%'/'%_timeStamp%'/'%XGFS_file
say 'OPENING 'full_file
'open ' full_file

full_file = _XGEMGLBDir%'/'%_timeStamp%'/'%XGEMGLB_file
say 'OPENING 'full_file
'open ' full_file

full_file = _NOGAPSDir%'/'20%_timeStamp%'/'%NOGAPS_file
say 'OPENING 'full_file
'open ' full_file

********************
* Set the forecast time for the printmap script call
********************
ftime = _tind - 1
ftime = ftime + 2*(_tind-1)
if (ftime<10)
   ftime = '00'ftime
endif
if (ftime>10 & ftime<100)
    ftime = '0'ftime
endif
say "PROCESSING FORECAST HOUR "ftime

********************
* This loads the color pallete used for the color bars
********************

'run '_scriptData'rgb'

* Set background to white
'set background 1'

********************
* Scripting for each product available is found below.
* Coarse map products start here.
********************

if (SURFACEWIND=1)
   name = "SURFACEWIND"
   rc = gsfallow("on")
   nrows = 2
   ncols = 2

  'reset'

   panels(2 2)

   ptot = nrows * ncols
  'set mproj latlon'

    set_maps(_map)
   'set map 0 1 5'

    _ioffset = -ftime    
        'set t '_tind
        'c'

    'set line 0 1 6'
    'draw line 0 4.15 11 4.15'
    'set line 0 1 6'
    'draw line 5.5 0 5.5 8.18'
    'set line 0 1 6'
    'draw rec 0 0 11 8.18'
    paneltitle(title)
    'set ylopts 0 5 0.15'
    'set xlopts 0 5 0.15'
    'set parea 0.25 10.7 1.15 7.55'
    'set grid off'
    
* fill in data and plot information

* plot 1
        p = 1
        _vpg.p

        'set mpdset hires'
        'set lat 40 60'
        'set lon -76 -40'

        'e=mag('gfssurfacewindu','gfssurfacewindv')'
        'set gxout shaded'
        'set rgb 60 78 244 91'
        'set rgb 17 186 60 0'
        'set rgb 46 0 186 0'
        'set rgb 38 0 138 0'
        'set rgb 61 116 154 255'
        'set rgb 21 0 160 147'
        'set rgb 70 138 161 255'
        'set rgb 80 198 231 255'
        'set rgb 47 187 110 255'
        'set rgb 34 193 0 227'
        'set rgb 36 107 0 147'
        'set rgb 45 194 0 0'
        'set rgb 55 238 64 24'
        'set rbcols 1 60 46 38 4 70 80 47 34 36'
        'set clevs 10 20 30 40 50 60 70 80 90'
        'set clopts 0 5 0.18'
        'set cthick 5'
        'set csmooth on'
        'd e'
        
        
        ptitle = "GFS Surface Wind"
        title(ptitle)
        dvalidvp(' ',' ',' ','GFS Surface Wind Speed (Knots)')
        'run '_scriptData'gfsprovinces 0 1 1'
        'run '_scriptData'drawtowns-panel2 '_model''
        _drawinitflag = 0
        
*plot #2
             p = 2 
             _vpg.p
        if (_tind<18)
        'set mpdset hires'
        'set lat 40 60'
        'set lon 284 320'

        'e=mag('regsurfacewindu','regsurfacewindv')'
        'set gxout shaded'
        'set rgb 60 78 244 91'
        'set rgb 17 186 60 0'
        'set rgb 46 0 186 0'
        'set rgb 38 0 138 0'
        'set rgb 61 116 154 255'
        'set rgb 21 0 160 147'
        'set rgb 70 138 161 255'
        'set rgb 80 198 231 255'
        'set rgb 47 187 110 255'
        'set rgb 34 193 0 227'
        'set rgb 36 107 0 147'
        'set rgb 45 194 0 0'
        'set rgb 55 238 64 24'
        'set rbcols 1 60 46 38 4 70 80 47 34 36'
        'set clevs 10 20 30 40 50 60 70 80 90'
        'set clopts 0 5 0.18'
        'set cthick 5'
        'set csmooth on'
        'd e'

        'set ccolor 0'
        ptitle = "GEM REGIONAL Surface Wind"
        title(ptitle)
        _drawinitflag = 0
        dvalidvp(' ',' ',' ','GEM Regional Surface Wind Speed (Knots)')
        'run '_scriptData'cmcprovinces 0 1 1'
        'run '_scriptData'drawtowns-panel '_model''
        endif
        
        if (_tind=18)
        'set string 0 c'
        'set strsiz 0.35'
        'draw string 5.4 4.7 Switching from GEMREG to NOGAPS'
        endif 
        
        if (_tind>18)
        'set mpdset hires'
        'set lat 40 60'
        'set lon -76 -40'

        'e=mag('nogapssfcwindu','nogapssfcwindv')'
        'set gxout shaded'
        'set rgb 60 78 244 91'
        'set rgb 17 186 60 0'
        'set rgb 46 0 186 0'
        'set rgb 38 0 138 0'
        'set rgb 61 116 154 255'
        'set rgb 21 0 160 147'
        'set rgb 70 138 161 255'
        'set rgb 80 198 231 255'
        'set rgb 47 187 110 255'
        'set rgb 34 193 0 227'
        'set rgb 36 107 0 147'
        'set rgb 45 194 0 0'
        'set rgb 55 238 64 24'
        'set rbcols 1 60 46 38 4 70 80 47 34 36'
        'set clevs 10 20 30 40 50 60 70 80 90'
        'set clopts 0 5 0.18'
        'set cthick 5'
        'set csmooth on'
        'd e'

        'set ccolor 0'
        ptitle = "NOGAPS Surface Wind"
        title(ptitle)
        dvalidvp(' ',' ',' ','NOGAPS Surface Wind Speed (Knots)')
        'run '_scriptData'gfsprovinces 0 1 1'
        'run '_scriptData'drawtowns-panel2 '_model''
        _drawinitflag = 0
        endif
        
*plot #3
             p = 3 
             _vpg.p
             
        'set mpdset hires'
        'set lat 40 60'
        'set lon 284 320'
             
        'e=mag('glbsurfacewindu','glbsurfacewindv')'
        'set gxout shaded'
        'set rgb 60 78 244 91'
        'set rgb 17 186 60 0'
        'set rgb 46 0 186 0'
        'set rgb 38 0 138 0'
        'set rgb 61 116 154 255'
        'set rgb 21 0 160 147'
        'set rgb 70 138 161 255'
        'set rgb 80 198 231 255'
        'set rgb 47 187 110 255'
        'set rgb 34 193 0 227'
        'set rgb 36 107 0 147'
        'set rgb 45 194 0 0'
        'set rgb 55 238 64 24'
        'set rbcols 1 60 46 38 4 70 80 47 34 36'
        'set clevs 10 20 30 40 50 60 70 80 90'
        'set clopts 0 5 0.18'
        'set cthick 5'
        'set csmooth on'
        'd e'
        
        'set ccolor 0'
        ptitle = "GEM GLOBAL Surface Wind"
        title(ptitle)
        dvalidvp(' ',' ',' ','GEM Global Surface Wind Speed (Knots)')
        'run '_scriptData'cmcprovinces 0 1 1'
        'run '_scriptData'drawtowns-panel '_model''
        _drawinitflag = 1
        
*plot #4
             p = 4 
             _vpg.p
             
        'set mpdset hires'
        'set lat 40 60'
        'set lon -76 -40'
             
        'e=mag('wrfsurfacewindu','wrfsurfacewindv')'
        'set gxout shaded'
        'set rgb 60 78 244 91'
        'set rgb 17 186 60 0'
        'set rgb 46 0 186 0'
        'set rgb 38 0 138 0'
        'set rgb 61 116 154 255'
        'set rgb 21 0 160 147'
        'set rgb 70 138 161 255'
        'set rgb 80 198 231 255'
        'set rgb 47 187 110 255'
        'set rgb 34 193 0 227'
        'set rgb 36 107 0 147'
        'set rgb 45 194 0 0'
        'set rgb 55 238 64 24'
        'set rbcols 1 60 46 38 4 70 80 47 34 36'
        'set clevs 10 20 30 40 50 60 70 80 90'
        'set clopts 0 5 0.18'
        'set cthick 5'
        'set csmooth on'
        'd e'
        
        if (_tind<22)
        'set ccolor 1'
        ptitle = "WRF-GEAST Surface Wind"
        title(ptitle)
        dvalidvp(' ',' ',' ','WRF Surface Wind Speed (Knots)')
        'run '_scriptData'gfsprovinces 0 1 1'
        'run '_scriptData'drawtowns-panel2 '_model''
        _drawinitflag = 0
        endif
        'run '_scriptData'printmap 1 0 'ftime' 'name' '_timeStamp' '_model
endif

********************
* Exit GrADS.  Comment out when testing individual products
********************

'quit'

********************
* Function: title
* Input: A string to be viewed at the top of the map as the title.
* Output: A title is drawn at the top of the map.
********************
function paneltitle(str)
'set string 0 c'
'set strsiz 0.12'
if (_product = 2)
'draw string 5.6 8.43 Ensemble Wind Forecast'
endif
'q dims'
b = sublin(result, 5)
word = subwrd(b, 6)
hour = substr(word,1,2)
day = substr(word,4,2)
month = substr(word,6,3)
year = substr(word,9,4)
'query time'
if (subwrd(result,6)= Sat)
daywk = Saturday
endif
if (subwrd(result,6)= Sun)
daywk = Sunday
endif
if (subwrd(result,6)= Mon)
daywk = Monday
endif
if (subwrd(result,6)= Tue)
daywk = Tuesday
endif
if (subwrd(result,6)= Wed)
daywk = Wednesday
endif
if (subwrd(result,6)= Thu)
daywk = Thursday
endif
if (subwrd(result,6)= Fri)
daywk = Friday
endif

* Check Daylight savings flag and display accordingly
ISDSTFILE="/home/amec/gradscripts/ISDST"
file = read(ISDSTFILE)
error = close(ISDSTFILE)
ISDST = sublin(file,2)
if (ISDST = 0)
   time_convert(hour,day,month,year,-4,nst)
   time_convert(hour,day,month,year,-4,ast)
   'set strsiz 0.10'
*   'draw string 5.5 1.10 Forecast valid '_hour.nst':30 NST '_month.nst' '_day.nst' 'year' ('daywk' 'month' 'day' 'hour'Z)'
   'draw string 5.8 8.28 Forecast valid '_hour.nst':30 NST '_month.nst' '_day.nst' 'year' ('daywk' 'month' 'day' 'hour'Z)'
else
   time_convert(hour,day,month,year,-3,ndt)
   time_convert(hour,day,month,year,-3,adt)
   'set strsiz 0.10'
*   'draw string 5.5 1.10 Forecast valid '_hour.ndt':30 NDT '_month.ndt' '_day.ndt' 'year' ('daywk' 'month' 'day' 'hour'Z)'
   'draw string 5.8 8.28 Forecast valid '_hour.ndt':30 NDT '_month.ndt' '_day.ndt' 'year' ('daywk' 'month' 'day' 'hour'Z)'
endif
time_convert(hour,day,month,year,_ioffset,init) 
'set strsiz 0.09'
'draw string 9.9 8.45 Initialized: '_hour.init'Z, '_month.init' '_day.init' '_year.init

'set string 0 tc'
'set strsiz 0.1'
'draw string 1.18 8.5 AMEC Earth & Environmental' 

return

function title(str)
'set strsiz  0.2 0.2'
'set string 0 tc'
'draw string 5.5 7.8 'str

return 

*function amec()
*'set strsiz 0.13'
*'set string 0 bl'
*'draw string 0.05 0.05 AMEC Earth & Environmental' 
*return

function website()
'set string 0 bl'
'set strsiz 0.21'
'draw string 0.05 0.08 http:/www.weather.na.amec.com' 
return

********************
* Function: MaxMin
* Plots the maximum and minimum of the variable and places it in the bottom left corner
********************
function minmax(var)
'set gxout stat'
'd 'var
ret = result
line = sublin(ret,8)
dmin = subwrd(line,4)
dmax = subwrd(line,5)
return dmin' 'dmax

return

********************
* Function: dvalid
* Input: A string to be viewed on the weather map, usually a title.  
*        Strings entered as individual arguments (ie. seperated by a comma)
*        will be placed on consecutive lines after another.
* Output: Draws the color bar, forcast time, copyright and 
*         a the title of the product.
********************

function dvalid(str.1,str.2,str.3,str.4,str.5)
'run cbarn 0.62 0 5.5 0.75'
'q dims'
b = sublin(result, 5)
word = subwrd(b, 6)
hour = substr(word,1,2)
day = substr(word,4,2)
month = substr(word,6,3)
year = substr(word,9,4)

* Check Daylight savings flag and display accordingly
ISDSTFILE="/home/amec/gradscripts/ISDST"
file = read(ISDSTFILE)
error = close(ISDSTFILE)
ISDST = sublin(file,2)
if (ISDST = 0)
   time_convert(hour,day,month,year,-4,nst)
   time_convert(hour,day,month,year,-4,ast)
   'set strsiz 0.12'
   'set string 0 tc'
   'draw string 5.5 1.10 Forecast valid '_hour.nst':30 NST '' ('_hour.ast':00 AST) ['hour'Z] 'month' 'day' 'year
   'draw string 5.5 8.25 Forecast valid '_hour.nst':30 NST '' ('_hour.ast':00 AST) ['hour'Z] 'month' 'day' 'year
else
   time_convert(hour,day,month,year,-3,ndt)
   time_convert(hour,day,month,year,-3,adt)
   'set strsiz 0.12'
   'set string 0 tc'
   'draw string 5.5 1.10 Forecast valid '_hour.ndt':30 NDT ' '('_hour.adt':00 ADT) ['hour'Z] 'month' 'day' 'year
   'draw string 5.5 8.25 Forecast valid '_hour.ndt':30 NDT ' '('_hour.adt':00 ADT) ['hour'Z] 'month' 'day' 'year
endif

* Display the string list
'set string 0 tc'
'set strsiz 0.12'
i=1
height = 0.65
while(str.i != 'str.'i)
   'draw string 5.5 'height' 'str.i
   height=height-0.14
   i = i + 1
endwhile

* Display initialization time
time_convert(hour,day,month,year,_ioffset,init) 
'set strsiz 0.09'
'set string 0 tr'
'draw string 10.95 8.45 Initialized: '_hour.init'Z, '_month.init' '_day.init' '_year.init

return

********************
* Function: arrlab
* Input: x,y - location on map to display
*        length - length of arrow to be drawn
*        scale - the scale label of the arrow
* Output: An arrow with a label below it is draw.  This is 
*         used in lieu of the grads arrow label which is 
*         drawn too close to the color bar.
********************

function arrlab(x,y,len,scale)
'set line 0 0 4'
'draw line 'x+len/2-0.05' 'y+0.025' 'x+len/2' 'y 
'draw line 'x+len/2-0.05' 'y-0.025' 'x+len/2' 'y 
'draw line 'x-len/2' 'y' 'x+len/2' 'y
'set string 0 c'
'set strsiz 0.1'
'draw string 'x' 'y' 'scale' km/hr'
return

********************
* Function: time_convert
* Input: hour - hour to convert
*        day - day to convert
*        year - year to convert 
*        offset - number of hours to shift time by
*        name - name of new time variables in script.  eg. _hour.pst
* Output: The function creates global variables with the name tacked
*         onto the end of the time unit.  The variables follow
*         the format _unit.name.  eg _day.mst
********************

function time_convert(hour,day,month,year,offset,name)

_hour.name = hour
*say "Hour: "_hour.name
_day.name = day
*say "Day: "_day.name
_month.name = month
*say "Month: "_month.name
_year.name = year
*say "Year: "_year.name

*say "Offset: "offset
*say "Name: "name


feb_date = 28

* Do leap year calculation
if (math_fmod(_year.name,4) =  0)
   feb_date = 29
endif
if (math_fmod(_year.name,100) = 0)
   feb_date = 28
endif
if (math_fmod(_year.name,400) = 0)
   feb_date = 29
endif

* Calculate the offset and adjust the day
_hour.name = hour + offset
if (_hour.name > 24)
   dayinc = math_int(_hour.name/24)
   _day.name = day + dayinc
   _hour.name = _hour.name - dayinc*24
endif
if (_hour.name < 0)
   dayinc = math_int(_hour.name/24)
   _day.name = day + dayinc
   _hour.name = dayinc*-24 + _hour.name
endif

* Adjust month and year if nessecary
if (_month.name=JAN & _day.name>31)
   _month.iname = FEB
   _day.name = _day.name - 31
endif
if (_month.name=JAN & _day.name<1)
   _month.name = DEC
   _day.name = 31 + _day.name 
   _year.name = year - 1
endif
if (_month.name=FEB & _day.name>feb_date)
   _month.name = MAR
   _day.name = _day.name - feb_date
endif
if (_month.name=FEB & _day.name<1)
   _month.name = JAN
   _day.name = 31 + _day.name
endif
if (_month.name=MAR & _day.name>31)
   _month.name = APR
   _day.name = _day.name - 31
endif
if (_month.name=MAR & _day.name<1)
   _month.name = FEB
   _day.name = feb_date + _day.name
endif
if (_month.name=APR & _day.name>30)
   _month.name = MAY
   _day.name = _day.name - 30
endif
if (_month.name=APR & _day.name<1)
   _month.name = MAR
   _day.name = 31 + _day.name
endif
if (_month.name=MAY & _day.name>31)
   _month.name = JUN
   _day.name = _day.name - 31
endif
if (_month.name=MAY & _day.name<1)
   _month.name = APR
   _day.name = 30 + _day.name
endif
if (_month.name=JUN & _day.name>30)
   _month.name = JUL
   _day.name = _day.name - 30
endif
if (_month.name=JUN & _day.name<1)
   _month.name = MAY
   _day.name = 31 + _day.name
endif
if (_month.name=JUL & _day.name>31)
   _month.name = AUG
   _day.name = _day.name - 31
endif
if (_month.name=JUL & _day.name<1)
   _month.name = JUN
   _day.name = 30 + _day.name
endif
if (_month.name=AUG & _day.name>31)
   _month.name = SEP
   _day.name = _day.name - 31
endif
if (_month.name=AUG & _day.name<1)
   _month.name = JUL
   _day.name = 31 + _day.name
endif
if (_month.name=SEP & _day.name>30)
   _month.name = OCT
   _day.name = _day.name - 30
endif
if (_month.name=SEP & _day.name<1)
   _month.name = AUG
   _day.name = 31 + _day.name
endif
if (_month.name=OCT & _day.name>31)
   _month.name = NOV
   _day.name = _day.name - 31
endif
if (_month.name=OCT & _day.name<1)
   _month.name = SEP
   _day.name = 30 + _day.name
endif
if (_month.name=NOV & _day.name>30)
   _month.name = DEC
   _day.name = _day.name - 30
endif
if (_month.name=NOV & _day.name<1)
   _month.name = OCT
   _day.name = 31 + _day.name 
endif
if (_month.name=DEC & _day.name>31)
   _month.name = JAN 
   _day.name = _day.name - 31
   _year.name = year + 1
endif
if (_month.name=DEC & _day.name<1)
   _month.name = NOV
   _day.name = 30 + _day.name
endif
return 

function set_maps(params)
*********************************************************
* This script sets the map properties and the lat and long 
* boundaries for amecmaps.gs.  The positions of the map on 
* the page is also set here.
*
* USAGE: set_maps map
*
*        map - the map number to display
*********************************************************

map = subwrd(params, 1)

**********************
* Set the drawing area and turn off the time labels.
**********************

*'set parea 0.5 10.7 1.15 7.55'
*'set timelab off'

**********************
* Below the lat and long domains as well as the map properties
* are set for each map and each model.
**********************


if (map=0)
'set mpdset hires'
*'set mpt 10 50 1 1'
*'set mpt 0 1 1 1'
*'set mpt 15 1 1'
*'set mpt 0 -1'
*'set map 1 1 10'
'set lat 0 90'
'set lon -90 -50'
endif

if (map=1)
'set mpdset hires'
*'set mpt 10 50 1 1'
*'set mpt 0 1 1 1'
*'set mpt 15 1 1'
*'set mpt 0 -1'
*'set map 1 1 10'
'set lat 20 50'
'set lon 275 330'
endif

if (map=2)
'set mpdset hires'
'set mpt 15 1 1'
'set mpt 0 -1'
'set lat 38.5 50'
'set lon 296 318'
endif

if (map=3)
'set mpdset hires'
'set mpt 0 1'
'set mpt 0 -1'
*'set mpt 1 1'
*'set mpt 15 1'
*'set mpt 16 15'
'set lat 41.5 50'
'set lon 300 316'
endif

********************
* Function: dvalidvp (panel plot)
* Input: A string to be viewed on the weather map, usually a title.  
*        Strings entered as individual arguments (ie. seperated by a comma)
*        will be placed on consecutive lines after another.
* Output: Draws the color bar, forecast time, copyright and 
*         a the title of the product.
********************

function dvalidvp(str.1,str.2,str.3,str.4,str.5)
'run '_scriptData'cbarn 1 0 5.5 0.6'
'q dims'
b = sublin(result, 5)
word = subwrd(b, 6)
hour = substr(word,1,2)
day = substr(word,4,2)
month = substr(word,6,3)
year = substr(word,9,4)

* Check Daylight savings flag and display accordingly
ISDSTFILE="/home/amec/gradscripts/ISDST"
file = read(ISDSTFILE)
error = close(ISDSTFILE)
ISDST = sublin(file,2)
if (ISDST = 0)
   time_convert(hour,day,month,year,-4,nst)
   time_convert(hour,day,month,year,-4,ast)
   'set strsiz 0.13'
   'set string 0 tc'
*   'draw string 5.5 1.10 Forecast valid '_hour.nst':30 NST '' ('_hour.ast':00 AST) ['hour'Z] 'month' 'day' 'year
*   'draw string 5.5 8.25 Forecast valid '_hour.nst':30 NST '' ('_hour.ast':00 AST) ['hour'Z] 'month' 'day' 'year
else
   time_convert(hour,day,month,year,-3,ndt)
   time_convert(hour,day,month,year,-3,adt)
   'set strsiz 0.12'
   'set string 0 tc'
*   'draw string 5.5 1.10 Forecast valid '_hour.ndt':30 NDT ' '('_hour.adt':00 ADT) ['hour'Z] 'month' 'day' 'year
*   'draw string 5.5 8.25 Forecast valid '_hour.ndt':30 NDT ' '('_hour.adt':00 ADT) ['hour'Z] 'month' 'day' 'year
endif

* Display the string list
'set string 0 tc'
'set strsiz 0.17'
i=1
height = 0.75
while(str.i != 'str.'i)
   'draw string 5.5 'height' 'str.i
   height=height-0.15
   i = i + 1
endwhile

*if (_drawinitflag = 1)
* Display initialization time
*time_convert(hour,day,month,year,_ioffset,init) 
*'set strsiz 0.15'
*'set string 0 bl'
*'draw string 0.05 0.05 Initialized: '_hour.init'Z, '_day.init' '_month.init' '_year.init
*endif

return

********************
* Function: dvalidncvp (panel plot)
* Input: A string to be viewed on the weather map, usually a title.
* Output: Identical to dvalid though no color bar is plotted.
********************

function dvalidncvp(str.1,str.2,str.3,str.4,str.5)
'q dims'
b = sublin(result, 5)
word = subwrd(b, 6)
hour = substr(word,1,2)
day = substr(word,4,2)
month = substr(word,6,3)
year = substr(word,9,4)

* Check Daylight savings flag and display accordingly
ISDSTFILE="/home/amec/gradscripts/ISDST"
file = read(ISDSTFILE)
error = close(ISDSTFILE)
ISDST = sublin(file,2)
if (ISDST = 0)
   time_convert(hour,day,month,year,-4,nst)
   time_convert(hour,day,month,year,-4,ast)
   'set strsiz 0.12'
   'set string 0 tc'
*   'draw string 5.5 1.10 Forecast valid '_hour.nst':30 NST '' ('_hour.ast':00 AST) ['hour'Z] 'month' 'day' 'year
*   'draw string 5.5 7.85 Forecast valid '_hour.nst':30 NST '' ('_hour.ast':00 AST) ['hour'Z] 'month' 'day' 'year
else
   time_convert(hour,day,month,year,-3,ndt)
   time_convert(hour,day,month,year,-3,adt)
   'set strsiz 0.12'
   'set string 0 tc'
*   'draw string 5.5 1.10 Forecast valid '_hour.ndt':30 NDT ' '('_hour.adt':00 ADT) ['hour'Z] 'month' 'day' 'year
*   'draw string 5.5 7.85 Forecast valid '_hour.ndt':30 NDT ' '('_hour.adt':00 ADT) ['hour'Z] 'month' 'day' 'year
endif

* Display the string list
'set string 0 tc'
'set strsiz 0.2'
i=1
*height = 1.13
height = 1.3
while(str.i != 'str.'i)
   'draw string 5.5 'height' 'str.i
   height=height-0.24
   i = i + 1
endwhile

* Display initialization time
*if (_drawinitflag = 1)
*time_convert(hour,day,month,year,_ioffset,init) 
*'set strsiz 0.15'
*'set string 0 bl'
*'draw string 0.05 0.05 Initialized: '_hour.init'Z, '_day.init' '_month.init' '_year.init
*'draw string 0.05 0.05 http:AMEC Website' 

*endif

return