function hplt (args)

expid   = subwrd(args,1)
field   = subwrd(args,2)
season  = subwrd(args,3)
output  = subwrd(args,4)
level   = subwrd(args,5)
nmod    = subwrd(args,6)
nobs    = subwrd(args,7)
expfile = subwrd(args,8)
anafile = subwrd(args,9)
anal    = subwrd(args,10)
obsname = subwrd(args,11)
debug   = subwrd(args,12)
expdsc  = subwrd(args,13)
 
blak   = 0

* Get Dates
* ---------
'getenv "BEGDATEO"'
         begdateo = result
'getenv "ENDDATEO"'
         enddateo = result

'getenv "BEGDATE"'
         begdate  = result
'getenv "ENDDATE"'
         enddate  = result
if( begdate = "NULL" )
   'set dfile 'expfile
   'set t    '1
   'getinfo date'
         begdate = result
endif
if( enddate = "NULL" )
   'set dfile 'expfile
   'getinfo tdim'
            tdim = result
   'set t  'tdim
   'getinfo date'
         enddate = result
endif

'getenv "CLIMATE"'
         climate = result
if( begdate = begdateo & enddate = enddateo ) 
         climate = 'Actual'
endif

* Check for Contour Level Type
* ----------------------------
'getenv "LEVTYPE"'
           CLEVS  = result

'getenv "STD_DEV"'
         std_dev  = result

* Set GC Components
* -----------------
if( field = 'SLP'   ) ; GC = DYN   ; name = slp   ; endif
if( field = 'U'     ) ; GC = DYN   ; name = uwnd  ; endif
if( field = 'V'     ) ; GC = DYN   ; name = vwnd  ; endif
if( field = 'T'     ) ; GC = DYN   ; name = tmpu  ; endif
if( field = 'Q'     ) ; GC = MOIST ; name = sphu  ; endif
if( field = 'RH2'   ) ; GC = MOIST ; name = rh    ; endif
if( field = 'OMEGA' ) ; GC = DYN   ; name = omega ; endif
if( field = 'ZLE'   ) ; GC = DYN   ; name = hght  ; endif
if( field = 'HE'    ) ; GC = DYN   ; name = hghte ; endif
if( field = 'CHI'   ) ; GC = DYN   ; name = chi   ; endif
if( field = 'PSI'   ) ; GC = DYN   ; name = psi   ; endif

* Get Plotting Values from Resource File
* --------------------------------------
'getenv "GEOSUTIL"'
         geosutil = result
PLOTRC = geosutil'/plots/grads_util/plot.rc'

                        'getresource 'PLOTRC' 'field'_'GC'_'level'_FACTOR'
if( result = 'NULL' ) ; 'getresource 'PLOTRC' 'field'_'GC'_FACTOR' ; endif
                                                           fact = result

                        'getresource 'PLOTRC' 'field'_'GC'_'level'_TITLE'
if( result = 'NULL' ) ; 'getresource 'PLOTRC' 'field'_'GC'_TITLE' ; endif
                                                           desc = result

                        'getresource 'PLOTRC' 'field'_'GC'_'level'_CINT'
if( result = 'NULL' ) ; 'getresource 'PLOTRC' 'field'_'GC'_CINT' ; endif
                                                           cint = result

                        'getresource 'PLOTRC' 'field'_'GC'_'level'_CCOLS'
if( result = 'NULL' ) ; 'getresource 'PLOTRC' 'field'_'GC'_CCOLS' ; endif
                                                           ccols = result

                        'getresource 'PLOTRC' 'field'_'GC'_'level'_CVALS'
if( result = 'NULL' ) ; 'getresource 'PLOTRC' 'field'_'GC'_CVALS' ; endif
                                                           cvals = result

                        'getresource 'PLOTRC' 'field'_'GC'_'level'_DCOLS'
if( result = 'NULL' ) ; 'getresource 'PLOTRC' 'field'_'GC'_DCOLS' ; endif
                                                           dcols = result

                        'getresource 'PLOTRC' 'field'_'GC'_'level'_'CLEVS
if( result = 'NULL' ) ; 'getresource 'PLOTRC' 'field'_'GC'_'CLEVS ; endif
                                                           dlevs = result

                        'getresource 'PLOTRC' 'field'_'GC'_'level'_CLAB'
