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

'setenv "ANTIALIAS" TRUE'

*******************************************************
****              Horizontal Means                 ****
*******************************************************

'set lev 1000'
'movie statdplt.gs "p 0.5 geos ncep" -name std.geos-ncep.p1000 -print -rotate'
'c'

fields = 'h u v t'
levels = '1000 850 700 500 400 300 250 200 150 100'

'numargs  'levels
 numlevs = result

'numargs  'fields
 numflds = result

* Loop over Fields
* ----------------
n = 1
while ( n<=numflds )
              field = subwrd(fields,n)
'set dfile 1'
'setlons'
'sety'
if( field = h ) ; cint = 5   ; endif
if( field = u ) ; cint = 0.8 ; endif
if( field = v ) ; cint = 0.8 ; endif
if( field = t ) ; cint = 0.2 ; endif

* Loop over Levels
* ----------------
z = 1
while ( z<=numlevs )
              level = subwrd(levels,z)
'set dfile 1'
'set lev 'level
'c'
'movie statdplt.gs "'field' 'cint' geos ncep" -name std.geos-ncep.'field''level' -print -rotate'
'c'
z = z + 1
endwhile


*******************************************************
****                Vertical Means                 ****
*******************************************************

if( field = h ) ; cint = 2   ; endif
if( field = u ) ; cint = 0.2 ; endif
if( field = v ) ; cint = 0.2 ; endif
if( field = t ) ; cint = 0.1 ; endif

'set dfile 1'
'set x 1'
'sety'
'setz'
'c'
'movie statdpltz.gs "'field' 'cint' geos ncep" -name std.geos-ncep.'field'z -print -rotate'
'c'

n = n + 1
endwhile


return