* +-======-+ 
*  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 hstat (args)
field = subwrd (args,1)

'numargs  'args
 numargs = result

        num = 0
while ( num < numargs )
        num = num + 1

if( subwrd(args,num) = '-desc'   ) ; DESC = subwrd(args,num+1) ; endif

endwhile

'getinfo pagex'
         pagex = result
     if( pagex = 8.5 ) ; 'setenv ORIENTATION PORTRAIT'  ; endif
     if( pagex = 11  ) ; 'setenv ORIENTATION LANDSCAPE' ; endif

'getinfo level'
         level = result

* Initialize Plot Values
* ----------------------
    CCOLS = '59   57   55   47   44   37   36   34   33    0     21   22   23   24   25   26   27   28   29'
    CLAB  =  off
    scale =  1

if( field = "p" )
    name  = "Sea-Level Pressure"
    unit  = "(mb)"
    label = "slp"
    cint  = 4
        rbrange = '970 1025'
        DLEVS = '-4.5 -4 -3.5 -3 -2.5 -2 -1.5 -1 -.5 .5 1 1.5 2 2.5 3 3.5 4 4.5'
        DCINT = 0.5
        CLEVS = '-13.5 -12.0 -10.5 -9.0 -7.5 -6.0 -4.5 -3.0 -1.5 1.5 3.0 4.5 6.0 7.5 9.0 10.5 12.0 13.5'
        CCINT = 1.5
        RCOLS = '0  50  42  44  46  48  39  37  36  34  32  31  21  22  24  25  26  27  28   29'
        RLEVS = '0.5 1 1.5 2 2.5 3 3.5  4 4.5  5 5.5  6  7  8   9  10 11 12'
        RCINT = 0.5
endif

if( field = "h" )
    name  = "Heights"
    unit  = "(m)"
    label = "hght"
    RCOLS = '0  50  42  44  46  48  39  37  36  34  32  31  21  22  24  25  26  27  28   29'
    RLEVS = '1  2  3  4  5  6  7  8  9  10 11   12  14  16  18  20   22   24'
    RCINT = 1
    RLEVS = '2  4  6  8  10  15  20  25  30  40  50  60  70  80  90  100  120  140'
    RCINT = 2
    CLEVS = '-45 -40 -35 -30 -25 -20 -15 -10 -5 5 10 15 20 25 30 35 40 45'
    CCINT = 5
    CLEVS = '-90 -80 -70 -60 -50 -40 -30 -20 -10 10 20 30 40 50 60 70 80 90'
    CCINT = 10
    CLEVS = '-135 -120 -105 -90 -75 -60 -45 -30 -15 15 30 45 60 75 90 105 120 135'
    CCINT = 15
    if( level = 1000 )
        cint  = 20
        rbrange = '-200 200'
        DLEVS = '-36 -32 -28 -24 -20 -16 -12 -8 -4 4 8 12 16 20 24 28 32 36'
        DCINT = 4
    endif
    if( level = 850 )
        cint  = 40
        rbrange = '1000 1600'
        DLEVS = '-36 -32 -28 -24 -20 -16 -12 -8 -4 4 8 12 16 20 24 28 32 36'
        DCINT = 4
    endif
    if( level = 700 )
        cint  = 50
        rbrange = '2800 3100'
        DLEVS = '-36 -32 -28 -24 -20 -16 -12 -8 -4 4 8 12 16 20 24 28 32 36'
        DCINT = 4
    endif
    if( level = 500 )
        cint  = 25
        rbrange = '-200 200'
        DLEVS = '-18 -16 -14 -12 -10 -8 -6 -4 -2 2 4 6 8 10 12 14 16 18'
        DCINT = 2
        DLEVS = '-36 -32 -28 -24 -20 -16 -12 -8 -4 4 8 12 16 20 24 28 32 36'
        DCINT = 4
        cint  = 60
        rbrange = '4860 5880'
        DLEVS = '-36 -32 -28 -24 -20 -16 -12 -8 -4 4 8 12 16 20 24 28 32 36'
        DCINT = 4
    endif
    if( level = 400 )
        cint  = 60
        rbrange = '6400 7600'
        DLEVS = '-36 -32 -28 -24 -20 -16 -12 -8 -4 4 8 12 16 20 24 28 32 36'
        DCINT = 4
    endif
    if( level = 300 )
        cint  = 60
        rbrange = '8200 9600'
        DLEVS = '-36 -32 -28 -24 -20 -16 -12 -8 -4 4 8 12 16 20 24 28 32 36'
        DCINT = 4
    endif
    if( level = 250 )
        cint  = 60
        rbrange = '9400 10800'
        DLEVS = '-45 -40 -35 -30 -25 -20 -15 -10 -5 5 10 15 20 25 30 35 40 45'
        DCINT = 5
    endif
    if( level = 200 )
        cint  = 60
        rbrange = '10800 12400'
        DLEVS = '-45 -40 -35 -30 -25 -20 -15 -10 -5 5 10 15 20 25 30 35 40 45'
        DCINT = 5
    endif
    if( level = 150 )
        cint  = 100
        rbrange = '12400 14200'
        DLEVS = '-45 -40 -35 -30 -25 -20 -15 -10 -5 5 10 15 20 25 30 35 40 45'
        DCINT = 5
    endif
    if( level = 100 )
        cint  = 100
        rbrange = '14800 16600'
        DLEVS = '-45 -40 -35 -30 -25 -20 -15 -10 -5 5 10 15 20 25 30 35 40 45'
        DCINT = 5
    endif