if( result = 'NULL' ) ; 'getresource 'PLOTRC' 'field'_'GC'_CLAB' ; endif
                                                           clab = result

                        'getresource 'PLOTRC' 'field'_'GC'_'level'_CBSCALE'
if( result = 'NULL' ) ; 'getresource 'PLOTRC' 'field'_'GC'_CBSCALE' ; endif
                                                           cbscale = result

if( cbscale = 'NULL' ) ; cbscale =  0.8      ; endif
if( clab    = 'NULL' ) ; clab    =  on       ; endif

say ''

* Set Default Plotting Variables
********************************

if( fact    = 'NULL' ) ; fact    = 1         ; endif
if( ccols   = 'NULL' ) ; ccols   = "84 0 85" ; endif


* Sea-Level Pressure
* ------------------
if( field = 'SLP' )
        cvals = "988 1020"
        ccols = "84 0 85"
         fact = 1
    std_cint  = 2  
    std_cvals = '4 8'
    std_dlevs = "-3 -2.6 -2.3 -2 -1.6 -1.3 -1 -.6  -.3 .3 .6 1 1.3 1.6 2 2.3 2.6 3"
endif

* Temperature
* -----------
if( field = 'T' )
    'getresource 'PLOTRC' T_DYN_'level'_CVALS' ; cvals = result
     std_cint  = 0.5  
     std_cvals = '1 2'
     std_dlevs = "-2.25 -2 -1.75 -1.5 -1.25 -1 -.75 -.5 -.25 .25 .5 .75 1 1.25 1.5 1.75 2 2.25"
endif

* OMEGA
* -----
if( field = 'OMEGA' )
     std_cint  = 0.5  
     std_cvals = '1 2'
     std_dlevs = "-2.25 -2 -1.75 -1.5 -1.25 -1 -.75 -.5 -.25 .25 .5 .75 1 1.25 1.5 1.75 2 2.25"
endif

* Winds
* -----
if( field = 'U' ) 
         cvals = "0"
     std_cint  = 1.0  
     std_cvals = '2 5'
     std_dlevs = "-2.25 -2 -1.75 -1.5 -1.25 -1 -.75 -.5 -.25 .25 .5 .75 1 1.25 1.5 1.75 2 2.25"
endif

if( field = 'V' ) 
         cvals = "0"
     std_cint  = 1.0  
     std_cvals = '2 5'
     std_dlevs = "-2.25 -2 -1.75 -1.5 -1.25 -1 -.75 -.5 -.25 .25 .5 .75 1 1.25 1.5 1.75 2 2.25"
endif


* Specific Humidity
* -----------------
if( field = 'Q' )

    'getresource 'PLOTRC' Q_MOIST_'level'_CVALS' ; cvals = result

     std_cint  = 0.2  
     std_cvals = '.4 1'
     std_dlevs = "-.9 -.8 -.7 -.6 -.5 -.4 -.3 -.2 -.1 .1 .2 .3 .4 .5 .6 .7 .8 .9"

    if( level = 100 )
     std_cint  = 0.1  
     std_cvals = '.2 .4'
     std_dlevs = "-.45 -.4 -.35 -.3 -.25 -.2 -.15 -.1 -.05 .05 .1 .15 .2 .25 .3 .35 .4 .45"
    endif
    if( level = 150 )
     std_cint  = 2 
     std_cvals = '4 8'
     std_dlevs = "-9 -8 -7 -6 -5 -4  -3  -2  -1 1 2 3 4 5 6 7 8 9"
    endif
    if( level = 200 )
     std_cint  = 10 
     std_cvals = '10 40'
     std_dlevs = "-18 -16 -14 -12 -10 -8 -6 -4 -2 2 4 6 8 10 12 14 16 18"
    endif
    if( level = 250 )
     std_cint  = 20 
     std_cvals = '40 80'
     std_dlevs = "-45 -40 -35 -30 -25 -20 -15 -10 -5 5 10 15 20 25 30 35 40 45"
    endif
    if( level = 300 )
     std_cint  = 20 
     std_cvals = '40 80'
     std_dlevs = "-45 -40 -35 -30 -25 -20 -15 -10 -5 5 10 15 20 25 30 35 40 45"
    endif
endif

