c
c------------------------------------------------------------------
c Begin file /home/sting/foster/timegcm/mksrc/glbm.h
c------------------------------------------------------------------
c
c Modified from ~/timegcm/mksrc/glbm.h to include an ixglbm index
c for each of the 27 fields of a timegcm history. The index is
c into glbmfld for each field. For example, igmt ixglbm(1) = 5,
c i.e., the first timegcm field is tn, and will be set to the 5th
c global mean field (glbmfld(k,5).
c
c (Note: nixglbm should = nfldsd, and kmxglbm should = kmx)
c
c If ixglbm(ip) = 0, then it is either a wind (u,v, or w) and is
c zeroed out, or it is hox, which is oh + h2o + h (see settime.f)
c
c     parameter(nglbmfld=96,kmxglbm=37,nixglbm=27,glbmzp1=-13.,
      parameter(nglbmfld=96,kmxglbm=45,nixglbm=27,glbmzp1=-17.,
     +  glbmdelz=0.5)
      character*24 glbmvol
      common /glbmch/ glbmvol
c
c There are nglbmfld fields in global mean history, each is a
c   column of kmxglbm values: (see also ~/glbmproc)
c (Note winds indices are 0 and will be zeroed out, and that
c  hox index is 0, and will be defined hox = oh + h2o + h
c
      common /glbmcom/ glbmfld(kmxglbm,nglbmfld),glbmzp(kmxglbm)
      common/glbmix/ 
     +  igmtn ,igmun ,igmvn ,igmo2 ,igmo1 ,igmn4s,igmnoz,igmco ,igmco2,
     +  igmh2o,igmh2 ,igmhox,igmop ,igmno2,igmno ,igmo3 ,igmno1,igmoh ,
     +  igmho2,igmh  ,igmn2d,igmti ,igmte ,igmne ,igmo2p,igmw  ,igmz  
      data 
     +  igmtn ,igmun ,igmvn ,igmo2 ,igmo1 ,igmn4s,igmnoz,igmco ,igmco2,
     +  igmh2o,igmh2 ,igmhox,igmop ,igmno2,igmno ,igmo3 ,igmno1,igmoh ,
     +  igmho2,igmh  ,igmn2d,igmti ,igmte ,igmne ,igmo2p,igmw  ,igmz  
     + /    5 ,    0 ,    0 ,   12 ,    7 ,    32,    36,   50 ,    51,
     +      41,   42 ,     0,   26 ,    35,   34 ,   60 ,     7,   44 ,
     +      45,  47  ,    31,   21 ,   20 ,   22 ,    30,   0  ,   3  / 
      dimension ixglbm(nixglbm)
      equivalence (ixglbm,igmtn)
c
c     data fldnm/
c    +'zp      ','zpms    ','zpht    ','zpmht   ','tn      ','tnms    ',
c    +'xno     ','xom     ','xno1d   ','xnn2    ','xn2m    ','xno2    ',
c    +'xo2m    ','sht     ','shtms   ','amas    ','amass   ','rho     ', 
c    +'rhoms   ','te      ','ti      ','xne     ','xneee   ','xnpi    ',
c    +'xnni    ','xiop    ','xinop   ','xinp    ','xin2p   ','xio2p   ',
c    +'xn2d    ','xn4s    ','xn4sm   ','xnno    ','XNNO2   ','xnnoz   ',
c    +'xnhe    ','xnhem   ','xnarg   ','xnargm  ','XNH2O   ','XNH2    ',
c    +'xnch4   ','XNOH    ','XNHO2   ','xnh2o2  ','XNH     ','xnhm    ',
c    +'xnhtot  ','XNCO    ','XNCO2   ','qntot   ','xircool ','dso2    ',
c    +'do3t    ','sped    ','shal    ','parcd   ','qtin    ','XNO3    ',
c    +'xnas    ','xnasp   ','edy     ','qnn     ','qnc     ','qsr     ', 
c    +'qsrb    ','qn1d    ','xlen    ','qo3pr   ','qian    ','qjoul   ',
c    +'qno3    ','qnlya   ','xnoc    ','xco2c   ','xo3pc   ','hcm     ', 
c    +'hce     ','dsichm  ','dsnchm  ','dso2lya ','dso2sr  ','dso2srb ',
c    +'qnight  ','qop     ','qnop    ','qo2p    ','qn2p    ','qinp    ',
c    +'qnspe   ','gwheat  ','qpro    ','qproh   ','qcr     ','qixray  '/
c