* +-======-+ 
*  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 gwd_2 (args)

expid  = subwrd(args,1)
season = subwrd(args,2)
output = subwrd(args,3)
num    = subwrd(args,4)
debug  = subwrd(args,5)

'run getinfo undef'
             undef = result

'setlons'
'set z 1'
'set clab off'

* Total Surface Stress
* --------------------
'vpage 1 1 2 2'

'run qminmax mag(tauorox'season'+taubkgx'season',tauoroy'season'+taubkgy'season')'
     qmin = subwrd(result,1)
     qmax = subwrd(result,2)
 if( qmin != undef | qmax != undef )
'set gxout shaded'
'black 0.0149'
'set cint 0.015'
'set rbrange 0 .15'
'd mag( tauorox'season'+taubkgx'season' , tauoroy'season'+taubkgy'season' )'
'cbarn'
'set gxout contour'
'black 0.0149'
'set clevs 0.015'
'set ccolor 1'
'd mag( tauorox'season'+taubkgx'season' , tauoroy'season'+taubkgy'season' )'
 else
   'set strsiz .25'
   'set string  1 c 6'
   'draw string 5.5 5 Not Exported'
 endif
'draw title Total Surface Stress from GWD (N/m`a2`n)'

* Orographic Surface Stress
* -------------------------
'vpage 2 1 2 2'
'run qminmax mag(tauorox'season',tauoroy'season')'
     qmin = subwrd(result,1)
     qmax = subwrd(result,2)
 if( qmin != undef | qmax != undef )
'set gxout shaded'
'black 0.0149'
'set cint 0.015'
'set rbrange 0 .15'
'd mag( tauorox'season' , tauoroy'season' )'
'cbarn'
'set gxout contour'
'black 0.0149'
'set clevs 0.015'
'set ccolor 1'
'd mag( tauorox'season' , tauoroy'season' )'
 else
   'set strsiz .25'
   'set string  1 c 6'
   'draw string 5.5 5 Not Exported'
 endif
'draw title Orographic Surface Stress from GWD (N/m`a2`n)'

* Moist-Processes Surface Stress
* ------------------------------
*vpage 1 2 2 2'
*set gxout shaded'
*black 0.0149'
*set cint 0.015'
*set rbrange 0 .15'
*d mag( taumstx'season' , taumsty'season' )*10'
*cbarn'
*set gxout contour'
*black 0.0149'
*set clevs 0.015'
*set ccolor 1'
*d mag( taumstx'season' , taumsty'season' )*10'
*draw title Moist-Processes Surface Stress from GWD (N/m`a2`n x 10`a-1`n)'

* Background Surface Stress
* -------------------------
'vpage 2 2 2 2'
'run qminmax mag(taubkgx'season',taubkgy'season')'
     qmin = subwrd(result,1)
     qmax = subwrd(result,2)
 if( qmin != undef | qmax != undef )
'set gxout shaded'
'black 0.0149'
'set cint 0.015'
'set rbrange 0 .15'
'd mag( taubkgx'season' , taubkgy'season' )*100'
'cbarn'
'set gxout contour'
'black 0.0149'
'set clevs 0.015'
'set ccolor 1'
'd mag( taubkgx'season' , taubkgy'season' )*100'
 else
   'set strsiz .25'
   'set string  1 c 6'
   'draw string 5.5 5 Not Exported'
 endif
'draw title Background Surface Stress from GWD (N/m`a2`n x 10`a-2`n)'

'set vpage off'
'set strsiz .13'
'set string 1 c 6'

'xlabel 1 5.5 8.3'
'draw string 5.5 8.1 'season' ('num')'

'myprint -name 'output'/gwds1.'season

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

'c'


* Total Surface Stress
* --------------------
'vpage 1 1 2 2'
'run qminmax mag(tauorox'season'+taubkgx'season',tauoroy'season'+taubkgy'season')'
     qmin = subwrd(result,1)
     qmax = subwrd(result,2)
 if( qmin != undef | qmax != undef )
