! +-======-+ 
!  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 ps1 (args)

expid  = subwrd(args,1)
expdsc = subwrd(args,2)
output = subwrd(args,3)
qlname = subwrd(args,4)
qiname = subwrd(args,5)

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

'set t 1 'tdim


* Compute Total, Dry, and Moist Instantaneous Values
* --------------------------------------------------
'define qtot = VAR1g/100'
'define qdry = VAR1g/100-(VAR2g+VAR3g+VAR4g)/9.8'
'define qwet =           (VAR2g+VAR3g+VAR4g)/9.8'


* Compute Total, Dry, and Moist Time-Averaged Values
* --------------------------------------------------
'define qtotave = ave(qtot,t=1,t='tdim')'
'define qdryave = ave(qdry,t=1,t='tdim')'
'define qwetave = ave(qwet,t=1,t='tdim')'

* -------------------------------------------------

'run minmax qtot'
       tmax = subwrd(result,1)
       tmin = subwrd(result,2)
'run minmax qtotave'
       tave = subwrd(result,1)
*      tave = 985.505
*      tmax = 985.562
*      tmin = 985.442

       tdif1  = tmax-tave
       tdif2  = tave-tmin
   if( tdif1 >= tdif2 )
       tdif   = 2*tdif1
   else
       tdif   = 2*tdif2
   endif

'run minmax qdry'
       dmax = subwrd(result,1)
       dmin = subwrd(result,2)
'run minmax qdryave'
       dave = subwrd(result,1)
*      dave = 983.099

       ddif1  = dmax-dave
       ddif2  = dave-dmin
   if( ddif1 >= ddif2 )
       ddif   = 2*ddif1
   else
       ddif   = 2*ddif2
   endif

'run minmax qwetave'
     wetave = subwrd(result,1)
*    wetave = 2.408

if( tdif > ddif ) ; ddif = tdif ; endif

* Plot Total Mass
* ---------------
'set vpage 0 11 3.9 8.4'

'set grads off'
'set ylab %.2f'

totave =   tave
totmax = totave + 0.6*ddif
totmin = totave - 0.6*ddif

'set t 1 'tdim
'set axlim 'totmin' 'totmax

'set grid on'
'set ylpos 0 l'
'set ylopts 1 3 0.11'

'set cmark  0'
'set cstyle 1'
'set ccolor 1'
'set cthick 6'
'd qtot'
'set cmark  0'
'set cstyle 1'
'set ccolor 1'
'set cthick 1'
'd qtotave'

* Plot Time-Average Values
* ------------------------
'set t 'tdim
'getinfo date'
         date = result
'd qtotave'
   qtotave = subwrd(result,4)
'q w2xy 'date' 'qtotave
    ytot = subwrd(result,6)
'd qdryave'
   qdryave = subwrd(result,4)

'set strsiz .08'
 'set string 1 l 6'
 'draw string 10.52 'ytot ' 'qtotave


* Plot Dry Mass and QV
* --------------------
'set vpage 0 11 0 4.5'
'set grads off'
'set ylab %.2f'

dryave =   dave
drymax = dryave + 0.6* ddif
drymin = dryave - 0.6* ddif

wetmax = wetave + 0.6* ddif
wetmin = wetave - 0.6* ddif

'set t 1 'tdim
'set axlim 'drymin' 'drymax

'set ylpos 0 l'
'set ylopts 4 3 0.11'

'set cmark  0'
'set cstyle 1'
'set ccolor 4'
'set cthick 6'
'd qdry'
'set cmark  0'
'set cstyle 1'
'set ccolor 4'
'set cthick 1'
'd qdryave'

'set axlim 'wetmin' 'wetmax
'set ylopts 2 3 0.11'
'set ylpos 0 r'
'set grid off'
'set cmark  0'
'set cstyle 1'
'set ccolor 2'
'set cthick 6'
'd qwet'


* Plot Time-Average Values
* ------------------------
'set strsiz .08'
 'set string 4 l 6'
 'draw string 0.78 'ytot ' 'qdryave

'set vpage off'
'set string 1 c 6'
'set strsiz .14'
'draw string 6.15 8.25 'expid': 'expdsc
'draw string 6.15 7.95 Globally Integrated Total (Dry + QV+QL+QI) Surface Pressure (mb)'

    wetname = "QV"
if( qlname != NULL ) 
    wetname = wetname"+QL"
endif
if( qiname != NULL ) 
    wetname = wetname"+QI"
endif

'draw string 6.15 4.00 Dry (Blue) and Wet ['wetname'] (Red) Surface Pressure (mb)'

'myprint -name 'output'/PS_TIME_SERIES'
return