endif

if( field = "u" )
    name  = "U-Wind"
    unit  = "(m/sec)"
    label = "uwnd"
    RCOLS = '0  50  42  44  46  48  39  37  36  34  32  31  21  22  24  25  26  27  28   29'
    RLEVS = '2  3  4   5   6   7   8   9  10  11  12  13  14  15   16   17   18 20'
    RCINT = 2
    CLEVS = '-18 -16 -14 -12 -10 -8 -6 -4 -2 2  4 6   8 10 12 14 16 18'
    CCINT = 2.0
    if( level = 1000 )
        cint  = 2
        rbrange = '-10 10'
        DLEVS = '-4.5 -4 -3.5 -3 -2.5 -2 -1.5 -1 -.5 .5 1 1.5 2 2.5 3 3.5 4 4.5'
        DCINT = 0.5
    endif
    if( level = 850 )
        cint  = 5
        rbrange = '-10 30'
        DLEVS = '-4.5 -4 -3.5 -3 -2.5 -2 -1.5 -1 -.5 .5 1 1.5 2 2.5 3 3.5 4 4.5'
        DCINT = 0.5
    endif
    if( level = 700 )
        cint  = 5
        rbrange = '-10 30'
        DLEVS = '-4.5 -4 -3.5 -3 -2.5 -2 -1.5 -1 -.5 .5 1 1.5 2 2.5 3 3.5 4 4.5'
        DCINT = 0.5
    endif
    if( level = 500 )
        cint  = 5
        rbrange = '-10 30'
        DLEVS = '-4.5 -4 -3.5 -3 -2.5 -2 -1.5 -1 -.5 .5 1 1.5 2 2.5 3 3.5 4 4.5'
        DCINT = 0.5
    endif
    if( level = 400 )
        cint  = 5
        rbrange = '-10 30'
        DLEVS = '-4.5 -4 -3.5 -3 -2.5 -2 -1.5 -1 -.5 .5 1 1.5 2 2.5 3 3.5 4 4.5'
        DCINT = 0.5
    endif
    if( level = 300 )
        cint  = 5
        rbrange = '-10 35'
        DLEVS = '-4.5 -4 -3.5 -3 -2.5 -2 -1.5 -1 -.5 .5 1 1.5 2 2.5 3 3.5 4 4.5'
        DCINT = 0.5
    endif
    if( level = 250 )
        cint  = 5
        rbrange = '-10 50'
        DLEVS = '-4.5 -4 -3.5 -3 -2.5 -2 -1.5 -1 -.5 .5 1 1.5 2 2.5 3 3.5 4 4.5'
        DCINT = 0.5
    endif
    if( level = 200 )
        cint  = 5
        rbrange = '-10 60'
        DLEVS = '-4.5 -4 -3.5 -3 -2.5 -2 -1.5 -1 -.5 .5 1 1.5 2 2.5 3 3.5 4 4.5'
        DCINT = 0.5
    endif
    if( level = 150 )
        cint  = 5
        rbrange = '-10 50'
        DLEVS = '-9 -8 -7 -6 -5 -4 -3 -2 -1 1 2 3 4 5 6 7 8 9'
        DCINT = 1.0
    endif
    if( level = 100 )
        cint  = 5
        rbrange = '-10 40'
        DLEVS = '-9 -8 -7 -6 -5 -4 -3 -2 -1 1 2 3 4 5 6 7 8 9'
        DCINT = 1.0
    endif
