c
c------------------------------------------------------------------
c
c nfhist = number of fields on history (tn through epot,ui,vi,wi,n2)
c nftot = number of fields available to plot: tn -> poten, n2 plus
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=37,nfhist=34)
      parameter(mxvols=12,mxtms=mxvols*24+1)
c
c Field indices:
      common/indices/ixt,ixu,ixv,ixo2,ixox,ixn4s,ixnoz,ixco,
     +               ixco2,ixh2o,ixh2,ixhox,ixop,ixch4,ixo21d,
     +               ixno2,ixno,ixo3,ixo1,ixoh,ixho2,ixh,
     +               ixn2d,ixti,ixte,ixne,ixo2p,ixw,ixz,ixpot,
     +               ixui,ixvi,ixwi,ixn2,ixe5577,ixe6300,ixeo200
      dimension ixgcm(nftot)
      equivalence(ixgcm,ixt)
c
c Grid sizes:
      parameter (jmx=36,imx=73,kmx=42)
      parameter (glat1=-87.5, glon1=-180., zp1=-15.5,
     +           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)
      character*16 locname(mxloc)
      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,locname,mutlines(3,mxtms),nutlines,
     +  modelhts,zphtslt(2,mxloc),izphtslt(2,mxloc),nzphtslt
      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 If doing ions, htions will be allocated if needed in height
c If doing e5577, hte5577 will be allocated if needed in height
c If doing e6300, hte6300 will be allocated if needed in height
c If doing doppler t,u,v,brightness, then idop5577=1 and pdop5577 will
c   be allocated for (ntms,nloc,4) and plotted from linezpht
c Emissions and doppler fields are found in mkfields
c
      character*8 labshort
      character*56 flab,dirascii,dircgm
      pointer(phtscale,htscale(1))
      pointer(phtions,htions(1))
      pointer(phte5577,hte5577(1)),(pdop5577,dop5577(1))
      pointer(phte6300,hte6300(1)),(pdop6300,dop6300(1))
      pointer(phteo200,hteo200(1)),(pdopo200,dopo200(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,nzprange,
     +  nhtscale,gmflag,flab(nftot),nfplt,logloc,logutlat,ibilin,
     +  iwrascii,dirascii,dircgm,labels,finc,ionvel,phtions,
     +  phte5577,pdop5577,idop5577, phte6300,pdop6300,idop6300,
     +  phteo200,pdopo200,idopo200,ifdep(nfhist,nftot-nfhist),locslt
c
c History input:
      character*24 histvol(mxvols)
      logical istimes
      common/chinp/ histvol
      common/input/ iden,mtimes(3,mxtms),ntms,ut(mxtms),nhvols,luhist,
     +  ivol,lumag,utinc(mxtms),luascii,irmvols
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)