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

*************************************************************
*****                                                   *****
*****  Usage:  set_axlim name field1 field2 ...         *****
*****                                                   *****
*****  This function sets the AXLIM                     *****
*****  based on an arbitrary number of input fields.    *****
*****  If NAME exists as a file, AXLIM will             *****
*****  be read in, otherwise AXLIM will be computed.    *****
*****                                                   *****
*************************************************************

k = 1
while( k > 0 )
field.k = subwrd(args,k)
if( field.k = '' )
k = 1-k
else
k = k+1
endif
endwhile
k = -k

name = ''
n = k
while( n > 1 )
name = name' 'field.n
n = n-1
endwhile

file = subwrd(args,1)
       ioflag = sublin( read(file),1 )
    if(ioflag = 0)
    dum =          close(file)
    min =   sublin( read(file),2 )
    max =   sublin( read(file),2 )
    else
    min = ''
    max = ''
    endif

if( min = '' )
'getinfo latmin'
         latmin = result
'getinfo latmax'
         latmax = result
'getinfo lonmin'
         lonmin = result
'getinfo lonmax'
         lonmax = result
say 'Lons: 'lonmin' 'lonmax'   Lats: 'latmin' 'latmax

'minmax 'name
max = subwrd( result,1 )
min = subwrd( result,2 )
write(file,min)
write(file,max)
endif

dif = max-min
max = max + 0.25*dif
min = min - 0.25*dif

say 'File: 'file' setting axlim for field: 'name' with 'min' 'max
'set axlim 'min' 'max
return min' 'max