function radg (args)

expid   = subwrd(args,1)
expfile = subwrd(args,2)
expdsc  = subwrd(args,3)
obsnam  = subwrd(args,4)
obsdsc  = subwrd(args,5)
obsfile = subwrd(args,6)
output  = subwrd(args,7)

'run getenv "GEOSUTIL"'
         geosutil = result

'run getenv "RADLWAVE"'
'define radlwave = 'result
'run getenv "RADSWAVE"'
'define radswave = 'result

'run getenv "ORADLWAVE"'
'define oradlwave = 'result
'run getenv "ORADSWAVE"'
'define oradswave = 'result

'getinfo tmin'
         tmin = result
'getinfo tmax'
         tmax = result

'set dfile 'expfile
'set x 1'
'set y 1'
'set z 1'
'define net = -( radlwg- radswg)'

'run minmax radlwg'
       radmax = subwrd(result,1)
       radmin = subwrd(result,2)

'run minmax radswg'
       dummax = subwrd(result,1)
       dummin = subwrd(result,2)
   if( dummax > radmax ) ; radmax = dummax ; endif
   if( dummin < radmin ) ; radmin = dummin ; endif

'set dfile 'obsfile
'set x 1'
'set y 1'
'set z 1'
'define onet = -(oradlwg-oradswg)'

'run minmax oradlwg'
       dummax = subwrd(result,1)
       dummin = subwrd(result,2)
   if( dummax > radmax ) ; radmax = dummax ; endif
   if( dummin < radmin ) ; radmin = dummin ; endif

'run minmax oradswg'
       dummax = subwrd(result,1)
       dummin = subwrd(result,2)
   if( dummax > radmax ) ; radmax = dummax ; endif
   if( dummin < radmin ) ; radmin = dummin ; endif

radave = (radmax+radmin)/2
raddif = 1.1*(radmax-radave)

radmax = radave + raddif
radmin = radave - raddif

* Plot Global Mean Values
* -----------------------
'set vpage 0 11 3.5 8'
'set grads off'
'set ylab %.0f'
'set axlim 'radmin' 'radmax

'set dfile 'expfile
'set x 1'
'set y 1'
'set z 1'
'set cmark  0'
'set cstyle 1'
'set ccolor 2'
'set cthick 6'
'd radlwg'
'set cmark  0'
'set cstyle 1'
'set ccolor 4'
'set cthick 6'
'd radswg'
'set cmark  0'
'set cstyle 1'
'set ccolor 2'
'set cthick 1'
'd radlwave'
'set cmark  0'
'set cstyle 1'
'set ccolor 4'
'set cthick 1'
'd radswave'

'set dfile 'obsfile
'set x 1'
'set y 1'
'set z 1'
'set cmark  0'
'set cstyle 3'
'set ccolor 2'
'set cthick 6'
'd oradlwg'
'set cmark  0'
'set cstyle 3'
'set ccolor 4'
'set cthick 6'
'd oradswg'

'set cmark  0'
'set cstyle 3'
'set ccolor 2'
'set cthick 1'
'd oradlwave'
'set cmark  0'
'set cstyle 3'
'set ccolor 4'
'set cthick 1'
'd oradswave'


* Plot Global Mean Time-Average Values
* ------------------------------------
'set t 'tmin
'getinfo date'
         date = result

'd radlwave'
   radlwave = subwrd(result,4)
'd radswave'
   radswave = subwrd(result,4)

'q w2xy 'date' 'radlwave
    ylwg = subwrd(result,6)
'q w2xy 'date' 'radswave
    yswg = subwrd(result,6)
say 'ylwg: 'ylwg
say 'yswg: 'yswg
    ydif = ( ylwg-yswg )
if( ydif < 0 ) ; ydif = -ydif ; endif

if( ydif < 0.1 )
if( ylwg <= yswg )
    yswg  = ylwg + 0.1
endif
if( yswg <= ylwg )
    ylwg  = yswg + 0.1
endif
endif

'set strsiz .08'

'set string 2 r 5'
'draw string 1.4 'ylwg ' 'radlwave
'set string 4 r 5'
'draw string 1.4 'yswg ' 'radswave


* Plot NET
* --------
'set t 'tmin' 'tmax

'run minmax net'
       netmax = subwrd(result,1)
       netmin = subwrd(result,2)

'run minmax onet'
       dummax = subwrd(result,1)
       dummin = subwrd(result,2)
   if( dummax > netmax ) ; netmax = dummax ; endif
   if( dummin < netmin ) ; netmin = dummin ; endif

       netmax = netmax * 1.05
       netmin = netmin * 1.05

'set vpage 0 11 0 4.5'
'set grads off'
'set axlim 'netmin' 'netmax
'set cmark  0'
'set cstyle 1'
'set ccolor 3'
'set cthick 6'
'd net'
'set cmark  0'
'set cstyle 3'
'set ccolor 3'
'set cthick 6'
'd onet'
'define radave = radswave-radlwave'
'set cmark  0'
'set cstyle 1'
'set ccolor 3'
'set cthick 1'
'd radave'
'define oradave = oradswave-oradlwave'
'set cmark  0'
'set cstyle 3'
'set ccolor 3'
'set cthick 1'
'd oradave'

'set t 'tmin
'getinfo date'
         date = result
'd radave'
   radave = subwrd(result,4)
'q w2xy 'date' 'radave
    yrad = subwrd(result,6)

'set strsiz .08'

'set string 3 r 5'
'draw string 1.4 'yrad ' 'radave
'set t 'tmin' 'tmax



'set vpage off'
'set string 1 c 6'
'set strsiz .14'
'draw string 6.15 8.35 'expid'  vs.  'obsnam
'draw string 6.15 7.80 Globally Averaged TOA Radiation Budget (W/m`a2`n)'
'set strsiz .12'
'draw string 6.15 8.12 'expdsc

'set string 1 l 6'
'set line 2 1 6'
'draw line   2.0 7.5 2.5 7.5'
'draw string 2.6 7.5 Outgoing LW Rad'
'set line 4 1 6'
'draw line   4.6 7.5 5.1 7.5'
'draw string 5.2 7.5 Net Downward SW Rad'
'set line 3 1 6'
'draw line   7.9 7.5 8.4 7.5'
'draw string 8.5 7.5 Net Downward Rad'

'myprint -name 'output'/RAD_TIME_SERIES.'obsnam
return