* Relative Humidity
* -----------------
if( field = 'RH2' )

     std_cint  = 4.0  
     std_cvals = '8 12'
     std_dlevs = "-9 -8 -7 -6 -5 -4  -3  -2  -1 1 2 3 4 5 6 7 8 9"

    if( level =  30 )
        dlevs = "-18 -16 -14 -12 -10 -8  -6  -4  -2 2 4 6 8 10 12 14 16 18"
    endif
    if( level =  50 )
        dlevs = "-27 -24 -21 -18 -15 -12  -9 -6  -3 3 6 9 12 15 18 21 24 27"
    endif
    if( level =  70 )
        dlevs = "-45 -40 -35 -30 -25 -20 -15 -10 -5 5 10 15 20 25 30 35 40 45"
    endif
    if( level = 100 )
        dlevs = "-45 -40 -35 -30 -25 -20 -15 -10 -5 5 10 15 20 25 30 35 40 45"
    endif
    if( level = 150 )
        dlevs = "-45 -40 -35 -30 -25 -20 -15 -10 -5 5 10 15 20 25 30 35 40 45"
    endif
    if( level = 200 )
        dlevs = "-90 -80 -70 -60 -50 -40 -30 -20 -10 10 20 30 40 50 60 70 80 90"
    endif
    if( level = 250 )
        dlevs = "-90 -80 -70 -60 -50 -40 -30 -20 -10 10 20 30 40 50 60 70 80 90"
    endif
    if( level = 300 )
        dlevs = "-90 -80 -70 -60 -50 -40 -30 -20 -10 10 20 30 40 50 60 70 80 90"
    endif
    if( level = 400 )
        dlevs = "-90 -80 -70 -60 -50 -40 -30 -20 -10 10 20 30 40 50 60 70 80 90"
    endif
    if( level = 500 )
        dlevs = "-45 -40 -35 -30 -25 -20 -15 -10 -5 5 10 15 20 25 30 35 40 45"
    endif
    if( level = 600 )
        dlevs = "-45 -40 -35 -30 -25 -20 -15 -10 -5 5 10 15 20 25 30 35 40 45"
    endif
    if( level = 700 )
        dlevs = "-45 -40 -35 -30 -25 -20 -15 -10 -5 5 10 15 20 25 30 35 40 45"
    endif
    if( level = 850 )
        dlevs = "-27 -24 -21 -18 -15 -12  -9 -6  -3 3 6 9 12 15 18 21 24 27"
    endif
    if( level = 925 )
        dlevs = "-18 -16 -14 -12 -10 -8  -6  -4  -2 2 4 6 8 10 12 14 16 18"
    endif
    if( level = 1000 )
        dlevs = "-18 -16 -14 -12 -10 -8  -6  -4  -2 2 4 6 8 10 12 14 16 18"
    endif
endif

