      SUBROUTINE CKWC  (P, T, C, ICKWRK, RCKWRK, WDOT)
C
C  START PROLOGUE
C
C SUBROUTINE CKWC  (T, C, ICKWRK, RCKWRK, WDOT)
C    RETURNS THE MOLAR PRODUCTION RATES OF THE
C    SPECIES GIVEN THE TEMPERATURE AND MOLAR
C    CONCENTRATIONS;  SEE EQ. (49).
C
C INPUT
C   T      - TEMPERATURE.
C               CGS UNITS - K
C                 DATA TYPE - REAL SCALAR
C   C      - MOLAR CONCENTRATIONS.
C               CGS UNITS - MOLE/CM**3
C               DATA TYPE - REAL ARRAY
C               DIMENSION C(*) AT LEAST KK, 
C               THE TOTAL NUMBER OF SPECIES.
C   ICKWRK - ARRAY OF INTEGER WORKSPACE.
C               DATA TYPE - INTEGER ARRAY
C               DIMENSION ICKWRK(*) 
C               AT LEAST LENIWK.
C   RCKWRK - ARRAY OF REAL WORK SPACE.
C               DATA TYPE - REAL ARRAY
C               DIMENSION RCKWRK(*) 
C               AT LEAST LENRWK.
C
C OUTPUT
C   WDOT   - CHEMICAL MOLAR PRODUCTION RATES.
C               CGS UNITS - MOLES/(CM**3*SEC)
C               DATA TYPE - REAL ARRAY
C               DIMENSION WDOT(*) AT LEAST KK,
C               THE TOTAL NUMBER OF SPECIES.
C
C  END PROLOGUE
C
C*****DOUBLE PRECISION
       IMPLICIT DOUBLE PRECISION (A-H, O-Z)
       IMPLICIT INTEGER (I-N)
C*****END DOUBLE PRECISION
C*****SINGLE PRECISION
C      IMPLICIT REAL (A-H, O-Z), INTEGER (I-N)
C*****END SINGLE PRECISION
C
      PARAMETER(NR=209)
C
C
C*****ANZAHL DER INNEREN ITERATIONEN
      PARAMETER(NIT=20)  
C
C*****KONVERGENZMASS
      PARAMETER(GDIFF=1.0E-4) 
C
C
      DIMENSION C(*), ICKWRK(*), RCKWRK(*)
      DIMENSION WDOT(*), W(NR)
      DIMENSION WGLOB(16)
      DIMENSION RF(NR), RB(NR), CM(NR)
      DATA SMALL/1.D-50/
C
C-----------------------------------------------------------------------
C --- SET LOCAL VALUES FOR THE CONCENTRATIONS OF THE RED. MECH.
C-----------------------------------------------------------------------

      CO2     = MAX(SMALL, C( 1))
      CH      = MAX(SMALL, C( 2))
      CH2     = MAX(SMALL, C( 3))
      CH2O    = MAX(SMALL, C( 4))
      CN2     = MAX(SMALL, C( 5))
      CCO2    = MAX(SMALL, C( 6))
      CCO     = MAX(SMALL, C( 7))
      CCH2O   = MAX(SMALL, C( 8))
      CCH4    = MAX(SMALL, C( 9))
      CCH3OH  = MAX(SMALL, C(10))
      CC2H6   = MAX(SMALL, C(11))
      CC2H4   = MAX(SMALL, C(12))
      CC2H2   = MAX(SMALL, C(13))
      CNO     = MAX(SMALL, C(14))
      CNH3    = MAX(SMALL, C(15))
      CN2O    = MAX(SMALL, C(16))
      CHCN    = MAX(SMALL, C(17))
      CHNCO   = MAX(SMALL, C(18))
      CC4H5N  = MAX(SMALL, C(19))
      CC6H11NO= MAX(SMALL, C(20))


C-----------------------------------------------------------------------
C --- CALCULATE THE EFFECTIVE THIRD BODY CONCENTRATIONS
C-----------------------------------------------------------------------

      CM(  5)= + CO2       + CCO       + CCO2      + CCH4      + CH2       
     &         + 5.0 * CH2O      + CC2H2     + CC2H4     + CNO       
     &         + CHCN      + CNH3      + CN2O      + CHNCO     + CC2H6     
     &         + CH        + CC4H5N    + CC6H11NO  + CCH2O     + CCH3OH    
     &         + CN2        
      CM( 15)= + CO2       + CCO       + CCO2      + CCH4      + CH2       
     &         + 5.0 * CH2O      + CC2H2     + CC2H4     + CNO       
     &         + CHCN      + CNH3      + CN2O      + CHNCO     + CC2H6     
     &         + CH        + CC4H5N    + CC6H11NO  + CCH2O     + CCH3OH    
     &         + CN2        
      CM( 23)= + CO2       + CCO       + CCO2      + CCH4      + CH2       
     &         + 5.0 * CH2O      + CC2H2     + CC2H4     + CNO       
     &         + CHCN      + CNH3      + CN2O      + CHNCO     + CC2H6     
     &         + CH        + CC4H5N    + CC6H11NO  + CCH2O     + CCH3OH    
     &         + CN2        
      CM( 26)= + CO2       + CCO       + CCO2      + CCH4      + CH2       
     &         + 8.6 * CH2O      + CC2H2     + CC2H4     + CNO       
     &         + CHCN      + CNH3      + CN2O      + CHNCO     + CC2H6     
     &         + CH        + CC4H5N    + CC6H11NO  + CCH2O     + CCH3OH    
     &         + 1.4 * CN2        
      CM( 34)= + CO2       + CCO       + CCO2      + CCH4      + CH2       
     &         + 8.6 * CH2O      + CC2H2     + CC2H4     + CNO       
     &         + CHCN      + CNH3      + CN2O      + CHNCO     + CC2H6     
     &         + CH        + CC4H5N    + CC6H11NO  + CCH2O     + CCH3OH    
     &         + 1.4 * CN2        
      CM( 38)= + CO2       + 2.0 * CCO       + 3.0 * CCO2      
     &         + CCH4      + 2.0 * CH2       + 8.6 * CH2O      
     &         + CC2H2     + CC2H4     + CNO       + CHCN      + CNH3      
     &         + CN2O      + CHNCO     + CC2H6     + CH        + CC4H5N    
     &         + CC6H11NO  + CCH2O     + CCH3OH    + 1.4 * CN2       
     &          
      CM( 43)= + CO2       + CCO       + CCO2      + CCH4      + CH2       
     &         + 8.6 * CH2O      + CC2H2     + CC2H4     + CNO       
     &         + CHCN      + CNH3      + CN2O      + CHNCO     + CC2H6     
     &         + CH        + CC4H5N    + CC6H11NO  + CCH2O     + CCH3OH    
     &         + 1.4 * CN2        
      CM( 44)= + CO2       + 2.0 * CCO       + 3.0 * CCO2      
     &         + CCH4      + 2.0 * CH2       + 8.6 * CH2O      
     &         + CC2H2     + CC2H4     + CNO       + CHCN      + CNH3      
     &         + CN2O      + CHNCO     + CC2H6     + CH        + CC4H5N    
     &         + CC6H11NO  + CCH2O     + CCH3OH    + 1.4 * CN2       
     &          
      CM( 50)= + CO2       + CCO       + CCO2      + CCH4      + CH2       
     &         + 5.0 * CH2O      + CC2H2     + CC2H4     + CNO       
     &         + CHCN      + CNH3      + CN2O      + CHNCO     + CC2H6     
     &         + CH        + CC4H5N    + CC6H11NO  + CCH2O     + CCH3OH    
     &         + CN2        
      CM( 57)= + CO2       + 2.0 * CCO       + 3.0 * CCO2      
     &         + CCH4      + 2.0 * CH2       + 5.0 * CH2O      
     &         + CC2H2     + CC2H4     + CNO       + CHCN      + CNH3      
     &         + CN2O      + CHNCO     + CC2H6     + CH        + CC4H5N    
     &         + CC6H11NO  + CCH2O     + CCH3OH    + CN2        
      CM( 67)= + 1.5 * CO2       + CCO       + 3.0 * CCO2      
     &         + CCH4      + 2.0 * CH2       +10.0 * CH2O      
     &         + CC2H2     + CC2H4     + CNO       + CHCN      + CNH3      
     &         + CN2O      + CHNCO     + CC2H6     + CH        + CC4H5N    
     &         + CC6H11NO  + CCH2O     + CCH3OH    + 0.0 * CN2       
     &          
      CM( 69)= + 1.5 * CO2       + CCO       + CCO2      + CCH4      
     &         + CH2       +10.0 * CH2O      + CC2H2     + CC2H4     
     &         + CNO       + CHCN      + CNH3      + CN2O      + CHNCO     
     &         + CC2H6     + CH        + CC4H5N    + CC6H11NO  + CCH2O     
     &         + CCH3OH    + 1.7 * CN2        
      CM( 70)= + CO2       + CCO       + CCO2      + CCH4      + CH2       
     &         + 5.0 * CH2O      + CC2H2     + CC2H4     + CNO       
     &         + CHCN      + CNH3      + CN2O      + CHNCO     + CC2H6     
     &         + CH        + CC4H5N    + CC6H11NO  + CCH2O     + CCH3OH    
     &         + CN2        
      CM( 78)= + CO2       + CCO       + CCO2      + CCH4      + CH2       
     &         + 5.0 * CH2O      + CC2H2     + CC2H4     + CNO       
     &         + CHCN      + CNH3      + CN2O      + CHNCO     + CC2H6     
     &         + CH        + CC4H5N    + CC6H11NO  + CCH2O     + CCH3OH    
     &         + 2.0 * CN2        
      CM( 87)= 0.012186/T
      CM(115)= + 2.0 * CO2       + CCO       + CCO2      + CCH4      
     &         + 2.0 * CH2       +15.0 * CH2O      + CC2H2     
     &         + CC2H4     + CNO       + CHCN      + CNH3      + CN2O      
     &         + CHNCO     + CC2H6     + CH        + CC4H5N + CC6H11NO  
     &         + CCH2O     + CCH3OH    + 2.0 * CN2        
      CM(123)= + 1.4 * CO2       + 1.5 * CCO       + 3.0 * CCO2      
     &         + CCH4      + CH2       +12.0 * CH2O      + CC2H2     
     &         + CC2H4     + CNO       + CHCN      + CNH3      + CN2O      
     &         + CHNCO     + CC2H6     + CH        + CC4H5N + CC6H11NO  
     &         + CCH2O     + CCH3OH    + 1.7 * CN2        
      CM(150)= 0.012186/T
      CM(185)= 0.012186/T


