c
c tracsat.h:
c
c nfhist = 31 -> tn,un,vn,o2,ox,n4s,noz,co,co2,h2o,h2,hox,o+,ch4,o21d,
c                no2,no,o3,o1,oh,ho2,h,n2d,ti,te,ne,o2+,w,z,epot,n2
c ifget(nfhist) -> fields to get from hist
c nfget -> number of fields to get from hist (sum of ifget(i) > 0)
c
c nftot = nfhist+4 = 37 -> nfhist + ui,vi,wi,ver6300,ver15m,ver5.3m
c ifplt(nftot) -> fields to plot
c nfplt -> number of fields to plot (sum of ifplt(i) > 0)
c
c
      parameter(jmx=36,imx=73,nfhist=31,glat1=-87.5,glon1=-180.,
     +  dlat=5.,dlon=5.,dzp=0.5,spval=1.e36,zp1times=-17.,zp1tigcm=-7.,
     +  ispval=9999,rminlog=1.e-20,nftot=nfhist+6)
      character*8 flab(nftot),model
c
c gcmzp = pressure array (will be allocated and defined after first
c         call to gettgcm, when we know kmx (ti(e)gcm or timesgcm)
c dipdec is used in ion drift calc
c
      pointer(pgcmzp,gcmzp(1)),(pdipdec,dipdec(imx,jmx,2))
      pointer(phtscale,htscale(1))
      character*80 unixdat,unixcgm
      common/gcm/ gcmlat(jmx),gcmlon(imx),pgcmzp,kmx,ifplt(nftot),
     +  nfplt,flab,logint(nftot),ionvel,ifget(nfhist),nfget,pdipdec,
     +  inttime,phtscale,nhtscale,iwrdat,unixdat,unixcgm,iydgcm,model,
     +  cmin(nftot),cmax(nftot),cint(nftot)
c
c History input:
c
      parameter(mxvols=10,mxtms=241)
      character*24 histvol(mxvols)
      logical istimes,isdyn
      common/input/ histvol,iden,mt(3,mxtms),ntms,ut(mxtms),
     +  nhvols,luhist,ivol,lumag,istimes,isdyn,ludipdec,ludat
      common/ixgcm/
     +  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   ,
     +  ixn2    ,ixui    ,ixvi    ,ixwi    ,ixver   ,ixver15m,
     +  ixver53m
c
c satfile = mss path to file containing sat data (attitude):
      parameter(mxpts=2000)
      character*80 satfile
      character*16 chsat
      common/sat/ satfile,satut(mxpts),satlat(mxpts),satlon(mxpts),
     +  satalt(mxpts),iydsat,lusat,iseclim(2),norbfr,mtsat(3,mxpts),
     +  chsat