* Heights
* -------
if( field = 'ZLE' )
    if( level >= 300 )
    desc  = 'Height (m)'
    else
    desc  = 'Height (km) '
    endif

    if( level >= 300 )
     std_cint  = 10.0  
     std_cvals = '30 60'
     std_dlevs = "-27 -24 -21 -18 -15 -12  -9  -6  -3 3 6 9 12 15 18 21 24 27"
    else
     std_cint  = 0.02  
     std_cvals = '0.04 0.1'
     std_dlevs = "-.09 -.08 -.07 -.06 -.05 -.04 -.03 -.02 -.01 .01 .02 .03 .04 .05 .06 .07 .08 .09"
    endif

    if( level = 30 )
        cvals = "23.0 23.80"
        dlevs = "-.9 -.8 -.7 -.6 -.5 -.4 -.3 -.2 -.1 .1 .2 .3 .4 .5 .6 .7 .8 .9"
    endif
    if( level =  50 )
        cvals = "19.8 20.60"
        dlevs = "-.9 -.8 -.7 -.6 -.5 -.4 -.3 -.2 -.1 .1 .2 .3 .4 .5 .6 .7 .8 .9"
    endif
    if( level =  70 )
        cvals = "17.8 18.60"
        dlevs = "-.9 -.8 -.7 -.6 -.5 -.4 -.3 -.2 -.1 .1 .2 .3 .4 .5 .6 .7 .8 .9"
    endif
    if( level = 100 )
        cvals = "15.4 16.60"
        dlevs = "-.9 -.8 -.7 -.6 -.5 -.4 -.3 -.2 -.1 .1 .2 .3 .4 .5 .6 .7 .8 .9"
    endif
    if( level = 150 )
        cvals = "13.0 14.20"
        dlevs = "-.450 -.400 -.350 -.300 -.250 -.200 -.150 -.100 -.050 .050 .100 .150 .200 .250 .300 .350 .400 .450"
    endif
    if( level = 200 )
        cvals = "11.0 12.40"
        dlevs = "-.270 -.240 -.210 -.180 -.150 -.120 -.090 -.060 -.030 .030 .060 .090 .120 .150 .180 .210 .240 .270"
    endif
    if( level = 250 )
        cvals = "9.8 10.80"
        dlevs = "-.270 -.240 -.210 -.180 -.150 -.120 -.090 -.060 -.030 .030 .060 .090 .120 .150 .180 .210 .240 .270"
    endif
    if( level = 300 )
        cvals = "8500 9600"
        dlevs = "-180 -160 -140 -120 -100 -80 -60 -40 -20 20 40 60 80 100 120 140 160 180"
    endif
    if( level = 400 )
        cvals = "6600 7500"
        dlevs = "-135 -120 -105 -90 -75 -60 -45 -30 -15 15 30 45 60 75 90 105 120 135"
    endif
    if( level = 500 )
        cvals = "5220 5820"
        dlevs = "-90 -80 -70 -60 -50 -40 -30 -20 -10 10 20 30 40 50 60 70 80 90"
    endif
    if( level = 600 )
        cvals = "4000 4400"
        dlevs = "-90 -80 -70 -60 -50 -40 -30 -20 -10 10 20 30 40 50 60 70 80 90"
    endif
    if( level = 700 )
        cvals = "2800 3100"
        dlevs = "-90 -80 -70 -60 -50 -40 -30 -20 -10 10 20 30 40 50 60 70 80 90"
    endif
    if( level = 850 )
        cvals = "1200 1500"
        dlevs = "-90 -80 -70 -60 -50 -40 -30 -20 -10 10 20 30 40 50 60 70 80 90"
    endif
    if( level = 925 )
        cvals = "630 810"
        dlevs = "-90 -80 -70 -60 -50 -40 -30 -20 -10 10 20 30 40 50 60 70 80 90"
    endif
    if( level = 1000 )
        cvals = "0 120"
        dlevs = "-90 -80 -70 -60 -50 -40 -30 -20 -10 10 20 30 40 50 60 70 80 90"
    endif
endif