C-----------------------------------------------------------------------
C --- CALCULATE THE REACTION RATE CONSTANTS
C-----------------------------------------------------------------------

      ALOGT = DLOG(T)

      RF(  1) =  0.20000D+15 * EXP( -0.40*ALOGT) 
      DG = +0.120229400E+00*(ALOGT-1.) +0.510315957E-03*T 
     &     -0.154403390E-06*T*T +0.256195405E-10*T*T*T 
     &     -0.173206767E-14*T*T*T*T +0.869103392E+04/T 
     &     -0.425448636E+01
      RB(  1) = RF(  1)/MAX(EXP(DG),SMALL) 

      RF(  2) =  0.50000D+05 * EXP(  2.67*ALOGT  -3166./T) 
      DG = -0.313812300E+00*(ALOGT-1.) +0.341903140E-03*T 
     &     -0.691499988E-07*T*T +0.846983446E-11*T*T*T 
     &     -0.475631821E-15*T*T*T*T -0.106489672E+04/T 
     &     +0.245069786E+01
      RB(  2) = RF(  2)/MAX(EXP(DG),SMALL) 

      RF(  3) =  0.21000D+09 * EXP(  1.52*ALOGT  -1736./T) 
      DG = -0.896174370E+00*(ALOGT-1.) +0.838945387E-03*T 
     &     -0.131672418E-06*T*T +0.142058798E-10*T*T*T 
     &     -0.747333775E-15*T*T*T*T +0.743913528E+04/T 
     &     +0.324841440E+01
      RB(  3) = RF(  3)/MAX(EXP(DG),SMALL) 

      RF(  4) =  0.43000D+04 * EXP(  2.70*ALOGT  +1251./T) 
      DG = -0.582362070E+00*(ALOGT-1.) +0.497042247E-03*T 
     &     -0.625224192E-07*T*T +0.573604536E-11*T*T*T 
     &     -0.271701955E-15*T*T*T*T +0.850403200E+04/T 
     &     +0.797716540E+00
      RB(  4) = RF(  4)/MAX(EXP(DG),SMALL) 

      RF(  5) =  0.16000D+23 * EXP( -2.00*ALOGT) 
      DG = -0.255889519E+01*(ALOGT-1.) +0.814244174E-03*T 
     &     -0.484296270E-07*T*T -0.757985555E-12*T*T*T 
     &     +0.253943055E-15*T*T*T*T +0.593366140E+05/T 
     &     +0.936756914E+00
      RB(  5) = RF(  5)/MAX(EXP(DG),SMALL) / (82.06*T) 

      RF(  6) =  0.67000D+20 * EXP( -1.42*ALOGT) 
      DG = -0.176532695E+01*(ALOGT-1.) +0.378366307E-03*T 
     &     +0.207180838E-07*T*T -0.793534836E-11*T*T*T 
     &     +0.544046270E-15*T*T*T*T +0.242733455E+05/T 
     &     -0.122144623E+01
      RB(  6) = RF(  6)/MAX(EXP(DG),SMALL) / (82.06*T) 

      RF(  7) =  0.36000D+20 * EXP( -1.00*ALOGT) 
      DG = -0.176532695E+01*(ALOGT-1.) +0.378366307E-03*T 
     &     +0.207180838E-07*T*T -0.793534836E-11*T*T*T 
     &     +0.544046270E-15*T*T*T*T +0.242733455E+05/T 
     &     -0.122144623E+01
      RB(  7) = RF(  7)/MAX(EXP(DG),SMALL) / (82.06*T) 

      RF(  8) =  0.20000D+21 * EXP( -1.42*ALOGT) 
      DG = -0.176532695E+01*(ALOGT-1.) +0.378366307E-03*T 
     &     +0.207180838E-07*T*T -0.793534836E-11*T*T*T 
     &     +0.544046270E-15*T*T*T*T +0.242733455E+05/T 
     &     -0.122144623E+01
      RB(  8) = RF(  8)/MAX(EXP(DG),SMALL) / (82.06*T) 

      RF(  9) =  0.43000D+14 * EXP(   -710./T) 
      DG = +0.102606130E+00*(ALOGT-1.) -0.403067520E-03*T 
     &     +0.625247071E-07*T*T -0.702851702E-11*T*T*T 
     &     +0.457230560E-15*T*T*T*T +0.276241333E+05/T 
     &     -0.109021126E+01
      RB(  9) = RF(  9)/MAX(EXP(DG),SMALL) 

      RF( 10) =  0.17000D+15 * EXP(   -440./T) 
      DG = -0.331435570E+00*(ALOGT-1.) -0.571480337E-03*T 
     &     +0.147778098E-06*T*T -0.241782230E-10*T*T*T 
     &     +0.171366641E-14*T*T*T*T +0.178682026E+05/T 
     &     +0.561497296E+01
      RB( 10) = RF( 10)/MAX(EXP(DG),SMALL) 

      RF( 11) =  0.30000D+14 * EXP(   -866./T) 
      DG = -0.913797640E+00*(ALOGT-1.) -0.744380903E-04*T 
     &     +0.852556791E-07*T*T -0.184421777E-10*T*T*T 
     &     +0.144196446E-14*T*T*T*T +0.263722346E+05/T 
     &     +0.641268950E+01
      RB( 11) = RF( 11)/MAX(EXP(DG),SMALL) 

      RF( 12) =  0.33000D+14 
      DG = -0.211206170E+00*(ALOGT-1.) -0.611643802E-04*T 
     &     -0.662529165E-08*T*T +0.144131744E-11*T*T*T 
     &     -0.184012605E-16*T*T*T*T +0.265592365E+05/T 
     &     +0.136048660E+01
      RB( 12) = RF( 12)/MAX(EXP(DG),SMALL) 

      RF( 13) =  0.19000D+17 * EXP( -1.00*ALOGT) 
      DG = -0.793568240E+00*(ALOGT-1.) +0.435877867E-03*T 
     &     -0.691477108E-07*T*T +0.717736280E-11*T*T*T 
     &     -0.290103215E-15*T*T*T*T +0.350632685E+05/T 
     &     +0.215820314E+01
      RB( 13) = RF( 13)/MAX(EXP(DG),SMALL) 

      RF( 14) =  0.42000D+15 * EXP(  -6030./T) 
      DG = -0.586881110E+00*(ALOGT-1.) +0.955882110E-03*T 
     &     -0.232007103E-06*T*T +0.293469834E-10*T*T*T 
     &     -0.144404514E-14*T*T*T*T +0.191739588E+05/T 
     &     +0.799183610E+00
      RB( 14) = RF( 14)/MAX(EXP(DG),SMALL) 

      RF( 15) =  0.13000D+18 * EXP( -22899./T) 
      DG = +0.202077249E+01*(ALOGT-1.) -0.190572875E-02*T 
     &     +0.359067118E-06*T*T -0.455898581E-10*T*T*T 
     &     +0.261366529E-14*T*T*T*T -0.255791017E+05/T 
     &     +0.603723558E+01
      RB( 15) = RF( 15)/MAX(EXP(DG),SMALL) * (82.06*T) 

      RF( 16) =  0.58000D+15 * EXP(  -4811./T) 
      DG = -0.206687130E+00*(ALOGT-1.) -0.520004243E-03*T 
     &     +0.162859392E-06*T*T -0.221696206E-10*T*T*T 
     &     +0.115394193E-14*T*T*T*T +0.158893097E+05/T 
     &     +0.135901953E+01
      RB( 16) = RF( 16)/MAX(EXP(DG),SMALL) 

      RF( 17) =  0.15000D+08 * EXP(  1.30*ALOGT   +381./T) 
      DG = +0.549387020E+00*(ALOGT-1.) +0.901706545E-03*T 
     &     -0.223748913E-06*T*T +0.317819442E-10*T*T*T 
     &     -0.192924410E-14*T*T*T*T +0.129922907E+05/T 
     &     -0.104704287E+02
      RB( 17) = RF( 17)/MAX(EXP(DG),SMALL) 

      RF( 18) =  0.58000D+14 * EXP( -11542./T) 
      DG = +0.217951450E+00*(ALOGT-1.) +0.330226208E-03*T 
     &     -0.759708152E-07*T*T +0.760372120E-11*T*T*T 
     &     -0.215577685E-15*T*T*T*T +0.308604933E+05/T 
     &     -0.485545571E+01
      RB( 18) = RF( 18)/MAX(EXP(DG),SMALL) 

      RF( 19) =  0.13000D+09 * EXP(  1.62*ALOGT  -1090./T) 
      DG = +0.184876349E+01*(ALOGT-1.) -0.214622401E-02*T 
     &     +0.406264794E-06*T*T -0.497347276E-10*T*T*T 
     &     +0.275301150E-14*T*T*T*T +0.841606837E+04/T 
     &     -0.661631848E+01
      RB( 19) = RF( 19)/MAX(EXP(DG),SMALL) 

      RF( 20) =  0.18000D+14 * EXP(  -1550./T) 
      DG = +0.153495119E+01*(ALOGT-1.) -0.180432087E-02*T 
     &     +0.337114795E-06*T*T -0.412648931E-10*T*T*T 
     &     +0.227737968E-14*T*T*T*T +0.735117165E+04/T 
     &     -0.416562062E+01
      RB( 20) = RF( 20)/MAX(EXP(DG),SMALL) 

      RF( 21) =  0.34000D+10 * EXP(  1.18*ALOGT   +225./T) 
      DG = +0.952589120E+00*(ALOGT-1.) -0.130727862E-02*T 
     &     +0.274592376E-06*T*T -0.355288478E-10*T*T*T 
     &     +0.200567772E-14*T*T*T*T +0.158552037E+05/T 
     &     -0.336790408E+01
      RB( 21) = RF( 21)/MAX(EXP(DG),SMALL) 

      RF( 22) =  0.60000D+14 * EXP( -20463./T) 
      DG = +0.174615736E+01*(ALOGT-1.) -0.174315649E-02*T 
     &     +0.343740086E-06*T*T -0.427062106E-10*T*T*T 
     &     +0.229578094E-14*T*T*T*T -0.192080649E+05/T 
     &     -0.552610722E+01
      RB( 22) = RF( 22)/MAX(EXP(DG),SMALL) 

      RF( 23) =  0.19000D+18 * EXP( -1.00*ALOGT  -8556./T) 
      DG = +0.244301124E+01*(ALOGT-1.) -0.144721393E-02*T 
     &     +0.247605063E-06*T*T -0.299257312E-10*T*T*T 
     &     +0.164930500E-14*T*T*T*T -0.730986935E+04/T 
     &     -0.242634011E+01
      RB( 23) = RF( 23)/MAX(EXP(DG),SMALL) * (82.06*T) 

      RF( 24) =  0.12000D+14 * EXP(  0.25*ALOGT) 
      DG = +0.780290420E+00*(ALOGT-1.) -0.147191514E-02*T 
     &     +0.330847853E-06*T*T -0.448895966E-10*T*T*T 
     &     +0.265058183E-14*T*T*T*T +0.445876094E+05/T 
     &     -0.473799760E+01
      RB( 24) = RF( 24)/MAX(EXP(DG),SMALL) 

      RF( 25) =  0.76000D+13 * EXP(   -201./T) 
      DG = +0.677684290E+00*(ALOGT-1.) -0.106884762E-02*T 
     &     +0.268323146E-06*T*T -0.378610796E-10*T*T*T 
     &     +0.219335127E-14*T*T*T*T +0.169634761E+05/T 
     &     -0.364778634E+01
      RB( 25) = RF( 25)/MAX(EXP(DG),SMALL) 

      RFH =  0.13000D+17 * T**(-0.63) * EXP(   -193./T)
      RFL =  0.17500D+34 * T**(-4.76) * EXP(  -1228./T)
      PR = CM( 26)*RFL/MAX(RFH,SMALL)
      F4 =   0.783000
      F5 =      74.00
      F6 =    2941.00
      F7 =    6964.00
      FC = (1.-F4)*EXP(-T/F5)+F4*EXP(-T/F6)+EXP(-F7/T)
      FCLOG = DLOG10(FC)
      PRLOG = DLOG10(PR)
      X = -0.4-0.67*FCLOG
      Y = 0.75-1.27*FCLOG
      A = (PRLOG+X)/(Y-0.14*(PRLOG+X))
      F = 10.**(FCLOG/(1.+A*A))
      RF( 26) = RFH * PR/(1.+PR)*F 
      DG = -0.471086624E+01*(ALOGT-1.) +0.307552318E-02*T 
     &     -0.457619104E-06*T*T +0.522700592E-10*T*T*T 
     &     -0.275498955E-14*T*T*T*T +0.517175888E+05/T 
     &     +0.104039291E+02
      RB( 26) = RF( 26)/MAX(EXP(DG),SMALL) / (82.06*T) 

      RF( 27) =  0.13000D+05 * EXP(  3.00*ALOGT  -4046./T) 
      DG = +0.304814542E+01*(ALOGT-1.) -0.310022439E-02*T 
     &     +0.540861895E-06*T*T -0.672339246E-10*T*T*T 
     &     +0.375626639E-14*T*T*T*T +0.179889932E+03/T 
     &     -0.127155866E+02
      RB( 27) = RF( 27)/MAX(EXP(DG),SMALL) 

      RF( 28) =  0.10000D+10 * EXP(  1.50*ALOGT  -4328./T) 
      DG = +0.273433312E+01*(ALOGT-1.) -0.275832125E-02*T 
     &     +0.471711897E-06*T*T -0.587640901E-10*T*T*T 
     &     +0.328063456E-14*T*T*T*T -0.885006790E+03/T 
     &     -0.102648888E+02
      RB( 28) = RF( 28)/MAX(EXP(DG),SMALL) 

      RF( 29) =  0.13600D+07 * EXP(  2.18*ALOGT  -1349./T) 
      DG = +0.215197105E+01*(ALOGT-1.) -0.226127900E-02*T 
     &     +0.409189477E-06*T*T -0.530280448E-10*T*T*T 
     &     +0.300893261E-14*T*T*T*T +0.761902521E+04/T 
     &     -0.946717221E+01
      RB( 29) = RF( 29)/MAX(EXP(DG),SMALL) 

      RF( 30) =  0.79000D+14 * EXP( -28183./T) 
      DG = +0.294553929E+01*(ALOGT-1.) -0.269715687E-02*T 
     &     +0.478337188E-06*T*T -0.602054076E-10*T*T*T 
     &     +0.329903582E-14*T*T*T*T -0.274442433E+05/T 
     &     -0.116253754E+02
      RB( 30) = RF( 30)/MAX(EXP(DG),SMALL) 

      RF( 31) =  0.84000D+14 
      DG = -0.594448410E+00*(ALOGT-1.) +0.102303727E-02*T 
     &     -0.246232182E-06*T*T +0.350621043E-10*T*T*T 
     &     -0.214492305E-14*T*T*T*T +0.345153358E+05/T 
     &     -0.547703440E-01
      RB( 31) = RF( 31)/MAX(EXP(DG),SMALL) 

      RF( 32) =  0.30000D+16 * EXP( -0.60*ALOGT) 
      DG = -0.173059770E+01*(ALOGT-1.) +0.201675002E-02*T 
     &     -0.431417407E-06*T*T +0.574477427E-10*T*T*T 
     &     -0.329137925E-14*T*T*T*T -0.631046400E+03/T 
     &     +0.953802939E+01
      RB( 32) = RF( 32)/MAX(EXP(DG),SMALL) 

      RF( 33) =  0.10000D+15 
      DG = -0.814060310E+00*(ALOGT-1.) -0.428718930E-03*T 
     &     +0.148395489E-06*T*T -0.232913452E-10*T*T*T 
     &     +0.149406089E-14*T*T*T*T +0.159691233E+04/T 
     &     +0.759301865E+01
      RB( 33) = RF( 33)/MAX(EXP(DG),SMALL) 

      RFH =  0.63000D+14
      RFL =  0.18900D+39 * T**(-6.30) * EXP(  -1560./T)
      PR = CM( 34)*RFL/MAX(RFH,SMALL)
      F4 =   0.210500
      F5 =      83.50
      F6 =    5398.00
      F7 =    8370.00
      FC = (1.-F4)*EXP(-T/F5)+F4*EXP(-T/F6)+EXP(-F7/T)
      FCLOG = DLOG10(FC)
      PRLOG = DLOG10(PR)
      X = -0.4-0.67*FCLOG
      Y = 0.75-1.27*FCLOG
      A = (PRLOG+X)/(Y-0.14*(PRLOG+X))
      F = 10.**(FCLOG/(1.+A*A))
      RF( 34) = RFH * PR/(1.+PR)*F 
      DG = -0.358889748E+01*(ALOGT-1.) +0.315274956E-02*T 
     &     -0.584061683E-06*T*T +0.728310301E-10*T*T*T 
     &     -0.410430091E-14*T*T*T*T +0.460091160E+05/T 
     &     +0.154559441E+01
      RB( 34) = RF( 34)/MAX(EXP(DG),SMALL) / (82.06*T) 

      RF( 35) =  0.80000D+13 
      DG = +0.560758050E+00*(ALOGT-1.) -0.529896892E-03*T 
     &     +0.181820476E-06*T*T -0.336195058E-10*T*T*T 
     &     +0.240606852E-14*T*T*T*T +0.129009515E+05/T 
     &     -0.485890284E+01
      RB( 35) = RF( 35)/MAX(EXP(DG),SMALL) 

      RF( 36) =  0.29000D+14 * EXP( -15340./T) 
      DG = +0.771964220E+00*(ALOGT-1.) -0.468732512E-03*T 
     &     +0.188445768E-06*T*T -0.350608232E-10*T*T*T 
     &     +0.242446979E-14*T*T*T*T -0.136582850E+05/T 
     &     -0.621938944E+01
      RB( 36) = RF( 36)/MAX(EXP(DG),SMALL) 

      RF( 37) =  0.19000D+13 * EXP( -10224./T) 
      DG = -0.714677810E+00*(ALOGT-1.) +0.512721313E-03*T 
     &     -0.918287922E-07*T*T +0.944256378E-11*T*T*T 
     &     -0.412855380E-15*T*T*T*T +0.258243019E+05/T 
     &     +0.419971602E+01
      RB( 37) = RF( 37)/MAX(EXP(DG),SMALL) 

      RFH =  0.21000D+17 * T**(-0.97) * EXP(   -312./T)
      RFL =  0.12600D+51 * T**(-9.67) * EXP(  -3130./T)
      PR = CM( 38)*RFL/MAX(RFH,SMALL)
      F4 =   0.532500
      F5 =     151.00
      F6 =    1038.00
      F7 =    4970.00
      FC = (1.-F4)*EXP(-T/F5)+F4*EXP(-T/F6)+EXP(-F7/T)
      FCLOG = DLOG10(FC)
      PRLOG = DLOG10(PR)
      X = -0.4-0.67*FCLOG
      Y = 0.75-1.27*FCLOG
      A = (PRLOG+X)/(Y-0.14*(PRLOG+X))
      F = 10.**(FCLOG/(1.+A*A))
      RF( 38) = RFH * PR/(1.+PR)*F 
      DG = -0.349955394E+01*(ALOGT-1.) +0.360273348E-02*T 
     &     -0.675219957E-06*T*T +0.852292268E-10*T*T*T 
     &     -0.482860056E-14*T*T*T*T +0.449775618E+05/T 
     &     -0.184453288E+01
      RB( 38) = RF( 38)/MAX(EXP(DG),SMALL) / (82.06*T) 

      RF( 39) =  0.78000D-07 * EXP(  6.10*ALOGT   -990./T) 
      DG = -0.119938193E+01*(ALOGT-1.) +0.954000385E-03*T 
     &     -0.134597102E-06*T*T +0.174991970E-10*T*T*T 
     &     -0.100325488E-14*T*T*T*T +0.823617844E+04/T 
     &     +0.609926813E+01
      RB( 39) = RF( 39)/MAX(EXP(DG),SMALL) 

      RF( 40) =  0.22000D+23 * EXP( -3.30*ALOGT  -1443./T) 
      DG = -0.150711174E+01*(ALOGT-1.) +0.573303865E-03*T 
     &     -0.495134050E-07*T*T +0.828365667E-12*T*T*T 
     &     +0.191295550E-15*T*T*T*T +0.886766848E+05/T 
     &     +0.725667680E+01
      RB( 40) = RF( 40)/MAX(EXP(DG),SMALL) 

      RF( 41) =  0.13000D+14 * EXP(   -216./T) 
      DG = -0.101458600E+00*(ALOGT-1.) +0.159191000E-04*T 
     &     +0.350166667E-10*T*T -0.150623750E-12*T*T*T 
     &     +0.133956000E-16*T*T*T*T +0.439851900E+04/T 
     &     +0.162358133E+01
      RB( 41) = RF( 41)/MAX(EXP(DG),SMALL) 

      RF( 42) =  0.70000D+14 
      DG = +0.950324850E+00*(ALOGT-1.) -0.225021209E-03*T 
     &     -0.104876131E-08*T*T +0.143572747E-11*T*T*T 
     &     -0.492519051E-16*T*T*T*T +0.337385898E+05/T 
     &     +0.794350122E+01
      RB( 42) = RF( 42)/MAX(EXP(DG),SMALL) * (82.06*T) 

      RFH =  0.54000D+12 * T**( 0.45) * EXP(  -1309./T)
      RFL =  0.15400D+31 * T**(-4.80) * EXP(  -2798./T)
      PR = CM( 43)*RFL/MAX(RFH,SMALL)
      F4 =   0.758000
      F5 =      94.00
      F6 =    1555.00
      F7 =    4200.00
      FC = (1.-F4)*EXP(-T/F5)+F4*EXP(-T/F6)+EXP(-F7/T)
      FCLOG = DLOG10(FC)
      PRLOG = DLOG10(PR)
      X = -0.4-0.67*FCLOG
      Y = 0.75-1.27*FCLOG
      A = (PRLOG+X)/(Y-0.14*(PRLOG+X))
      F = 10.**(FCLOG/(1.+A*A))
      RF( 43) = RFH * PR/(1.+PR)*F 
      DG = -0.489891090E+00*(ALOGT-1.) -0.664251898E-03*T 
     &     +0.294367352E-06*T*T -0.509974179E-10*T*T*T 
     &     +0.336297018E-14*T*T*T*T +0.113499951E+05/T 
     &     -0.102800651E+02
      RB( 43) = RF( 43)/MAX(EXP(DG),SMALL) / (82.06*T) 

      RFH =  0.54000D+12 * T**( 0.45) * EXP(  -1812./T)
      RFL =  0.91000D+32 * T**(-4.82) * EXP(  -3286./T)
      PR = CM( 44)*RFL/MAX(RFH,SMALL)
      F4 =   0.718700
      F5 =     103.00
      F6 =    1291.00
      F7 =    4160.00
      FC = (1.-F4)*EXP(-T/F5)+F4*EXP(-T/F6)+EXP(-F7/T)
      FCLOG = DLOG10(FC)
      PRLOG = DLOG10(PR)
      X = -0.4-0.67*FCLOG
      Y = 0.75-1.27*FCLOG
      A = (PRLOG+X)/(Y-0.14*(PRLOG+X))
      F = 10.**(FCLOG/(1.+A*A))
      RF( 44) = RFH * PR/(1.+PR)*F 
      DG = -0.568024400E+00*(ALOGT-1.) -0.277116413E-03*T 
     &     +0.111929486E-06*T*T -0.182647899E-10*T*T*T 
     &     +0.117650717E-14*T*T*T*T +0.147203339E+05/T 
     &     -0.739920794E+01
      RB( 44) = RF( 44)/MAX(EXP(DG),SMALL) / (82.06*T) 

      RF( 45) =  0.72000D+14 * EXP(  -1800./T) 
      DG = -0.119730255E+01*(ALOGT-1.) +0.655482720E-03*T 
     &     -0.912114018E-07*T*T +0.103294416E-10*T*T*T 
     &     -0.632460900E-15*T*T*T*T +0.955301160E+04/T 
     &     +0.617776171E+01
      RB( 45) = RF( 45)/MAX(EXP(DG),SMALL) 

      RF( 46) =  0.54000D+03 * EXP(  3.50*ALOGT  -2622./T) 
      DG = +0.850031339E+01*(ALOGT-1.) -0.814045123E-02*T 
     &     +0.169143195E-05*T*T -0.225555467E-09*T*T*T 
     &     +0.131436214E-13*T*T*T*T +0.473215662E+04/T 
     &     -0.409780371E+02
      RB( 46) = RF( 46)/MAX(EXP(DG),SMALL) 

      RF( 47) =  0.30000D+08 * EXP(  2.00*ALOGT  -2574./T) 
      DG = +0.818650109E+01*(ALOGT-1.) -0.779854809E-02*T 
     &     +0.162228195E-05*T*T -0.217085632E-09*T*T*T 
     &     +0.126679895E-13*T*T*T*T +0.366725990E+04/T 
     &     -0.385273392E+02
      RB( 47) = RF( 47)/MAX(EXP(DG),SMALL) 

      RF( 48) =  0.72000D+07 * EXP(  2.00*ALOGT   -435./T) 
      DG = +0.760413902E+01*(ALOGT-1.) -0.730150584E-02*T 
     &     +0.155975953E-05*T*T -0.211349587E-09*T*T*T 
     &     +0.123962876E-13*T*T*T*T +0.121712919E+05/T 
     &     -0.377296227E+02
      RB( 48) = RF( 48)/MAX(EXP(DG),SMALL) 

      RF( 49) =  0.55000D+00 * EXP(  4.00*ALOGT  -4177./T) 
      DG = +0.545216798E+01*(ALOGT-1.) -0.504022684E-02*T 
     &     +0.115057005E-05*T*T -0.158321542E-09*T*T*T 
     &     +0.938735497E-14*T*T*T*T +0.455226669E+04/T 
     &     -0.282624505E+02
      RB( 49) = RF( 49)/MAX(EXP(DG),SMALL) 

      RFH =  0.11000D+13 * T**( 0.45) * EXP(   -917./T)
      RFL =  0.11100D+35 * T**(-5.00) * EXP(  -2239./T)
      PR = CM( 50)*RFL/MAX(RFH,SMALL)
      F4 =   0.500000
      F5 =      95.00
      F6 =      95.00
      F7 =     200.00
      FC = (1.-F4)*EXP(-T/F5)+F4*EXP(-T/F6)+EXP(-F7/T)
      FCLOG = DLOG10(FC)
      PRLOG = DLOG10(PR)
      X = -0.4-0.67*FCLOG
      Y = 0.75-1.27*FCLOG
      A = (PRLOG+X)/(Y-0.14*(PRLOG+X))
      F = 10.**(FCLOG/(1.+A*A))
      RF( 50) = RFH * PR/(1.+PR)*F 
      DG = +0.419880453E+01*(ALOGT-1.) -0.459582369E-02*T 
     &     +0.105571789E-05*T*T -0.148767369E-09*T*T*T 
     &     +0.892750303E-14*T*T*T*T +0.201482770E+05/T 
     &     -0.329627724E+02
      RB( 50) = RF( 50)/MAX(EXP(DG),SMALL) / (82.06*T) 

      RF( 51) =  0.49000D+13 * EXP(  0.35*ALOGT) 
      DG = -0.666348027E+01*(ALOGT-1.) +0.451301653E-02*T 
     &     -0.932969201E-06*T*T +0.125362374E-09*T*T*T 
     &     -0.731374396E-14*T*T*T*T +0.218776030E+04/T 
     &     +0.405109125E+02
      RB( 51) = RF( 51)/MAX(EXP(DG),SMALL) 

      RF( 52) =  0.10000D+11 * EXP(   1102./T) 
      DG = -0.596413148E+01*(ALOGT-1.) +0.497418999E-02*T 
     &     -0.103499981E-05*T*T +0.140832020E-09*T*T*T 
     &     -0.838345676E-14*T*T*T*T +0.412506843E+04/T 
     &     +0.317413262E+02
      RB( 52) = RF( 52)/MAX(EXP(DG),SMALL) 

      RF( 53) =  0.54000D+15 * EXP(  -7499./T) 
      DG = +0.598733573E+01*(ALOGT-1.) -0.561977468E-02*T 
     &     +0.115264703E-05*T*T -0.154884587E-09*T*T*T 
     &     +0.912646861E-14*T*T*T*T -0.865922038E+03/T 
     &     -0.290412827E+02
      RB( 53) = RF( 53)/MAX(EXP(DG),SMALL) 

      RF( 54) =  0.47000D+07 * EXP(  1.88*ALOGT    -91./T) 
      DG = +0.387013807E+01*(ALOGT-1.) -0.321442500E-02*T 
     &     +0.654201118E-06*T*T -0.879317546E-10*T*T*T 
     &     +0.513587773E-14*T*T*T*T +0.819348354E+04/T 
     &     -0.205816419E+02
      RB( 54) = RF( 54)/MAX(EXP(DG),SMALL) 

      RF( 55) =  0.81000D+07 * EXP(  1.88*ALOGT    -91./T) 
      DG = +0.452360160E+00*(ALOGT-1.) -0.118129268E-02*T 
     &     +0.199538514E-06*T*T -0.231137523E-10*T*T*T 
     &     +0.122057068E-14*T*T*T*T +0.133699628E+05/T 
     &     +0.318870877E+01
      RB( 55) = RF( 55)/MAX(EXP(DG),SMALL) 

      RF( 56) =  0.20000D+14 * EXP(  -2989./T) 
      DG = +0.509116136E+01*(ALOGT-1.) -0.478082929E-02*T 
     &     +0.102097461E-05*T*T -0.140678707E-09*T*T*T 
     &     +0.837913483E-14*T*T*T*T +0.657321324E+04/T 
     &     -0.257928683E+02
      RB( 56) = RF( 56)/MAX(EXP(DG),SMALL) 

      RFH =  0.31000D+12 * T**( 0.58) * EXP(  -1303./T)
      RFL =  0.22500D+41 * T**(-7.27) * EXP(  -3310./T)
      PR = CM( 57)*RFL/MAX(RFH,SMALL)
      F4 =   0.500000
      F5 =     675.00
      F6 =     675.00
      FC = (1.-F4)*EXP(-T/F5)+F4*EXP(-T/F6)
      FCLOG = DLOG10(FC)
      PRLOG = DLOG10(PR)
      X = -0.4-0.67*FCLOG
      Y = 0.75-1.27*FCLOG
      A = (PRLOG+X)/(Y-0.14*(PRLOG+X))
      F = 10.**(FCLOG/(1.+A*A))
      RF( 57) = RFH * PR/(1.+PR)*F 
      DG = +0.538598050E+00*(ALOGT-1.) -0.125319921E-02*T 
     &     +0.346432462E-06*T*T -0.561859655E-10*T*T*T 
     &     +0.364606474E-14*T*T*T*T +0.191800321E+05/T 
     &     -0.143006089E+02
      RB( 57) = RF( 57)/MAX(EXP(DG),SMALL) / (82.06*T) 

      RF( 58) =  0.11000D+24 * EXP( -3.29*ALOGT  -1958./T) 
      DG = -0.593584108E+01*(ALOGT-1.) +0.460930017E-02*T 
     &     -0.975787310E-06*T*T +0.132743564E-09*T*T*T 
     &     -0.784312148E-14*T*T*T*T +0.411250834E+05/T 
     &     +0.339790096E+02
      RB( 58) = RF( 58)/MAX(EXP(DG),SMALL) 

      RF( 59) =  0.25000D+16 * EXP( -0.78*ALOGT  -1578./T) 
      DG = -0.192361476E+01*(ALOGT-1.) +0.155313057E-02*T 
     &     -0.274892524E-06*T*T +0.328634571E-10*T*T*T 
     &     -0.178289138E-14*T*T*T*T +0.143326838E+04/T 
     &     +0.102634293E+02
      RB( 59) = RF( 59)/MAX(EXP(DG),SMALL) 

      RF( 60) =  0.61000D+07 * EXP(  2.00*ALOGT   -956./T) 
      DG = -0.280528100E-01*(ALOGT-1.) -0.110170745E-02*T 
     &     +0.204053850E-06*T*T -0.261451088E-10*T*T*T 
     &     +0.144158044E-14*T*T*T*T +0.236964777E+05/T 
     &     +0.434018002E+01
      RB( 60) = RF( 60)/MAX(EXP(DG),SMALL) 

      RF( 61) =  0.14000D+08 * EXP(  2.00*ALOGT   -956./T) 
      DG = +0.141121539E+01*(ALOGT-1.) -0.895176225E-03*T 
     &     +0.122924230E-06*T*T -0.142611563E-10*T*T*T 
     &     +0.774367644E-15*T*T*T*T +0.103527034E+05/T 
     &     -0.793099984E+01
      RB( 61) = RF( 61)/MAX(EXP(DG),SMALL) 

      RF( 62) =  0.20000D+09 * EXP(  1.50*ALOGT -15148./T) 
      DG = -0.188575872E+01*(ALOGT-1.) +0.123457817E-02*T 
     &     -0.306332845E-06*T*T +0.417867358E-10*T*T*T 
     &     -0.244532207E-14*T*T*T*T +0.168237052E+05/T 
     &     +0.153737398E+02
      RB( 62) = RF( 62)/MAX(EXP(DG),SMALL) 

      RF( 63) =  0.10000D+14 * EXP( -21137./T) 
      DG = -0.974766670E+00*(ALOGT-1.) +0.585918400E-03*T 
     &     -0.207057542E-06*T*T +0.348922710E-10*T*T*T 
     &     -0.226600205E-14*T*T*T*T -0.213339010E+04/T 
     &     +0.213440105E+02
      RB( 63) = RF( 63)/MAX(EXP(DG),SMALL) * (82.06*T) 

      RF( 64) =  0.22000D+12 * EXP(   -755./T) 
      DG = -0.168944448E+01*(ALOGT-1.) +0.109863971E-02*T 
     &     -0.298886334E-06*T*T +0.443348348E-10*T*T*T 
     &     -0.267885743E-14*T*T*T*T +0.236909118E+05/T 
     &     +0.255437265E+02
      RB( 64) = RF( 64)/MAX(EXP(DG),SMALL) * (82.06*T) 

      RF( 65) =  0.14000D+08 * EXP(  1.70*ALOGT   -503./T) 
      DG = +0.161902270E+00*(ALOGT-1.) +0.454235013E-03*T 
     &     -0.143703071E-06*T*T +0.214843430E-10*T*T*T 
     &     -0.132467881E-14*T*T*T*T +0.556761358E+05/T 
     &     +0.812067108E+01
      RB( 65) = RF( 65)/MAX(EXP(DG),SMALL) * (82.06*T) 

      RF( 66) =  0.29000D+08 * EXP(  1.70*ALOGT   -503./T) 
      DG = -0.387484750E+00*(ALOGT-1.) -0.447471532E-03*T 
     &     +0.800458425E-07*T*T -0.102976012E-10*T*T*T 
     &     +0.604565290E-15*T*T*T*T +0.426838451E+05/T 
     &     +0.185910997E+02
      RB( 66) = RF( 66)/MAX(EXP(DG),SMALL) * (82.06*T) 

      RF( 67) =  0.27000D+16 * EXP(    302./T) 
      DG = -0.214546161E+01*(ALOGT-1.) +0.101069086E-02*T 
     &     -0.138527756E-06*T*T +0.131060113E-10*T*T*T 
     &     -0.566437530E-15*T*T*T*T +0.236255265E+05/T 
     &     -0.111235579E+01
      RB( 67) = RF( 67)/MAX(EXP(DG),SMALL) / (82.06*T) 

      RF( 68) =  0.70000D+20 * EXP( -1.50*ALOGT) 
      DG = -0.214546161E+01*(ALOGT-1.) +0.101069086E-02*T 
     &     -0.138527756E-06*T*T +0.131060113E-10*T*T*T 
     &     -0.566437530E-15*T*T*T*T +0.236255265E+05/T 
     &     -0.111235579E+01
      RB( 68) = RF( 68)/MAX(EXP(DG),SMALL) / (82.06*T) 

      RF( 69) =  0.75000D+20 * EXP( -1.41*ALOGT) 
      DG = -0.945272180E+00*(ALOGT-1.) +0.533632707E-03*T 
     &     -0.734745065E-07*T*T +0.816960091E-11*T*T*T 
     &     -0.385281100E-15*T*T*T*T +0.368220554E+05/T 
     &     -0.112710583E+02
      RB( 69) = RF( 69)/MAX(EXP(DG),SMALL) / (82.06*T) 

      RF( 70) =  0.51000D+24 * EXP( -2.51*ALOGT    +34./T) 
      DG = -0.866601270E+00*(ALOGT-1.) +0.123926499E-02*T 
     &     -0.224412891E-06*T*T +0.263063405E-10*T*T*T 
     &     -0.139595538E-14*T*T*T*T +0.250482776E+05/T 
     &     -0.138430008E+02
      RB( 70) = RF( 70)/MAX(EXP(DG),SMALL) / (82.06*T) 

      RF( 71) =  0.21000D+13 * EXP(    241./T) 
      DG = +0.699825370E+00*(ALOGT-1.) -0.355049557E-03*T 
     &     +0.602107997E-07*T*T -0.951459122E-11*T*T*T 
     &     +0.802740300E-15*T*T*T*T +0.385767601E+04/T 
     &     -0.579512570E+01
      RB( 71) = RF( 71)/MAX(EXP(DG),SMALL) 

      RF( 72) =  0.84000D+14 
      DG = -0.103126094E+01*(ALOGT-1.) -0.216430780E-03*T 
     &     +0.875672986E-07*T*T -0.146636318E-10*T*T*T 
     &     +0.910926110E-15*T*T*T*T +0.140105266E+05/T 
     &     +0.114100987E+02
      RB( 72) = RF( 72)/MAX(EXP(DG),SMALL) 

      RF( 73) =  0.39000D+13 * EXP(    120./T) 
      DG = -0.911031540E+00*(ALOGT-1.) +0.293885177E-03*T 
     &     -0.668360913E-07*T*T +0.109559087E-10*T*T*T 
     &     -0.821141560E-15*T*T*T*T +0.227015605E+05/T 
     &     +0.715561230E+01
      RB( 73) = RF( 73)/MAX(EXP(DG),SMALL) 

      RF( 74) =  0.45000D+12 * EXP(  0.72*ALOGT   -330./T) 
      DG = +0.482740790E+00*(ALOGT-1.) -0.103539208E-02*T 
     &     +0.221770547E-06*T*T -0.280698767E-10*T*T*T 
     &     +0.156771436E-14*T*T*T*T +0.282719522E+05/T 
     &     -0.119930170E+01
      RB( 74) = RF( 74)/MAX(EXP(DG),SMALL) 

      RF( 75) =  0.10000D+14 
      DG = +0.168928490E+00*(ALOGT-1.) -0.693488935E-03*T 
     &     +0.152620548E-06*T*T -0.196000422E-10*T*T*T 
     &     +0.109208254E-14*T*T*T*T +0.272070555E+05/T 
     &     +0.125139616E+01
      RB( 75) = RF( 75)/MAX(EXP(DG),SMALL) 

      RF( 76) =  0.36000D+14 
      DG = -0.413433580E+00*(ALOGT-1.) -0.196446688E-03*T 
     &     +0.900981290E-07*T*T -0.138639969E-10*T*T*T 
     &     +0.820380585E-15*T*T*T*T +0.357110875E+05/T 
     &     +0.204911270E+01
      RB( 76) = RF( 76)/MAX(EXP(DG),SMALL) 

      RF( 77) =  0.10000D+14 * EXP( -12582./T) 
      DG = +0.380134660E+00*(ALOGT-1.) -0.632324555E-03*T 
     &     +0.159245840E-06*T*T -0.210413597E-10*T*T*T 
     &     +0.111048380E-14*T*T*T*T +0.647818967E+03/T 
     &     -0.109090440E+00
      RB( 77) = RF( 77)/MAX(EXP(DG),SMALL) 

      RF( 78) =  0.25000D+16 * EXP( -25164./T) 
      DG = +0.185291462E+01*(ALOGT-1.) -0.694112927E-03*T 
     &     +0.444218010E-07*T*T -0.445985728E-12*T*T*T 
     &     -0.105078440E-15*T*T*T*T -0.309733885E+05/T 
     &     +0.246834227E+01
      RB( 78) = RF( 78)/MAX(EXP(DG),SMALL) * (82.06*T) 

      RF( 79) =  0.30000D+08 * EXP(  2.00*ALOGT  -1007./T) 
      DG = +0.190193800E+00*(ALOGT-1.) -0.718814140E-03*T 
     &     +0.127664592E-06*T*T -0.154098511E-10*T*T*T 
     &     +0.896198390E-15*T*T*T*T +0.209240902E+05/T 
     &     +0.156684784E+00
      RB( 79) = RF( 79)/MAX(EXP(DG),SMALL) 

      RF( 80) =  0.50000D+14 
      DG = -0.834599630E+00*(ALOGT-1.) -0.422486945E-03*T 
     &     +0.119910827E-06*T*T -0.152498537E-10*T*T*T 
     &     +0.854070670E-15*T*T*T*T +0.571242531E+04/T 
     &     +0.954855649E+01
      RB( 80) = RF( 80)/MAX(EXP(DG),SMALL) 

      RF( 81) =  0.30000D+08 * EXP(  2.00*ALOGT  -1007./T) 
      DG = -0.123618500E+00*(ALOGT-1.) -0.376911000E-03*T 
     &     +0.585145932E-07*T*T -0.694001664E-11*T*T*T 
     &     +0.420566569E-15*T*T*T*T +0.198591935E+05/T 
     &     +0.260738264E+01
      RB( 81) = RF( 81)/MAX(EXP(DG),SMALL) 

      RF( 82) =  0.60000D+14 * EXP( -15098./T) 
      DG = +0.714370230E+00*(ALOGT-1.) -0.878290119E-04*T 
     &     +0.344925630E-07*T*T -0.103696867E-10*T*T*T 
     &     +0.877997000E-15*T*T*T*T -0.144034592E+05/T 
     &     -0.529407013E+01
      RB( 82) = RF( 82)/MAX(EXP(DG),SMALL) 

      RF( 83) =  0.24000D+07 * EXP(  2.00*ALOGT   +600./T) 
      DG = -0.705980570E+00*(ALOGT-1.) +0.120131247E-03*T 
     &     -0.400782600E-08*T*T -0.120397128E-11*T*T*T 
     &     +0.148864615E-15*T*T*T*T +0.283632255E+05/T 
     &     +0.340509918E+01
      RB( 83) = RF( 83)/MAX(EXP(DG),SMALL) 

      RF( 84) =  0.60000D+12 * EXP(  -1007./T) 
      DG = -0.449475900E-01*(ALOGT-1.) +0.328721285E-03*T 
     &     -0.924237917E-07*T*T +0.111967230E-10*T*T*T 
     &     -0.590107715E-15*T*T*T*T +0.808541571E+04/T 
     &     +0.354401300E-01
      RB( 84) = RF( 84)/MAX(EXP(DG),SMALL) 

      RF( 85) =  0.12000D+14 * EXP(  -3700./T) 
      DG = +0.235141390E+00*(ALOGT-1.) -0.104753543E-02*T 
     &     +0.220088384E-06*T*T -0.266065741E-10*T*T*T 
     &     +0.148630611E-14*T*T*T*T +0.128386745E+05/T 
     &     +0.121244654E+00
      RB( 85) = RF( 85)/MAX(EXP(DG),SMALL) 

      RF( 86) =  0.40000D+13 
      DG = -0.661032980E+00*(ALOGT-1.) -0.208590038E-03*T 
     &     +0.884159657E-07*T*T -0.124006943E-10*T*T*T 
     &     +0.738972330E-15*T*T*T*T +0.202778098E+05/T 
     &     +0.336965905E+01
      RB( 86) = RF( 86)/MAX(EXP(DG),SMALL) 

      RF( 87) =  0.22000D+17 * EXP( -47041./T) 
      DG = +0.270029001E+01*(ALOGT-1.) -0.122947391E-02*T 
     &     +0.132326596E-06*T*T -0.847017373E-11*T*T*T 
     &     +0.232908605E-15*T*T*T*T -0.541903127E+05/T 
     &     -0.492559414E+00
      RB( 87) = RF( 87)/MAX(EXP(DG),SMALL) * (82.06*T) 

      RF( 88) =  0.64000D+06 * EXP(  2.39*ALOGT  -5119./T) 
      DG = +0.103756919E+01*(ALOGT-1.) -0.125417513E-02*T 
     &     +0.215569387E-06*T*T -0.234340391E-10*T*T*T 
     &     +0.123418544E-14*T*T*T*T -0.229283398E+04/T 
     &     -0.280421690E+01
      RB( 88) = RF( 88)/MAX(EXP(DG),SMALL) 

      RF( 89) =  0.94000D+07 * EXP(  1.94*ALOGT  -3251./T) 
      DG = +0.723756890E+00*(ALOGT-1.) -0.912271985E-03*T 
     &     +0.146419388E-06*T*T -0.149642046E-10*T*T*T 
     &     +0.758553614E-15*T*T*T*T -0.335773070E+04/T 
     &     -0.353519040E+00
      RB( 89) = RF( 89)/MAX(EXP(DG),SMALL) 

      RF( 90) =  0.20000D+07 * EXP(  2.04*ALOGT   -285./T) 
      DG = +0.141394820E+00*(ALOGT-1.) -0.415229738E-03*T 
     &     +0.838969690E-07*T*T -0.922815928E-11*T*T*T 
     &     +0.486851660E-15*T*T*T*T +0.514630130E+04/T 
     &     +0.444197500E+00
      RB( 90) = RF( 90)/MAX(EXP(DG),SMALL) 

      RF( 91) =  0.40000D+14 * EXP(  -1837./T) 
      DG = +0.786229890E+00*(ALOGT-1.) -0.963433825E-03*T 
     &     +0.168097389E-06*T*T -0.198539515E-10*T*T*T 
     &     +0.112208523E-14*T*T*T*T +0.647492782E+04/T 
     &     -0.353797680E+01
      RB( 91) = RF( 91)/MAX(EXP(DG),SMALL) 

      RF( 92) =  0.46000D+14 
      DG = +0.710981130E+00*(ALOGT-1.) +0.455759453E-04*T 
     &     -0.613962338E-07*T*T +0.830983710E-11*T*T*T 
     &     -0.433504101E-15*T*T*T*T +0.141467682E+05/T 
     &     -0.694117385E+01
      RB( 92) = RF( 92)/MAX(EXP(DG),SMALL) 

      RF( 93) =  0.68000D+13 
      DG = +0.472417590E+00*(ALOGT-1.) -0.621530685E-03*T 
     &     +0.989473898E-07*T*T -0.113841171E-10*T*T*T 
     &     +0.646453405E-15*T*T*T*T +0.541003110E+04/T 
     &     -0.108727894E+01
      RB( 93) = RF( 93)/MAX(EXP(DG),SMALL) 

      RF( 94) =  0.24000D+07 * EXP(  2.00*ALOGT    -25./T) 
      DG = -0.109944480E+00*(ALOGT-1.) -0.124488438E-03*T 
     &     +0.364249707E-07*T*T -0.564807170E-11*T*T*T 
     &     +0.374751450E-15*T*T*T*T +0.139140631E+05/T 
     &     -0.289562400E+00
      RB( 94) = RF( 94)/MAX(EXP(DG),SMALL) 

      RF( 95) =  0.20000D+14 
      DG = +0.503164060E+00*(ALOGT-1.) -0.148993392E-03*T 
     &     +0.278672713E-07*T*T -0.892836930E-11*T*T*T 
     &     +0.859595740E-15*T*T*T*T +0.121557773E+05/T 
     &     -0.393358353E+01
      RB( 95) = RF( 95)/MAX(EXP(DG),SMALL) 

      RF( 96) =  0.10000D+14 
      DG = -0.934963060E+00*(ALOGT-1.) +0.851107605E-03*T 
     &     -0.153044680E-06*T*T +0.164055221E-10*T*T*T 
     &     -0.776954875E-15*T*T*T*T +0.299169672E+05/T 
     &     +0.171400564E+01
      RB( 96) = RF( 96)/MAX(EXP(DG),SMALL) 

      RF( 97) =  0.25000D+16 * EXP( -1.20*ALOGT   +107./T) 
      DG = +0.764294370E+00*(ALOGT-1.) -0.479237292E-03*T 
     &     +0.746536247E-07*T*T -0.105006179E-10*T*T*T 
     &     +0.680178295E-15*T*T*T*T -0.416422400E+03/T 
     &     -0.394251620E+01
      RB( 97) = RF( 97)/MAX(EXP(DG),SMALL) 

      RF( 98) =  0.43400D+16 * EXP( -1.20*ALOGT   +107./T) 
      DG = -0.134715210E+00*(ALOGT-1.) -0.366758830E-03*T 
     &     +0.105088334E-06*T*T -0.168799005E-10*T*T*T 
     &     +0.110104863E-14*T*T*T*T +0.630200264E+05/T 
     &     -0.194242090E+01
      RB( 98) = RF( 98)/MAX(EXP(DG),SMALL) 

      RF( 99) =  0.16600D+15 * EXP( -0.74*ALOGT) 
      DG = +0.137569090E+00*(ALOGT-1.) -0.287880330E-03*T 
     &     +0.106565974E-06*T*T -0.192949581E-10*T*T*T 
     &     +0.127381891E-14*T*T*T*T +0.464193476E+05/T 
     &     -0.363794995E+01
      RB( 99) = RF( 99)/MAX(EXP(DG),SMALL) 

      RF(100) =  0.35000D+13 
      DG = -0.196661310E+00*(ALOGT-1.) +0.206056165E-03*T 
     &     -0.323435283E-07*T*T +0.586221917E-12*T*T*T 
     &     +0.568554400E-16*T*T*T*T +0.829810129E+04/T 
     &     +0.186154217E+01
      RB(100) = RF(100)/MAX(EXP(DG),SMALL) 

      RF(101) =  0.30000D+13 * EXP(   -503./T) 
      DG = -0.847375390E+00*(ALOGT-1.) +0.535360985E-03*T 
     &     -0.879047950E-07*T*T +0.802418800E-11*T*T*T 
     &     -0.337987045E-15*T*T*T*T +0.232169242E+05/T 
     &     +0.296090168E+01
      RB(101) = RF(101)/MAX(EXP(DG),SMALL) 

      RF(102) =  0.85000D+12 
      DG = +0.861484300E+00*(ALOGT-1.) -0.622176865E-04*T 
     &     -0.138138570E-07*T*T -0.146239367E-11*T*T*T 
     &     +0.328703820E-15*T*T*T*T +0.209404559E+05/T 
     &     -0.101601301E+02
      RB(102) = RF(102)/MAX(EXP(DG),SMALL) 

      RF(103) =  0.50000D+14 
      DG = +0.752544100E-01*(ALOGT-1.) +0.901216138E-03*T 
     &     -0.181911246E-06*T*T +0.183915579E-10*T*T*T 
     &     -0.793381405E-15*T*T*T*T +0.144655281E+05/T 
     &     -0.662215331E+01
      RB(103) = RF(103)/MAX(EXP(DG),SMALL) 

      RF(104) =  0.30000D+14 
      DG = +0.469529290E+00*(ALOGT-1.) -0.602177400E-03*T 
     &     +0.134202257E-06*T*T -0.189740538E-10*T*T*T 
     &     +0.117469430E-14*T*T*T*T +0.124108938E+05/T 
     &     -0.384951070E+01
      RB(104) = RF(104)/MAX(EXP(DG),SMALL) 

      RF(105) =  0.92000D+14 
      DG = +0.407492030E+00*(ALOGT-1.) -0.263823047E-04*T 
     &     -0.772307546E-08*T*T +0.939119304E-13*T*T*T 
     &     +0.121250344E-16*T*T*T*T +0.359437926E+05/T 
     &     -0.460249875E+01
      RB(105) = RF(105)/MAX(EXP(DG),SMALL) 

      RF(106) =  0.20000D+14 
      DG = +0.238563540E+00*(ALOGT-1.) +0.667106630E-03*T 
     &     -0.160343624E-06*T*T +0.196939542E-10*T*T*T 
     &     -0.107995751E-14*T*T*T*T +0.873673710E+04/T 
     &     -0.585389491E+01
      RB(106) = RF(106)/MAX(EXP(DG),SMALL) 

      RF(107) =  0.50000D+12 * EXP(  0.50*ALOGT  -1007./T) 
      DG = -0.426645080E+00*(ALOGT-1.) +0.236767987E-03*T 
     &     +0.252983900E-08*T*T -0.476817395E-11*T*T*T 
     &     +0.427360520E-15*T*T*T*T +0.198500291E+05/T 
     &     -0.601096300E+00
      RB(107) = RF(107)/MAX(EXP(DG),SMALL) 

      RF(108) =  0.46000D+06 * EXP(  2.00*ALOGT  -3271./T) 
      DG = +0.118334140E+00*(ALOGT-1.) +0.156790673E-03*T 
     &     -0.594023368E-08*T*T -0.592558633E-11*T*T*T 
     &     +0.652110165E-15*T*T*T*T +0.457031800E+02/T 
     &     -0.159940855E+01
      RB(108) = RF(108)/MAX(EXP(DG),SMALL) 

      RF(109) =  0.13000D+07 * EXP(  1.50*ALOGT    -50./T) 
      DG = +0.287262630E+00*(ALOGT-1.) -0.536698262E-03*T 
     &     +0.146680314E-06*T*T -0.255256286E-10*T*T*T 
     &     +0.174419271E-14*T*T*T*T +0.272527587E+05/T 
     &     -0.348012390E+00
      RB(109) = RF(109)/MAX(EXP(DG),SMALL) 

      RF(110) =  0.27400D+16 * EXP( -0.78*ALOGT    -40./T) 
      DG = +0.127877451E+01*(ALOGT-1.) +0.530388885E-04*T 
     &     -0.174262956E-07*T*T +0.101674544E-11*T*T*T 
     &     -0.118586451E-16*T*T*T*T +0.184947579E+05/T 
     &     -0.147584862E+02
      RB(110) = RF(110)/MAX(EXP(DG),SMALL) 

      RF(111) =  0.68400D+15 * EXP( -0.78*ALOGT    -40./T) 
      DG = -0.247707300E-01*(ALOGT-1.) -0.242270392E-03*T 
     &     +0.686633630E-07*T*T -0.112318288E-10*T*T*T 
     &     +0.726297175E-15*T*T*T*T +0.491059633E+05/T 
     &     -0.165285850E+01
      RB(111) = RF(111)/MAX(EXP(DG),SMALL) 

      RF(112) =  0.38000D+14 
      DG = +0.251775040E+00*(ALOGT-1.) +0.233891955E-03*T 
     &     -0.727753336E-07*T*T +0.105981312E-10*T*T*T 
     &     -0.686937440E-15*T*T*T*T +0.245977955E+05/T 
     &     -0.320368591E+01
      RB(112) = RF(112)/MAX(EXP(DG),SMALL) 

      RF(113) =  0.64000D+10 * EXP(  1.00*ALOGT  -3161./T) 
      DG = +0.131545640E+00*(ALOGT-1.) -0.276424002E-03*T 
     &     +0.816280563E-07*T*T -0.150214092E-10*T*T*T 
     &     +0.104513023E-14*T*T*T*T +0.159067616E+05/T 
     &     +0.105080045E+01
      RB(113) = RF(113)/MAX(EXP(DG),SMALL) 

      RF(114) =  0.33000D+13 * EXP(  0.30*ALOGT) 
      DG = -0.180487720E+00*(ALOGT-1.) +0.180038681E-04*T 
     &     +0.361110482E-08*T*T -0.727609492E-12*T*T*T 
     &     +0.272347005E-16*T*T*T*T +0.377599662E+05/T 
     &     -0.254045660E+00
      RB(114) = RF(114)/MAX(EXP(DG),SMALL) 

      RF(115) =  0.50000D+17 * EXP( -25164./T) 
      DG = +0.307306511E+01*(ALOGT-1.) -0.172403766E-02*T 
     &     +0.234749003E-06*T*T -0.223042329E-10*T*T*T 
     &     +0.960039725E-15*T*T*T*T -0.297707399E+05/T 
     &     -0.206190201E+01
      RB(115) = RF(115)/MAX(EXP(DG),SMALL) * (82.06*T) 

      RF(116) =  0.30000D+13 
      DG = +0.120352010E+01*(ALOGT-1.) -0.848177250E-03*T 
     &     +0.164484950E-06*T*T -0.173748124E-10*T*T*T 
     &     +0.781522760E-15*T*T*T*T +0.402922980E+04/T 
     &     -0.813633290E+01
      RB(116) = RF(116)/MAX(EXP(DG),SMALL) 

      RF(117) =  0.50000D+14 * EXP(   -503./T) 
      DG = +0.141034429E+01*(ALOGT-1.) -0.174873888E-02*T 
     &     +0.317991794E-06*T*T -0.372680983E-10*T*T*T 
     &     +0.196131656E-14*T*T*T*T +0.221267388E+05/T 
     &     -0.437355950E+01
      RB(117) = RF(117)/MAX(EXP(DG),SMALL) 

      RF(118) =  0.10000D+14 * EXP(   -503./T) 
      DG = +0.514169920E+00*(ALOGT-1.) -0.909793488E-03*T 
     &     +0.186319376E-06*T*T -0.230622185E-10*T*T*T 
     &     +0.121398278E-14*T*T*T*T +0.295658741E+05/T 
     &     -0.112514510E+01
      RB(118) = RF(118)/MAX(EXP(DG),SMALL) 

      RF(119) =  0.10000D+14 * EXP(   -503./T) 
      DG = +0.372775100E+00*(ALOGT-1.) -0.494563750E-03*T 
     &     +0.102422407E-06*T*T -0.138340592E-10*T*T*T 
     &     +0.727131120E-15*T*T*T*T +0.244195728E+05/T 
     &     -0.156934260E+01
      RB(119) = RF(119)/MAX(EXP(DG),SMALL) 

      RF(120) =  0.10000D+08 
      DG = +0.165988561E+01*(ALOGT-1.) -0.701765712E-03*T 
     &     +0.788643360E-07*T*T -0.562129706E-11*T*T*T 
     &     +0.166927275E-15*T*T*T*T +0.409983480E+04/T 
     &     +0.106333839E+01
      RB(120) = RF(120)/MAX(EXP(DG),SMALL) * (82.06*T) 

      RF(121) =  0.20000D+15 
      DG = -0.105441340E+00*(ALOGT-1.) -0.323399405E-03*T 
     &     +0.995824198E-07*T*T -0.135566454E-10*T*T*T 
     &     +0.710973545E-15*T*T*T*T +0.283731803E+05/T 
     &     -0.158107840E+00
      RB(121) = RF(121)/MAX(EXP(DG),SMALL) 

      RF(122) =  0.50000D+14 
      DG = +0.165988561E+01*(ALOGT-1.) -0.701765712E-03*T 
     &     +0.788643360E-07*T*T -0.562129706E-11*T*T*T 
     &     +0.166927275E-15*T*T*T*T +0.409983480E+04/T 
     &     +0.106333839E+01
      RB(122) = RF(122)/MAX(EXP(DG),SMALL) * (82.06*T) 

      RF(123) =  0.40000D+15 * EXP( -28234./T) 
      DG = +0.672987880E+00*(ALOGT-1.) -0.612511207E-03*T 
     &     +0.719968665E-07*T*T -0.575454333E-11*T*T*T 
     &     +0.212510810E-15*T*T*T*T -0.202213766E+05/T 
     &     +0.129665873E+02
      RB(123) = RF(123)/MAX(EXP(DG),SMALL) * (82.06*T) 

      RF(124) =  0.78000D+15 * EXP(  -9814./T) 
      DG = -0.130354524E+01*(ALOGT-1.) -0.295309280E-03*T 
     &     +0.860896586E-07*T*T -0.122485742E-10*T*T*T 
     &     +0.738155820E-15*T*T*T*T +0.306112054E+05/T 
     &     +0.131056277E+02
      RB(124) = RF(124)/MAX(EXP(DG),SMALL) 

      RF(125) =  0.30000D+06 * EXP(  2.45*ALOGT  -1126./T) 
      DG = -0.115680269E+01*(ALOGT-1.) +0.196507084E-02*T 
     &     -0.399609100E-06*T*T +0.470345663E-10*T*T*T 
     &     -0.243837272E-14*T*T*T*T +0.102568492E+05/T 
     &     +0.357937530E+01
      RB(125) = RF(125)/MAX(EXP(DG),SMALL) 

      RF(126) =  0.14000D+05 * EXP(  2.64*ALOGT  -2506./T) 
      DG = +0.257646823E+01*(ALOGT-1.) -0.145771650E-02*T 
     &     +0.243455636E-06*T*T -0.294850046E-10*T*T*T 
     &     +0.160953043E-14*T*T*T*T +0.500123920E+04/T 
     &     -0.173461096E+02
      RB(126) = RF(126)/MAX(EXP(DG),SMALL) 

      RF(127) =  0.35000D+04 * EXP(  2.64*ALOGT  -2506./T) 
      DG = -0.496999370E+00*(ALOGT-1.) -0.222922728E-03*T 
     &     +0.225972167E-07*T*T +0.185224075E-12*T*T*T 
     &     -0.138452776E-15*T*T*T*T +0.161041235E+05/T 
     &     +0.516733398E+01
      RB(127) = RF(127)/MAX(EXP(DG),SMALL) 

      RF(128) =  0.39000D+07 * EXP(  1.83*ALOGT  -5184./T) 
      DG = +0.260628320E+00*(ALOGT-1.) -0.112612545E-02*T 
     &     +0.267936682E-06*T*T -0.328286864E-10*T*T*T 
     &     +0.169103895E-14*T*T*T*T -0.281771390E+04/T 
     &     -0.330960900E+00
      RB(128) = RF(128)/MAX(EXP(DG),SMALL) 

      RF(129) =  0.59000D+05 * EXP(  2.40*ALOGT  -6291./T) 
      DG = +0.200276734E+01*(ALOGT-1.) -0.127154487E-02*T 
     &     +0.221521648E-06*T*T -0.277359624E-10*T*T*T 
     &     +0.152677014E-14*T*T*T*T -0.271916190E+04/T 
     &     -0.143976070E+02
      RB(129) = RF(129)/MAX(EXP(DG),SMALL) 

      RF(130) =  0.20000D-02 * EXP(  4.00*ALOGT   -503./T) 
      DG = +0.252596234E+01*(ALOGT-1.) -0.125342987E-02*T 
     &     +0.219727731E-06*T*T -0.279735208E-10*T*T*T 
     &     +0.156670314E-14*T*T*T*T +0.988824710E+04/T 
     &     -0.185699740E+02
      RB(130) = RF(130)/MAX(EXP(DG),SMALL) 

      RF(131) =  0.40000D+14 
      DG = +0.173347784E+01*(ALOGT-1.) +0.165451193E-03*T 
     &     -0.870034653E-07*T*T +0.907972724E-11*T*T*T 
     &     -0.353210470E-15*T*T*T*T +0.163229851E+05/T 
     &     -0.162174321E+02
      RB(131) = RF(131)/MAX(EXP(DG),SMALL) 

      RF(132) =  0.75000D+13 * EXP(    196./T) 
      DG = +0.161324844E+01*(ALOGT-1.) -0.344864764E-03*T 
     &     +0.673999247E-07*T*T -0.165398132E-10*T*T*T 
     &     +0.137885720E-14*T*T*T*T +0.763195118E+04/T 
     &     -0.119629458E+02
      RB(132) = RF(132)/MAX(EXP(DG),SMALL) 

      RF(133) =  0.22000D+08 * EXP(  1.70*ALOGT  -1912./T) 
      DG = -0.349537930E+01*(ALOGT-1.) +0.165203783E-02*T 
     &     -0.296077905E-06*T*T +0.395428619E-10*T*T*T 
     &     -0.235160932E-14*T*T*T*T +0.805845300E+03/T 
     &     +0.248245869E+02
      RB(133) = RF(133)/MAX(EXP(DG),SMALL) 

      RF(134) =  0.98000D+08 * EXP(  1.41*ALOGT  -4290./T) 
      DG = -0.247357469E+01*(ALOGT-1.) +0.193221369E-02*T 
     &     -0.420879428E-06*T*T +0.599406891E-10*T*T*T 
     &     -0.363440002E-14*T*T*T*T +0.192081671E+05/T 
     &     +0.132668793E+02
      RB(134) = RF(134)/MAX(EXP(DG),SMALL) 

      RF(135) =  0.32700D+08 * EXP(  1.50*ALOGT  -1812./T) 
      DG = -0.531856180E+00*(ALOGT-1.) +0.292755612E-03*T 
     &     -0.387945143E-07*T*T +0.422456155E-11*T*T*T 
     &     -0.228874670E-15*T*T*T*T +0.361702410E+04/T 
     &     +0.202158100E+01
      RB(135) = RF(135)/MAX(EXP(DG),SMALL) 

      RF(136) =  0.36000D+07 * EXP(  1.50*ALOGT  -1812./T) 
      DG = -0.294599228E+01*(ALOGT-1.) +0.255374437E-02*T 
     &     -0.519826818E-06*T*T +0.713248061E-10*T*T*T 
     &     -0.428085342E-14*T*T*T*T +0.137981360E+05/T 
     &     +0.143541583E+02
      RB(136) = RF(136)/MAX(EXP(DG),SMALL) 

      RF(137) =  0.50000D+13 * EXP(  -3120./T) 
      DG = -0.673251000E+00*(ALOGT-1.) +0.707985350E-03*T 
     &     -0.122691483E-06*T*T +0.134527208E-10*T*T*T 
     &     -0.715726330E-15*T*T*T*T -0.152927720E+04/T 
     &     +0.157738350E+01
      RB(137) = RF(137)/MAX(EXP(DG),SMALL) 

      RF(138) =  0.20000D+08 * EXP(  2.00*ALOGT  -1007./T) 
      DG = +0.887513190E+00*(ALOGT-1.) -0.528074775E-03*T 
     &     +0.910839870E-07*T*T -0.102188766E-10*T*T*T 
     &     +0.558392105E-15*T*T*T*T +0.878529782E+04/T 
     &     -0.539920040E+01
      RB(138) = RF(138)/MAX(EXP(DG),SMALL) 

      RF(139) =  0.64000D+06 * EXP(  2.00*ALOGT  -1290./T) 
      DG = -0.866118000E-02*(ALOGT-1.) +0.310870612E-03*T 
     &     -0.405884310E-07*T*T +0.398700322E-11*T*T*T 
     &     -0.188941670E-15*T*T*T*T +0.162244331E+05/T 
     &     -0.215078600E+01
      RB(139) = RF(139)/MAX(EXP(DG),SMALL) 

      RF(140) =  0.20000D+15 
      DG = -0.313524536E+01*(ALOGT-1.) +0.160312371E-02*T 
     &     -0.313088601E-06*T*T +0.414142249E-10*T*T*T 
     &     -0.242344394E-14*T*T*T*T +0.332508203E+05/T 
     &     +0.217139689E+02
      RB(140) = RF(140)/MAX(EXP(DG),SMALL) 

      RF(141) =  0.40000D+14 
      DG = -0.467019655E+01*(ALOGT-1.) +0.340744457E-02*T 
     &     -0.650203396E-06*T*T +0.826791180E-10*T*T*T 
     &     -0.470082362E-14*T*T*T*T +0.258996486E+05/T 
     &     +0.258795895E+02
      RB(141) = RF(141)/MAX(EXP(DG),SMALL) 

      RF(142) =  0.50000D+14 
      DG = -0.307346760E+01*(ALOGT-1.) +0.123479378E-02*T 
     &     -0.220858420E-06*T*T +0.296702286E-10*T*T*T 
     &     -0.174798320E-14*T*T*T*T +0.111028843E+05/T 
     &     +0.225134435E+02
      RB(142) = RF(142)/MAX(EXP(DG),SMALL) 

      RF(143) =  0.43200D+14 
      DG = -0.266597557E+01*(ALOGT-1.) +0.120841147E-02*T 
     &     -0.228581495E-06*T*T +0.297641406E-10*T*T*T 
     &     -0.173585817E-14*T*T*T*T +0.470466769E+05/T 
     &     +0.179109448E+02
      RB(143) = RF(143)/MAX(EXP(DG),SMALL) 

      RF(144) =  0.76000D+03 * EXP(  3.00*ALOGT  -2013./T) 
      DG = -0.364318190E+00*(ALOGT-1.) +0.546189775E-03*T 
     &     -0.928779036E-07*T*T +0.998131827E-11*T*T*T 
     &     -0.518459105E-15*T*T*T*T +0.382211118E+04/T 
     &     +0.122683340E+01
      RB(144) = RF(144)/MAX(EXP(DG),SMALL) 

      RF(145) =  0.80000D+18 * EXP( -1.73*ALOGT   -380./T) 
      DG = -0.179469309E+01*(ALOGT-1.) +0.128783267E-02*T 
     &     -0.238284715E-06*T*T +0.306869741E-10*T*T*T 
     &     -0.175984185E-14*T*T*T*T +0.295976422E+05/T 
     &     +0.775495732E+01
      RB(145) = RF(145)/MAX(EXP(DG),SMALL) 

      RF(146) =  0.58000D+18 * EXP( -1.73*ALOGT   -380./T) 
      DG = -0.254885131E+01*(ALOGT-1.) +0.189422993E-02*T 
     &     -0.375943970E-06*T*T +0.502203441E-10*T*T*T 
     &     -0.295093013E-14*T*T*T*T +0.732011383E+05/T 
     &     +0.103901564E+02
      RB(146) = RF(146)/MAX(EXP(DG),SMALL) 

      RF(147) =  0.10000D+15 * EXP(  -4902./T) 
      DG = +0.819770110E+00*(ALOGT-1.) -0.222739130E-02*T 
     &     +0.553961367E-06*T*T -0.854460043E-10*T*T*T 
     &     +0.575378101E-14*T*T*T*T +0.234249330E+04/T 
     &     +0.138931370E+01
      RB(147) = RF(147)/MAX(EXP(DG),SMALL) 

      RF(148) =  0.10000D+15 
      DG = -0.166856480E+01*(ALOGT-1.) -0.595451457E-04*T 
     &     +0.906926919E-07*T*T -0.159348005E-10*T*T*T 
     &     +0.115444682E-14*T*T*T*T +0.208487094E+05/T 
     &     +0.106503030E+02
      RB(148) = RF(148)/MAX(EXP(DG),SMALL) 

      RF(149) =  0.20000D+08 * EXP(  2.00*ALOGT  -1007./T) 
      DG = +0.215451482E+01*(ALOGT-1.) -0.227093434E-02*T 
     &     +0.484973176E-06*T*T -0.672777860E-10*T*T*T 
     &     +0.414545008E-14*T*T*T*T +0.131914058E+05/T 
     &     -0.964672050E+01
      RB(149) = RF(149)/MAX(EXP(DG),SMALL) 

      RF(150) =  0.30000D+15 * EXP( -11072./T) 
      DG = +0.716754010E+00*(ALOGT-1.) +0.477449438E-03*T 
     &     -0.219263179E-06*T*T +0.399792996E-10*T*T*T 
     &     -0.273805443E-14*T*T*T*T -0.126520709E+05/T 
     &     +0.760559009E+01
      RB(150) = RF(150)/MAX(EXP(DG),SMALL) * (82.06*T) 

      RF(151) =  0.32000D+12 * EXP( -10185./T) 
      DG = -0.754158220E+00*(ALOGT-1.) +0.606397265E-03*T 
     &     -0.137659255E-06*T*T +0.195333700E-10*T*T*T 
     &     -0.119108828E-14*T*T*T*T +0.436034961E+05/T 
     &     +0.263519904E+01
      RB(151) = RF(151)/MAX(EXP(DG),SMALL) 

      RF(152) =  0.19000D+12 * EXP(  -1711./T) 
      DG = -0.297611980E+00*(ALOGT-1.) -0.667814590E-03*T 
     &     +0.150973580E-06*T*T -0.211838130E-10*T*T*T 
     &     +0.124230666E-14*T*T*T*T +0.116055048E+05/T 
     &     +0.726674276E+01
      RB(152) = RF(152)/MAX(EXP(DG),SMALL) 

      RF(153) =  0.60000D+13 
      DG = +0.148444530E+01*(ALOGT-1.) -0.160003423E-02*T 
     &     +0.313386890E-06*T*T -0.397534093E-10*T*T*T 
     &     +0.223455240E-14*T*T*T*T +0.122381796E+05/T 
     &     -0.538948508E+01
      RB(153) = RF(153)/MAX(EXP(DG),SMALL) 

      RF(154) =  0.72000D+13 
      DG = +0.297549630E+00*(ALOGT-1.) -0.436523065E-03*T 
     &     +0.109077307E-06*T*T -0.168197199E-10*T*T*T 
     &     +0.108286747E-14*T*T*T*T +0.163156572E+05/T 
     &     -0.353869590E+01
      RB(154) = RF(154)/MAX(EXP(DG),SMALL) 

      RF(155) =  0.62000D+05 * EXP(  2.64*ALOGT   +220./T) 
      DG = +0.189134273E+01*(ALOGT-1.) -0.113515355E-02*T 
     &     +0.141252795E-06*T*T -0.201993583E-10*T*T*T 
     &     +0.131789367E-14*T*T*T*T +0.104367391E+05/T 
     &     -0.913621131E+01
      RB(155) = RF(155)/MAX(EXP(DG),SMALL) 

      RF(156) =  0.98000D+13 * EXP(  -4087./T) 
      DG = +0.268382723E+01*(ALOGT-1.) -0.255403462E-02*T 
     &     +0.447983992E-06*T*T -0.572526063E-10*T*T*T 
     &     +0.323780728E-14*T*T*T*T +0.400200111E+04/T 
     &     -0.114887532E+02
      RB(156) = RF(156)/MAX(EXP(DG),SMALL) 

      RF(157) =  0.15000D+00 * EXP(  3.52*ALOGT  -1988./T) 
      DG = +0.914126170E+00*(ALOGT-1.) -0.116494832E-02*T 
     &     +0.275183907E-06*T*T -0.371656417E-10*T*T*T 
     &     +0.216203243E-14*T*T*T*T +0.418453360E+05/T 
     &     -0.627480939E+01
      RB(157) = RF(157)/MAX(EXP(DG),SMALL) 

      RF(158) =  0.15000D+00 * EXP(  3.52*ALOGT  -1988./T) 
      DG = +0.275626735E+01*(ALOGT-1.) -0.245664193E-02*T 
     &     +0.542876713E-06*T*T -0.763869557E-10*T*T*T 
     &     +0.464614380E-14*T*T*T*T -0.483920710E+04/T 
     &     -0.148171564E+02
      RB(158) = RF(158)/MAX(EXP(DG),SMALL) 

      RF(159) =  0.14000D+14 
      DG = -0.139067320E+00*(ALOGT-1.) -0.174847335E-03*T 
     &     +0.121609677E-06*T*T -0.241049146E-10*T*T*T 
     &     +0.160332823E-14*T*T*T*T +0.904327548E+04/T 
     &     +0.936222860E+00
      RB(159) = RF(159)/MAX(EXP(DG),SMALL) 

      RF(160) =  0.71000D+14 
      DG = +0.300804239E+01*(ALOGT-1.) -0.222274997E-02*T 
     &     +0.470101379E-06*T*T -0.657888245E-10*T*T*T 
     &     +0.395920635E-14*T*T*T*T +0.197585884E+05/T 
     &     -0.180208423E+02
      RB(160) = RF(160)/MAX(EXP(DG),SMALL) 

      RF(161) =  0.50000D+13 * EXP(  -1007./T) 
      DG = -0.282253745E+01*(ALOGT-1.) +0.349893617E-02*T 
     &     -0.775483015E-06*T*T +0.113181967E-09*T*T*T 
     &     -0.728055115E-14*T*T*T*T +0.376668600E+03/T 
     &     +0.130082933E+02
      RB(161) = RF(161)/MAX(EXP(DG),SMALL) 

      RF(162) =  0.22000D+13 * EXP(    190./T) 
      DG = -0.715012930E+00*(ALOGT-1.) +0.835882608E-03*T 
     &     -0.156268517E-06*T*T +0.204327248E-10*T*T*T 
     &     -0.114274242E-14*T*T*T*T +0.368157706E+05/T 
     &     +0.163963867E+01
      RB(162) = RF(162)/MAX(EXP(DG),SMALL) 

      RF(163) =  0.13000D+13 * EXP(    190./T) 
      DG = +0.186424697E+01*(ALOGT-1.) +0.113465868E-03*T 
     &     -0.900036290E-07*T*T +0.157173980E-10*T*T*T 
     &     -0.102057624E-14*T*T*T*T +0.120901538E+05/T 
     &     -0.172221952E+02
      RB(163) = RF(163)/MAX(EXP(DG),SMALL) 

      RF(164) =  0.15000D-08 * EXP(  6.89*ALOGT  +1465./T) 
      DG = +0.813599520E+01*(ALOGT-1.) -0.759426145E-02*T 
     &     +0.159855405E-05*T*T -0.215574148E-09*T*T*T 
     &     +0.126251622E-13*T*T*T*T +0.855426780E+04/T 
     &     -0.397512037E+02
      RB(164) = RF(164)/MAX(EXP(DG),SMALL) 

      RF(165) =  0.59000D+15 * EXP( -0.24*ALOGT) 
      DG = +0.483053304E+01*(ALOGT-1.) -0.365470384E-02*T 
     &     +0.753037930E-06*T*T -0.107850020E-09*T*T*T 
     &     +0.668809588E-14*T*T*T*T +0.939092714E+04/T 
     &     -0.254619074E+02
      RB(165) = RF(165)/MAX(EXP(DG),SMALL) 

      RF(166) =  0.10000D+13 * EXP(   -503./T) 
      DG = -0.268405966E+01*(ALOGT-1.) +0.226389007E-02*T 
     &     -0.484960218E-06*T*T +0.692919768E-10*T*T*T 
     &     -0.421250227E-14*T*T*T*T +0.444549440E+04/T 
     &     +0.131882531E+02
      RB(166) = RF(166)/MAX(EXP(DG),SMALL) 

      RF(167) =  0.72000D+13 
      DG = +0.424978770E+00*(ALOGT-1.) -0.930653550E-04*T 
     &     -0.887400850E-08*T*T +0.383344558E-11*T*T*T 
     &     -0.353363440E-15*T*T*T*T +0.254339281E+05/T 
     &     -0.495101538E+01
      RB(167) = RF(167)/MAX(EXP(DG),SMALL) 

      RF(168) =  0.16000D+14 
      DG = -0.160489411E+01*(ALOGT-1.) +0.153105793E-02*T 
     &     -0.298887810E-06*T*T +0.403307166E-10*T*T*T 
     &     -0.240477373E-14*T*T*T*T +0.631518356E+05/T 
     &     +0.344038986E+01
      RB(168) = RF(168)/MAX(EXP(DG),SMALL) 

      RF(169) =  0.40000D+14 
      DG = +0.197796744E+01*(ALOGT-1.) -0.696266957E-03*T 
     &     +0.155428496E-06*T*T -0.238617332E-10*T*T*T 
     &     +0.148042897E-14*T*T*T*T +0.247822072E+05/T 
     &     -0.158344089E+02
      RB(169) = RF(169)/MAX(EXP(DG),SMALL) 

      RF(170) =  0.64000D+53 * EXP(-11.99*ALOGT  -8450./T) 
      DG = +0.113196360E+00*(ALOGT-1.) +0.202669665E-03*T 
     &     +0.448331267E-07*T*T -0.143468026E-10*T*T*T 
     &     +0.111903547E-14*T*T*T*T +0.444758371E+05/T 
     &     -0.202155617E+02
      RB(170) = RF(170)/MAX(EXP(DG),SMALL) / (82.06*T) 

      RF(171) =  0.61000D+29 * EXP( -5.69*ALOGT -12215./T) 
      DG = +0.156275868E+01*(ALOGT-1.) -0.489844142E-03*T 
     &     +0.218447094E-07*T*T +0.336672319E-11*T*T*T 
     &     -0.401001470E-15*T*T*T*T -0.886719940E+04/T 
     &     +0.193753104E+01
      RB(171) = RF(171)/MAX(EXP(DG),SMALL) * (82.06*T) 

      RF(172) =  0.24000D+09 * EXP(  1.50*ALOGT  -3685./T) 
      DG = +0.149960424E+01*(ALOGT-1.) -0.212866042E-02*T 
     &     +0.448875140E-06*T*T -0.609011450E-10*T*T*T 
     &     +0.365347168E-14*T*T*T*T +0.738727498E+04/T 
     &     -0.603594413E+01
      RB(172) = RF(172)/MAX(EXP(DG),SMALL) 

      RF(173) =  0.12000D+07 * EXP(  2.00*ALOGT    +45./T) 
      DG = +0.603429870E+00*(ALOGT-1.) -0.128971503E-02*T 
     &     +0.317202722E-06*T*T -0.466952652E-10*T*T*T 
     &     +0.290613790E-14*T*T*T*T +0.148264103E+05/T 
     &     -0.278752973E+01
      RB(173) = RF(173)/MAX(EXP(DG),SMALL) 

      RF(174) =  0.82000D+06 * EXP(  1.87*ALOGT  -3585./T) 
      DG = -0.154854117E+01*(ALOGT-1.) +0.971563975E-03*T 
     &     -0.919867550E-07*T*T +0.633277958E-11*T*T*T 
     &     -0.102794705E-15*T*T*T*T +0.720738505E+04/T 
     &     +0.667964248E+01
      RB(174) = RF(174)/MAX(EXP(DG),SMALL) 

      RF(175) =  0.30000D+09 * EXP(  1.50*ALOGT  -3085./T) 
      DG = +0.653599570E+00*(ALOGT-1.) -0.116136684E-02*T 
     &     +0.207767252E-06*T*T -0.242885686E-10*T*T*T 
     &     +0.131641873E-14*T*T*T*T +0.360240348E+04/T 
     &     -0.367885076E+00
      RB(175) = RF(175)/MAX(EXP(DG),SMALL) 

      RF(176) =  0.24000D+07 * EXP(  2.00*ALOGT   -230./T) 
      DG = -0.242574800E+00*(ALOGT-1.) -0.322421448E-03*T 
     &     +0.760948340E-07*T*T -0.100826888E-10*T*T*T 
     &     +0.569084950E-15*T*T*T*T +0.110415388E+05/T 
     &     +0.288052932E+01
      RB(176) = RF(176)/MAX(EXP(DG),SMALL) 

      RF(177) =  0.17000D+07 * EXP(  2.08*ALOGT) 
      DG = -0.257241585E+01*(ALOGT-1.) +0.171930423E-02*T 
     &     -0.401660678E-06*T*T +0.589238375E-10*T*T*T 
     &     -0.362535202E-14*T*T*T*T +0.973312856E+04/T 
     &     +0.157796385E+02
      RB(177) = RF(177)/MAX(EXP(DG),SMALL) 

      RF(178) =  0.13000D+43 * EXP( -9.24*ALOGT -20805./T) 
      DG = +0.214209383E+01*(ALOGT-1.) -0.663323117E-03*T 
     &     +0.573686060E-07*T*T -0.346462414E-11*T*T*T 
     &     +0.959240799E-16*T*T*T*T -0.146600748E+05/T 
     &     -0.260655814E+00
      RB(178) = RF(178)/MAX(EXP(DG),SMALL) * (82.06*T) 

      RF(179) =  0.24000D+49 * EXP(-10.82*ALOGT -26190./T) 
      DG = +0.933055080E+00*(ALOGT-1.) -0.257026617E-03*T 
     &     +0.184840527E-07*T*T -0.835102395E-12*T*T*T 
     &     +0.138510099E-16*T*T*T*T -0.185307113E+05/T 
     &     +0.731698062E+01
      RB(179) = RF(179)/MAX(EXP(DG),SMALL) * (82.06*T) 

      RF(180) =  0.10000D+23 * EXP( -3.09*ALOGT  -3400./T) 
      DG = -0.832271870E+00*(ALOGT-1.) +0.121339690E-03*T 
     &     +0.392021365E-07*T*T -0.877045075E-11*T*T*T 
     &     +0.557897280E-15*T*T*T*T +0.574263421E+04/T 
     &     +0.609553439E+01
      RB(180) = RF(180)/MAX(EXP(DG),SMALL) 

      RF(181) =  0.56000D+09 * EXP(  1.50*ALOGT  -2750./T) 
      DG = +0.171794589E+01*(ALOGT-1.) -0.160534383E-02*T 
     &     +0.260208705E-06*T*T -0.285337798E-10*T*T*T 
     &     +0.146962772E-14*T*T*T*T +0.531200926E+04/T 
     &     -0.647380461E+01
      RB(181) = RF(181)/MAX(EXP(DG),SMALL) 

      RF(182) =  0.36000D+07 * EXP(  2.00*ALOGT   -120./T) 
      DG = +0.821771520E+00*(ALOGT-1.) -0.766398443E-03*T 
     &     +0.128536287E-06*T*T -0.143278999E-10*T*T*T 
     &     +0.722293940E-15*T*T*T*T +0.127511445E+05/T 
     &     -0.322539021E+01
      RB(182) = RF(182)/MAX(EXP(DG),SMALL) 

      RF(183) =  0.48000D+09 * EXP(  1.50*ALOGT  -4885./T) 
      DG = +0.508907140E+00*(ALOGT-1.) -0.119904733E-02*T 
     &     +0.221324152E-06*T*T -0.259042580E-10*T*T*T 
     &     +0.138755465E-14*T*T*T*T +0.144137276E+04/T 
     &     +0.110383182E+01
      RB(183) = RF(183)/MAX(EXP(DG),SMALL) 

      RF(184) =  0.24000D+07 * EXP(  2.00*ALOGT   -225./T) 
      DG = -0.387267230E+00*(ALOGT-1.) -0.360101943E-03*T 
     &     +0.896517340E-07*T*T -0.116983782E-10*T*T*T 
     &     +0.640220870E-15*T*T*T*T +0.888050804E+04/T 
     &     +0.435224622E+01
      RB(184) = RF(184)/MAX(EXP(DG),SMALL) 

      RF(185) =  0.12000D+37 * EXP( -5.50*ALOGT -28789./T) 
      DG = -0.170428100E+01*(ALOGT-1.) +0.183737000E-02*T 
     &     -0.434118333E-06*T*T +0.660092500E-10*T*T*T 
     &     -0.442935000E-14*T*T*T*T -0.741745000E+04/T 
     &     +0.976190600E+01
      RB(185) = RF(185)/MAX(EXP(DG),SMALL) 

      RF(186) =  0.31000D+16 * EXP( -38826./T) 
      DG = +0.652711000E+00*(ALOGT-1.) -0.411700000E-03*T 
     &     +0.939750000E-07*T*T -0.156601667E-10*T*T*T 
     &     +0.116842000E-14*T*T*T*T -0.216360500E+05/T 
     &     +0.207240700E+01
      RB(186) = RF(186)/MAX(EXP(DG),SMALL) 

      RF(187) =  0.14000D+82 * EXP(-18.80*ALOGT -57037./T) 
      DG = +0.536498000E+00*(ALOGT-1.) -0.105824500E-02*T 
     &     +0.254161667E-06*T*T -0.394101667E-10*T*T*T 
     &     +0.268130000E-14*T*T*T*T +0.371321000E+04/T 
     &     +0.659206000E+00
      RB(187) = RF(187)/MAX(EXP(DG),SMALL) 

      RF(188) =  0.31000D+15 * EXP( -18235./T) 
      DG = +0.123392200E+01*(ALOGT-1.) -0.234960500E-02*T 
     &     +0.517151167E-06*T*T -0.749092833E-10*T*T*T 
     &     +0.484538500E-14*T*T*T*T -0.362803000E+04/T 
     &     +0.894732200E+01
      RB(188) = RF(188)/MAX(EXP(DG),SMALL) * (82.06*T) 

      RF(189) =  0.72000D+15 * EXP( -29619./T) 
      DG = -0.479149000E+00*(ALOGT-1.) -0.458195000E-03*T 
     &     +0.185185000E-06*T*T -0.379572500E-10*T*T*T 
     &     +0.310054500E-14*T*T*T*T +0.301344000E+04/T 
     &     +0.133464200E+01
      RB(189) = RF(189)/MAX(EXP(DG),SMALL) 

      RF(190) =  0.40000D+14 
      DG = -0.463470972E+01*(ALOGT-1.) +0.482680981E-02*T 
     &     -0.896991087E-06*T*T +0.101823530E-09*T*T*T 
     &     -0.469413803E-14*T*T*T*T +0.490863443E+05/T 
     &     +0.245583421E+01
      RB(190) = RF(190)/MAX(EXP(DG),SMALL) / (82.06*T) 

      RF(191) =  0.40000D+14 
      DG = -0.493464972E+01*(ALOGT-1.) +0.507403981E-02*T 
     &     -0.960974420E-06*T*T +0.113933363E-09*T*T*T 
     &     -0.572037803E-14*T*T*T*T +0.498290843E+05/T 
     &     +0.401976021E+01
      RB(191) = RF(191)/MAX(EXP(DG),SMALL) / (82.06*T) 

      RF(192) =  0.50000D+14 * EXP( -26157./T) 
      DG = -0.491758400E+01*(ALOGT-1.) +0.400480500E-02*T 
     &     -0.899190000E-06*T*T +0.140134250E-09*T*T*T 
     &     -0.993645500E-14*T*T*T*T -0.137740200E+05/T 
     &     +0.309136260E+02
      RB(192) = RF(192)/MAX(EXP(DG),SMALL) 

      RF(193) =  0.50000D+14 * EXP( -14088./T) 
      DG = +0.725953324E+01*(ALOGT-1.) -0.736594868E-02*T 
     &     +0.166207976E-05*T*T -0.253330094E-09*T*T*T 
     &     +0.174415829E-13*T*T*T*T -0.101309052E+05/T 
     &     -0.249299297E+02
      RB(193) = RF(193)/MAX(EXP(DG),SMALL) * (82.06*T) 

      RF(194) =  0.23000D+14 * EXP( -30279./T) 
      DG = -0.211335299E+01*(ALOGT-1.) +0.257518999E-02*T 
     &     -0.720349831E-06*T*T +0.115796108E-09*T*T*T 
     &     -0.772188748E-14*T*T*T*T -0.299486799E+05/T 
     &     +0.229457521E+02
      RB(194) = RF(194)/MAX(EXP(DG),SMALL) * (82.06*T) 

      RF(195) =  0.80000D+15 * EXP( -47805./T) 
      DG = +0.471341001E+01*(ALOGT-1.) -0.308517851E-02*T 
     &     +0.533402803E-06*T*T -0.665198004E-10*T*T*T 
     &     +0.389150702E-14*T*T*T*T -0.461452082E+05/T 
     &     -0.105834774E+02
      RB(195) = RF(195)/MAX(EXP(DG),SMALL) * (82.06*T) 

      RF(196) =  0.23000D+14 
      DG = -0.282451920E+01*(ALOGT-1.) +0.329682400E-02*T 
     &     -0.681299033E-06*T*T +0.967821333E-10*T*T*T 
     &     -0.623510800E-14*T*T*T*T +0.375796520E+05/T 
     &     -0.634762500E+01
      RB(196) = RF(196)/MAX(EXP(DG),SMALL) / (82.06*T) 

      RF(197) =  0.13000D+06 * EXP(  2.50*ALOGT   -505./T) 
      DG = +0.850623350E+00*(ALOGT-1.) -0.227480648E-02*T 
     &     +0.578765497E-06*T*T -0.897434317E-10*T*T*T 
     &     +0.607345694E-14*T*T*T*T +0.318693640E+04/T 
     &     +0.283347849E+01
      RB(197) = RF(197)/MAX(EXP(DG),SMALL) 

      RF(198) =  0.50000D+12 * EXP(  -4424./T) 
      DG = +0.279809499E+01*(ALOGT-1.) -0.393142431E-02*T 
     &     +0.124593918E-05*T*T -0.228553144E-09*T*T*T 
     &     +0.178012078E-13*T*T*T*T +0.112309959E+05/T 
     &     -0.131102166E+02
      RB(198) = RF(198)/MAX(EXP(DG),SMALL) 

      RF(199) =  0.50000D+14 * EXP(  -4027./T) 
      DG = +0.294605319E+01*(ALOGT-1.) -0.319184123E-02*T 
     &     +0.558142627E-06*T*T -0.528111158E-10*T*T*T 
     &     +0.156109186E-14*T*T*T*T -0.310891178E+03/T 
     &     -0.122232644E+02
      RB(199) = RF(199)/MAX(EXP(DG),SMALL) 

      RF(200) =  0.22000D+15 * EXP( -27659./T) 
      DG = +0.602941050E+01*(ALOGT-1.) -0.708981461E-02*T 
     &     +0.160437181E-05*T*T -0.228258984E-09*T*T*T 
     &     +0.144405659E-13*T*T*T*T -0.328860680E+05/T 
     &     -0.146127449E+02
      RB(200) = RF(200)/MAX(EXP(DG),SMALL) * (82.06*T) 

      RF(201) =  0.20000D+13 * EXP(  -2512./T) 
      DG = -0.688654323E+01*(ALOGT-1.) +0.629298816E-02*T 
     &     -0.130675910E-05*T*T +0.184319892E-09*T*T*T 
     &     -0.116081845E-13*T*T*T*T +0.831779889E+04/T 
     &     +0.365247582E+02
      RB(201) = RF(201)/MAX(EXP(DG),SMALL) 

      RF(202) =  0.10000D+15 * EXP( -32787./T) 

      RF(203) =  0.11000D+14 * EXP( -18235./T) 

      RF(204) =  0.24500D+17 * EXP( -36739./T) 

      RF(205) =  0.10600D+48 * EXP( -9.30*ALOGT -52618./T) 
      DG = +0.773240801E+01*(ALOGT-1.) -0.468974501E-02*T 
     &     +0.743530003E-06*T*T -0.821601921E-10*T*T*T 
     &     +0.421205752E-14*T*T*T*T -0.406325575E+05/T 
     &     -0.307623529E+02
      RB(205) = RF(205)/MAX(EXP(DG),SMALL) * (82.06*T) 

      RF(206) =  0.77000D+06 * EXP(  2.21*ALOGT   -313./T) 
      DG = +0.517351282E+01*(ALOGT-1.) -0.387550084E-02*T 
     &     +0.695100376E-06*T*T -0.829181776E-10*T*T*T 
     &     +0.446600058E-14*T*T*T*T +0.187040565E+05/T 
     &     -0.298255960E+02
      RB(206) = RF(206)/MAX(EXP(DG),SMALL) 

      RF(207) =  0.72300D+13 * EXP(   -655./T) 
      DG = -0.641780130E+00*(ALOGT-1.) -0.474872535E-04*T 
     &     -0.121026108E-06*T*T +0.276272399E-10*T*T*T 
     &     -0.191550004E-14*T*T*T*T +0.454913186E+04/T 
     &     +0.744589400E+01
      RB(207) = RF(207)/MAX(EXP(DG),SMALL) 

      RF(208) =  0.17300D+06 * EXP(  2.50*ALOGT  -1254./T) 
      DG = +0.606968719E+01*(ALOGT-1.) -0.471444623E-02*T 
     &     +0.826772794E-06*T*T -0.971240574E-10*T*T*T 
     &     +0.521333436E-14*T*T*T*T +0.112649212E+05/T 
     &     -0.330740104E+02
      RB(208) = RF(208)/MAX(EXP(DG),SMALL) 

      RF(209) =  0.93300D+13 * EXP( -14701./T) 
      DG = -0.467896558E+01*(ALOGT-1.) +0.365417268E-02*T 
     &     -0.779287197E-06*T*T +0.105821887E-09*T*T*T 
     &     -0.626231742E-14*T*T*T*T -0.139364344E+05/T 
     &     +0.426861018E+02
      RB(209) = RF(209)/MAX(EXP(DG),SMALL) * (82.06*T) 



