! +-======-+ 
!  Copyright (c) 2003-2007 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)

expdsc = subwrd(args,1)
output = subwrd(args,2)

'run getenv "GEOSUTIL"'
         geosutil = result
'sett'
'getinfo tdim'
         tdim = result

'set grads off'
'set ylab %.1f'
'set axlim -3 4'

'set cmark  0'
'set cstyle 1'
'set ccolor 2'
'set cthick 6'
'd kegeng'
'define kegenave = ave(kegeng,t=1,t='tdim')'
'set cmark  0'
'set cstyle 1'
'set ccolor 2'
'set cthick 1'
'd kegenave'
'set cmark  0'
'set cstyle 3'
'set ccolor 6'
'set cthick 6'
'd kepgg + tedyng'

'set cmark  0'
'set cstyle 1'
'set ccolor 4'
'set cthick 6'
'd kedyng + tedyng'
'define kedynave = ave(kedyng + tedyng,t=1,t='tdim')'
'set cmark  0'
'set cstyle 1'
'set ccolor 4'
'set cthick 1'
'd kedynave'

'set cmark  0'
'set cstyle 2'
'set ccolor 4'
'set cthick 6'
'd kegeng - (kedyng + tedyng)'
'define kedynlos = ave(kegeng-(kedyng+tedyng),t=1,t='tdim')'
'set cmark  0'
'set cstyle 2'
'set ccolor 4'
'set cthick 1'
'd kedynlos'

'set cmark  0'
'set cstyle 1'
'set ccolor 3'
'set cthick 6'
'd kephyg + tephyg'
'define kephyave = ave(kephyg + tephyg,t=1,t='tdim')'
'set cmark  0'
'set cstyle 1'
'set ccolor 3'
'set cthick 1'
'd kephyave'

'set cmark  0'
'set cstyle 1'
'set ccolor 5'
'set cthick 6'
'd keanag + teanag'
'define keanaave = ave(keanag + teanag,t=1,t='tdim')'
'set cmark  0'
'set cstyle 1'
'set ccolor 5'
'set cthick 1'
'd keanaave'

'set cmark  0'
'set cstyle 1'
'set ccolor 1'
'set cthick 6'
'd kedyng+tedyng + kephyg+tephyg + keanag+teanag'

'set cmark  0'
'set cstyle 1'
'set ccolor 1'
'set cthick 3'
'define ketotave = ave( KEDYNg+TEDYNg + KEPHYg+TEPHYg + KEANAg+TEANAg, t=1,t='tdim')'
'd ketotave'

* Plot Time-Average Values
* ------------------------
'set t 'tdim
'getinfo date'
         date = result
'd kegenave'
   kegenave = subwrd(result,4)
'd kedynave'
   kedynave = subwrd(result,4)
'd kedynlos'
   kedynlos = subwrd(result,4)
'd kephyave'
   kephyave = subwrd(result,4)
'd keanaave'
   keanaave = subwrd(result,4)
'd ketotave'
   ketotave = subwrd(result,4)

'q w2xy 'date' 'kegenave
    ygen = subwrd(result,6)
'q w2xy 'date' 'kedynave
    ydyn = subwrd(result,6)
'q w2xy 'date' 'kedynlos
    ylos = subwrd(result,6)
'q w2xy 'date' 'kephyave
    yphy = subwrd(result,6)
'q w2xy 'date' 'keanaave
    yana = subwrd(result,6)
'q w2xy 'date' 'ketotave
    ytot = subwrd(result,6)

'set strsiz .08'

 'set string 1 l 5'
 'draw string 10.52 'ytot ' 'ketotave

 'set string 2 r 5'
 'draw string 1.58 'ygen ' 'kegenave
 'set string 4 r 5'
 'draw string 1.58 'ydyn ' 'kedynave
 'set string 4 r 5'
 'draw string 1.58 'ylos ' 'kedynlos
 'set string 3 r 5'
 'draw string 1.58 'yphy ' 'kephyave
 'set string 5 r 5'
 'draw string 1.58 'yana ' 'keanaave

'set vpage off'
'set string 1 c 6'
'set strsiz .14'
'draw string 6.15 8.25 'expdsc
'draw string 6.15 7.95 Globally Integrated Kinetic Energy Generation (W/m`a2`n)'
'lines 'geosutil'/plots/tseries/ke.stk 1'

'myprint -name 'output'/KE_TIME_SERIES'
return