endif

if( field = "v" )
    name  = "V-Wind"
    unit  = "(m/sec)"
    label = "vwnd"
    RCOLS = '0  50  42  44  46  48  39  37  36  34  32  31  21  22  24  25  26  27  28   29'
    RLEVS = '2  3  4   5   6   7   8   9  10  11  12  13  14  15   16   17   18 20'
    RCINT = 2
    if( level = 1000 )
        cint  = 2
        rbrange = '-10 10'
        DLEVS = '-4.5 -4 -3.5 -3 -2.5 -2 -1.5 -1 -.5 .5 1 1.5 2 2.5 3 3.5 4 4.5'
        DCINT = 0.5
        CLEVS = '-18 -16 -14 -12 -10 -8 -6 -4 -2 2  4 6   8 10 12 14 16 18'
        CCINT = 2.0
    endif
    if( level = 850 )
        cint  = 2
        rbrange = '-10 10'
        DLEVS = '-4.5 -4 -3.5 -3 -2.5 -2 -1.5 -1 -.5 .5 1 1.5 2 2.5 3 3.5 4 4.5'
        DCINT = 0.5
        CLEVS = '-18 -16 -14 -12 -10 -8 -6 -4 -2 2  4 6   8 10 12 14 16 18'
        CCINT = 2.0
    endif
    if( level = 700 )
        cint  = 2
        rbrange = '-10 10'
        DLEVS = '-4.5 -4 -3.5 -3 -2.5 -2 -1.5 -1 -.5 .5 1 1.5 2 2.5 3 3.5 4 4.5'
        DCINT = 0.5
        CLEVS = '-18 -16 -14 -12 -10 -8 -6 -4 -2 2  4 6   8 10 12 14 16 18'
        CCINT = 2.0
    endif
    if( level = 500 )
        cint  = 2
        rbrange = '-10 10'
        DLEVS = '-4.5 -4 -3.5 -3 -2.5 -2 -1.5 -1 -.5 .5 1 1.5 2 2.5 3 3.5 4 4.5'
        DCINT = 0.5
        CLEVS = '-18 -16 -14 -12 -10 -8 -6 -4 -2 2  4 6   8 10 12 14 16 18'
        CCINT = 2.0
    endif
    if( level = 400 )
        cint  = 2
        rbrange = '-10 10'
        DLEVS = '-4.5 -4 -3.5 -3 -2.5 -2 -1.5 -1 -.5 .5 1 1.5 2 2.5 3 3.5 4 4.5'
        DCINT = 0.5
        CLEVS = '-18 -16 -14 -12 -10 -8 -6 -4 -2 2  4 6   8 10 12 14 16 18'
        CCINT = 2.0
    endif
    if( level = 300 )
        cint  = 2
        rbrange = '-10 10'
        DLEVS = '-4.5 -4 -3.5 -3 -2.5 -2 -1.5 -1 -.5 .5 1 1.5 2 2.5 3 3.5 4 4.5'
        DCINT = 0.5
        CLEVS = '-18 -16 -14 -12 -10 -8 -6 -4 -2 2  4 6   8 10 12 14 16 18'
        CCINT = 2.0
    endif
    if( level = 250 )
        cint  = 2
        rbrange = '-10 10'
        DLEVS = '-4.5 -4 -3.5 -3 -2.5 -2 -1.5 -1 -.5 .5 1 1.5 2 2.5 3 3.5 4 4.5'
        DCINT = 0.5
        CLEVS = '-18 -16 -14 -12 -10 -8 -6 -4 -2 2  4 6   8 10 12 14 16 18'
        CCINT = 2.0
    endif
    if( level = 200 )
        cint  = 2
        rbrange = '-10 10'
        DLEVS = '-4.5 -4 -3.5 -3 -2.5 -2 -1.5 -1 -.5 .5 1 1.5 2 2.5 3 3.5 4 4.5'
        DCINT = 0.5
        CLEVS = '-18 -16 -14 -12 -10 -8 -6 -4 -2 2  4 6   8 10 12 14 16 18'
        CCINT = 2.0
    endif
    if( level = 150 )
        cint  = 3
        rbrange = '-10 10'
        DLEVS = '-4.5 -4 -3.5 -3 -2.5 -2 -1.5 -1 -.5 .5 1 1.5 2 2.5 3 3.5 4 4.5'
        DCINT = 0.5
        CLEVS = '-18 -16 -14 -12 -10 -8 -6 -4 -2 2  4 6   8 10 12 14 16 18'
        CCINT = 2.0
    endif
    if( level = 100 )
        cint  = 3
        rbrange = '-10 10'
        DLEVS = '-4.5 -4 -3.5 -3 -2.5 -2 -1.5 -1 -.5 .5 1 1.5 2 2.5 3 3.5 4 4.5'
        DCINT = 0.5
        CLEVS = '-18 -16 -14 -12 -10 -8 -6 -4 -2 2  4 6   8 10 12 14 16 18'
        CCINT = 2.0
    endif