* Eddy Heights
* ------------
if( field = 'HE' )
    blak  = 1
    desc  = 'Eddy Height (m)'
    cvals = "0"

     std_cint  = 10.0  
     std_cvals = '30 60'
     std_dlevs = "-27 -24 -21 -18 -15 -12  -9  -6  -3 3 6 9 12 15 18 21 24 27"

    if( level = 30 )
        cint  = 60
        cvals = "-60 60"
        dlevs = "-90 -80 -70 -60 -50 -40 -30 -20 -10 10 20 30 40 50 60 70 80 90"
    endif
    if( level =  50 )
        cint  = 60
        cvals = "-60 60"
        dlevs = "-90 -80 -70 -60 -50 -40 -30 -20 -10 10 20 30 40 50 60 70 80 90"
    endif
    if( level =  70 )
        cint  = 60
        cvals = "-60 60"
        dlevs = "-90 -80 -70 -60 -50 -40 -30 -20 -10 10 20 30 40 50 60 70 80 90"
    endif
    if( level = 100 )
        cint  = 60
        cvals = "-60 60"
        dlevs = "-90 -80 -70 -60 -50 -40 -30 -20 -10 10 20 30 40 50 60 70 80 90"
    endif
    if( level = 150 )
        cint  = 60
        cvals = "-60 60"
        dlevs = "-90 -80 -70 -60 -50 -40 -30 -20 -10 10 20 30 40 50 60 70 80 90"
    endif
    if( level = 200 )
        cint  = 60
        cvals = "-60 60"
        dlevs = "-90 -80 -70 -60 -50 -40 -30 -20 -10 10 20 30 40 50 60 70 80 90"
    endif
    if( level = 250 )
        cint  = 60
        cvals = "-60 60"
        dlevs = "-90 -80 -70 -60 -50 -40 -30 -20 -10 10 20 30 40 50 60 70 80 90"
    endif
    if( level = 300 )
        cint  = 60
        cvals = "-60 60"
        dlevs = "-90 -80 -70 -60 -50 -40 -30 -20 -10 10 20 30 40 50 60 70 80 90"
    endif
    if( level = 400 )
        cint  = 50
        cvals = "-50 50"
        dlevs = "-90 -80 -70 -60 -50 -40 -30 -20 -10 10 20 30 40 50 60 70 80 90"
    endif
    if( level = 500 )
        cint  = 40
        cvals = "-40 40"
        dlevs = "-90 -80 -70 -60 -50 -40 -30 -20 -10 10 20 30 40 50 60 70 80 90"
    endif
    if( level = 600 )
        cint  = 30
        cvals = "-30 30"
        dlevs = "-90 -80 -70 -60 -50 -40 -30 -20 -10 10 20 30 40 50 60 70 80 90"
    endif
    if( level = 700 )
        cint  = 30
        cvals = "-30 30"
        dlevs = "-90 -80 -70 -60 -50 -40 -30 -20 -10 10 20 30 40 50 60 70 80 90"
    endif
    if( level = 850 )
        cint  = 20
        cvals = "-20 20"
        dlevs = "-90 -80 -70 -60 -50 -40 -30 -20 -10 10 20 30 40 50 60 70 80 90"
    endif
    if( level = 925 )
        cint  = 20
        cvals = "-20 20"
        dlevs = "-90 -80 -70 -60 -50 -40 -30 -20 -10 10 20 30 40 50 60 70 80 90"
    endif
    if( level = 1000 )
        cint  = 20
        cvals = "-20 20"
        dlevs = "-45 -40 -35 -30 -25 -20 -15 -10 -5 5 10 15 20 25 30 35 40 45"
    endif
endif

    dcols = "55  49  47  45  44  36  34  33  32  0  21  22  23  24  25  26  27  28 69"
std_dcols = "55  49  47  45  44  36  34  33  32  0  21  22  23  24  25  26  27  28 69"

* Plot Mean Field
* ---------------
'c'
'set display color white'
'set vpage off'
'set parea off'
'set grid  off'
'set mproj scaled'
'set frame on'
'set xlopts 1 3 .11'
'set ylopts 1 3 .11'
'rgbset'
'set rgb 84 204 204 204'
'set rgb 85 137 137 137'

'set dfile 'expfile
'setlons'
'set lat -90 90'
'set lev 'level
'set t 1'
'q dims'
say 'EXP DIMS Environment: 'result

'define qmod  = mod'season'*'fact
'define maskm = 1 + qmod-qmod'

* Get Dimension of Environment
* ----------------------------
'getinfo lonmin'
         lonbeg = result
'getinfo lonmax'
         lonend = result
'getinfo latmin'
         latbeg = result
'getinfo latmax'
         latend = result

say 'Environment Dimension: 'lonbeg' 'lonend' 'latbeg' 'latend

* Determine DLAT & DLON of Analysis
* ---------------------------------
'set gxout stat'
'set dfile 'anafile
'set z 1'
'set t 1'
'd lev'
dlon = sublin(result,3)
dlat = sublin(result,4)
dlon = subwrd(dlon,9 )
dlat = subwrd(dlat,9 )
'set gxout shaded'

say 'Analysis DLAT: 'dlat
say 'Analysis DLON: 'dlon

'set lon 'lonbeg' 'lonend
'set lat 'latbeg' 'latend
'define qobs  = obs'season'*'fact

say 'define masko = regrid2( maskm,'dlon','dlat',bs_p1,'lonbeg','latbeg' )'
    'define masko = regrid2( maskm,'dlon','dlat',bs_p1,'lonbeg','latbeg' )'
    'define qobso = regrid2( qobs ,'dlon','dlat',bs_p1,'lonbeg','latbeg' )'

'set dfile 'expfile
'set lon 'lonbeg' 'lonend
'set lat 'latbeg' 'latend


