! +-======-+ 
!  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
!  
! +-======-+ 
 reinit
 set display color white
 set clab off
 set csmooth on
 c

* Note:  &0 => experiment id
* Note:  &1 => period to be processed  (eg DJF, DEC, JJA)
* Note:  &2 => path for data xdf_file
* Note:  &3 => output location
* Note:  &4 => GEOS_Util location


xdfopen &2

setx
sety
setz
sett

chckname rice
chckname rliq

makez   qv z
makez   ql z
makez   qi z
makez rliq z
makez rice z
makez   qr z
makez  uim z
makez  vim z
makez thim z
makez fcld z

set x 1
seasonal   qvz
seasonal   qlz
seasonal   qiz
seasonal   qrz
seasonal rliqz
seasonal ricez
seasonal  uimz
seasonal  vimz
seasonal thimz
seasonal fcldz

set lev 1000 100
define kappa = 2/7
define pk = pow(1000/lev,kappa)
set xlopts 1 4 .14
set ylopts 1 4 .14

exec &4/plots/moist/moist_1 &0 &1 &3
exec &4/plots/moist/moist_2 &0 &1 &3
exec &4/plots/moist/moist_3 &0 &1 &3