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/n2
c ifget(nfhist) = fields to allocate for fields array, and get from history
c nfget = number of fields to get from history (sum of ifget(ip) > 0)
c ixplt(nfhist) = indices to fields array 
c   (ixplt(ip)=0 if not getting field ip) 
c ifplt(nftot) = fields to plot
c ixgcm(nftot) = indices of fields to plot (not indices to any array)
c
      parameter (nftot=29,nfhist=17)
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,ixunvn,ixuivi,ixvpar,ixwpar,ixon2,ixn2o
      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)
c (vars with 0 at end are for cntr, others for pert)
      pointer(pglbui,glbui(imx,kmx,jmx)),(pglbui0,glbui0(imx,kmx,jmx)), 
     +       (pglbvi,glbvi(imx,kmx,jmx)),(pglbvi0,glbvi0(imx,kmx,jmx)),
     +       (pglbwi,glbwi(imx,kmx,jmx)),(pglbwi0,glbwi0(imx,kmx,jmx))
      common/drifts/ pglbui,pglbui0,pglbvi,pglbvi0,pglbwi,pglbwi0
c
c Selected grid:
      parameter(mxloc=15)
      character*16 locname(mxloc)
      common/sgrid/ spls(kmx),npls,gloc(2,mxloc),nloc,slon(imx),
     +  nlon,slat(jmx),nlat,shts(kmx),nhts,sslt(imx),nslt,locname,
     +  zphtlon(2,mxloc),izphtlon(2,mxloc),nzphtlon
      common/gcmgrid/ gcmlat(jmx),gcmlon(imx),gcmzp(kmx)
c
c Plot flags:
c (zmflag (special slon) means find zonal means in pltlon and pltlat)
c (r12flag means use local noon longitude for censatv (sat view proj)
c (gmflag means find global means for location line plots)
c pointer phtscale is allocated in getinp when nhtscale is known
c
      parameter(mxmap=25)
      character*8 labshort
      character*56 flab
      pointer(phtscale,htscale(1))
      common/pltcom/ nfget,ifget(nfhist),icont,ituv,icolor,
     +  cint(nftot),cmin(nftot),cmax(nftot),labshort(nftot),ilog(nftot),
     +  iframe,ifplt(nftot),ixfld(nfhist),zmflag,zprange(2),
     +  izprange(2),censatv(2,mxmap),nsatv,r12flag,phtscale,
     +  nhtscale,perimlat(mxmap),npollat,log10map,ipuv,gmflag,
     +  cenglb(2,mxmap),nglb,cenmolw(2,mxmap),nmolw,flab(nftot),idif,
     +  iperc(nftot),imap,islice,iutloc,lineloc,nfplt,iutline,
     +  nzprange,dipdec(imx,jmx,2),ibilin,uvmax,iutlat,lupert,lucntr,
     +  ludipdec
c
c History input:
      parameter(mxvols=15,mxtms=241)
      character*24 pertvol(mxvols), cntrvol(mxvols)
      character*56 dirascii,dircgm,mssascii
      logical istiegcm
      common/chinp/ pertvol,cntrvol
      common/input/ iden,mtpert(3,mxtms),mtcntr(3,mxtms),ntms,
     +  ut(mxtms),npvols,ncvols,ipvol,icvol,lumag,istiegcm,
     +  utinc(mxtms),luascii,iwrascii,dirascii,dircgm,mssascii
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)