C-----------------------------------------------------------------------
C --- CALCULATE THE STEADY STATE AND THE PART. EQUIL. CONDITIONS
C-----------------------------------------------------------------------


      CNO2       = SMALL 
      CO         = SMALL 
      COH        = SMALL 
      CHO2       = SMALL 
      CCH3CN     = SMALL 
      CH2O2      = SMALL 
      CHCO       = SMALL 
      CCH3       = SMALL 
      CCH2       = SMALL 
      CCH2S      = SMALL 
      CCH3O      = SMALL 
      CCH2OH     = SMALL 
      CC2H5      = SMALL 
      CC2H3      = SMALL 
      CCH2HCO    = SMALL 
      CHCCO      = SMALL 
      CHNO       = SMALL 
      CHONO      = SMALL 
      CH2NO      = SMALL 
      CNH2       = SMALL 
      CNH        = SMALL 
      CN         = SMALL 
      CN2H2      = SMALL 
      CNNH       = SMALL 
      CCN        = SMALL 
      CNCO       = SMALL 
      CHOCN      = SMALL 
      CHCNO      = SMALL 
      CCH2CN     = SMALL 
      CH2CN      = SMALL 
      CCH3NH2    = SMALL 
      CCH2NH2    = SMALL 
      CCH3NH     = SMALL 
      CH2CNH     = SMALL 
      CHCNH      = SMALL 
      CC3H3      = SMALL 
      CC3H4P     = SMALL 
      CC3H5      = SMALL 
      CHCCHCN    = SMALL 
      CH2CCHCN   = SMALL 
      CAC3H4CN   = SMALL 
      CCC3H4CN   = SMALL 
      CPYRLNE    = SMALL 
      CALLYLCN   = SMALL 
      CHNCPROP   = SMALL 
      CCIC3H5CN  = SMALL 
      CTRC3H5CN  = SMALL 
      CNCCHCH2CN = SMALL 
      CSUCC      = SMALL 
      CHNC6H10O  = SMALL 
      CC4H9CHO   = SMALL 
      CSC4H9     = SMALL 
      CC3H6      = SMALL 

      DO 100 I=1,NIT   !BEGINN DER ITERATION

      DIFFMAX = 0.


