! +-======-+ 
!  Copyright (c) 2003-2018 United States Government as represented by 
!  the Admistrator of the National Aeronautics and Space Administration.  
!  All Rights Reserved.
!  
!  THIS OPEN  SOURCE  AGREEMENT  ("AGREEMENT") DEFINES  THE  RIGHTS  OF USE,
!  REPRODUCTION,  DISTRIBUTION,  MODIFICATION AND REDISTRIBUTION OF CERTAIN 
!  COMPUTER SOFTWARE ORIGINALLY RELEASED BY THE UNITED STATES GOVERNMENT AS 
!  REPRESENTED BY THE GOVERNMENT AGENCY LISTED BELOW ("GOVERNMENT AGENCY").  
!  THE UNITED STATES GOVERNMENT, AS REPRESENTED BY GOVERNMENT AGENCY, IS AN 
!  INTENDED  THIRD-PARTY  BENEFICIARY  OF  ALL  SUBSEQUENT DISTRIBUTIONS OR 
!  REDISTRIBUTIONS  OF THE  SUBJECT  SOFTWARE.  ANYONE WHO USES, REPRODUCES, 
!  DISTRIBUTES, MODIFIES  OR REDISTRIBUTES THE SUBJECT SOFTWARE, AS DEFINED 
!  HEREIN, OR ANY PART THEREOF,  IS,  BY THAT ACTION, ACCEPTING IN FULL THE 
!  RESPONSIBILITIES AND OBLIGATIONS CONTAINED IN THIS AGREEMENT.
!  
!  Government Agency: National Aeronautics and Space Administration
!  Government Agency Original Software Designation: GSC-15354-1
!  Government Agency Original Software Title:  GEOS-5 GCM Modeling Software
!  User Registration Requested.  Please Visit http://opensource.gsfc.nasa.gov
!  Government Agency Point of Contact for Original Software:  
!  			Dale Hithon, SRA Assistant, (301) 286-2691
!  
! +-======-+ 
function keg (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

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

'set t 'tmin
'getinfo year'
         year = result
'set time JAN'year' DEC'year

'getinfo tmin'
         tmin = result
'getinfo tmax'
         tmax = 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

'set dfile 'expfile
'set x 1'
'set y 1'
'set z 1'
'define net = -( radlwgclim- radswgclim)'

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

'run minmax radswgclim'
       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 = -(oradlwgclim-oradswgclim)'

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

'run minmax oradswgclim'
       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 radlwgclim'
'set cmark  0'
'set cstyle 1'
'set ccolor 4'
'set cthick 6'
'd radswgclim'

'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 oradlwgclim'
'set cmark  0'
'set cstyle 3'
'set ccolor 4'
'set cthick 6'
'd oradswgclim'

'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)
    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 time JAN'year' DEC'year

'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.1
     netmin = netmin * 1.1

'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_CLIM.'obsnam
return