endif

if( field = "t" )
    name  = "Temperature"
    unit  = "(K)"
    label = "tmpu"
    RCOLS = '0  50  42  44  46  48  39  37  36  34  32  31  21  22  24  25  26  27  28   29'
    RLEVS = '.2  .4  .6 .8  1  1.2 1.4 1.6 1.8 2 2.3 2.5 2.7  3   3.5  4   4.5  5.0  '
    RCINT = 0.2
    CLEVS = '-4.5 -4 -3.5 -3 -2.5 -2 -1.5 -1 -.5 .5 1 1.5 2 2.5 3 3.5 4 4.5'
    CCINT = 0.5
    if( level = 1000 )
        cint  = 2
        rbrange = '270 300'
        DLEVS = '-1.8 -1.6 -1.4 -1.2 -1.0 -.8 -.6 -.4 -.2 .2 .4 .6 .8 1.0 1.2 1.4 1.6 1.8'
        DCINT = 0.2
    endif
    if( level = 850 )
        cint  = 2
        rbrange = '250 290'
        DLEVS = '-3.6 -3.2 -2.8 -2.4 -2.0 -1.6 -1.2 -.8 -.4 .4 .8 1.2 1.6 2.0 2.4 2.8 3.2 3.6'
        DCINT = 0.4
        CLEVS = '-4.5 -4 -3.5 -3 -2.5 -2 -1.5 -1 -.5 .5 1 1.5 2 2.5 3 3.5 4 4.5'
        CCINT = 0.5
    endif
    if( level = 700 )
        cint  = 2
        rbrange = '250 290'
        DLEVS = '-3.6 -3.2 -2.8 -2.4 -2.0 -1.6 -1.2 -.8 -.4 .4 .8 1.2 1.6 2.0 2.4 2.8 3.2 3.6'
        DCINT = 0.4
        CLEVS = '-4.5 -4 -3.5 -3 -2.5 -2 -1.5 -1 -.5 .5 1 1.5 2 2.5 3 3.5 4 4.5'
        CCINT = 0.5
    endif
    if( level = 500 )
        cint  = 2
        rbrange = '230 270'
        DLEVS = '-1.8 -1.6 -1.4 -1.2 -1.0 -.8 -.6 -.4 -.2 .2 .4 .6 .8 1.0 1.2 1.4 1.6 1.8'
        DCINT = 0.2
        CLEVS = '-4.5 -4 -3.5 -3 -2.5 -2 -1.5 -1 -.5 .5 1 1.5 2 2.5 3 3.5 4 4.5'
        CCINT = 0.5
    endif
    if( level = 400 )
        cint  = 2
        rbrange = '220 260'
        DLEVS = '-1.8 -1.6 -1.4 -1.2 -1.0 -.8 -.6 -.4 -.2 .2 .4 .6 .8 1.0 1.2 1.4 1.6 1.8'
        DCINT = 0.2
        CLEVS = '-4.5 -4 -3.5 -3 -2.5 -2 -1.5 -1 -.5 .5 1 1.5 2 2.5 3 3.5 4 4.5'
        CCINT = 0.5
    endif
    if( level = 300 )
        cint  = 2
        rbrange = '210 244'
        DLEVS = '-1.8 -1.6 -1.4 -1.2 -1.0 -.8 -.6 -.4 -.2 .2 .4 .6 .8 1.0 1.2 1.4 1.6 1.8'
        DCINT = 0.2
        CLEVS = '-4.5 -4 -3.5 -3 -2.5 -2 -1.5 -1 -.5 .5 1 1.5 2 2.5 3 3.5 4 4.5'
        CCINT = 0.5
    endif
    if( level = 250 )
        cint  = 2
        rbrange = '206 232'
        DLEVS = '-1.8 -1.6 -1.4 -1.2 -1.0 -.8 -.6 -.4 -.2 .2 .4 .6 .8 1.0 1.2 1.4 1.6 1.8'
        DCINT = 0.2
        CLEVS = '-4.5 -4 -3.5 -3 -2.5 -2 -1.5 -1 -.5 .5 1 1.5 2 2.5 3 3.5 4 4.5'
        CCINT = 0.5
    endif
    if( level = 200 )
        cint  = 2
        rbrange = '200 220'
        DLEVS = '-2.7 -2.4 -2.1 -1.8 -1.5 -1.2 -.9 -.6 -.3 .3 .6 .9 1.2 1.5 1.8 2.1 2.4 2.7'
        DCINT = 0.3
        CLEVS = '-4.5 -4 -3.5 -3 -2.5 -2 -1.5 -1 -.5 .5 1 1.5 2 2.5 3 3.5 4 4.5'
        CCINT = 0.5
    endif
    if( level = 150 )
        cint  = 2
        rbrange = '200 220'
        DLEVS = '-2.7 -2.4 -2.1 -1.8 -1.5 -1.2 -.9 -.6 -.3 .3 .6 .9 1.2 1.5 1.8 2.1 2.4 2.7'
        DCINT = 0.3
        CLEVS = '-4.5 -4 -3.5 -3 -2.5 -2 -1.5 -1 -.5 .5 1 1.5 2 2.5 3 3.5 4 4.5'
        CCINT = 0.5
    endif
    if( level = 100 )
        cint  = 2
        rbrange = '195 225'
        DLEVS = '-2.7 -2.4 -2.1 -1.8 -1.5 -1.2 -.9 -.6 -.3 .3 .6 .9 1.2 1.5 1.8 2.1 2.4 2.7'
        DCINT = 0.3
        CLEVS = '-4.5 -4 -3.5 -3 -2.5 -2 -1.5 -1 -.5 .5 1 1.5 2 2.5 3 3.5 4 4.5'
        CCINT = 0.5
    endif
