!
! Grid dimensions and derived parameters:
! NLEV, ZTOP, and ZBOT are from #defines in dims.h:
! 
      integer,parameter ::
     | ZIMX=72,        ! number of longitudinal grid points
     | ZJMX=36,        ! number of latitudinal grid points
     | ZKMX=NLEV       ! number of vertical layers (see dims.h)
!
      real,parameter ::
     | ZST= ZTOP,      ! height of top of model (see dims.h)
     | ZSB= ZBOT,      ! height of bottom of model (see dims.h)
     | DZ =(ZTOP-ZBOT)/NLEV ! delta zp (0.5 or 0.25)
!
! 3/8/01: set zfldx,zphyx according to vertical resolution (dims.h)
! ZKBM:  Pressure level at which upward dynamo integrals begin (ZP -8.)
!
#if (NLEV==44)               
! NLEV==44 for -17 to +5 by 0.50
      integer,parameter ::
     | ZFLDX = 3045,
     | ZPHYX = 2926,
     | ZKBM  = 19
#elif (NLEV==88)             
! NLEV==88 for -17 to +5 by 0.25
      integer,parameter ::
     | ZFLDX = 6037,
     | ZPHYX = 5786,
     | ZKBM  = 37
#elif (NLEV==96)             
! NLEV==96 for -17 to +7 by 0.25
      integer,parameter ::
     | ZFLDX = 6581,
     | ZPHYX = 6306,
     | ZKBM  = 37
#endif
!
! Derived parameters based on above grid dimensions:
      integer,parameter ::
     | ZIMXP=ZIMX+4,
     | ZJMX2=ZJMX+ZJMX,
     | ZKMXP=ZKMX+1,
     | ZIXKX=ZIMXP*ZKMXP,
     | ZIMXP1=ZIMX+1,
     | ZJMXP1=ZJMX+1,
c
c ntndown = number of pressure levels for tndown (tn lower boundary)
c (see LOWBND and DLOWBND)
c (ntndown = 18 corresponds to zp -21.5 to -17.25 by .25)
c
     | ntndown=18,
!
! Derived parameters:
     | IMAXM=80,      ! number of apex longitude grid points
     | JMAXM=97,      ! number of apex latitude grid points
     | IMAXMP=IMAXM+1,
     | IMX0=IMAXM+1,
     | JMX0=(JMAXM+1)/2,
     | IMX1=(IMX0+1)/2,
     | JMX1=(JMX0+1)/2,
     | IMX2=(IMX1+1)/2,
     | JMX2=(JMX1+1)/2,
     | IMX3=(IMX2+1)/2,
     | JMX3=(JMX2+1)/2,
     | IMX4=(IMX3+1)/2,
     | JMX4=(JMX3+1)/2,
C      ****           SPACE NEEDED FOR DESCRETIZED COEFFICIENTS OF
C      ****             OF DYNAMO PDE AT ALL 5 LEVELS OF RESOLUTION
     | NCEE=10*IMX0*JMX0+9*(IMX1*JMX1+IMX2*JMX2+IMX3*JMX3+IMX4*JMX4),
     | JMAXMH=(JMAXM+1)/2,
C      ****           COEFFICIENTS ARE STORED IN 1-D ARRAY CEE(NCEE).
C      ****             FOLLOWING PARAMETERS ARE POINTERS TO THE
C      ****             BEGINNING OF THE COEFFICIENTS FOR EACH LEVEL
C      ****             OF RESOLUTION.
     | NC0=1,
     | NC1=NC0+10*IMX0*JMX0,
     | NC2=NC1+9*IMX1*JMX1,
     | NC3=NC2+9*IMX2*JMX2,
     | NC4=NC3+9*IMX3*JMX3
!
      integer,parameter :: IMAXG=ZIMX,JMAXG=ZJMX,IMAXGP=IMAXG+1,
     |                     JMAXGP=JMAXG+1
!
! Special values:
      real,parameter :: spval=1.e36
      integer,parameter  :: ispval = 9999
!
! For modules hist_mod and input_mod:
      character(len=8),parameter :: tgcm_version = 'tgcm24  '
!
! tgcm_name is either "tiegcm" or "time-gcm"
      character(len=8),parameter :: tgcm_name    = 'time-gcm'
!
      integer,parameter :: 
     |  nflds=69,  ! number of prognostics for current and previous iteration
     |  nphys=66,  ! number of physics fields
     |  ntape=51,  ! number of fields written to tape
     |  ndisk=51,  ! number of fields written to disk
     |  mxday=366,         ! max allowed calendar day
     |  mxhvols = 100,     ! max number of output history files
     |  mxseries = 10,     ! max number of time series for primary histories
     |  mxseries_sech = 50,! max number of secondary history time series
     |  mxfsech = 50       ! max number of fields on secondary histories
!