* Make Plot
* ---------------------------------
'set vpage 0 8.5 0.0 11'
'set parea 1.5 7.0 7.70 10.50'
'set grads off'

   'set gxout shaded'
   'set clevs 'cvals
   'set ccols 'ccols
   'set clab  'clab
   'd qmod'

if( field != "OMEGA" & field != "RH2" )
   'set gxout contour'
   'set ccolor 1'
   'set cint 'cint
   if( blak = 1 ) ; 'black' ;endif
   'd qmod'
endif

'set parea off'
'set vpage 0 8.5 0.0 11'
'set parea 1.5 7.0 4.30 7.10'
'set grads off'

   'set gxout shaded'
   'set clevs 'cvals
   'set ccols 'ccols
   'd maskout(qobso,masko)'

if( field != "OMEGA" & field != "RH2" )
   'set gxout contour'
   'set ccolor 1'
   'set cint 'cint
   if( blak = 1 ) ; 'black' ;endif
   'd maskout(qobso,masko)'
else
   'cbarn -vert -snum 0.8 -ymid 6.4 -scaley 0.9 '
endif

'set parea off'
'set vpage 0 8.5 0.0 11'
'set parea 1.5 7.0 0.90 3.70'
'set grads off'
'getinfo lon'
         lon = result
'set dfile 'expfile
'set lev 'level
'set t 1'
'q dims'
say 'DIMS before Final Regridding: 'result
'define obsg = regrid2(  qobs,1,1,bs_p1,'lon',-90)'
'define modg = regrid2(  qmod,1,1,bs_p1,'lon',-90)'
'define mask = regrid2( maskm,1,1,bs_p1,'lon',-90)'
'define obsg = maskout( obsg,mask )'

'set gxout shaded'
'set clevs 'dlevs
'set ccols 'dcols
'd modg-obsg'
'cbarn -snum 'cbscale

'stats maskout(modg,abs(obsg))'
 avgmod = subwrd(result,1)
 stdmod = subwrd(result,2)
'stats maskout(obsg,abs(obsg))'
 avgobs = subwrd(result,1)
 stdobs = subwrd(result,2)
'stats maskout(modg-obsg,abs(obsg))'
 avgdif = subwrd(result,1)
 stddif = subwrd(result,2)

'set vpage off'
'set string 1 c 6'
'set strsiz .11'
'draw string 4.25 10.9 EXPID: 'expid'  'expdsc
'draw string 4.25 10.6 'level'-mb 'desc' 'season' ('nmod')'
'draw string 4.25 7.2 'obsname'  'season' ('nobs')  ('climate')'
'draw string 4.25 3.80 'level'-mb  Difference (Top-Middle)'

'set string 1 c 4'
'set strsiz .08'
'draw string 0.52 10  Mean: 'avgmod
'draw string 0.52 9.8  Std: 'stdmod
'draw string 0.52 6.6 Mean: 'avgobs
'draw string 0.52 6.4  Std: 'stdobs
'draw string 0.52 3.2 Mean: 'avgdif
'draw string 0.52 3.0  Std: 'stddif

                date = getdate (begdate)
bmnthm = subwrd(date,1)
byearm = subwrd(date,2)
                date = getdate (enddate)
emnthm = subwrd(date,1)
eyearm = subwrd(date,2)
                date = getdate (begdateo)
bmntho = subwrd(date,1)
byearo = subwrd(date,2)
                date = getdate (enddateo)
emntho = subwrd(date,1)
eyearo = subwrd(date,2)

if( ccols = "84 0 85" )
'set string 1 l 4'
'set strsiz .08'
   'draw string 7.10 10  Beg: 'bmnthm' 'byearm
   'draw string 7.10 9.8 End: 'emnthm' 'eyearm
   'draw string 7.10 6.6 Beg: 'bmntho' 'byearo
   'draw string 7.10 6.4 End: 'emntho' 'eyearo
else
'set string 1 r 4'
'set strsiz .08'
   'draw string 1.05 9.5 Beg: 'bmnthm' 'byearm
   'draw string 1.05 9.3 End: 'emnthm' 'eyearm
   'draw string 1.05 6.1 Beg: 'bmntho' 'byearo
   'draw string 1.05 5.9 End: 'emntho' 'eyearo
endif

'myprint -name 'output'/horiz_'anal'_'name'_'level'.'season
'set clab on'