endif

if( field = "q" )
    name  = "Specific Humidity"
    unit  = "(g/g)"
    label = "sphu"
    scale = 1000
    RCOLS = '0  50  42  44  46  48  39  37  36  34  32  31  21  22  24  25  26  27  28   29'
    RLEVS = '.2  .4  .6 .8  1  1.2 1.4 1.6 1.8 2 2.3 2.5 2.7  3   3.5  4   4.5  5.0  '
    RCINT = 0.2
    if( level = 1000 )
        cint  = 1
        rbrange = '1 16'
        CLEVS = '-4.5 -4 -3.5 -3 -2.5 -2 -1.5 -1 -.5 .5 1 1.5 2 2.5 3 3.5 4 4.5'
        CCINT = 0.5
        DLEVS = '-3.6 -3.2 -2.8 -2.4 -2.0 -1.6 -1.2 -.8 -.4 .4 .8 1.2 1.6 2.0 2.4 2.8 3.2 3.6'
        DCINT = 0.4
    endif
    if( level = 850 )
        cint  = 1
        rbrange = '1 13'
        CLEVS = '-4.5 -4 -3.5 -3 -2.5 -2 -1.5 -1 -.5 .5 1 1.5 2 2.5 3 3.5 4 4.5'
        CCINT = 0.5
        DLEVS = '-3.6 -3.2 -2.8 -2.4 -2.0 -1.6 -1.2 -.8 -.4 .4 .8 1.2 1.6 2.0 2.4 2.8 3.2 3.6'
        DCINT = 0.4
    endif
    if( level = 700 )
        cint  = 1
        rbrange = '1 10'
        CLEVS = '-4.5 -4 -3.5 -3 -2.5 -2 -1.5 -1 -.5 .5 1 1.5 2 2.5 3 3.5 4 4.5'
        CCINT = 0.5
        DLEVS = '-3.6 -3.2 -2.8 -2.4 -2.0 -1.6 -1.2 -.8 -.4 .4 .8 1.2 1.6 2.0 2.4 2.8 3.2 3.6'
        DCINT = 0.4
    endif
    if( level = 500 )
        cint  = 0.5
        rbrange = '.5 4'
        CLEVS = '-2.7 -2.4 -2.1 -1.8 -1.5 -1.2 -.9 -.6 -.3 .3 .6 .9 1.2 1.5 1.8 2.1 2.4 2.7'
        CCINT = 0.3
        DLEVS = '-2.7 -2.4 -2.1 -1.8 -1.5 -1.2 -.9 -.6 -.3 .3 .6 .9 1.2 1.5 1.8 2.1 2.4 2.7'
        DCINT = 0.3
        RLEVS = '.1  .2  .3  .4  .5  .6  .7  .8  .9  1.0 1.1   1.2  1.3  1.4  1.5  1.6   1.8   2.0'
        RCINT = 0.1
    endif
    if( level = 400 )
        cint  = 0.2
        rbrange = '.2 2'
        CLEVS = '-1.8 -1.6 -1.4 -1.2 -1.0 -.8 -.6 -.4 -.2 .2 .4 .6 .8 1.0 1.2 1.4 1.6 1.8'
        CCINT = 0.2
        DLEVS = '-.9 -.8 -.7 -.6 -.5 -.4 -.3 -.2 -.1 .1 .2 .3 .4 .5 .6 .7 .8 .9'
        DCINT = 0.1
        RLEVS = '.1  .2  .3  .4  .5  .6  .7  .8  .9  1.0 1.1   1.2  1.3  1.4  1.5  1.6   1.8   2.0'
        RCINT = 0.1
    endif
    if( level = 300 )
        cint  = 0.1
        rbrange = '.1 1'
        CLEVS = '-.9 -.8 -.7 -.6 -.5 -.4 -.3 -.2 -.1 .1 .2 .3 .4 .5 .6 .7 .8 .9'
        CCINT = 0.1
        DLEVS = '-.9 -.8 -.7 -.6 -.5 -.4 -.3 -.2 -.1 .1 .2 .3 .4 .5 .6 .7 .8 .9'
        DCINT = 0.1
        RLEVS = '.1  .2  .3  .4  .5  .6  .7  .8  .9  1.0 1.1   1.2  1.3  1.4  1.5  1.6   1.8   2.0'
        RCINT = 0.1
    endif
    if( level = 250 )
        cint  = 0.1
        rbrange = '.1 1'
        CLEVS = '-.9 -.8 -.7 -.6 -.5 -.4 -.3 -.2 -.1 .1 .2 .3 .4 .5 .6 .7 .8 .9'
        CCINT = 0.1
        DLEVS = '-.9 -.8 -.7 -.6 -.5 -.4 -.3 -.2 -.1 .1 .2 .3 .4 .5 .6 .7 .8 .9'
        DCINT = 0.1
        RLEVS = '.1  .2  .3  .4  .5  .6  .7  .8  .9  1.0 1.1   1.2  1.3  1.4  1.5  1.6   1.8   2.0'
        RCINT = 0.1
    endif
    if( level = 200 )
        cint  = 0.1
        rbrange = '.1 1'
        CLEVS = '-.9 -.8 -.7 -.6 -.5 -.4 -.3 -.2 -.1 .1 .2 .3 .4 .5 .6 .7 .8 .9'
        CCINT = 0.1
        DLEVS = '-.9 -.8 -.7 -.6 -.5 -.4 -.3 -.2 -.1 .1 .2 .3 .4 .5 .6 .7 .8 .9'
        DCINT = 0.1
        RLEVS = '.1  .2  .3  .4  .5  .6  .7  .8  .9  1.0 1.1   1.2  1.3  1.4  1.5  1.6   1.8   2.0'
        RCINT = 0.1
    endif
    if( level = 150 )
        cint  = 0.1
        rbrange = '.1 1'
        CLEVS = '-.9 -.8 -.7 -.6 -.5 -.4 -.3 -.2 -.1 .1 .2 .3 .4 .5 .6 .7 .8 .9'
        CCINT = 0.1
        DLEVS = '-.9 -.8 -.7 -.6 -.5 -.4 -.3 -.2 -.1 .1 .2 .3 .4 .5 .6 .7 .8 .9'
        DCINT = 0.1
        RLEVS = '.1  .2  .3  .4  .5  .6  .7  .8  .9  1.0 1.1   1.2  1.3  1.4  1.5  1.6   1.8   2.0'
        RCINT = 0.1
    endif
    if( level = 100 )
        cint  = 0.1
        rbrange = '.1 1'
        CLEVS = '-.9 -.8 -.7 -.6 -.5 -.4 -.3 -.2 -.1 .1 .2 .3 .4 .5 .6 .7 .8 .9'
        CCINT = 0.1
        DLEVS = '-.9 -.8 -.7 -.6 -.5 -.4 -.3 -.2 -.1 .1 .2 .3 .4 .5 .6 .7 .8 .9'
        DCINT = 0.1
        RLEVS = '.1  .2  .3  .4  .5  .6  .7  .8  .9  1.0 1.1   1.2  1.3  1.4  1.5  1.6   1.8   2.0'
        RCINT = 0.1
    endif