C --- STEADY STATE ---

      COLD = CNO2      
      ABV=+RF( 69)*CM( 69)*CNO       *CO        +RF( 71)*CNO *CHO2      
     &    +RB( 72)*CNO       *COH       +RB( 73)*CO2       *CNO       
     &    +RB( 84)*CHNO      *CHONO     +RF( 85)*CH        *CHONO     
     &    +RF( 86)*COH       *CHONO     +RB( 99)*CH2O      *CN2O      
     &    +RB(100)*CNO       *CH2NO     +RB(159)*CNO       *CCH3O     
     &    
      DEN=+RB( 69)*CM( 69)+RB( 71)*COH       
     &    +RF( 72)*CH        +RF( 73)*CO        
     &    +RF( 84)*CH2NO     +RB( 85)*CH2       
     &    +RB( 86)*CH2O      +RF( 99)*CNH2      
     &    +RF(100)*CNH2      +RF(159)*CCH3      
     &    
      CNO2      = ABV/MAX(DEN,SMALL)
      DIFF = DABS((COLD - CNO2      )/MAX(SMALL,CNO2      ,COLD))
      DIFFMAX = MAX(DIFF,DIFFMAX)

      COLD = CO        
      ABV=+RB(  1)*CO2       *CH        +RB(  2)*CH        *COH       
     &    +RF(  4)*COH       *COH       +RF( 11)*CH        *CHO2      
     &    +RB( 12)*CO2       *COH       +RB( 20)*COH       *CHCO      
     &    +RB( 28)*COH       *CCH3      +RB( 31)*CH        *CCH2O     
     &    +RF( 36)*CO2       *CCH3      +RB( 47)*COH       *CC2H5     
     &    +RB( 54)*CH        *CCH2HCO   +RB( 55)*CHCO      *CCH3      
     &    +RF( 59)*CO2       *CC2H3     +RB( 60)*CCO       *CCH2      
     &    +RB( 61)*CH        *CHCCO     +RB( 69)*CM( 69)*CNO2      
     &    +RB( 73)*CO2       *CNO       +RB( 75)*CNO       *COH       
     &    +RB( 81)*COH       *CHNO      +RF( 82)*CO2       *CNH2      
     &    +RB( 89)*COH       *CNH2      +RB( 92)*CH        *CHNO      
     &    +RB( 93)*COH       *CNH       +RB(105)*CNO       *CH        
     &    +RF(108)*CO2       *CNH       +RF(113)*CO2       *CN        
     &    +RF(114)*CNO       *CN        +RF(123)*CM(123)*CN2O      
     &    +RB(126)*CH        *CNCO      +RB(127)*CCO       *CNH       
     &    +RF(132)*CO2       *CCN       +RB(134)*CCO2      *CNH       
     &    +RB(140)*CNO       *CHCO      +RB(143)*CCO       *CNO       
     &    +RB(148)*CCH2O     *CCN       +RB(177)*CCH2O     *CNH       
     &    
      DEN=+RF(  1)*COH       +RF(  2)*CH2       
     &    +RB(  4)*CH2O      +RB( 11)*CH2O      
     &    +RF( 12)*CHO2      +RF( 20)*CCH2O     
     &    +RF( 28)*CCH4      +RF( 31)*CCH3      
     &    +RB( 36)*CCH3O     +RF( 47)*CC2H6     
     &    +RF( 54)*CC2H4     +RF( 55)*CC2H4     
     &    +RB( 59)*CCH2HCO   +RF( 60)*CC2H2     
     &    +RF( 61)*CC2H2     +RF( 69)*CM( 69)*CNO       
     &    +RF( 73)*CNO2      +RF( 75)*CHNO      
     &    +RF( 81)*CH2NO     +RB( 82)*CH2NO     
     &    +RF( 89)*CNH3      +RF( 92)*CNH2      
     &    +RF( 93)*CNH2      +RF(105)*CNH       
     &    +RB(108)*CHNO      +RB(113)*CNO       
     &    +RB(114)*CN2       +RB(123)*CM(123)*CN2       
     &    +RF(126)*CHCN      +RF(127)*CHCN      
     &    +RB(132)*CNCO      +RF(134)*CHNCO     
     &    +RF(140)*CHCNO     +RF(143)*CNCO      
     &    +RF(148)*CCH2CN    +RF(177)*CH2CNH    
     &    
      CO        = ABV/MAX(DEN,SMALL)
      DIFF = DABS((COLD - CO        )/MAX(SMALL,CO        ,COLD))
      DIFFMAX = MAX(DIFF,DIFFMAX)

      COLD = COH       
      ABV=+RB(  1)*CO2       *CH        +RF(  2)*CH2       *CO        
     &    +RB(  3)*CH2O      *CH        +RB(  4)*2.0*CH2O      *CO        
     &    +RB(  5)*CM(  5)*CH2O      +RF( 10)*2.0*CH        *CHO2      
     &    +RF( 12)*CO        *CHO2      +RB( 13)*CO2       *CH2O      
     &    +RF( 15)*2.0*CM( 15)*CH2O2     +RB( 16)*CH2O      *CHO2      
     &    +RB( 17)*CCO2      *CH        +RF( 18)*CCO       *CHO2      
     &    +RF( 20)*CO        *CCH2O     +RB( 21)*CH2O      *CHCO      
     &    +RF( 28)*CCH4      *CO        +RB( 29)*CH2O      *CCH3      
     &    +RF( 32)*CH2O      *CCH2S     +RF( 33)*CH        *CCH2OH    
     &    +RB( 34)*CCH3OH    +RF( 35)*CHO2      *CCH3      
     &    +RF( 37)*CO2       *CCH3      +RF( 42)*CO2       *CCH2S     
     &    +RF( 47)*CC2H6     *CO        +RB( 48)*CH2O      *CC2H5     
     &    +RB( 56)*CH2O      *CC2H3     +RF( 64)*CO2       *CCH2HCO   
     &    +RF( 66)*CO2       *CHCCO     +RB( 70)*CM( 70)*CHONO     
     &    +RF( 71)*CNO       *CHO2      +RF( 72)*CNO2      *CH        
     &    +RF( 75)*CO        *CHNO      +RB( 76)*CH2O      *CNO       
     &    +RF( 80)*CH        *CH2NO     +RF( 81)*CO        *CH2NO     
     &    +RB( 83)*CH2O      *CHNO      +RB( 86)*CH2O      *CNO2      
     &    +RF( 89)*CNH3      *CO        +RB( 90)*CH2O      *CNH2      
     &    +RF( 93)*CO        *CNH2      +RB( 94)*CH2O      *CNH       
     &    +RF( 95)*CHO2      *CNH2      +RF( 97)*CNO       *CNH2      
     &    +RB(106)*CH        *CHNO      +RB(107)*CH2O      *CN        
     &    +RF(109)*CO2       *CNH       +RF(111)*CNO       *CNH       
     &    +RB(112)*CNO       *CH        +RB(118)*CH2O      *CNNH      
     &    +RF(124)*CN2O      *CH        +RB(128)*CH2O      *CCN       
     &    +RB(129)*CH        *CHOCN     +RB(130)*CHNCO     *CH        
     &    +RB(131)*CH        *CNCO      +RB(135)*CH2O      *CNCO      
     &    +RB(136)*CCO2      *CNH2      +RB(139)*CH2O      *CNCO      
     &    +RB(141)*CNO       *CCH2O     +RB(149)*CH2O      *CCH2CN    
     &    +RF(158)*CNO       *CCH3      +RF(161)*CCH3      *CHOCN     
     &    +RF(162)*CNO       *CCH2      +RB(173)*CH2O      *CH2CN     
     &    +RB(176)*CH2O      *CHCNH     +RB(182)*CH2O      *CCH2NH2   
     &    +RB(184)*CH2O      *CCH3NH    +RB(206)*CH2O      *CC3H5     
     &    
      DEN=+RF(  1)*CO        +RB(  2)*CH        
     &    +RF(  3)*CH2       +RF(  4)*2.0*COH       
     &    +RF(  5)*CM(  5)*CH        +RB( 10)*2.0*COH       
     &    +RB( 12)*CO2       +RF( 13)*CHO2      
     &    +RB( 15)*2.0*CM( 15)*COH       +RF( 16)*CH2O2     
     &    +RF( 17)*CCO       +RB( 18)*CCO2      
     &    +RB( 20)*CHCO      +RF( 21)*CCH2O     
     &    +RB( 28)*CCH3      +RF( 29)*CCH4      
     &    +RB( 32)*CCH3      +RB( 33)*CCH3      
     &    +RF( 34)*CCH3      +RB( 35)*CCH3O     
     &    +RB( 37)*CCH2O     +RB( 42)*CCO       *CH        
     &    +RB( 47)*CC2H5     +RF( 48)*CC2H6     
     &    +RF( 56)*CC2H4     +RB( 64)*CCO       *CCH2O     
     &    +RB( 66)*CCO       *CCO       +RF( 70)*CM( 70)*CNO       
     &    +RB( 71)*CNO2      +RB( 72)*CNO       
     &    +RB( 75)*CNO       +RF( 76)*CHNO      
     &    +RB( 80)*CNH2      +RB( 81)*CHNO      
     &    +RF( 83)*CH2NO     +RF( 86)*CHONO     
     &    +RB( 89)*CNH2      +RF( 90)*CNH3      
     &    +RB( 93)*CNH       +RF( 94)*CNH2      
     &    +RB( 95)*CH2NO     +RB( 97)*CNNH      
     &    +RF(106)*CNH       +RF(107)*CNH       
     &    +RB(109)*CNO       +RB(111)*CN2       
     &    +RF(112)*CN        +RF(118)*CN2H2     
     &    +RB(124)*CN2       +RF(128)*CHCN      
     &    +RF(129)*CHCN      +RF(130)*CHCN      
     &    +RF(131)*CCN       +RF(135)*CHNCO     
     &    +RF(136)*CHNCO     +RF(139)*CHOCN     
     &    +RF(141)*CHCNO     +RF(149)*CCH3CN    
     &    +RB(158)*CH2CN     +RB(161)*CCH3CN    
     &    +RB(162)*CHCN      +RF(173)*CH2CNH    
     &    +RF(176)*CH2CNH    +RF(182)*CCH3NH2   
     &    +RF(184)*CCH3NH2   +RF(206)*CC3H6     
     &    
      COH       = ABV/MAX(DEN,SMALL)
      DIFF = DABS((COLD - COH       )/MAX(SMALL,COH       ,COLD))
      DIFFMAX = MAX(DIFF,DIFFMAX)

      COLD = CHO2      
      ABV=+RF(  6)*CO2       *CH        *CN2    +RF(  7)*CO2 *CH2O*CH        
     &    +RF(  8)*CO2       *CCO2      *CH        +RB(  9)*CO2 *CH2       
     &    +RB( 10)*COH       *COH       +RB( 11)*CH2O      *CO        
     &    +RB( 12)*CO2       *COH       +RB( 13)*CO2       *CH2O      
     &    +RB( 14)*2.0*CO2       *CH2O2     +RF( 16)*COH       *CH2O2     
     &    +RB( 18)*CCO2      *COH       +RF( 22)*CO2       *CCH2O     
     &    +RF( 25)*CO2       *CHCO      +RF( 30)*CO2       *CCH4      
     &    +RB( 35)*COH       *CCH3O     +RF( 45)*CO2       *CCH2OH    
     &    +RF( 52)*CO2       *CC2H5     +RB( 71)*CNO2      *COH       
     &    +RF( 77)*CO2       *CHNO      +RB( 95)*COH       *CH2NO     
     &    +RB( 96)*CO2       *CNH3      +RF(121)*CO2       *CNNH      
     &    +RF(180)*CO2       *CCH2NH2   
      DEN=+RB(  6)*CN2       +RB(  7)*CH2O      
     &    +RB(  8)*CCO2      +RF(  9)*CH        
     &    +RF( 10)*CH        +RF( 11)*CH        
     &    +RF( 12)*CO        +RF( 13)*COH       
     &    +RF( 14)*2.0*CHO2      +RB( 16)*CH2O      
     &    +RF( 18)*CCO       +RB( 22)*CHCO      
     &    +RB( 25)*CCO       +RB( 30)*CCH3      
     &    +RF( 35)*CCH3      +RB( 45)*CCH2O     
     &    +RB( 52)*CC2H4     +RF( 71)*CNO       
     &    +RB( 77)*CNO       +RF( 95)*CNH2      
     &    +RF( 96)*CNH2      +RB(121)*CN2       
     &    +RB(180)*CH2CNH    
      CHO2      = ABV/MAX(DEN,SMALL)
      DIFF = DABS((COLD - CHO2      )/MAX(SMALL,CHO2      ,COLD))
      DIFFMAX = MAX(DIFF,DIFFMAX)

      COLD = CCH3CN    
      ABV=+RB(147)*CHCN      *CCH3      +RB(149)*CH2O      *CCH2CN    
     &    +RF(161)*CCH3      *CHOCN     +RB(195)*CH        *CCH2CN    
     &    
      DEN=+RF(147)*CH        +RF(149)*COH       
     &    +RB(161)*COH       +RF(195)
     &    
      CCH3CN    = ABV/MAX(DEN,SMALL)
      DIFF = DABS((COLD - CCH3CN    )/MAX(SMALL,CCH3CN    ,COLD))
      DIFFMAX = MAX(DIFF,DIFFMAX)

      COLD = CH2O2     
      ABV=+RF( 14)*CHO2      *CHO2      +RB( 15)*CM( 15)*COH    *COH       
     &    +RB( 16)*CH2O      *CHO2      
      DEN=+RB( 14)*CO2       +RF( 15)*CM( 15)
     &    +RF( 16)*COH       
      CH2O2     = ABV/MAX(DEN,SMALL)
      DIFF = DABS((COLD - CH2O2     )/MAX(SMALL,CH2O2     ,COLD))
      DIFFMAX = MAX(DIFF,DIFFMAX)

      COLD = CHCO      
      ABV=+RF( 19)*CH        *CCH2O     +RF( 20)*CO        *CCH2O     
     &    +RF( 21)*COH       *CCH2O     +RF( 22)*CO2       *CCH2O     
     &    +RB( 23)*CM( 23)*CCO       *CH        +RB( 24)*CCO    *CH2       
     &    +RB( 25)*CCO       *CHO2      +RF( 39)*CCH2O     *CCH3      
     &    +RF( 55)*CC2H4     *CO        +RF( 58)*CO2       *CC2H3     
     &    +RF( 62)*2.0*CO2       *CC2H2     +RF(140)*CO        *CHCNO     
     &    +RF(153)*CCH2O     *CNCO      +RB(154)*CCO       *CHNO      
     &    +RF(204)*CC4H9CHO  
      DEN=+RB( 19)*CH2       +RB( 20)*COH       
     &    +RB( 21)*CH2O      +RB( 22)*CHO2      
     &    +RF( 23)*CM( 23)+RF( 24)*CH        
     &    +RF( 25)*CO2       +RB( 39)*CCH4      
     &    +RB( 55)*CCH3      +RB( 58)*CCH2O     
     &    +RB( 62)*2.0*CHCO      +RB(140)*CNO       
     &    +RB(153)*CHNCO     +RF(154)*CNO       
     &    
      CHCO      = ABV/MAX(DEN,SMALL)
      DIFF = DABS((COLD - CHCO      )/MAX(SMALL,CHCO      ,COLD))
      DIFFMAX = MAX(DIFF,DIFFMAX)

      COLD = CCH3      
      ABV=+RB( 26)*CCH4      +RF( 27)*CCH4      *CH        
     &    +RF( 28)*CCH4      *CO        +RF( 29)*CCH4      *COH       
     &    +RF( 30)*CO2       *CCH4      +RB( 31)*CH        *CCH2O     
     &    +RF( 32)*CH2O      *CCH2S     +RF( 33)*CH        *CCH2OH    
     &    +RB( 34)*CCH3OH    +RB( 35)*COH       *CCH3O     
     &    +RB( 36)*CO        *CCH3O     +RB( 37)*COH       *CCH2O     
     &    +RB( 38)*2.0*CC2H6     +RB( 39)*CCH4      *CHCO      
     &    +RB( 49)*CCH4      *CC2H5     +RF( 51)*2.0*CH        *CC2H5     
     &    +RF( 55)*CC2H4     *CO        +RF( 63)*CCH2HCO   
     &    +RF(147)*CH        *CCH3CN    +RF(155)*CCH4      *CCN       
     &    +RF(156)*CCH4      *CNCO      +RB(157)*CH2O      *CHCN      
     &    +RB(158)*COH       *CH2CN     +RB(159)*CNO       *CCH3O     
     &    +RB(160)*CH        *CH2CN     +RB(161)*COH       *CCH3CN    
     &    +RB(169)*CH        *CH2CNH    +RB(170)*CCH3NH2   
     &    +RB(174)*CCH4      *CH2CN     +RB(190)*CTRC3H5CN 
     &    +RB(191)*CCIC3H5CN +RF(197)*CH        *CC3H4P    
     &    +RB(198)*CCH4      *CC3H3     +RB(201)*CCH4      *CNCCHCH2CN
     &    +RF(207)*CH        *CC3H6     +RF(209)*CSC4H9    
     &    
      DEN=+RF( 26)*CH        +RB( 27)*CH2       
     &    +RB( 28)*COH       +RB( 29)*CH2O      
     &    +RB( 30)*CHO2      +RF( 31)*CO        
     &    +RB( 32)*COH       +RB( 33)*COH       
     &    +RF( 34)*COH       +RF( 35)*CHO2      
     &    +RF( 36)*CO2       +RF( 37)*CO2       
     &    +RF( 38)*2.0*CCH3      +RF( 39)*CCH2O     
     &    +RF( 49)*CC2H6     +RB( 51)*2.0*CCH3      
     &    +RB( 55)*CHCO      +RB( 63)*CCO       
     &    +RB(147)*CHCN      +RB(155)*CHCN      
     &    +RB(156)*CHNCO     +RF(157)*CNO       
     &    +RF(158)*CNO       +RF(159)*CNO2      
     &    +RF(160)*CN        +RF(161)*CHOCN     
     &    +RF(169)*CNH       +RF(170)*CNH2      
     &    +RF(174)*CH2CNH    +RF(190)*CHCCHCN   
     &    +RF(191)*CHCCHCN   +RB(197)*CC2H2     
     &    +RF(198)*CC3H4P    +RF(201)*CSUCC     
     &    +RB(207)*CC2H4     +RB(209)*CC3H6     
     &    
      CCH3      = ABV/MAX(DEN,SMALL)
      DIFF = DABS((COLD - CCH3      )/MAX(SMALL,CCH3      ,COLD))
      DIFFMAX = MAX(DIFF,DIFFMAX)

      COLD = CCH2      
      ABV=+RB( 40)*CCO       *CH2O      +RF( 41)*CCH2S     *CN2       
     &    +RF( 60)*CC2H2     *CO        +RB(162)*CHCN      *COH       
     &    +RB(163)*CH        *CHCNO     
      DEN=+RF( 40)*CO2       +RB( 41)*CN2       
     &    +RB( 60)*CCO       +RF(162)*CNO       
     &    +RF(163)*CNO       
      CCH2      = ABV/MAX(DEN,SMALL)
      DIFF = DABS((COLD - CCH2      )/MAX(SMALL,CCH2      ,COLD))
      DIFFMAX = MAX(DIFF,DIFFMAX)

      COLD = CCH2S     
      ABV=+RB( 32)*COH       *CCH3      +RB( 41)*CCH2      *CN2       
     &    +RB( 42)*CCO       *CH        *COH       
      DEN=+RF( 32)*CH2O      +RF( 41)*CN2       
     &    +RF( 42)*CO2       
      CCH2S     = ABV/MAX(DEN,SMALL)
      DIFF = DABS((COLD - CCH2S     )/MAX(SMALL,CCH2S     ,COLD))
      DIFFMAX = MAX(DIFF,DIFFMAX)

      COLD = CCH3O     
      ABV=+RF( 35)*CHO2      *CCH3      +RF( 36)*CO2       *CCH3      
     &    +RF( 43)*CH        *CCH2O     +RF(159)*CNO2      *CCH3      
     &    
      DEN=+RB( 35)*COH       +RB( 36)*CO        
     &    +RB( 43)+RB(159)*CNO       
     &    
      CCH3O     = ABV/MAX(DEN,SMALL)
      DIFF = DABS((COLD - CCH3O     )/MAX(SMALL,CCH3O     ,COLD))
      DIFFMAX = MAX(DIFF,DIFFMAX)

      COLD = CCH2OH    
      ABV=+RB( 33)*COH       *CCH3      +RF( 44)*CH        *CCH2O     
     &    +RB( 45)*CHO2      *CCH2O     
      DEN=+RF( 33)*CH        +RB( 44)
     &    +RF( 45)*CO2       
      CCH2OH    = ABV/MAX(DEN,SMALL)
      DIFF = DABS((COLD - CCH2OH    )/MAX(SMALL,CCH2OH    ,COLD))
      DIFFMAX = MAX(DIFF,DIFFMAX)

      COLD = CC2H5     
      ABV=+RF( 46)*CC2H6     *CH        +RF( 47)*CC2H6     *CO        
     &    +RF( 48)*CC2H6     *COH       +RF( 49)*CC2H6     *CCH3      
     &    +RF( 50)*CC2H4     *CH        +RB( 51)*CCH3      *CCH3      
     &    +RB( 52)*CC2H4     *CHO2      +RF(164)*CC2H6     *CNCO      
     &    
      DEN=+RB( 46)*CH2       +RB( 47)*COH       
     &    +RB( 48)*CH2O      +RB( 49)*CCH4      
     &    +RB( 50)+RF( 51)*CH        
     &    +RF( 52)*CO2       +RB(164)*CHNCO     
     &    
      CC2H5     = ABV/MAX(DEN,SMALL)
      DIFF = DABS((COLD - CC2H5     )/MAX(SMALL,CC2H5     ,COLD))
      DIFFMAX = MAX(DIFF,DIFFMAX)

      COLD = CC2H3     
      ABV=+RF( 53)*CC2H4     *CH        +RF( 56)*CC2H4     *COH       
     &    +RF( 57)*CC2H2     *CH        +RB( 58)*CCH2O     *CHCO      
     &    +RB( 59)*CO        *CCH2HCO   +RF(165)*CC2H4     *CCN       
     &    +RB(166)*CC2H2     *CHNO      
      DEN=+RB( 53)*CH2       +RB( 56)*CH2O      
     &    +RB( 57)+RF( 58)*CO2       
     &    +RF( 59)*CO2       +RB(165)*CHCN      
     &    +RF(166)*CNO       
      CC2H3     = ABV/MAX(DEN,SMALL)
      DIFF = DABS((COLD - CC2H3     )/MAX(SMALL,CC2H3     ,COLD))
      DIFFMAX = MAX(DIFF,DIFFMAX)

      COLD = CCH2HCO   
      ABV=+RF( 54)*CC2H4     *CO        +RF( 59)*CO2       *CC2H3     
     &    +RB( 63)*CCO       *CCH3      +RB( 64)*CCO     *COH *CCH2O     
     &    
      DEN=+RB( 54)*CH        +RB( 59)*CO        
     &    +RF( 63)+RF( 64)*CO2       
     &    
      CCH2HCO   = ABV/MAX(DEN,SMALL)
      DIFF = DABS((COLD - CCH2HCO   )/MAX(SMALL,CCH2HCO   ,COLD))
      DIFFMAX = MAX(DIFF,DIFFMAX)

      COLD = CHCCO     
      ABV=+RF( 61)*CC2H2     *CO        +RB( 65)*CCO       *CCO2 *CH        
     &    +RB( 66)*CCO       *CCO       *COH    +RB(167)*CCO  *CHCNO     
     &    +RB(168)*CCO2      *CHCN      
      DEN=+RB( 61)*CH        +RF( 65)*CO2       
     &    +RF( 66)*CO2       +RF(167)*CNO       
     &    +RF(168)*CNO       
      CHCCO     = ABV/MAX(DEN,SMALL)
      DIFF = DABS((COLD - CHCCO     )/MAX(SMALL,CHCCO     ,COLD))
      DIFFMAX = MAX(DIFF,DIFFMAX)

      COLD = CHNO      
      ABV=+RF( 67)*CM( 67)*CNO       *CH        +RF( 68)*CNO *CH *CN2       
     &    +RB( 74)*CH2       *CNO       +RB( 75)*CNO       *COH       
     &    +RB( 76)*CH2O      *CNO       +RB( 77)*CNO       *CHO2      
     &    +RF( 78)*CM( 78)*CH2NO     +RF( 79)*CH        *CH2NO     
     &    +RF( 81)*CO        *CH2NO     +RF( 83)*COH       *CH2NO     
     &    +RF( 84)*CNO2      *CH2NO     +RF( 92)*CO        *CNH2      
     &    +RF(101)*CH2NO     *CNH2      +RF(106)*COH       *CNH       
     &    +RF(108)*CO2       *CNH       +RF(154)*CNO       *CHCO      
     &    +RF(166)*CNO       *CC2H3     
      DEN=+RB( 67)*CM( 67)+RB( 68)*CN2       
     &    +RF( 74)*CH        +RF( 75)*CO        
     &    +RF( 76)*COH       +RF( 77)*CO2       
     &    +RB( 78)*CM( 78)*CH        +RB( 79)*CH2       
     &    +RB( 81)*COH       +RB( 83)*CH2O      
     &    +RB( 84)*CHONO     +RB( 92)*CH        
     &    +RB(101)*CNH3      +RB(106)*CH        
     &    +RB(108)*CO        +RB(154)*CCO       
     &    +RB(166)*CC2H2     
      CHNO      = ABV/MAX(DEN,SMALL)
      DIFF = DABS((COLD - CHNO      )/MAX(SMALL,CHNO      ,COLD))
      DIFFMAX = MAX(DIFF,DIFFMAX)

      COLD = CHONO     
      ABV=+RF( 70)*CM( 70)*CNO       *COH       +RF( 84)*CNO2*CH2NO     
     &    +RB( 85)*CH2       *CNO2      +RB( 86)*CH2O      *CNO2      
     &    
      DEN=+RB( 70)*CM( 70)+RB( 84)*CHNO      
     &    +RF( 85)*CH        +RF( 86)*COH       
     &    
      CHONO     = ABV/MAX(DEN,SMALL)
      DIFF = DABS((COLD - CHONO     )/MAX(SMALL,CHONO     ,COLD))
      DIFFMAX = MAX(DIFF,DIFFMAX)

      COLD = CH2NO     
      ABV=+RB( 78)*CM( 78)*CH        *CHNO      +RB( 79)*CH2 *CHNO      
     &    +RB( 80)*COH       *CNH2      +RB( 81)*COH       *CHNO      
     &    +RF( 82)*CO2       *CNH2      +RB( 83)*CH2O      *CHNO      
     &    +RB( 84)*CHNO      *CHONO     +RF( 95)*CHO2      *CNH2      
     &    +RF(100)*CNO2      *CNH2      +RB(101)*CNH3      *CHNO      
     &    
      DEN=+RF( 78)*CM( 78)+RF( 79)*CH        
     &    +RF( 80)*CH        +RF( 81)*CO        
     &    +RB( 82)*CO        +RF( 83)*COH       
     &    +RF( 84)*CNO2      +RB( 95)*COH       
     &    +RB(100)*CNO       +RF(101)*CNH2      
     &    
      CH2NO     = ABV/MAX(DEN,SMALL)
      DIFF = DABS((COLD - CH2NO     )/MAX(SMALL,CH2NO     ,COLD))
      DIFFMAX = MAX(DIFF,DIFFMAX)

      COLD = CNH2      
      ABV=+RF( 80)*CH        *CH2NO     +RB( 82)*CO        *CH2NO     
     &    +RF( 87)*CM( 87)*CNH3      +RF( 88)*CNH3      *CH        
     &    +RF( 89)*CNH3      *CO        +RF( 90)*CNH3      *COH       
     &    +RB( 91)*CH2       *CNH       +RB( 92)*CH        *CHNO      
     &    +RB( 93)*COH       *CNH       +RB( 94)*CH2O      *CNH       
     &    +RB( 95)*COH       *CH2NO     +RB( 96)*CO2       *CNH3      
     &    +RB( 97)*COH       *CNNH      +RB( 98)*CH2O      *CN2       
     &    +RB( 99)*CH2O      *CN2O      +RB(100)*CNO       *CH2NO     
     &    +RB(101)*CNH3      *CHNO      +RB(102)*2.0*CH2       *CN2H2     
     &    +RB(103)*CH        *CN2H2     +RF(116)*CNO       *CN2H2     
     &    +RB(119)*CNH3      *CNNH      +RF(133)*CHNCO     *CH        
     &    +RF(136)*CHNCO     *COH       +RB(137)*CNH3      *CNCO      
     &    +RB(170)*CCH3NH2   
      DEN=+RB( 80)*COH       +RF( 82)*CO2       
     &    +RB( 87)*CM( 87)*CH        +RB( 88)*CH2       
     &    +RB( 89)*COH       +RB( 90)*CH2O      
     &    +RF( 91)*CH        +RF( 92)*CO        
     &    +RF( 93)*CO        +RF( 94)*COH       
     &    +RF( 95)*CHO2      +RF( 96)*CHO2      
     &    +RF( 97)*CNO       +RF( 98)*CNO       
     &    +RF( 99)*CNO2      +RF(100)*CNO2      
     &    +RF(101)*CH2NO     +RF(102)*2.0*CNH2      
     &    +RF(103)*CNH       +RB(116)*CN2O      
     &    +RF(119)*CN2H2     +RB(133)*CCO       
     &    +RB(136)*CCO2      +RF(137)*CHNCO     
     &    +RF(170)*CCH3      
      CNH2      = ABV/MAX(DEN,SMALL)
      DIFF = DABS((COLD - CNH2      )/MAX(SMALL,CNH2      ,COLD))
      DIFFMAX = MAX(DIFF,DIFFMAX)

      COLD = CNH       
      ABV=+RF( 91)*CH        *CNH2      +RF( 93)*CO        *CNH2      
     &    +RF( 94)*COH       *CNH2      +RB(103)*CH        *CN2H2     
     &    +RB(104)*CH2       *CN        +RB(105)*CNO       *CH        
     &    +RB(106)*CH        *CHNO      +RB(107)*CH2O      *CN        
     &    +RB(108)*CO        *CHNO      +RB(109)*CNO       *COH       
     &    +RB(110)*CN2O      *CH        +RB(111)*COH       *CN2       
     &    +RF(127)*CHCN      *CO        +RF(134)*CHNCO     *CO        
     &    +RF(142)*CH        *CNCO      +RB(169)*CH        *CH2CNH    
     &    +RF(177)*CO        *CH2CNH    
      DEN=+RB( 91)*CH2       +RB( 93)*COH       
     &    +RB( 94)*CH2O      +RF(103)*CNH2      
     &    +RF(104)*CH        +RF(105)*CO        
     &    +RF(106)*COH       +RF(107)*COH       
     &    +RF(108)*CO2       +RF(109)*CO2       
     &    +RF(110)*CNO       +RF(111)*CNO       
     &    +RB(127)*CCO       +RB(134)*CCO2      
     &    +RB(142)*CCO       +RF(169)*CCH3      
     &    +RB(177)*CCH2O     
      CNH       = ABV/MAX(DEN,SMALL)
      DIFF = DABS((COLD - CNH       )/MAX(SMALL,CNH       ,COLD))
      DIFFMAX = MAX(DIFF,DIFFMAX)

      COLD = CN        
      ABV=+RF(104)*CH        *CNH       +RF(107)*COH       *CNH       
     &    +RB(112)*CNO       *CH        +RB(113)*CNO       *CO        
     &    +RB(114)*CO        *CN2       +RB(152)*CCO       *CNO       
     &    +RB(160)*CH        *CH2CN     
      DEN=+RB(104)*CH2       +RB(107)*CH2O      
     &    +RF(112)*COH       +RF(113)*CO2       
     &    +RF(114)*CNO       +RF(152)*CCO2      
     &    +RF(160)*CCH3      
      CN        = ABV/MAX(DEN,SMALL)
      DIFF = DABS((COLD - CN        )/MAX(SMALL,CN        ,COLD))
      DIFFMAX = MAX(DIFF,DIFFMAX)

      COLD = CN2H2     
      ABV=+RF(102)*CNH2      *CNH2      +RF(103)*CNH2      *CNH       
     &    +RB(115)*CM(115)*CH        *CNNH      +RB(116)*CN2O *CNH2      
     &    +RB(117)*CH2       *CNNH      +RB(118)*CH2O      *CNNH      
     &    +RB(119)*CNH3      *CNNH      
      DEN=+RB(102)*CH2       +RB(103)*CH        
     &    +RF(115)*CM(115)+RF(116)*CNO       
     &    +RF(117)*CH        +RF(118)*COH       
     &    +RF(119)*CNH2      
      CN2H2     = ABV/MAX(DEN,SMALL)
      DIFF = DABS((COLD - CN2H2     )/MAX(SMALL,CN2H2     ,COLD))
      DIFFMAX = MAX(DIFF,DIFFMAX)

      COLD = CNNH      
      ABV=+RF( 97)*CNO       *CNH2      +RF(115)*CM(115)*CN2H2     
     &    +RF(117)*CH        *CN2H2     +RF(118)*COH       *CN2H2     
     &    +RF(119)*CNH2      *CN2H2     +RB(120)*CH        *CN2       
     &    +RB(121)*CHO2      *CN2       +RB(122)*CO2       *CH *CN2       
     &    
      DEN=+RB( 97)*COH       +RB(115)*CM(115)*CH        
     &    +RB(117)*CH2       +RB(118)*CH2O      
     &    +RB(119)*CNH3      +RF(120)
     &    +RF(121)*CO2       +RF(122)*CO2       
     &    
      CNNH      = ABV/MAX(DEN,SMALL)
      DIFF = DABS((COLD - CNNH      )/MAX(SMALL,CNNH      ,COLD))
      DIFFMAX = MAX(DIFF,DIFFMAX)

      COLD = CCN       
      ABV=+RB(125)*CHCN      *CH        +RF(128)*CHCN      *COH       
     &    +RB(131)*CH        *CNCO      +RB(132)*CO        *CNCO      
     &    +RF(148)*CO        *CCH2CN    +RB(155)*CHCN      *CCH3      
     &    +RB(165)*CHCN      *CC2H3     +RF(200)*CNCCHCH2CN
     &    
      DEN=+RF(125)*CH2       +RB(128)*CH2O      
     &    +RF(131)*COH       +RF(132)*CO2       
     &    +RB(148)*CCH2O     +RF(155)*CCH4      
     &    +RF(165)*CC2H4     +RB(200)*CH2CCHCN  
     &    
      CCN       = ABV/MAX(DEN,SMALL)
      DIFF = DABS((COLD - CCN       )/MAX(SMALL,CCN       ,COLD))
      DIFFMAX = MAX(DIFF,DIFFMAX)

      COLD = CNCO      
      ABV=+RF(126)*CHCN      *CO        +RF(131)*COH       *CCN       
     &    +RF(132)*CO2       *CCN       +RF(135)*CHNCO     *COH       
     &    +RF(137)*CHNCO     *CNH2      +RF(138)*CH        *CHOCN     
     &    +RF(139)*COH       *CHOCN     +RB(142)*CCO       *CNH       
     &    +RB(143)*CCO       *CNO       +RB(144)*CHNCO     *CH        
     &    +RB(145)*CCO       *CN2O      +RB(146)*CCO2      *CN2       
     &    +RB(153)*CHNCO     *CHCO      +RB(156)*CHNCO     *CCH3      
     &    +RB(164)*CHNCO     *CC2H5     
      DEN=+RB(126)*CH        +RB(131)*CH        
     &    +RB(132)*CO        +RB(135)*CH2O      
     &    +RB(137)*CNH3      +RB(138)*CH2       
     &    +RB(139)*CH2O      +RF(142)*CH        
     &    +RF(143)*CO        +RF(144)*CH2       
     &    +RF(145)*CNO       +RF(146)*CNO       
     &    +RF(153)*CCH2O     +RF(156)*CCH4      
     &    +RF(164)*CC2H6     
      CNCO      = ABV/MAX(DEN,SMALL)
      DIFF = DABS((COLD - CNCO      )/MAX(SMALL,CNCO      ,COLD))
      DIFFMAX = MAX(DIFF,DIFFMAX)

      COLD = CHOCN     
      ABV=+RF(129)*CHCN      *COH       +RB(138)*CH2       *CNCO      
     &    +RB(139)*CH2O      *CNCO      +RB(161)*COH       *CCH3CN    
     &    
      DEN=+RB(129)*CH        +RF(138)*CH        
     &    +RF(139)*COH       +RF(161)*CCH3      
     &    
      CHOCN     = ABV/MAX(DEN,SMALL)
      DIFF = DABS((COLD - CHOCN     )/MAX(SMALL,CHOCN     ,COLD))
      DIFFMAX = MAX(DIFF,DIFFMAX)

      COLD = CHCNO     
      ABV=+RB(140)*CNO       *CHCO      +RB(141)*CNO       *CCH2O     
     &    +RF(163)*CNO       *CCH2      +RF(167)*CNO       *CHCCO     
     &    
      DEN=+RF(140)*CO        +RF(141)*COH       
     &    +RB(163)*CH        +RB(167)*CCO       
     &    
      CHCNO     = ABV/MAX(DEN,SMALL)
      DIFF = DABS((COLD - CHCNO     )/MAX(SMALL,CHCNO     ,COLD))
      DIFFMAX = MAX(DIFF,DIFFMAX)

      COLD = CCH2CN    
      ABV=+RB(148)*CCH2O     *CCN       +RF(149)*COH       *CCH3CN    
     &    +RF(193)*CCC3H4CN  +RF(195)*CCH3CN    
     &    +RB(196)*2.0*CSUCC     
      DEN=+RF(148)*CO        +RB(149)*CH2O      
     &    +RB(193)*CC2H2     +RB(195)*CH        
     &    +RF(196)*2.0*CCH2CN    
      CCH2CN    = ABV/MAX(DEN,SMALL)
      DIFF = DABS((COLD - CCH2CN    )/MAX(SMALL,CCH2CN    ,COLD))
      DIFFMAX = MAX(DIFF,DIFFMAX)

      COLD = CH2CN     
      ABV=+RB(150)*CM(150)*CHCN      *CH        +RF(158)*CNO    *CCH3      
     &    +RF(160)*CCH3      *CN        +RF(172)*CH        *CH2CNH    
     &    +RF(173)*COH       *CH2CNH    +RF(174)*CCH3      *CH2CNH    
     &    
      DEN=+RF(150)*CM(150)+RB(158)*COH       
     &    +RB(160)*CH        +RB(172)*CH2       
     &    +RB(173)*CH2O      +RB(174)*CCH4      
     &    
      CH2CN     = ABV/MAX(DEN,SMALL)
      DIFF = DABS((COLD - CH2CN     )/MAX(SMALL,CH2CN     ,COLD))
      DIFFMAX = MAX(DIFF,DIFFMAX)

      COLD = CCH3NH2   
      ABV=+RF(170)*CCH3      *CNH2      +RB(181)*CH2       *CCH2NH2   
     &    +RB(182)*CH2O      *CCH2NH2   +RB(183)*CH2       *CCH3NH    
     &    +RB(184)*CH2O      *CCH3NH    
      DEN=+RB(170)+RF(181)*CH        
     &    +RF(182)*COH       +RF(183)*CH        
     &    +RF(184)*COH       
      CCH3NH2   = ABV/MAX(DEN,SMALL)
      DIFF = DABS((COLD - CCH3NH2   )/MAX(SMALL,CCH3NH2   ,COLD))
      DIFFMAX = MAX(DIFF,DIFFMAX)

      COLD = CCH2NH2   
      ABV=+RB(179)*CH        *CH2CNH    +RB(180)*CHO2      *CH2CNH    
     &    +RF(181)*CH        *CCH3NH2   +RF(182)*COH       *CCH3NH2   
     &    
      DEN=+RF(179)+RF(180)*CO2       
     &    +RB(181)*CH2       +RB(182)*CH2O      
     &    
      CCH2NH2   = ABV/MAX(DEN,SMALL)
      DIFF = DABS((COLD - CCH2NH2   )/MAX(SMALL,CCH2NH2   ,COLD))
      DIFFMAX = MAX(DIFF,DIFFMAX)

      COLD = CCH3NH    
      ABV=+RB(178)*CH        *CH2CNH    +RF(183)*CH        *CCH3NH2   
     &    +RF(184)*COH       *CCH3NH2   
      DEN=+RF(178)+RB(183)*CH2       
     &    +RB(184)*CH2O      
      CCH3NH    = ABV/MAX(DEN,SMALL)
      DIFF = DABS((COLD - CCH3NH    )/MAX(SMALL,CCH3NH    ,COLD))
      DIFFMAX = MAX(DIFF,DIFFMAX)

      COLD = CH2CNH    
      ABV=+RF(169)*CCH3      *CNH       +RB(172)*CH2       *CH2CN     
     &    +RB(173)*CH2O      *CH2CN     +RB(174)*CCH4      *CH2CN     
     &    +RB(175)*CH2       *CHCNH     +RB(176)*CH2O      *CHCNH     
     &    +RB(177)*CCH2O     *CNH       +RF(178)*CCH3NH    
     &    +RF(179)*CCH2NH2   +RF(180)*CO2       *CCH2NH2   
     &    
      DEN=+RB(169)*CH        +RF(172)*CH        
     &    +RF(173)*COH       +RF(174)*CCH3      
     &    +RF(175)*CH        +RF(176)*COH       
     &    +RF(177)*CO        +RB(178)*CH        
     &    +RB(179)*CH        +RB(180)*CHO2      
     &    
      CH2CNH    = ABV/MAX(DEN,SMALL)
      DIFF = DABS((COLD - CH2CNH    )/MAX(SMALL,CH2CNH    ,COLD))
      DIFFMAX = MAX(DIFF,DIFFMAX)

      COLD = CHCNH     
      ABV=+RB(171)*CHCN      *CH        +RF(175)*CH        *CH2CNH    
     &    +RF(176)*COH       *CH2CNH    
      DEN=+RF(171)+RB(175)*CH2       
     &    +RB(176)*CH2O      
      CHCNH     = ABV/MAX(DEN,SMALL)
      DIFF = DABS((COLD - CHCNH     )/MAX(SMALL,CHCNH     ,COLD))
      DIFFMAX = MAX(DIFF,DIFFMAX)

      COLD = CC3H3     
      ABV=+RF(198)*CCH3      *CC3H4P    
      DEN=+RB(198)*CCH4      
      CC3H3     = ABV/MAX(DEN,SMALL)
      DIFF = DABS((COLD - CC3H3     )/MAX(SMALL,CC3H3     ,COLD))
      DIFFMAX = MAX(DIFF,DIFFMAX)

      COLD = CC3H4P    
      ABV=+RF(188)*CHNCPROP  +RF(194)*CC3H5     
     &    +RB(197)*CC2H2     *CCH3      +RB(198)*CCH4      *CC3H3     
     &    
      DEN=+RB(188)*CHCN      +RB(194)*CH        
     &    +RF(197)*CH        +RF(198)*CCH3      
     &    
      CC3H4P    = ABV/MAX(DEN,SMALL)
      DIFF = DABS((COLD - CC3H4P    )/MAX(SMALL,CC3H4P    ,COLD))
      DIFFMAX = MAX(DIFF,DIFFMAX)

      COLD = CC3H5     
      ABV=+RB(194)*CH        *CC3H4P    +RF(205)*CC3H6     
     &    +RF(206)*COH       *CC3H6     +RF(208)*CH        *CC3H6     
     &    
      DEN=+RF(194)+RB(205)*CH        
     &    +RB(206)*CH2O      +RB(208)*CH2       
     &    
      CC3H5     = ABV/MAX(DEN,SMALL)
      DIFF = DABS((COLD - CC3H5     )/MAX(SMALL,CC3H5     ,COLD))
      DIFFMAX = MAX(DIFF,DIFFMAX)

      COLD = CHCCHCN   
      ABV=+RB(190)*CTRC3H5CN +RB(191)*CCIC3H5CN 
     &    +RF(199)*CH        *CH2CCHCN  
      DEN=+RF(190)*CCH3      +RF(191)*CCH3      
     &    +RB(199)*CH2       
      CHCCHCN   = ABV/MAX(DEN,SMALL)
      DIFF = DABS((COLD - CHCCHCN   )/MAX(SMALL,CHCCHCN   ,COLD))
      DIFFMAX = MAX(DIFF,DIFFMAX)

      COLD = CH2CCHCN  
      ABV=+RB(199)*CH2       *CHCCHCN   +RF(200)*CNCCHCH2CN
     &    
      DEN=+RF(199)*CH        +RB(200)*CCN       
     &    
      CH2CCHCN  = ABV/MAX(DEN,SMALL)
      DIFF = DABS((COLD - CH2CCHCN  )/MAX(SMALL,CH2CCHCN  ,COLD))
      DIFFMAX = MAX(DIFF,DIFFMAX)

      COLD = CAC3H4CN  
      ABV=+RB(192)*CCC3H4CN  
      DEN=+RF(192)
      CAC3H4CN  = ABV/MAX(DEN,SMALL)
      DIFF = DABS((COLD - CAC3H4CN  )/MAX(SMALL,CAC3H4CN  ,COLD))
      DIFFMAX = MAX(DIFF,DIFFMAX)

      COLD = CCC3H4CN  
      ABV=+RF(192)*CAC3H4CN  +RB(193)*CC2H2     *CCH2CN    
     &    
      DEN=+RB(192)+RF(193)
     &    
      CCC3H4CN  = ABV/MAX(DEN,SMALL)
      DIFF = DABS((COLD - CCC3H4CN  )/MAX(SMALL,CCC3H4CN  ,COLD))
      DIFFMAX = MAX(DIFF,DIFFMAX)

      COLD = CPYRLNE   
      ABV=+RF(185)*CM(185)*CC4H5N    +RB(187)*CCIC3H5CN 
     &    
      DEN=+RB(185)*CM(185)+RF(187)
     &    
      CPYRLNE   = ABV/MAX(DEN,SMALL)
      DIFF = DABS((COLD - CPYRLNE   )/MAX(SMALL,CPYRLNE   ,COLD))
      DIFFMAX = MAX(DIFF,DIFFMAX)

      COLD = CALLYLCN  
      ABV=+RB(189)*CCIC3H5CN 
      DEN=+RF(189)
      CALLYLCN  = ABV/MAX(DEN,SMALL)
      DIFF = DABS((COLD - CALLYLCN  )/MAX(SMALL,CALLYLCN  ,COLD))
      DIFFMAX = MAX(DIFF,DIFFMAX)

      COLD = CHNCPROP  
      ABV=+RF(186)*CC4H5N    +RB(188)*CHCN      *CC3H4P    
     &    
      DEN=+RB(186)+RF(188)
     &    
      CHNCPROP  = ABV/MAX(DEN,SMALL)
      DIFF = DABS((COLD - CHNCPROP  )/MAX(SMALL,CHNCPROP  ,COLD))
      DIFFMAX = MAX(DIFF,DIFFMAX)

      COLD = CCIC3H5CN 
      ABV=+RF(187)*CPYRLNE   +RF(189)*CALLYLCN  
     &    +RF(191)*CCH3      *CHCCHCN   
      DEN=+RB(187)+RB(189)
     &    +RB(191)
      CCIC3H5CN = ABV/MAX(DEN,SMALL)
      DIFF = DABS((COLD - CCIC3H5CN )/MAX(SMALL,CCIC3H5CN ,COLD))
      DIFFMAX = MAX(DIFF,DIFFMAX)

      COLD = CTRC3H5CN 
      ABV=+RF(190)*CCH3      *CHCCHCN   
      DEN=+RB(190)
      CTRC3H5CN = ABV/MAX(DEN,SMALL)
      DIFF = DABS((COLD - CTRC3H5CN )/MAX(SMALL,CTRC3H5CN ,COLD))
      DIFFMAX = MAX(DIFF,DIFFMAX)

      COLD = CNCCHCH2CN
      ABV=+RB(200)*CCN       *CH2CCHCN  +RF(201)*CCH3      *CSUCC     
     &    
      DEN=+RF(200)+RB(201)*CCH4      
     &    
      CNCCHCH2CN= ABV/MAX(DEN,SMALL)
      DIFF = DABS((COLD - CNCCHCH2CN)/MAX(SMALL,CNCCHCH2CN,COLD))
      DIFFMAX = MAX(DIFF,DIFFMAX)

      COLD = CSUCC     
      ABV=+RF(196)*CCH2CN    *CCH2CN    +RB(201)*CCH4      *CNCCHCH2CN
     &    
      DEN=+RB(196)+RF(201)*CCH3      
     &    
      CSUCC     = ABV/MAX(DEN,SMALL)
      DIFF = DABS((COLD - CSUCC     )/MAX(SMALL,CSUCC     ,COLD))
      DIFFMAX = MAX(DIFF,DIFFMAX)

      COLD = CHNC6H10O 
      ABV=+RF(202)*CC6H11NO  
      DEN=+RF(203)
      CHNC6H10O = ABV/MAX(DEN,SMALL)
      DIFF = DABS((COLD - CHNC6H10O )/MAX(SMALL,CHNC6H10O ,COLD))
      DIFFMAX = MAX(DIFF,DIFFMAX)

      COLD = CC4H9CHO  
      ABV=+RF(203)*CHNC6H10O 
      DEN=+RF(204)
      CC4H9CHO  = ABV/MAX(DEN,SMALL)
      DIFF = DABS((COLD - CC4H9CHO  )/MAX(SMALL,CC4H9CHO  ,COLD))
      DIFFMAX = MAX(DIFF,DIFFMAX)

      COLD = CSC4H9    
      ABV=+RF(204)*CC4H9CHO  +RB(209)*CCH3      *CC3H6     
     &    
      DEN=+RF(209)
      CSC4H9    = ABV/MAX(DEN,SMALL)
      DIFF = DABS((COLD - CSC4H9    )/MAX(SMALL,CSC4H9    ,COLD))
      DIFFMAX = MAX(DIFF,DIFFMAX)

      COLD = CC3H6     
      ABV=+RB(205)*CH        *CC3H5     +RB(206)*CH2O      *CC3H5     
     &    +RB(207)*CC2H4     *CCH3      +RB(208)*CH2       *CC3H5     
     &    +RF(209)*CSC4H9    
      DEN=+RF(205)+RF(206)*COH       
     &    +RF(207)*CH        +RF(208)*CH        
     &    +RB(209)*CCH3      
      CC3H6     = ABV/MAX(DEN,SMALL)
      DIFF = DABS((COLD - CC3H6     )/MAX(SMALL,CC3H6     ,COLD))
      DIFFMAX = MAX(DIFF,DIFFMAX)

