c
c------------------------------------------------------------------
c Begin file /home/sting/foster/tigcmloc/tigcmloc.h
c------------------------------------------------------------------
c
c nfhist = number of fields on history (tn through n2)
c          (actually, n2 is not on hist, but is returned by getgcm)
c nftot = number of fields available to plot: tn -> poten, n2, plus
c         ui, vi, wi, fof2, hmf2, rho, o/(o2+n2)
c ifget(nfhist) = fields to get from history, and indices to fields array
c nfget = number of fields to get from history (sum of ifget(ip) > 0)
c ifplt(nftot) = fields to plot, and indices to plot arrays
c nfplt = number of fields to plot (sum of ifplt(ip) > 0)
c ixgcm(nftot) = indices of fields to plot (not indices to any array)
c
      parameter (nftot=25,nfhist=17)
      parameter(mxvols=10,mxtms=241,mxdays=10)
c
c Field indices:
      common /indices/ ixt,ixu,ixv,ixo2,ixo1,ixn4s,ixno,ixop,ixn2d,ixti,
     +  ixte,ixne,ixo2p,ixw,ixz,ixpot,ixn2,ixui,ixvi,ixwi,ixfof2,ixhmf2,
     +  ixrho,ixamr,ixon2r
      dimension ixgcm(nftot)
      equivalence(ixgcm,ixt)
c
c Grid sizes:
      parameter (jmx=36,imx=73,kmx=25)
      parameter (glat1=-87.5, glon1=-180., zp1=-7.,
     +           dlat=5.,dlon=5.,dzp=0.5)
c
c Ion drifts (pointers allocated space only if drifts are to be plotted)
      pointer(pglbui,glbui(imx,kmx,jmx)), (pglbvi,glbvi(imx,kmx,jmx)),
     +       (pglbwi,glbwi(imx,kmx,jmx))
      common/drifts/ pglbui,pglbvi,pglbwi
c
c Selected grid:
      parameter(mxloc=15)
      common/sgrid/ gloc(2,mxloc),ixloc(2,mxloc),nloc,zphtlon(2,mxloc),
     +  izphtlon(2,mxloc),nzphtlon,zphtline(kmx),izphtlin(kmx),nzphtlin,
     +  mutline(3,mxtms),nutline
      common/gcmgrid/ gcmlat(jmx),gcmlon(imx),gcmzp(kmx)
c
c Plot flags:
c (gmflag means find global means, zmflag means find zonal means)
c pointer phtscale is allocated when nhtscale is known
c
      character*8 labshort
      character*56 flab,dirascii,dircgm
      pointer(phtscale,htscale(1))
      common/pltcom/ nfget,ifget(nfhist),icolor,cint(nftot),cmin(nftot),
     +  cmax(nftot),labshort(nftot),logterp(nftot),iframe,ifplt(nftot),
     +  zmflag,zprange(2),izprange(2),r12flag,phtscale,iutvert,
     +  nhtscale,gmflag,flab(nftot),nfplt,logloc,logutlat,ibilin,
     +  iwrascii,dirascii,dircgm,idop,iyd(mxdays)
c
c History input:
      character*24 histvol(mxvols)
      logical istiegcm
      common/chinp/ histvol
      common/input/ iden,mtimes(3,mxtms),ntms,ut(mxtms),nhvols,luhist,
     +  ivol,lumag,istiegcm,utinc(mxtms),luascii
c
c Special value for conpack:
      parameter (cpspval=1.e36)
c
c degrees to radians:
      parameter (pi = 3.14159)
      parameter (dtr = pi/180.)
c
c radians to degrees:
      parameter (rtd = 180./pi)