c
c------------------------------------------------------------------
c Begin file /home/sting/foster/timegcm/mksrc/blktime.f
c------------------------------------------------------------------
c
      block data blktime
      include 'mksrc.h'
c
c For timegcm:
      data          itxt,itxu,itxv,itxo2,itxo1,itxn4s,itxnoz,itxco,
     +              itxco2,itxh2o,itxh2,itxhox,itxop,itxno2,itxno,
     +              itxo3,itxnpo1,itxoh,itxho2,itxh,itxn2d,itxti,
     +              itxte,itxne,itxo2p,itxw,itxz
     +  /1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23,
     +   24,25,26,27/
c
c gcmwtime(1->ntigcmf) molecular weights
c
      data gcmwtime 
c       tn   un   vn   o2    o  n4s   noz   co  co2  h2o  h2 hox 
     +/  0.,  0.,  0., 32., 16., 14.,  30., 28., 44., 18., 2., 1.,
c       o+   no2  no   o3 npo1   oh   ho2    h  n2d   ti  te  ne 
     +   0.,  46.,30., 48., 16., 17.,  33.,  1., 14.,  0., 0., 0.,
c       o2+    w   z 
     +   0.,   0., 0. /
c
c If itimelog(1->ntigcmf) = 0, do not plot log10, otherwise do:
c
      data itimelog 
c       tn   un   vn   o2    o  n4s   noz   co  co2  h2o  h2 hox 
     +/  0,   0,   0,   1,   1,   1,    1,   1,   1,   1,  1,  1,
c       o+   no2  no   o3 npo1   oh   ho2    h  n2d   ti  te  ne 
     +   1,    1,  1,   1,   1,   1,    1,   1,   1,   0,  0,  1,
c       o2+    w   z 
     +    1,   0,  0 /
c
c Short field labels for timegcm:
c
      data timelab_short
     + /'TN      ','UN      ','VN      ','O2      ','O1      ',
     +  'N4S     ','NOZ     ','CO      ','CO2     ','H2O     ',
     +  'H2      ','HOX     ','O+      ','NO2     ','NO      ',
     +  'O3      ','NPO1    ','OH      ','HO2     ','H       ',
     +  'N2D     ','TI      ','TE      ','NE      ','O2+     ',
     +  'W       ','Z       '/
c
      data lutime/10/, iframe/0/, ludatv/3/, luacd/4/, lucira/22/
c
c Species available from sol9.f acd 2d processor (see ~/acd2d)
c
c -------------------------------------------------------------------------
c      1=n2o     2=ccl4    3=cfc13    4=cf2c12   5=ch3cc13   6=ch3cl   7=ch4
c      8=ccl2o   9=cclfo  10=cf2o    11=hcn     12=ch3cn    13=h2o    14=h2
c     15=co     16=co2	  17=tracer  18=noy     19=hno3     20=n2o5   21=clx
c     22=h      23=ox	  24=hf	     25=cfc113  26=cfc22    27=h-1211 28=h-1301
c     29=h2o2   30=ch3br  31=bry     32=cfc114  33=ocs      34=c2h6s  35=sox
c     36=hso3   37=h2so4  38=cfc115
c
c     39=old    40=oh	  41=cl	     42=o	43=o3	    44=ho2    45=no2
c     46=clo    47=hcl	  48=clono2  49=no	50=hno4	    51=hocl   52=ch3o2
c     53=no3    54=n	  55=h2co    56=oclo	57=cl2o2    58=br     59=bro
c     60=brono2 61=hbr	  62=hobr    63=temp.
c --------------------------------------------------------------------------
c
c For acd2d model lower boundary conditions:
c
c luacd (above) = unit for acd 2d model history
c                 (mss history path is input via getinp)
c hgtlb = height of desired lower boundaries (zp -17 = 30 km)
c ispecno = species numbers for acd2d processor (see ~/acd2d)
c ispecpsi(ip)=1 if species is in mass mix ratio on history
c                (will need to convert number density species to psi's)
c
      data hgtlb/30./
c                   h2o, h2, co, co2,  h, oh,  o, o3, ho2, no2, no
      data ispecno / 13, 14, 15,  16, 22, 40, 42, 43,  44,  45, 49/
      data ispecpsi/  1,  1,  1,   1,  1,  0,  0,  0,   0,   0,  0/
      data acdname /'H2O     ','H2      ','CO      ','CO2     ',
     +              'H       ','OH      ','O1      ','O3      ',
     +              'HO2     ','NO2     ','NO      '/
      data lbname /'H2O     ','H2      ','CO      ','CO2     ',
     +             'H       ','OH      ','O1      ','O3      ',
     +             'HO2     ','NO2     ','NO      ','TN(GND) ',
     +             'TN      ','HT      ','UN      ','OX      ',
     +             'NOZ     ','HOX     '/
c
      end