C
      IF (DIFFMAX .LT. GDIFF) GOTO 200
C
  100 CONTINUE   !ENDE DER ITERATION
  200 CONTINUE


C-----------------------------------------------------------------------
C --- CALCULATE THE NET REACTION RATE OF THE SCELETAL MECHANISM
C-----------------------------------------------------------------------

      W(  1)=RF(  1)*CO        *COH       -RB(  1)*CO2       *CH         
      W(  2)=RF(  2)*CH2       *CO        -RB(  2)*CH        *COH        
      W(  3)=RF(  3)*CH2       *COH       -RB(  3)*CH2O      *CH         
      W(  4)=RF(  4)*COH       *COH       -RB(  4)*CH2O      *CO         
      W(  5)=RF(  5)*CH        *COH       *CM(  5)-RB( 5)*CH2O *CM(5) 
      W(  6)=RF(  6)*CO2       *CH        *CN2  -RB(  6)*CHO2  *CN2        
      W(  7)=RF(  7)*CO2       *CH2O      *CH   -RB(  7)*CH2O  *CHO2       
      W(  8)=RF(  8)*CO2       *CCO2      *CH   -RB(  8)*CCO2  *CHO2       
      W(  9)=RF(  9)*CH        *CHO2      -RB(  9)*CO2       *CH2        
      W( 10)=RF( 10)*CH        *CHO2      -RB( 10)*COH       *COH        
      W( 11)=RF( 11)*CH        *CHO2      -RB( 11)*CH2O      *CO         
      W( 12)=RF( 12)*CO        *CHO2      -RB( 12)*CO2       *COH        
      W( 13)=RF( 13)*COH       *CHO2      -RB( 13)*CO2       *CH2O       
      W( 14)=RF( 14)*CHO2      *CHO2      -RB( 14)*CO2       *CH2O2      
      W( 15)=RF( 15)*CH2O2     *CM( 15)-RB( 15)*COH  *COH    *CM( 15) 
      W( 16)=RF( 16)*COH       *CH2O2     -RB( 16)*CH2O      *CHO2       
      W( 17)=RF( 17)*CCO       *COH       -RB( 17)*CCO2      *CH         
      W( 18)=RF( 18)*CCO       *CHO2      -RB( 18)*CCO2      *COH        
      W( 19)=RF( 19)*CH        *CCH2O     -RB( 19)*CH2       *CHCO       
      W( 20)=RF( 20)*CO        *CCH2O     -RB( 20)*COH       *CHCO       
      W( 21)=RF( 21)*COH       *CCH2O     -RB( 21)*CH2O      *CHCO       
      W( 22)=RF( 22)*CO2       *CCH2O     -RB( 22)*CHO2      *CHCO       
      W( 23)=RF( 23)*CHCO      *CM( 23)-RB( 23)*CCO    *CH  *CM( 23) 
      W( 24)=RF( 24)*CH        *CHCO      -RB( 24)*CCO       *CH2        
      W( 25)=RF( 25)*CO2       *CHCO      -RB( 25)*CCO       *CHO2       
      W( 26)=RF( 26)*CH        *CCH3      -RB( 26)*CCH4       
      W( 27)=RF( 27)*CCH4      *CH        -RB( 27)*CH2       *CCH3       
      W( 28)=RF( 28)*CCH4      *CO        -RB( 28)*COH       *CCH3       
      W( 29)=RF( 29)*CCH4      *COH       -RB( 29)*CH2O      *CCH3       
      W( 30)=RF( 30)*CO2       *CCH4      -RB( 30)*CHO2      *CCH3       
      W( 31)=RF( 31)*CO        *CCH3      -RB( 31)*CH        *CCH2O      
      W( 32)=RF( 32)*CH2O      *CCH2S     -RB( 32)*COH       *CCH3       
      W( 33)=RF( 33)*CH        *CCH2OH    -RB( 33)*COH       *CCH3       
      W( 34)=RF( 34)*COH       *CCH3      -RB( 34)*CCH3OH     
      W( 35)=RF( 35)*CHO2      *CCH3      -RB( 35)*COH       *CCH3O      
      W( 36)=RF( 36)*CO2       *CCH3      -RB( 36)*CO        *CCH3O      
      W( 37)=RF( 37)*CO2       *CCH3      -RB( 37)*COH       *CCH2O      
      W( 38)=RF( 38)*CCH3      *CCH3      -RB( 38)*CC2H6      
      W( 39)=RF( 39)*CCH2O     *CCH3      -RB( 39)*CCH4      *CHCO       
      W( 40)=RF( 40)*CO2       *CCH2      -RB( 40)*CCO       *CH2O       
      W( 41)=RF( 41)*CCH2S     *CN2       -RB( 41)*CCH2      *CN2        
      W( 42)=RF( 42)*CO2       *CCH2S     -RB( 42)*CCO   *CH *COH        
      W( 43)=RF( 43)*CH        *CCH2O     -RB( 43)*CCH3O      
      W( 44)=RF( 44)*CH        *CCH2O     -RB( 44)*CCH2OH     
      W( 45)=RF( 45)*CO2       *CCH2OH    -RB( 45)*CHO2      *CCH2O      
      W( 46)=RF( 46)*CC2H6     *CH        -RB( 46)*CH2       *CC2H5      
      W( 47)=RF( 47)*CC2H6     *CO        -RB( 47)*COH       *CC2H5      
      W( 48)=RF( 48)*CC2H6     *COH       -RB( 48)*CH2O      *CC2H5      
      W( 49)=RF( 49)*CC2H6     *CCH3      -RB( 49)*CCH4      *CC2H5      
      W( 50)=RF( 50)*CC2H4     *CH        -RB( 50)*CC2H5      
      W( 51)=RF( 51)*CH        *CC2H5     -RB( 51)*CCH3      *CCH3       
      W( 52)=RF( 52)*CO2       *CC2H5     -RB( 52)*CC2H4     *CHO2       
      W( 53)=RF( 53)*CC2H4     *CH        -RB( 53)*CH2       *CC2H3      
      W( 54)=RF( 54)*CC2H4     *CO        -RB( 54)*CH        *CCH2HCO    
      W( 55)=RF( 55)*CC2H4     *CO        -RB( 55)*CHCO      *CCH3       
      W( 56)=RF( 56)*CC2H4     *COH       -RB( 56)*CH2O      *CC2H3      
      W( 57)=RF( 57)*CC2H2     *CH        -RB( 57)*CC2H3      
      W( 58)=RF( 58)*CO2       *CC2H3     -RB( 58)*CCH2O     *CHCO       
      W( 59)=RF( 59)*CO2       *CC2H3     -RB( 59)*CO        *CCH2HCO    
      W( 60)=RF( 60)*CC2H2     *CO        -RB( 60)*CCO       *CCH2       
      W( 61)=RF( 61)*CC2H2     *CO        -RB( 61)*CH        *CHCCO      
      W( 62)=RF( 62)*CO2       *CC2H2     -RB( 62)*CHCO      *CHCO       
      W( 63)=RF( 63)*CCH2HCO   -RB( 63)*CCO       *CCH3       
      W( 64)=RF( 64)*CO2       *CCH2HCO   -RB( 64)*CCO   *COH   *CCH2O      
      W( 65)=RF( 65)*CO2       *CHCCO     -RB( 65)*CCO   *CCO2  *CH         
      W( 66)=RF( 66)*CO2       *CHCCO     -RB( 66)*CCO   *CCO   *COH        
      W( 67)=RF( 67)*CNO       *CH    *CM( 67)-RB( 67)*CHNO  *CM( 67) 
      W( 68)=RF( 68)*CNO       *CH    *CN2       -RB( 68)*CHNO *CN2        
      W( 69)=RF( 69)*CNO       *CO    *CM( 69)-RB( 69)*CNO2  *CM( 69) 
      W( 70)=RF( 70)*CNO       *COH   *CM( 70)-RB( 70)*CHONO *CM( 70) 
      W( 71)=RF( 71)*CNO       *CHO2      -RB( 71)*CNO2      *COH        
      W( 72)=RF( 72)*CNO2      *CH        -RB( 72)*CNO       *COH        
      W( 73)=RF( 73)*CNO2      *CO        -RB( 73)*CO2       *CNO        
      W( 74)=RF( 74)*CH        *CHNO      -RB( 74)*CH2       *CNO        
      W( 75)=RF( 75)*CO        *CHNO      -RB( 75)*CNO       *COH        
      W( 76)=RF( 76)*COH       *CHNO      -RB( 76)*CH2O      *CNO        
      W( 77)=RF( 77)*CO2       *CHNO      -RB( 77)*CNO       *CHO2       
      W( 78)=RF( 78)*CH2NO     *CM( 78)-RB( 78)*CH  *CHNO  *CM( 78) 
      W( 79)=RF( 79)*CH        *CH2NO     -RB( 79)*CH2       *CHNO       
      W( 80)=RF( 80)*CH        *CH2NO     -RB( 80)*COH       *CNH2       
      W( 81)=RF( 81)*CO        *CH2NO     -RB( 81)*COH       *CHNO       
      W( 82)=RF( 82)*CO2       *CNH2      -RB( 82)*CO        *CH2NO      
      W( 83)=RF( 83)*COH       *CH2NO     -RB( 83)*CH2O      *CHNO       
      W( 84)=RF( 84)*CNO2      *CH2NO     -RB( 84)*CHNO      *CHONO      
      W( 85)=RF( 85)*CH        *CHONO     -RB( 85)*CH2       *CNO2       
      W( 86)=RF( 86)*COH       *CHONO     -RB( 86)*CH2O      *CNO2       
      W( 87)=RF( 87)*CNH3      *CM( 87)-RB( 87)*CH  *CNH2   *CM( 87) 
      W( 88)=RF( 88)*CNH3      *CH        -RB( 88)*CH2       *CNH2       
      W( 89)=RF( 89)*CNH3      *CO        -RB( 89)*COH       *CNH2       
      W( 90)=RF( 90)*CNH3      *COH       -RB( 90)*CH2O      *CNH2       
      W( 91)=RF( 91)*CH        *CNH2      -RB( 91)*CH2       *CNH        
      W( 92)=RF( 92)*CO        *CNH2      -RB( 92)*CH        *CHNO       
      W( 93)=RF( 93)*CO        *CNH2      -RB( 93)*COH       *CNH        
      W( 94)=RF( 94)*COH       *CNH2      -RB( 94)*CH2O      *CNH        
      W( 95)=RF( 95)*CHO2      *CNH2      -RB( 95)*COH       *CH2NO      
      W( 96)=RF( 96)*CHO2      *CNH2      -RB( 96)*CO2       *CNH3       
      W( 97)=RF( 97)*CNO       *CNH2      -RB( 97)*COH       *CNNH       
      W( 98)=RF( 98)*CNO       *CNH2      -RB( 98)*CH2O      *CN2        
      W( 99)=RF( 99)*CNO2      *CNH2      -RB( 99)*CH2O      *CN2O       
      W(100)=RF(100)*CNO2      *CNH2      -RB(100)*CNO       *CH2NO      
      W(101)=RF(101)*CH2NO     *CNH2      -RB(101)*CNH3      *CHNO       
      W(102)=RF(102)*CNH2      *CNH2      -RB(102)*CH2       *CN2H2      
      W(103)=RF(103)*CNH2      *CNH       -RB(103)*CH        *CN2H2      
      W(104)=RF(104)*CH        *CNH       -RB(104)*CH2       *CN         
      W(105)=RF(105)*CO        *CNH       -RB(105)*CNO       *CH         
      W(106)=RF(106)*COH       *CNH       -RB(106)*CH        *CHNO       
      W(107)=RF(107)*COH       *CNH       -RB(107)*CH2O      *CN         
      W(108)=RF(108)*CO2       *CNH       -RB(108)*CO        *CHNO       
      W(109)=RF(109)*CO2       *CNH       -RB(109)*CNO       *COH        
      W(110)=RF(110)*CNO       *CNH       -RB(110)*CN2O      *CH         
      W(111)=RF(111)*CNO       *CNH       -RB(111)*COH       *CN2        
      W(112)=RF(112)*COH       *CN        -RB(112)*CNO       *CH         
      W(113)=RF(113)*CO2       *CN        -RB(113)*CNO       *CO         
      W(114)=RF(114)*CNO       *CN        -RB(114)*CO        *CN2        
      W(115)=RF(115)*CN2H2     *CM(115)-RB(115)*CH  *CNNH  *CM(115) 
      W(116)=RF(116)*CNO       *CN2H2     -RB(116)*CN2O      *CNH2       
      W(117)=RF(117)*CH        *CN2H2     -RB(117)*CH2       *CNNH       
      W(118)=RF(118)*COH       *CN2H2     -RB(118)*CH2O      *CNNH       
      W(119)=RF(119)*CNH2      *CN2H2     -RB(119)*CNH3      *CNNH       
      W(120)=RF(120)*CNNH      -RB(120)*CH        *CN2        
      W(121)=RF(121)*CO2       *CNNH      -RB(121)*CHO2      *CN2        
      W(122)=RF(122)*CO2       *CNNH      -RB(122)*CO2   *CH *CN2        
      W(123)=RF(123)*CN2O      *CM(123)-RB(123)*CO     *CN2 *CM(123) 
      W(124)=RF(124)*CN2O      *CH        -RB(124)*COH       *CN2        
      W(125)=RF(125)*CH2       *CCN       -RB(125)*CHCN      *CH         
      W(126)=RF(126)*CHCN      *CO        -RB(126)*CH        *CNCO       
      W(127)=RF(127)*CHCN      *CO        -RB(127)*CCO       *CNH        
      W(128)=RF(128)*CHCN      *COH       -RB(128)*CH2O      *CCN        
      W(129)=RF(129)*CHCN      *COH       -RB(129)*CH        *CHOCN      
      W(130)=RF(130)*CHCN      *COH       -RB(130)*CHNCO     *CH         
      W(131)=RF(131)*COH       *CCN       -RB(131)*CH        *CNCO       
      W(132)=RF(132)*CO2       *CCN       -RB(132)*CO        *CNCO       
      W(133)=RF(133)*CHNCO     *CH        -RB(133)*CCO       *CNH2       
      W(134)=RF(134)*CHNCO     *CO        -RB(134)*CCO2      *CNH        
      W(135)=RF(135)*CHNCO     *COH       -RB(135)*CH2O      *CNCO       
      W(136)=RF(136)*CHNCO     *COH       -RB(136)*CCO2      *CNH2       
      W(137)=RF(137)*CHNCO     *CNH2      -RB(137)*CNH3      *CNCO       
      W(138)=RF(138)*CH        *CHOCN     -RB(138)*CH2       *CNCO       
      W(139)=RF(139)*COH       *CHOCN     -RB(139)*CH2O      *CNCO       
      W(140)=RF(140)*CO        *CHCNO     -RB(140)*CNO       *CHCO       
      W(141)=RF(141)*COH       *CHCNO     -RB(141)*CNO       *CCH2O      
      W(142)=RF(142)*CH        *CNCO      -RB(142)*CCO       *CNH        
      W(143)=RF(143)*CO        *CNCO      -RB(143)*CCO       *CNO        
      W(144)=RF(144)*CH2       *CNCO      -RB(144)*CHNCO     *CH         
      W(145)=RF(145)*CNO       *CNCO      -RB(145)*CCO       *CN2O       
      W(146)=RF(146)*CNO       *CNCO      -RB(146)*CCO2      *CN2        
      W(147)=RF(147)*CH        *CCH3CN    -RB(147)*CHCN      *CCH3       
      W(148)=RF(148)*CO        *CCH2CN    -RB(148)*CCH2O     *CCN        
      W(149)=RF(149)*COH       *CCH3CN    -RB(149)*CH2O      *CCH2CN     
      W(150)=RF(150)*CH2CN     *CM(150)-RB(150)*CHCN    *CH  *CM(150) 
      W(151)=RF(151)*CCO       *CN2O      -RB(151)*CCO2      *CN2        
      W(152)=RF(152)*CCO2      *CN        -RB(152)*CCO       *CNO        
      W(153)=RF(153)*CCH2O     *CNCO      -RB(153)*CHNCO     *CHCO       
      W(154)=RF(154)*CNO       *CHCO      -RB(154)*CCO       *CHNO       
      W(155)=RF(155)*CCH4      *CCN       -RB(155)*CHCN      *CCH3       
      W(156)=RF(156)*CCH4      *CNCO      -RB(156)*CHNCO     *CCH3       
      W(157)=RF(157)*CNO       *CCH3      -RB(157)*CH2O      *CHCN       
      W(158)=RF(158)*CNO       *CCH3      -RB(158)*COH       *CH2CN      
      W(159)=RF(159)*CNO2      *CCH3      -RB(159)*CNO       *CCH3O      
      W(160)=RF(160)*CCH3      *CN        -RB(160)*CH        *CH2CN      
      W(161)=RF(161)*CCH3      *CHOCN     -RB(161)*COH       *CCH3CN     
      W(162)=RF(162)*CNO       *CCH2      -RB(162)*CHCN      *COH        
      W(163)=RF(163)*CNO       *CCH2      -RB(163)*CH        *CHCNO      
      W(164)=RF(164)*CC2H6     *CNCO      -RB(164)*CHNCO     *CC2H5      
      W(165)=RF(165)*CC2H4     *CCN       -RB(165)*CHCN      *CC2H3      
      W(166)=RF(166)*CNO       *CC2H3     -RB(166)*CC2H2     *CHNO       
      W(167)=RF(167)*CNO       *CHCCO     -RB(167)*CCO       *CHCNO      
      W(168)=RF(168)*CNO       *CHCCO     -RB(168)*CCO2      *CHCN       
      W(169)=RF(169)*CCH3      *CNH       -RB(169)*CH        *CH2CNH     
      W(170)=RF(170)*CCH3      *CNH2      -RB(170)*CCH3NH2    
      W(171)=RF(171)*CHCNH     -RB(171)*CHCN      *CH         
      W(172)=RF(172)*CH        *CH2CNH    -RB(172)*CH2       *CH2CN      
      W(173)=RF(173)*COH       *CH2CNH    -RB(173)*CH2O      *CH2CN      
      W(174)=RF(174)*CCH3      *CH2CNH    -RB(174)*CCH4      *CH2CN      
      W(175)=RF(175)*CH        *CH2CNH    -RB(175)*CH2       *CHCNH      
      W(176)=RF(176)*COH       *CH2CNH    -RB(176)*CH2O      *CHCNH      
      W(177)=RF(177)*CO        *CH2CNH    -RB(177)*CCH2O     *CNH        
      W(178)=RF(178)*CCH3NH    -RB(178)*CH        *CH2CNH     
      W(179)=RF(179)*CCH2NH2   -RB(179)*CH        *CH2CNH     
      W(180)=RF(180)*CO2       *CCH2NH2   -RB(180)*CHO2      *CH2CNH     
      W(181)=RF(181)*CH        *CCH3NH2   -RB(181)*CH2       *CCH2NH2    
      W(182)=RF(182)*COH       *CCH3NH2   -RB(182)*CH2O      *CCH2NH2    
      W(183)=RF(183)*CH        *CCH3NH2   -RB(183)*CH2       *CCH3NH     
      W(184)=RF(184)*COH       *CCH3NH2   -RB(184)*CH2O      *CCH3NH     
      W(185)=RF(185)*CC4H5N    *CM(185)-RB(185)*CPYRLNE   *CM(185) 
      W(186)=RF(186)*CC4H5N    -RB(186)*CHNCPROP   
      W(187)=RF(187)*CPYRLNE   -RB(187)*CCIC3H5CN  
      W(188)=RF(188)*CHNCPROP  -RB(188)*CHCN      *CC3H4P     
      W(189)=RF(189)*CALLYLCN  -RB(189)*CCIC3H5CN  
      W(190)=RF(190)*CCH3      *CHCCHCN   -RB(190)*CTRC3H5CN  
      W(191)=RF(191)*CCH3      *CHCCHCN   -RB(191)*CCIC3H5CN  
      W(192)=RF(192)*CAC3H4CN  -RB(192)*CCC3H4CN   
      W(193)=RF(193)*CCC3H4CN  -RB(193)*CC2H2     *CCH2CN     
      W(194)=RF(194)*CC3H5     -RB(194)*CH        *CC3H4P     
      W(195)=RF(195)*CCH3CN    -RB(195)*CH        *CCH2CN     
      W(196)=RF(196)*CCH2CN    *CCH2CN    -RB(196)*CSUCC      
      W(197)=RF(197)*CH        *CC3H4P    -RB(197)*CC2H2     *CCH3       
      W(198)=RF(198)*CCH3      *CC3H4P    -RB(198)*CCH4      *CC3H3      
      W(199)=RF(199)*CH        *CH2CCHCN  -RB(199)*CH2       *CHCCHCN    
      W(200)=RF(200)*CNCCHCH2CN-RB(200)*CCN       *CH2CCHCN   
      W(201)=RF(201)*CCH3      *CSUCC     -RB(201)*CCH4  *CNCCHCH2CN 
      W(202)=RF(202)*CC6H11NO   
      W(203)=RF(203)*CHNC6H10O  
      W(204)=RF(204)*CC4H9CHO   
      W(205)=RF(205)*CC3H6     -RB(205)*CH        *CC3H5      
      W(206)=RF(206)*COH       *CC3H6     -RB(206)*CH2O      *CC3H5      
      W(207)=RF(207)*CH        *CC3H6     -RB(207)*CC2H4     *CCH3       
      W(208)=RF(208)*CH        *CC3H6     -RB(208)*CH2       *CC3H5      
      W(209)=RF(209)*CSC4H9    -RB(209)*CCH3      *CC3H6      