endif

'getinfo lat'
         lat = result
'getinfo lon'
         lon = result

'set mproj latlon'


* Forecast (MEAN)
* ---------------
'set vpage 0 11 0 8.5'
'set parea 0.4 5.4 4.0 8.5'
'set grads off'
'set grid  off'
'set gxout shaded'
'set ccols 92 0'
'set clevs 0.5'
'd mask'
'set clab 'CLAB
'set gxout contour'
'set ccolor rainbow'
'set cint 'cint
'set rbrange 'rbrange
'd 'field'fmean*'scale


* Forecast-Climatology (MEAN)
* ---------------------------
'define diff = regrid2( 'field'fmean-'field'cmean, .25, .25, bs_p1, 'lon', 'lat' )*'scale
'set parea off'
'set vpage 0 11 0 8.5'
'set parea 5.9 10.9  4.0 8.5'
'set grads off'
'set grid  off'
'set gxout shaded'
'set ccols 92 0'
'set clevs 0.5'
'd mask'
'set gxout shaded'
'set CCOLS 'CCOLS
'set CLEVS 'CLEVS
'd maskout( diff,abs(diff)-'CCINT')'
'cbarn -scalex 0.55 -scaley 0.4 -xmid 8.45 -ymid 4.4 -snum 0.5'


