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

'run getenv ZFLIP'
            zflip = result

*************************************
* Note: To Write Data in Reverse Order
*       run setenv ZFLIP ON
*************************************

'setx'
'sety'
'set gxout fwrite'

'getinfo tdim'
         tdim  = result
'getinfo nvars'
         nvars = result

t=1
while(t<=tdim)
  'set t 't

n=1
while(n<=nvars)

'run getvarz 'n
    var  = result
   name  = subwrd(result,1)
   zdim  = subwrd(result,2)

   if( zdim = 0 )
       say 'Writing Variable: 'name
      'd 'name
   else

      if(zflip != 'ON' )
         z=1
         while(z<=zdim)
            'set z 'z
            'getinfo level'
             lev = result
             say 'Writing Variable: 'name' for Level: 'lev
            'd 'name
         z=z+1
         endwhile
      else
         z=zdim
         while(z>=1)
            'set z 'z
            'getinfo level'
             lev = result
             say 'Writing Variable: 'name' for Level: 'lev
            'd 'name
         z=z-1
         endwhile
      endif

   endif
       say '  '


n=n+1
endwhile
t=t+1
endwhile

return