C-----------------------------------------------------------------------
C --- CALCULATE THE GLOBAL REACTION RATES
C-----------------------------------------------------------------------

      WGLOB( 1) = - W(  2) + W(  4) + W(  5) + W(  9) + W( 11)  
     &           + W( 13) + W( 14) - W( 15) - W( 20) - W( 22)  
     &           + W( 24) + W( 26) - W( 28) - W( 30) + W( 38)  
     &           + W( 41) - W( 47) - W( 54) - W( 55) - W( 57)  
     &           - W( 60) -2.0*W( 62) -2.0*W( 65)  
     &           -2.0*W( 66) + W( 67) + W( 68) - W( 75)  
     &           - W( 77) - W( 78) - W( 81) - W( 84) + W( 85)  
     &           + W( 86) - W( 87) - W( 89) - W( 92) - W( 93)  
     &           + W( 96) + W( 98) + W( 99) - W(105) - W(109)  
     &           + W(110) + W(111) +2.0*W(114) - W(115)  
     &           + W(116) + W(131) + W(132) + W(133) + W(135)  
     &           + W(136) + W(137) + W(138) + W(139) - W(140)  
     &           - W(143) - W(144) + W(145) + W(146) - W(153)  
     &           + W(154) - W(156) + W(157) - W(164) +2.0*W(166)  
     &           - W(167) - W(168) -3.0*W(187) - W(195)  
     &           + W(206) + W(208)  
      WGLOB( 2) = + W(  2) - W(  4) + W( 10) + W( 12) + W( 15)  
     &           + W( 18) + W( 20) - W( 23) - W( 24) - W( 25)  
     &           + W( 28) - W( 36) - W( 38) - W( 41) - W( 42)  
     &           + W( 44) - W( 45) + W( 47) + W( 51) + W( 55)  
     &           - W( 59) - W( 60) - W( 65) - W( 66) + W( 69)  
     &           + W( 71) + W( 75) + W( 78) + W( 79) + W( 80)  
     &           +2.0*W( 81) - W( 82) + W( 83) + W( 84)  
     &           + W( 89) + W( 92) + W( 93) - W(100) + W(101)  
     &           + W(105) + W(109) - W(114) - W(123) - W(127)  
     &           - W(131) - W(132) -2.0*W(133) - W(134)  
     &           - W(135) -2.0*W(136) - W(137) - W(138)  
     &           - W(139) + W(140) - W(142) + W(144) - W(145)  
     &           - W(146) + W(153) - W(154) + W(156) - W(159)  
     &           + W(162) + W(164) - W(167) + W(187)  
      WGLOB( 3) = - W(  3) - W(  4) - W(  5) - W( 11) - W( 13)  
     &           - W( 14) + W( 15) + W( 19) + W( 20) + W( 22)  
     &           - W( 26) + W( 27) + W( 28) + W( 30) - W( 41)  
     &           - W( 48) - W( 49) + W( 53) + W( 54) + W( 55)  
     &           + W( 57) + W( 60) +2.0*W( 62) +2.0*W( 65)  
     &           +2.0*W( 66) - W( 67) - W( 68) + W( 74)  
     &           + W( 75) + W( 77) - W( 83) - W( 86) + W( 87)  
     &           + W( 88) + W( 89) + W( 91) + W( 93) - W( 96)  
     &           - W( 98) - W( 99) - W(101) + W(102) + W(104)  
     &           - W(106) - W(108) - W(110) - W(111) - W(112)  
     &           - W(113) -2.0*W(114) - W(116) - W(118)  
     &           - W(119) - W(125) + W(127) - W(131) - W(132)  
     &           + W(134) - W(135) - W(137) - W(139) + W(140)  
     &           + W(142) + W(143) - W(152) + W(153) - W(154)  
     &           + W(156) - W(157) - W(160) -2.0*W(166)  
     &           + W(167) + W(168) - W(169) + W(172) + W(175)  
     &           + W(177) - W(179) - W(180) + W(181) - W(184)  
     &           +2.0*W(187) + W(195) - W(206)  
      WGLOB( 4) = + W( 31) + W( 35) + W( 36) + W( 37) +2.0*W( 38)  
     &           + W( 41) + W( 42) - W( 44) + W( 45) -2.0*W( 51)  
     &           - W( 54) - W( 55) - W( 59) + W( 64) + W(127)  
     &           + W(133) + W(134) + W(136) + W(142) + W(143)  
     &           + W(145) + W(146) - W(147) + W(159) + W(161)  
     &           - W(162) - W(168) + W(177)  
      WGLOB( 5) = + W( 34)  
      WGLOB( 6) = - W( 38) + W( 46) + W( 47) + W( 48) + W( 49)  
     &           + W(164)  
      WGLOB( 7) = - W( 38) + W( 51) + W( 54) + W( 55) + W( 58)  
     &           + W( 59) + W( 60) + W( 62) + W( 65) + W( 66)  
     &           + W(167) + W(168) + W(187)  
      WGLOB( 8) = + W( 97) + W( 98) + W( 99) + W(110) + W(111)  
     &           + W(114) + W(115) + W(116) + W(117) + W(118)  
     &           + W(119) + W(145) + W(146)  
      WGLOB( 9) = - W( 99) - W(110) - W(116) + W(123) + W(124)  
     &           - W(145) + W(151)  
      WGLOB(10) = + W(130) - W(133) - W(134) - W(135) - W(136)  
     &           - W(137) + W(144) + W(153) + W(156) + W(164)  
     &            
      WGLOB(11) = - W( 23) - W( 24) - W( 25) - W( 41) - W( 42)  
     &           - W( 54) - W( 59) -2.0*W( 60) -2.0*W( 65)  
     &           -2.0*W( 66) - W(127) - W(133) - W(134)  
     &           - W(136) + W(140) + W(141) - W(142) - W(143)  
     &           - W(145) - W(146) - W(154) + W(162) -2.0*W(167)  
     &           - W(168)  
      WGLOB(12) = - W( 78) - W( 79) - W( 81) - W( 83) - W( 84)  
     &           - W( 92) - W(101) - W(105) - W(106) - W(108)  
     &           - W(109) - W(112) - W(113) - W(115) - W(117)  
     &           - W(118) - W(119) - W(143) - W(152) + W(157)  
     &           + W(158) + W(162) + W(168)  
      WGLOB(13) = + W(186) + W(187)  
      WGLOB(14) = - W( 57) - W( 60) - W( 62) - W( 65) - W( 66)  
     &           + W(166) - W(167) - W(168) - W(187) + W(194)  
     &            
      WGLOB(15) = + W(202)  
      WGLOB(16) = + W( 78) + W( 79) + W( 81) + W( 83) + W( 84)  
     &           + W( 92) + W(101) + W(105) + W(106) + W(108)  
     &           + W(109) + W(112) + W(113) + W(115) + W(117)  
     &           + W(118) + W(119) - W(127) - W(133) - W(134)  
     &           - W(136) - W(142) - W(145) - W(146) + W(152)  
     &           + W(160) + W(169) - W(177) + W(179) + W(180)  
     &           + W(183) + W(184)  