* Forecast-Analysis (MEAN)
* ------------------------
'define diff = regrid2( 'field'fmean-'field'amean, .25, .25, bs_p1, 'lon', 'lat' )*'scale

if(1>2)
        'stats 'diff
         MEAN  = subwrd(result,1)
         STDV  = subwrd(result,2)
         aGLOB = (MEAN + 2*STDV)/6
'd abs(' aGLOB ')'
         aGLOB = subwrd(result,4)
'getint 'aGLOB
      cint = result
 if( aGLOB = 0 ) ; cint = 1 ; endif
         n = 0
     while( cint = 0 )
     say 'N: 'n'  CINT: 'cint'  aGLOB: 'aGLOB
*    say 'Hit ENTER to continue'
*    pull flag
            n = n + 1
        aGLOB = 10*aGLOB
       'getint 'aGLOB
         cint = result
      endwhile
      'd pow(10,'n')'
      fact = subwrd(result,4)
      cint = cint/fact
      say '(F-A Mean):  MEAN: 'MEAN'  STDV: 'STDV'  CINT: 'cint'  FACT: 'fact
clevs = ''
n = -9
while ( n <= 9 )
  val = n * cint
  val = n
if( n != 0 ) ; clevs = clevs' 'val ; endif
    n = n + 1
endwhile