if( std_dev = 'true' )
if( debug = "debug" )
    say "Hit ENTER to continue..."
    pull flag
endif

* Plot STD Field
* --------------
'c'
'set display color white'
'set vpage off'
'set parea off'
'set grid  off'
'set mproj scaled'
'set frame on'
'set xlopts 1 3 .11'
'set ylopts 1 3 .11'
'rgbset'
'set rgb 84 204 204 204'
'set rgb 85 137 137 137'

'q dims'
say result

'define qmod = mod'season'std'
'define qobs = obs'season'std'

'set vpage 0 8.5 0.0 11'
'set parea 1.5 7.0 7.70 10.50'
'set grads off'
'set gxout shaded'
'set clevs 'std_cvals
'set ccols 0 84 85'
'd qmod'
'set gxout contour'
'set ccolor 1'
'set cint 'std_cint
if( blak = 1 ) ; 'black' ;endif
'd qmod'
'set parea 0 8.5 7.0 11'
'set parea off'

'set vpage 0 8.5 0.0 11'
'set parea 1.5 7.0 4.30 7.10'
'set grads off'
'set gxout shaded'
'set clevs 'std_cvals
'set ccols 0 84 85'
'd qobs'
'set gxout contour'
'set ccolor 1'
'set cint 'std_cint
if( blak = 1 ) ; 'black' ;endif
'd qobs'
'set parea off'

'set vpage 0 8.5 0.0 11'
'set parea 1.5 7.0 0.90 3.70'
'set grads off'
'getinfo lon'
         lon = result
'define obsg = regrid2( qobs,1,1,bs_p1,'lon',-90)'
'define modg = regrid2( qmod,1,1,bs_p1,'lon',-90)'

'set gxout shaded'
'set clevs 'std_dlevs
'set ccols 'std_dcols
'd modg-obsg'
'cbarn -snum 0.75'

'stats maskout(modg,abs(obsg))'
 avgmod = subwrd(result,1)
 stdmod = subwrd(result,2)
'stats maskout(obsg,abs(obsg))'
 avgobs = subwrd(result,1)
 stdobs = subwrd(result,2)
'stats maskout(modg-obsg,abs(obsg))'
 avgdif = subwrd(result,1)
 stddif = subwrd(result,2)

'set vpage off'
'set string 1 c 6'
'set strsiz .11'
'draw string 4.25 10.9 'expdsc
'draw string 4.25 10.6 'level'-mb 'desc' 'season' ('nmod') STD.DEV.'
'draw string 4.25 7.2 'obsname'  'season' ('climate')'
'draw string 4.25 3.80 'level'-mb  Difference (Top-Middle)'

'set string 1 c 4'
'set strsiz .08'
'draw string 0.52 10  Mean: 'avgmod
'draw string 0.52 9.8  Std: 'stdmod
'draw string 0.52 6.6 Mean: 'avgobs
'draw string 0.52 6.4  Std: 'stdobs
'draw string 0.52 3.2 Mean: 'avgdif
'draw string 0.52 3.0  Std: 'stddif

                date = getdate (begdate)
bmnthm = subwrd(date,1)
byearm = subwrd(date,2)
                date = getdate (enddate)
emnthm = subwrd(date,1)
eyearm = subwrd(date,2)
                date = getdate (begdateo)
bmntho = subwrd(date,1)
byearo = subwrd(date,2)
                date = getdate (enddateo)
emntho = subwrd(date,1)
eyearo = subwrd(date,2)

'set string 1 l 4'
'set strsiz .08'
'draw string 7.10 10  Beg: 'bmnthm' 'byearm
'draw string 7.10 9.8 End: 'emnthm' 'eyearm
'draw string 7.10 6.6 Beg: 'bmntho' 'byearo
'draw string 7.10 6.4 End: 'emntho' 'eyearo

'myprint -name 'output'/horiz_'anal'_'name'_'level'_std.'season

endif


'set mproj latlon'
return

function getdate (date,month,year)
       num = 1
       bit = substr(date,num,1)
while( bit != '' )
       num = num+1
       bit = substr(date,num,1)
endwhile
       loc = num-7
     month = substr(date,loc  ,3)
      year = substr(date,loc+3,4)
return month' 'year