C-----------------------------------------------------------------------
C --- CALCULATE THE PRODUCTION RATES OF THE GLOBAL SPECIES
C-----------------------------------------------------------------------

C  O2    
      WDOT( 1) = - WGLOB( 2) - WGLOB( 4) - WGLOB( 7)  
     &           + WGLOB(11)  
C  H     
      WDOT( 2) = -2.0*WGLOB( 1) -2.0*WGLOB( 5) + WGLOB( 8)  
     &           + WGLOB(12) +2.0*WGLOB(14) +2.0*WGLOB(15)  
     &           -2.0*WGLOB(16)  
C  H2    
      WDOT( 3) = + WGLOB( 1) + WGLOB( 3) + WGLOB( 6)  
     &           + WGLOB(16)  
C  H2O   
      WDOT( 4) = - WGLOB( 3) + WGLOB( 4) + WGLOB( 5)  
     &           + WGLOB( 8) - WGLOB(11) + WGLOB(12)  
     &           - WGLOB(13) -3.0*WGLOB(15) +3.0*WGLOB(16)  
     &            
C  N2    
      WDOT( 5) = + WGLOB( 8) + WGLOB( 9)  
C  CO2   
      WDOT( 6) = +2.0*WGLOB( 2) + WGLOB( 3) -2.0*WGLOB( 5)  
     &           + WGLOB( 9) - WGLOB(10) - WGLOB(11)  
     &           + WGLOB(13) +3.0*WGLOB(15) -3.0*WGLOB(16)  
     &            