'getinfo time'
        ntime = result
if( ntime = 1 ) ; 'setenv "CINT" 'cint ; endif
'getenv CINT'
     oldcint = result
if( oldcint != "NULL" )
    if( cint < oldcint )
        cint = oldcint
    endif
endif
DCINT = cint
DLEVS = clevs
endif

'set parea off'
'set vpage 0 11 0 8.5'
'set parea 0.4 5.4  0.2 4.7'
'set grads off'
'set grid  off'
'set gxout shaded'
'set ccols 92 0'
'set clevs 0.5'
'd  mask'
'set gxout shaded'
'set CCOLS 'CCOLS
'set CLEVS 'DLEVS
'd maskout( diff,abs(diff)-'DCINT')'
*'d maskout( diff/'DCINT',abs(diff/'DCINT')-1)'
'cbarn -scale 0.55 -xmid 2.95'


* Forecast-Analysis (STD)
* -----------------------
'define diff = regrid2( 'field'std, .25, .25, bs_p1, 'lon', 'lat' )*'scale
        'stats 'diff
         MEAN  = subwrd(result,1)
         STDV  = subwrd(result,2)
         aGLOB = (MEAN + 3*STDV)/18
'getint 'aGLOB
      cint = result
 if( aGLOB = 0 ) ; cint = 1 ; endif
         n = 0
     while( cint = 0 )
     say 'N: 'n'  CINT: 'cint'  aGLOB: 'aGLOB
*    say 'Hit ENTER to continue'
*    pull flag
            n = n + 1
        aGLOB = 10*aGLOB
       'getint 'aGLOB
         cint = result
      endwhile
      'd pow(10,'n')'
      fact = subwrd(result,4)
      cint = cint/fact
      say '(F-A STD):  MEAN: 'MEAN'  STDV: 'STDV'  CINT: 'cint'  FACT: 'fact
clevs = ''
n = 1
while ( n <= 19 )
val.n = n * cint
    n = n + 1
endwhile
n = 1
while ( n <= 19 )
        rat = val.n / val.19
'd pow('rat',2.0)'
        rat = subwrd(result,4)
val.n = fact * rat * val.19
'getint 'val.n
         val.n = result / fact
    n = n + 1
endwhile
n = 1
while ( n <= 19 )
clevs = clevs' 'val.n
    n = n + 1
endwhile
*RCINT = cint
*RLEVS = clevs

'set parea off'
'set vpage 0 11 0 8.5'
'set parea 5.9 10.9  0.2 4.7'
'set grads off'
'set grid  off'
'set gxout shaded'
'set ccols 92 0'
'set clevs 0.5'
'd mask'
'set CCOLS 'RCOLS
'set CLEVS 'RLEVS
'd maskout( diff,abs(diff)-'RCINT')'
'cbarn -sbar 0.55 -snum 0.45 -xmid 8.45'


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

'getinfo month'
         month = result
'getinfo year'
         year  = result
'getinfo time'
         time  = result
'getinfo tinc'
         tinc  = result
         hour  = (time-1)*tinc

'draw string 5.5  8.4 'DESC'   'month' 'year
'draw string 5.5  8.12 Field: 'field'  Level: 'level' mb   Hour: 'hour

'set strsiz .10'
'draw string 2.8  7.85 Forecast'
'draw string 8.2  7.85 Forecast-Climatology'
'draw string 2.8  4.10 Mean(Forecast-Analysis)   CINT: 'DCINT
'draw string 8.2  4.10 Std.Dev(Forecast-Analysis)'