'define us = tauorox'season'+taubkgx'season''
'define vs = tauoroy'season'+taubkgy'season''
'define sp = mag(us,vs)'
'define mask = sp-0.01'
'define us = -maskout(us,mask)'
'define vs = -maskout(vs,mask)'
'd skip(us,2) ; skip(vs,2)'
 else
   'set strsiz .25'
   'set string  1 c 6'
   'draw string 5.5 5 Not Exported'
 endif
'draw title Total Surface Stress from GWD (N/m`a2`n)'

* Orographic Surface Stress
* -------------------------
'vpage 2 1 2 2'
'run qminmax mag(tauorox'season',tauoroy'season')'
     qmin = subwrd(result,1)
     qmax = subwrd(result,2)
 if( qmin != undef | qmax != undef )
'define us = tauorox'season''
'define vs = tauoroy'season''
'define sp = mag(us,vs)'
'define mask = sp-0.01'
'define us = -maskout(us,mask)'
'define vs = -maskout(vs,mask)'
'd skip(us,2) ; skip(vs,2)'
 else
   'set strsiz .25'
   'set string  1 c 6'
   'draw string 5.5 5 Not Exported'
 endif
'draw title Orographic Surface Stress from GWD (N/m`a2`n)'

* Moist-Processes Surface Stress
* ------------------------------
*vpage 1 2 2 2'
*define us = taumstx'season'*10'
*define vs = taumsty'season'*10'
*define sp = mag(us,vs)'
*define mask = sp-0.01'
*define us = -maskout(us,mask)'
*define vs = -maskout(vs,mask)'
*d skip(us,2) ; skip(vs,2)'
*draw title Moist-Processes Surface Stress from GWD (N/m`a2`n x 10`a-1`n)'

* Background Surface Stress
* -------------------------
'vpage 2 2 2 2'
'run qminmax mag(taubkgx'season',taubkgy'season')'
     qmin = subwrd(result,1)
     qmax = subwrd(result,2)
 if( qmin != undef | qmax != undef )
'define us = taubkgx'season'*100'
'define vs = taubkgy'season'*100'
'define sp = mag(us,vs)'
'define mask = sp-0.01'
'define us = -maskout(us,mask)'
'define vs = -maskout(vs,mask)'
'd skip(us,2) ; skip(vs,2)'
 else
   'set strsiz .25'
   'set string  1 c 6'
   'draw string 5.5 5 Not Exported'
 endif
'draw title Background Surface Stress from GWD (N/m`a2`n x 10`a-2`n)'

'set vpage off'
'set strsiz .13'
'set string 1 c 6'

'xlabel 1 5.5 8.3'
'draw string 5.5 8.1 'season' ('num')'

'myprint -name 'output'/gwds2.'season

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

'c'


* Orographic STD'
* --------------'
'vpage 1 1 2 2'
'run qminmax sgh'season
     qmin = subwrd(result,1)
     qmax = subwrd(result,2)
 if( qmin != undef | qmax != undef )
'set gxout shaded'
'black 49.9'
'set cint 50'
'd sgh'season''
'cbarn'
'set gxout contour'
'black 49.9'
'set clevs 50'
'set ccolor 1'
'd sgh'season''
 else
   'set strsiz .25'
   'set string  1 c 6'
   'draw string 5.5 5 Not Exported'
 endif
'draw title Orographic Stnd.Dev. for GWD (m)'

* Cloud STD'
* ---------'
*vpage 2 1 2 2'
*set gxout shaded'
*black 9.9'
*set cint 10'
*d cldstd'season''
*cbarn'
*set gxout contour'
*black 9.9'
*set clevs 10'
*set ccolor 1'
*d cldstd'season''
*draw title Effective Moist-Processes Stnd.Dev. for GWD (m)'

'set vpage off'
'set strsiz .13'
'set string 1 c 6'

'xlabel 1 5.5 8.3'
'draw string 5.5 8.1 'season' ('num')'

'myprint -name 'output'/gwds3.'season

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

'c'

'set clab on'