C  CO    
      WDOT( 7) = -2.0*WGLOB( 2) - WGLOB( 3) +2.0*WGLOB( 5)  
     &           - WGLOB( 9) + WGLOB(10) - WGLOB(13)  
     &           -3.0*WGLOB(15) +3.0*WGLOB(16)  
C  CH2O  
      WDOT( 8) = + WGLOB( 4) +2.0*WGLOB( 7) + WGLOB(11)  
     &           + WGLOB(15)  
C  CH4   
      WDOT( 9) = - WGLOB( 4) - WGLOB( 5) - WGLOB(12)  
     &           + WGLOB(13) +2.0*WGLOB(15) - WGLOB(16)  
     &            
C  CH3OH 
      WDOT(10) = + WGLOB( 5)  
C  C2H6  
      WDOT(11) = - WGLOB( 6)  
C  C2H4  
      WDOT(12) = + WGLOB( 6) - WGLOB( 7) - WGLOB(14)  
     &           + WGLOB(15)  
C  C2H2  
      WDOT(13) = + WGLOB(13) + WGLOB(14)  
C  NO    
      WDOT(14) = - WGLOB( 8) - WGLOB(12)  
C  NH3   
      WDOT(15) = - WGLOB( 8) - WGLOB(16)  
C  N2O   
      WDOT(16) = - WGLOB( 9)  
C  HCN   
      WDOT(17) = - WGLOB(10) + WGLOB(12) + WGLOB(13)  
     &           + WGLOB(15) + WGLOB(16)  
C  HNCO  
      WDOT(18) = + WGLOB(10)  
C  C4H5N 
      WDOT(19) = - WGLOB(13)  
C  C6H11N
      WDOT(20) = - WGLOB(15)  

      RETURN
      END
