REACTION RATES
H1(E-,NU)N1
Q1 = -0.782
SEE AP.J.SUPPL. 48,279(1982) & AP.J. 293,1(1985)
H1(P,E+NU)H2
Q2 = 1.442
T9max2 = 3
f2 = 4.01e-15/T923*exp(-3.380/T913)*(1.0+0.123*T913+1.09*T923
| +0.938*T9)
H1(E-P,NU)H2
Q3 = 1.442
T9max3 = 3
f3 = 1.36e-20/T976*exp(-3.380/T913)*(1.0-0.729*T913+9.82*T923)
H2(P,G)HE3
Q4 = 5.494
f4 = 2.24e+03/T923*exp(-3.720/T913)*(1.0+0.112*T913+3.38*T923
| +2.65*T9)
r4 = f4*1.63e+10*T932*exp(-63.750/T9)
H2(P,N)2H1
Q5 = -2.225
f5 = 3.35e+07*exp(-3.720/T913-25.815/T9)*(1.0+0.784*T913+0.346
| *T923+0.690*T9)
r5 = 3.35e+07*exp(-3.720/T913)*(1.0+0.784*T913+0.346
| *T923+0.690*T9)*4.24e-10/T932
H2(D,G)HE4
Q6 = 23.847
f6 = 4.84e+01/T923*exp(-4.258/T913)*(1.0+0.098*T913-0.203*T923
| -0.139*T9+0.106*T943+0.185*T953)
r6 = f6*4.53e+10*T932*exp(-276.729/T9)
H2(D,N)HE3
Q7 = 3.269
f7 = 3.88e+08/T923*exp(-4.258/T913)*(1.0+0.098*T913+0.418*T923
| +0.287*T9+0.638*T943+1.112*T953)
r7 = f7*1.73e+00*exp(-37.935/T9)
H2(D,P)H3
Q8 = 4.033
f8 = 4.13e+08/T923*exp(-4.258/T913)*(1.0+0.098*T913+4.39e-02
| *T923+3.01e-02*T9+0.543*T943+0.946*T953)
r8 = f8*1.73e+00*exp(-46.798/T9)
H3(P,G)HE4
Q9 = 19.814
f9 = 2.20e+04/T923*exp(-3.869/T913)*(1.0+0.108*T913+1.68*T923
| +1.26*T9+0.551*T943+1.06*T953)
r9 = f9*2.61e+10*T932*exp(-229.932/T9)
H3(P,N)HE3
Q10 = -0.764
f10 = 7.07e+08*(1.0-0.150*T912+0.098*T9)*exp(-8.863/T9)
r10 = 7.07e+08*(1.0-0.150*T912+0.098*T9)*9.98e-01
H3(D,N)HE4
Q11 = 17.589
f11 = 8.09e+10/T923*exp(-4.524/T913-(T9/0.120)**2)*(1.0+0.092
| *T913+1.80*T923+1.16*T9+10.52*T943+17.24*T953)
| +8.73e+08/T923*exp(-0.523/T9)
r11 = f11*5.54e+00*exp(-204.117/T9)
H3(T,2N)HE4
Q12 = 11.332
f12 = 1.67e+09/T923*exp(-4.872/T913)*(1.0+0.086*T913-0.455*T923
| -0.272*T9+0.148*T943+0.225*T953)
r12 = f12*3.38e-10/T932*exp(-131.504/T9)
HE3(E-,NU)H3
Q13 = -0.019
T9max13 = 3
f13 = 7.71e-12*T932*(1.0+6.48*T9+7.48*T9**2+2.91*T9**3)
| *exp(-0.2158/T9)
HE3(P,E+NU)HE4
Q14 = 19.795
T9max14 = 3
f14 = 8.78e-13/T923*exp(-6.141/T913)
Q14 = 10.155 ! EXCLUSIVE OF NU-ENERGY
HE3(D,P)HE4
Q15 = 18.353
f15 = 5.86e+10/T923*exp(-7.181/T913-(T9/0.315)**2)*(1.0+0.058
| *T913+0.142*T923+5.78e-02*T9+2.25*T943+2.32*T953)+4.36e+08
| /T912*exp(-1.720/T9)
r15 = f15*5.55e+00*exp(-212.980/T9)
HE3(T,D)HE4
Q16 = 14.320
T9A = T9/(1.0+0.128*T9)
T9A13 = T9A**(1.0/3.0)
T9A56 = T9A**(5.0/6.0)
f16 = 5.46e+09*T9A56/T932*exp(-7.733/T9A13)
r16 = f16*1.6e+00*exp(-166.182/T9)
HE3(T,NP)HE4
Q17 = 12.096
T9A = T9/(1.0+0.115*T9)
T9A13 = T9A**(1.0/3.0)
T9A56 = T9A**(5.0/6.0)
f17 = 7.71e+09*T9A56/T932*exp(-7.733/T9A13)
r17 = f17*3.39e-10/T932*exp(-140.367/T9)
HE3(HE3,2P)HE4
Q18 = 12.860
f18 = 6.04e+10/T923*exp(-12.276/T913)*(1.0+0.034*T913-0.522*T923
| -0.124*T9+0.353*T943+0.213*T953)
r18 = f18*3.39e-10/T932*exp(-149.230/T9)
HE4(NN,G)HE6
Q19 = 0.975
f19 = rc19*4.04e-11/T9**2*(1.0+0.138*T9)*exp(-9.585/T9)
r19 = f19*1.08e+20*T9**3*exp(-11.319/T9)
HE4(NP,G)LI6
Q20 = 3.700
f20 = 4.62e-06/T9**2*(1.0+0.075*T9)*exp(-19.353/T9)
r20 = f20*7.22e+19*T9**3*exp(-42.933/T9)
HE4(D,G)LI6
Q21 = 1.475
f21 = 3.01e+01/T923*exp(-7.423/T913)*(1.0+0.056*T913-4.85*T923
| +8.85*T9-0.585*T943-0.584*T953)+8.55e+01/T932*exp(-8.228/T9)
r21 = f21*1.53e+10*T932*exp(-17.118/T9)
HE4(T,G)LI7
Q22 = 2.468
f22 = 8.67e+05/T923*exp(-8.080/T913)*(1.0+0.052*T913-0.448*T923
| -0.165*T9+0.144*T943+0.134*T953)
r22 = f22*1.11e+10*T932*exp(-28.640/T9)
HE4(T,N)LI6
Q23 = -4.782
T9A = T9/(1.0+49.18*T9)
T9A13 = T9A**(1.0/3.0)
T9A56 = T9A**(5.0/6.0)
T9A32 = T9A**(3.0/2.0)
f23 = 1.80e+08*exp(-55.494/T9)*(1.0-0.261*T9A32/T932)+2.72e+09
| /T932*exp(-57.884/T9)
r23 = 1.80e+08*(1.0-0.261*T9A32/T932)+2.72e+09/T932
| *exp(-57.884/T9)*9.35e-01
HE4(HE3,G)BE7
Q24 = 1.588
T9A = T9/(1.0+4.95e-02*T9)
T9A13 = T9A**(1.0/3.0)
T9A56 = T9A**(5.0/6.0)
f24 = 5.61e+06*T9A56/T932*exp(-12.826/T9A13)
r24 = f24*1.11e+10*T932*exp(-18.423/T9)
HE4(AN,G)BE9
Q25 = 1.574
f25 = 2.59e-06/(T9**2*(1.0+0.344*T9))*exp(-1.062/T9)
r25 = f25*5.84e+19*T9**3*exp(-18.260/T9)
HE4(A)BE8
Q26 = -0.092
T9min26 = 0.01
f26 = 7.40e+05/T932*exp(-1.0663/T9)+4.164e+09/T923
| *exp(-13.490/T913-(T9/0.098)**2)*(1.0+0.031*T913+8.009*T923
| +1.732*T9+49.883*T943+27.426*T953)
r26 = 7.40e+05/T932*1.40e+10*T932+4.164e+09/T923
| *exp(-13.490/T913-(T9/0.098)**2+1.0663/T9)
| *(1.0+0.031*T913+8.009*T923
| +1.732*T9+49.883*T943+27.426*T953)*1.40e+10*T932
BE8(A,G)C12
Q27 = 7.367
f27 = 1.30e+02/T932*exp(-3.3364/T9)+2.510e+07/T923
| *exp(-23.570/T913-(T9/0.235)**2)*(1.0+0.018*T913+5.249*T923
| +0.650*T9+19.176*T943+6.034*T953)
r27 = f27*4.30e+10*T932*exp(-85.486/T9)
HE4(2A,G)C12
Q28 = 7.275
T9<0.08 or T9=0.08
f28 = 2.90e-16*(7.40e+05/T932*exp(-1.0663/T9)+4.16e+09/T923
| *exp(-13.490/T913-(T9/0.098)**2)*(1.0+0.031*T913+8.009*T923
| +1.732*T9+49.883*T943+27.426*T953))*(1.30e+02/T932
| *exp(-3.3364/T9)+2.510e+07/T923
| *exp(-23.570/T913-(T9/0.235)**2)*(1.0+0.018*T913+5.249*T923
| +0.650*T9+19.176*T943+6.034*T953))*(0.01+0.2*(1.0+4.0
| *exp(-(0.025/T9)**3.263))/(1.0+4.0
| *exp(-(T9/0.025)**9.227)))
| +rc28*1.35e-07/T932*exp(-24.811/T9)
T9>0.08
f28 = 2.79e-08/T9**3*exp(-4.4027/T9)+rc28
| *1.35e-07/T932*exp(-24.811/T9)
Q = 14.437 IF Cl2(A,G)Ol6 ALWAYS FOLLOWS
r28 = f28*2.00e+20*T9**3*exp(-84.420/T9)
LI6(P,G)BE7
Q29 = 5.606
T9A = T9/(1.0-9.69e-02*T9+2.84e-02*T953
| /(1.0-9.69e-02*T9)**(2.0/3.0))
T9A13 = T9A**(1.0/3.0)
T9A56 = T9A**(5.0/6.0)
f29 = 6.69e+05*T9A56/T932*exp(-8.413/T9A13)
r29 = f29*1.19e+10*T932*exp(-65.054/T9)
LI6(P,HE3)HE4
Q30 = 4.018
f30 = 3.73e+10/T923*exp(-8.413/T913-(T9/5.50)**2)*(1.0+0.050*T913
| -0.061*T923-0.021*T9+0.006*T943+0.005*T953)+1.33e+10
| /T932*exp(-17.763/T9)+1.29e+09/T9*exp(-21.820/T9)
r30 = f30*1.07e+00*exp(-46.631/T9)
LI6(A,G)B10
Q31 = 4.460
f31 = 4.06e+06/T923*exp(-18.790/T913-(T9/1.326)**2)*(1.0+0.022
| *T913+1.54*T923+0.239*T9+2.20*T943+0.869*T953)+1.91e+03
| /T932*exp(-3.484/T9)+1.01e+04/T9*exp(-7.269/T9)
r31 = f31*1.58e+10*T932*exp(-51.753/T9)
LI7(P,N)BE7
Q32 = -1.644
f32 = 5.15e+09*exp(-1.167*T913-19.081/T9)+7.84e+09/T932
| *exp(-22.832/T9)
r32 = 5.15e+09*exp(-1.167*T913)+7.84e+09/T932
| *exp(-3.751/T9)*9.98e-01
LI7(P,G)BE8
Q33 = 17.254
f33 = 1.56e+05/T923*exp(-8.472/T913-(T9/1.696)**2)
| *(1.0+0.049*T913+2.498*T923
| +0.860*T9+3.518*T943+3.080*T953)+1.55e+06/T932
| *exp(-4.478/T9)
r33 = f33*6.55e+10*T932*exp(-200.225/T9)
LI7(P,A)HE4
Q34 = 17.346
T9A = T9/(1.0+0.759*T9)
T9A13 = T9A**(1.0/3.0)
T9A56 = T9A**(5.0/6.0)
f34 = 1.096e+09/T923*exp(-8.472/T913)-4.830e+08*T9A56/T932
| *exp(-8.472/T9A13)+1.06e+10/T932*exp(-30.442/T9)
r34 = f34*4.69e+00*exp(-201.291/T9)
r34 APPLIES TO BOTH LI7(P,A)HE4 AND LI7(P,AG)HE4
LI7(P,AG)HE4
Q35 = 17.346
T9A = T9/(1.0+0.759*T9)
T9A13 = T9A**(1.0/3.0)
T9A56 = T9A**(5.0/6.0)
Note: f35=f11+f14
f35 = (1.56e+05/T923*exp(-8.472/T913-(T9/1.696)**2)*(1.0+0.049
| *T913+2.498*T923+0.860*T9+3.518*T943+3.080*T953)+1.55e+06
| /T932*exp(-4.478/T9))+(1.096e+09/T923
| *exp(-8.472/T913)-4.830e+08*T9A56/T932*exp(-8.472/T9A13)
| +1.06e+10/T932*exp(-30.442/T9))
r35 = r34
LI7(D,N)2HE4
Q36 = 15.121
f36 = 2.92e+11/T923*exp(-10.259/T913)
r36 = f36*9.95e-10/T932*exp(-175.476/T9)
LI7(T,2N)2HE4
Q37= 8.864
f37 = 8.81e+11/T923*exp(-11.333/T913)
r37 = f37*1.22e-19/T9**3*exp(-102.864/T9)
LI7(HE3,NP)2HE4
Q38 = 9.628
f38 = 1.11e+13/T923*exp(-17.989/T913)
r38 = f38*6.09e-20/T9**3*exp(-111.727/T9)
LI7(A,G)B11
Q39 = 8.664
f39 = 3.55e+07/T923*exp(-19.161/T913-(T9/4.195)**2)*(1.0+0.022
| *T913+0.775*T923+0.118*T9+0.884*T943+0.342*T953)
| +3.33e+02/T932*exp(-2.977/T9)+4.10e+04/T9*exp(-6.227/T9)
r39 = f39*4.02e+10*T932*exp(-100.538/T9)
LI7(A,N)B10
Q40 = -2.790
f40 = 3.84e+08*exp(-32.382/T9)
r40 = 3.84e+08*1.32e+00
BE7(E-,NU+G)LI7
Q41 = 0.862
T9max41 = 3
f41 = 1.34e-10/T912*(1.0-0.537*T913+3.86*T923+0.0027/T9
| *exp(2.515e-03/T9))
Q = 0.049 EXCLUSIVE OF NU ENERGY
RATE MUST NOT EXCEED l.51e-07/(RHO*(1.0+X)/2.0)
F0R T9 LESS THAN 0.001.
BE7(P,G)B8
Q42 = 0.137
f42 = 3.11e+05/T923*exp(-10.262/T913)+2.53e+03/T932
| *exp(-7.306/T9)
r42 = f42*1.30e+10*T932*exp(-1.595/T9)
BE7(D,P)2HE4
Q43 = 16.766
f43 = 1.07e+12/T923*exp(-12.428/T913)
r43 = f43*9.97e-10/T932*exp(-194.557/T9)
BE7(T,NP)2HE4
Q44 = 10.508
f44 = 2.91e+12/T923*exp(-13.729/T913)
r44 = f44*6.09e-20/T9**3*exp(-121.944/T9)
BE7(HE3,2P)2HE4
Q45 = 11.272
f45 = 6.11e+13/T923*exp(-21.793/T913)
r45 = f45*1.22e-19/T9**3*exp(-130.807/T9)
BE7(A,G)C11
Q46 = 7.544
f46 = 8.45e+07/T923*exp(-23.212/T913-(T9/4.769)**2)*(1.0+0.018
| *T913+0.488*T923+0.061*T9+0.296*T943+0.095*T953)+1.25e+04
| /T932*exp(-6.510/T9)+1.29e+05/T954*exp(-10.039/T9)
r46 = f46*4.02e+10*T932*exp(-87.539/T9)
BE9(P,G)B10
Q47 = 6.586
f47 = 1.33e+07/T923*exp(-10.359/T913-(T9/0.846)**2)
| *(1.0+0.040*T913+1.52*T923+0.428*T9+2.15*T943+1.54*T953)
| +9.64e+04/T932*exp(-3.445/T9)+2.72e+06/T932*exp(-10.620/T9)
r47 = f47*9.73e+09*T932*exp(-76.427/T9)
BE9(P,N)B9
Q48 = -1.850
f48 = 5.58e+07*(1.0+0.042*T912+0.985*T9)*exp(-21.473/T9)
| +1.02e+09/T932*exp(-26.725/T9)
r48 = 5.58e+07*(1.0+0.042*T912+0.985*T9)*9.98e-01
| +1.02e+09/T932*exp(-5.252/T9)*9.98e-01
BE9(P,D)2HE4
Q49 = 0.651
f49 = 2.11e+11/T923*exp(-10.359/T913-(T9/0.520)**2)*(1.0+0.040
| *T913+1.09*T923+0.307*T9+3.21*T943+2.30*T953)+5.79e+08/T9
| *exp(-3.046/T9)+8.50e+08/T934*exp(-5.800/T9)
r49 = f49*8.07e-11/T932*exp(-7.555/T9)
BE9(P,A)LI6
Q50 = 2.126
f50 = 2.11e+11/T923*exp(-10.359/T913-(T9/0.520)**2)*(1.0+0.040
| *T913+1.09*T923+0.307*T9+3.21*T943+2.30*T953)+4.51e+08/T9
| *exp(-3.046/T9)+6.70e+08/T934*exp(-5.160/T9)
r50 = f50*6.18e-01*exp(-24.674/T9)
BE9(A,N)C12
Q51 = 5.701
f51 = 4.62e+13/T923*exp(-23.870/T913-(T9/0.049)**2)*(1.0+0.017
| *T913+8.57*T923+1.05*T9+74.51*T943+23.15*T953)+7.34e-05
| /T932*exp(-1.184/T9)+2.27e-01/T932*exp(-1.834/T9)
| +1.26e+05/T932*exp(-4.179/T9)+2.40e+08*exp(-12.732/T9)
r51 = f51*1.03e+01*exp(-66.160/T9)
B10(P,G)C11
Q52 = 8.690
f52 = 4.61e+05/T923*exp(-12.062/T913-(T9/4.402)**2)*(1.0+0.035
| *T913+0.426*T923+0.103*T9+0.281*T943+0.173*T953)+1.93e+05
| /T932*exp(-12.041/T9)+1.14e+04/T932*exp(-16.164/T9)
r52 = f52*3.03e+10*T932*exp(-100.840/T9)
B10(P,A)BE7
Q53 = 1.146
f53 = 1.26e+11/T923*exp(-12.062/T913-(T9/4.402)**2)*(1.0+0.035
| *T913-0.498*T923-0.121*T9+0.300*T943+0.184*T953)+2.59e+09
| /T9*exp(-12.260/T9)
r53 =f53*7.54e-01*exp(-13.301/T9)
B10(A,N)N13
Q54 = 1.059
f54 = 1.20e+13/T923*exp(-27.989/T913-(T9/9.589)**2)
r54 = f54*9.34e+00*exp(-12.287/T9)
B11(P,G)C12
Q55 = 15.957
f55 = 4.62e+07/T923*exp(-12.095/T913-(T9/0.239)**2)*(1.0+0.035
| *T913+3.00*T923+0.723*T9+9.91*T943+6.07*T953)+7.89e+03
| /T932*exp(-1.733/T9)+9.68e+04/T915*exp(-5.617/T9)
r55 = f55*7.01e+10*T932*exp(-185.173/T9)
B11(P,N)C11
Q56 = -2.764
f56 = 1.69e+08*(1.0-0.048*T912+0.010*T9)*exp(-32.080/T9)
r56 = 1.69e+08*(1.0-0.048*T912+0.010*T9)*9.98e-01
B11(P,A)2HE4
Q57 = 8.682
f57 = 2.20e+12/T923*exp(-12.095/T913-(T9/1.644)**2)*(1.0+0.034
| *T913+0.140*T923+0.034*T9+0.190*T943+0.116*T953)+4.03e+06
| /T932*exp(-1.734/T9)+6.73e+09/T932*exp(-6.262/T9)+3.88e+09/T9
| *exp(-14.154/T9)
r57 = f57*3.50e-10/T932*exp(-100.753/T9)
B11(A,N)N14
Q58 = 0.158
f58 = 6.97e+12/T923*exp(-28.234/T913-(T9/0.140)**2)*(1.0+0.015
| *T913+8.115*T923+0.838*T9+39.804*T943+10.456*T953)+1.79e+00
| /T932*exp(-2.827/T9)+1.71e+03/T932*exp(-5.178/T9)+4.49e+06
| *T935*exp(-8.596/T9)
r58 = f58*3.67e+00*exp(-1.835/T9)
B11(A,P)C14
Q59 = 0.784
f59 = 5.37e+11/T923*exp(-28.234/T913-(T9/0.347)**2)*(1.00+0.015
| *T913+5.575*T923+0.576*T9+15.888*T943+4.174*T953)+5.44e-03
| /T932*exp(-2.827/T9)+3.36e+02/T932*exp(-5.178/T9)+5.32e+06
| /T938*exp(-11.617/T9)
r59 = f59*1.10e+01*exp(-9.098/T9)
C11(P,G)N12
Q60 = 0.601
f60 = 4.24e+04/T923*exp(-13.658/T913-(T9/1.627)**2)*(1.00+0.031
| *T913+3.11*T923+0.665*T9+4.61*T943+2.50*T953)+8.84e+03
| /T932*exp(-7.021/T9)
r60 = f60*2.33e+10*T932*exp(-6.975/T9)
C12(P,G)N13
Q61 = 1.944
f61 = 2.04e+07/T923*exp(-13.690/T913-(T9/1.500)**2)*(1.00+0.030
| *T913+1.19*T923+0.254*T9+2.06*T943+1.12*T953)+1.08e+05
| /T932*exp(-4.925/T9)+2.15e+05/T932*exp(-18.179/T9)
r61 = f61*8.84e+09*T932*exp(-22.553/T9)
C12(A,G)O16
Q62 = 7.162
f62 = 1.04e+08/T9**2/(1.00+0.0489/T923)**2*exp(-32.120/T913-(T9
| /3.496)**2)+1.76e+08/T9**2/(1.00+0.2654/T923)**2
| *exp(-32.120/T913)+1.25e+03/T932*exp(-27.499/T9)
| +1.43e-02*T9**5*exp(-15.541/T9)
r62 = f62*5.13e+10*T932*exp(-83.111/T9)
C12(A,N)O15
Q63 = -8.502
f63 = 2.48e+07*(1.00*0.188*T912+0.015*T9)*exp(-98.661/T9)
r63 = 2.48e+07*(1.00*0.188*T912+0.015*T9)*1.41e+00
C13(P,G)N14
Q64 = 7.551
f64 = 8.01e+07/T923*exp(-13.717/T913-(T9/2.000)**2)*(1.0+0.030
| *T913+0.958*T923+0.204*T9+1.39*T943+0.753*T953)+1.21e+06
| /T965*exp(-5.701/T9)
r64 = f64*1.19e+10*T932*exp(-87.621/T9)
C13(P,N)N13
Q65 = -3.003
f65 = 1.88e+08*(1.0-0.167*T912+0.037*T9)*exp(-34.846/T9)
r65 = 1.88e+08*(1.0-0.167*T912+0.037*T9)*9.98e-01
C13(A,N)O16
Q66 = 2.216
f66 = 6.77e+15/T923*exp(-32.329/T913-(T9/1.284)**2)*(1.0+0.013
| *T913+2.04*T923*0.184*T9)+3.82e+05/T932*exp(-9.373/T9)
| +1.41e+06/T932*exp(-11.873/T9)+2.00e+09/T932*exp(-20.409/T9)
| +2.92e+09/T932*exp(-29.283/T9)
r66 = f66*5.79e+00*exp(-25.711/T9)
C14(P,G)N15
Q67 = 10.207
f67 = 6.80e+06/T923*exp(-13.741/T913-(T9/5.721)**2)*(1.0+0.030
| *T913+0.503*T923+0.107*T9+0.213*T943+0.115*T953)+5.36e+03
| /T932*exp(-3.811/T9)+9.82e+04/T913*exp(-4.739/T9)
r67 = f67*9.00e+09*T932*exp(-118.452/T9)
C14(P,N)N14
Q68 = -0.626
f68 = 7.19e+05*(1.0+0.361*T912+0.502*T9)*exp(-7.263/T9)
| +3.34e+08/T912*exp(-12.246/T9)
r68 = 7.19e+05*(1.0+0.361*T912+0.502*T9)*3.33e-01
| +3.34e+08/T912*exp(-4.983/T9)*3.33e-01
C14(A,G)O18
Q69 = 6.227
f69 = 3.375e+08/T9**2*exp(-32.513/T913)+1.528e+09/T923
| *exp(-32.513/T913-(T9/2.662)**2)*(1.0+0.0128*T913-0.869
| *T923-0.0779*T9+0.321*T943+0.0732*T953)+9.29e-08/T932
| *exp(-2.048/T9)+2.77e+03/T945*exp(-9.876/T9)
r69 = f69*5.42e+10*T932*exp(-72.262/T9)
N13(P,G)O14
Q70 = 4.628
f70 = 4.04e+07/T923*exp(-15.202/T913-(T9/1.191)**2)*(1.0+0.027
| *T913-0.803*T923-0.154*T9+5.00*T943+2.44*T953)+2.43e+05/T932
| *exp(-6.348/T9)
r70 = f70*3.57e+10*T932*exp(-53.706/T9)
N14(P,G)O15
Q71 = 7.297
f71 = 4.90e+07/T923*exp(-15.228/T913-(T9/3.294)**2)*(1.0+0.027
| *T913-0.778*T923-0.149*T9+0.261*T943+0.127*T953)+2.37e+03
| /T932*exp(-3.011/T9)+2.19e+04*exp(-12.530/T9)
r71 = f71*2.70e+10*T932*exp(-84.678/T9)
N14(P,N)O14
Q72 = -5.925
f72 = 6.74e+07*(1.0+0.658*T912+0.379*T9)*exp(-68.762/T9)
r72 = 6.74e+07*(1.0+0.658*T912+0.379*T9)*2.99e+00
N14(P,A)C11
Q73 = -2.923
T9A = T9/(1.0+4.78e-02*T9+7.56e-03*T953
| /(1.0+4.78e-02*T9)**(2.0/3.0))
T9A13 = T9A**(1.0/3.0)
T9A56 = T9A**(5.0/6.0)
f73 = 2.63e+16*T9A56/T932*exp(-31.883/T9A13-33.915/T9)
r73 = 2.63e+16*T9A56/T932*exp(-31.883/T9A13)*2.72e-01
N14(A,G)F18
Q74 = 4.415
f74 = 7.78e+09/T923*exp(-36.031/T913-(T9/0.881)**2)*(1.00+0.012
| *T913+1.45*T923+0.117*T9+1.97*T943+0.406*T953)+2.36e-10
| /T932*exp(-2.798/T9)+2.03e+00/T932*exp(-5.054/T9)+1.15e+04
| /T923*exp(-12.310/T9)
r74 = f74*5.42e+10*T932*exp(-51.236/T9)
N14(A,N)F17
Q75 = -4.735
f75 = 5.24e+09*(1.0-1.15*T912+0.365*T9)
| *exp(-54.942/T9-(T9/2.798)**2)+3.28e+10/T932*exp(-70.708/T9)
r75 = 5.24e+09*(1.0-1.15*T912+0.365*T9)
| *exp(-(T9/2.798)**2)*1.48e+00+3.28e+10/T932*exp(-15.766/T9)
| *1.48e+00
N15(P,G)O16
Q76 = 12.128
f76 = 9.78e+08/T923*exp(-15.251/T913-(T9/0.450)**2)*(1.00+0.027
| *T913+0.219*T923+0.042*T9+6.83*T943+3.32*T953)+1.11e+04
| /T932*exp(-3.328/T9)+1.49e+04/T932*exp(-4.665/T9)+3.80e+06
| /T932*exp(-11.048/T9)
r76 = f76*3.62e+10*T932*exp(-140.734/T9)
N15(P,N)O15
Q77 = -3.536
f77 = 3.51e+08*(1.0+0.452*T912-0.191*T9)*exp(-41.037/T9)
r77 = 3.51e+08*(1.0+0.452*T912-0.191*T9)*9.98e-01
N15(P,A)C12
Q78 = 4.966
f78 = 1.08e+12/T923*exp(-15.251/T913-(T9/0.522)**2)*(1.00+0.027
| *T913+2.62*T923+0.501*T9+5.36*T943+2.60*T953)+1.19e+08
| /T932*exp(-3.676/T9)+5.41e+08/T912*exp(-8.926/T9)+rc78
| *4.72e+08/T932*exp(-7.721/T9)+2.20e+09/T932*exp(-11.418/T9)
r78 = f78*7.06e-01*exp(-57.623/T9)
N15(A,G)F19
Q79 = 4.014
f79 = 2.54e+10/T923*exp(-36.211/T913-(T9/0.616)**2)*(1.0+0.012
| *T913+1.69*T923+0.136*T9+1.91*T943+0.391*T953)+9.83e-03
| /T932*exp(-4.232/T9)+1.52e+03*T9*exp(-9.747/T9)
r79 = f79*5.54e+10*T932*exp(-46.578/T9)
N15(A,N)F18
Q80 = -6.418
f80 = 3.14e+08*(1.0-0.641*T912+0.108*T9)*exp(-74.479/T9)
r80 = 3.14e+08*(1.0-0.641*T912+0.108*T9)*2.00e+00
O14(A,G)NE18
Q81 = 5.112
f81 = 9.47e+08/T923*exp(-39.388/T913-(T9/0.717)**2)*(1.0+0.011
| *T913+1.974*T923+0.146*T9+3.036*T943+0.572*T953)+1.16e-01
| /T932*exp(-11.733/T9)+3.39e+01/T932*exp(-22.609/T9)
| +9.10e-03*T9**5.0*exp(-12.159/T9)
r81 = f81*5.42e+10*T932*exp(-59.328/T9)
O14(A,P)F17
Q82 = 1.191
f82 = 1.68e+13/T923*exp(-39.388/T913-(T9/0.717)**2)*(1.0+0.011
| *T913+13.117*T923+0.971*T9+85.295*T943+16.061*T953)+3.31e+04
| /T932*exp(-11.733/T9)+1.79e+07/T932*exp(-22.609/T9)+9.00e+03
| *T9**(11.0/3.0)*exp(-12.517/T9)
r82 = f82*4.93e-01*exp(-13.820/T9)
O15(A,G)NE19
Q83 = 3.529
f83 = 3.57e+11/T923*exp(-39.584/T913-(T9/3.000)**2)*(1.0+0.011
| *T913-0.273*T923-0.020*T9)+5.10e+10/T923*exp(-39.584/T913
| -(T9/1.937)**2)*(1.0+0.011*T913+1.59*T923+0.117*T9+1.81
| *T943+0.338*T953)+3.95e-01/T932*exp(-5.849/T9)+1.90e+01
| *(T9**2.85)*exp(-7.356/T9-(T9/8.000)**2)
r83 = f83*5.54e+10*T932*exp(-40.957/T9)
O16(P,G)F17
Q84 = 0.600
f84 = 1.50e+08/(T923*(1.0+2.13*(1.0-exp(-0.728*T923))))
| *exp(-16.692/T913)
r84 = f84*3.03e+09*T932*exp(-6.968/T9)
O16(P,A)N13
Q85 = -5.218
T9A = T9/(1.0+7.76e-02*T9+2.64e-02*T953
| /(1.0+7.76e-02*T9)**(2.0/3.0))
T9A13 = T9A**(1.0/3.0)
T9A56 = T9A**(5.0/6.0)
f85 = 1.88e+18*T9A56/T932*exp(-35.829/T9A13-60.557/T9)
r85 = 1.88e+18*T9A56/T932*exp(-35.829/T9A13)*1.72e-01
O16(A,G)NE20
Q86 = 4.734
f86 = 9.37e+09/T923*exp(-39.757/T913-(T9/1.586)**2)+6.21e+01/T932
| *exp(-10.297/T9)+5.38e+02/T932*exp(-12.226/T9)+1.30e+01
| *T9**2*exp(-20.093/T9)
r86 = f86*5.65e+10*T932*exp(-54.937/T9)
O17(P,G)F18
Q87 = 5.607
T9A = T9/(1.0+2.69*T9)
T9A13 = T9A**(1.0/3.0)
T9A56 = T9A**(5.0/6.0)
f87 = 7.97e+07*T9A56/T932*exp(-16.712/T9A13)+1.51e+08/T923
| *exp(-16.712/T913)*(1.0+0.025*T913-0.051*T923-8.82e-03*T9)
| +1.56e+05/T9*exp(-6.272/T9)+rc87*1.31e+01/T932
| *exp(-1.961/T9)
r87 = f87*3.66e+10*T932*exp(-65.061/T9)
O17(P,A)N14
Q88 = 1.191
rc88=rc87
f88 = 1.53e+07/T923*exp(-16.712/T913-(T9/0.565)**2)*(1.0+0.025
| *T913+5.39*T923+0.940*T9+13.5*T943+5.98*T953)+2.92e+06
| *T9*exp(-4.247/T9)+rc88*(4.81e+10*T9*exp(-16.712/T913-(T9
| /0.040)**2)+5.05e-05/T932*exp(-0.723/T9))+rc88*1.31e+01
| /T932*exp(-1.961/T9)
r88 = f88*6.76e-01*exp(-13.825/T9)
O17(A,G)NE21
Q89 = 7.351
T9A = T9/(1.0+0.1646*T9)
T9A13 = T9A**(1.0/3.0)
T9A56 = T9A**(5.0/6.0)
GT9 = 1.0+exp(-10.106/T9)/3.0
FT9A = exp(-(0.786/T9A)**3.51)
FPT9A = exp(-(T9A/1.084)**1.69)
f89 = 1.73e+17*FPT9A/GT9*T9A56/T932*exp(-39.914/T9A13)+3.50e+15
| *FT9A/GT9*T9A56/T932*exp(-39.914/T9A13)
r89 = f89*8.63e+10*T932*exp(-85.305/T9)
O17(A,N)NE20
Q90 = 0.590
T9A = T9/(1.0+0.0268*T9+0.0232*T953
| /(1.0+0.0268*T9)**(2.0/3.0))
T9A13 = T9A**(1.0/3.0)
T9A56 = T9A**(5.0/6.0)
GT9 = 1.0+exp(-10.106/T9)/3.0
f90 = 1.03e+18/GT9*T9A56/T932*exp(-39.914/T9A13)
r90 = f90*1.86e+01*exp(-6.852/T9)
O18(P,G)F19
Q91 = 7.994
f91 = 3.45e+08/T923*exp(-16.729/T913-(T9/0.139)**2)*(1.0+0.025
| *T913+2.26*T923+0.394*T9+30.56*T943+13.55*T953)+1.25e-15
| /T932*exp(-0.231/T9)+1.64e+02/T932*exp(-1.670/T9)+1.28e+04
| *T912*exp(-5.098/T9)
r91 = f91*9.20e+09*T932*exp(-92.769/T9)
O18(P,A)N15
Q92 = 3.980
f92 = 3.63e+11/T923*exp(-16.729/T913-(T9/1.361)**2)*(1.0+0.025
| *T913+1.88*T923+0.327*T9+4.66*T943+2.06*T953)+9.90e-14
| /T932*exp(-0.231/T9)+2.66e+04/T932*exp(-1.670/T9)+2.41e+09
| /T932*exp(-7.638/T9)+1.46e+09/T9*exp(-8.310/T9)
r92 = f92*1.66e-01*exp(-46.191/T9)
O18(A,G)NE22
Q93 = 9.669
f93 = 1.82e+12/T923*exp(-40.057/T913-(T9/0.343)**2)*(1.00+0.010
| *T913+0.988*T923+0.072*T9+3.17*T943+0.586*T953)+7.54e+00
| /T932*exp(-6.228/T9)+3.48e+01/T932*exp(-7.301/T9)+6.23e+03
| *T9*exp(-16.987/T9)+rc93*1.0e-11/T932*exp(-1.994/T9)
r93 = f93*5.85e+10*T932*exp(-112.208/T9)
O18(A,N)NE21
Q94 = -0.693
T9A = T9/(1.0+0.0483*T9+0.00569*T953
| /(1.0+0.0483*T9)**(2.0/3.0))
T9A13 = T9A**(1.0/3.0)
T9A56 = T9A**(5.0/6.0)
GT9 = 1.0+5.0*exp(-23.002/T9)
FT9A = exp(-(0.431/T9A)**3.89)
T9min94 = 0.04
f94 = 7.22e+17*FT9A/GT9*T9A56/T932*exp(-40.056/T9A13)+150.31
| /GT9*exp(-8.045/T9)
r94 = 7.84e-01*7.22e+17/GT9*T9A56/T932
| *exp(-(0.431/T9A)**3.89 + -40.056/T9A13 + 8.045/T9)
| +150.31/GT9*7.84e-01
F19(P,G)NE20
Q95 = 12.848
f95 = (6.04e+07/T923*exp(-18.113/T913-(T9/0.416)**2)*(1.0+0.023
| *T913+2.06*T923+0.332*T9+3.16*T943+1.30*T953)+6.32e+02
| /T932*exp(-3.752/T9)+7.56e+04/T927*exp(-5.722/T9))
| /(1.0+4.0*exp(-2.090/T9)+7.00*exp(-16.440/T9))
r95 = f95*3.70e+10*T932*exp(-149.093/T9)
F19(P,N)NE19
Q96 = -4.021
f96 = 1.27e+08*(1.0-0.147*T912+0.069*T9)*exp(-46.659/T9)
r96 = 1.27e+08*(1.0-0.147*T912+0.069*T9)*9.98e-01
F19(P,A)O16
Q97 = 8.114
f97 = (3.55e+11/T923*exp(-18.113/T913-(T9/0.845)**2)*(1.0+0.023
| *T913+1.96*T923+0.316*T9+2.86*T943+1.17*T953)+3.67e+06
| /T932*exp(-3.752/T9)+3.07e+08*exp(-6.019/T9))/(1.0+4.0
| *exp(-2.090/T9)+7.0*exp(-16.440/T9))
r97 = f97*6.54e-01*exp(-94.156/T9)
F19(A,P)NE22
Q98 = 1.675
f98 = 4.50e+18/T923*exp(-43.467/T913-(T9/0.637)**2)+7.98e+04*T932
| *exp(-12.760/T9)
r98 = f98*6.36e+00*exp(-19.439/T9)
NE19(P,G)NA20
Q99 = 2.199
f99 = 1.71e+06/T923*exp(-19.431/T913)*(1.0+0.021*T913+0.130*T923
| +1.95e-02*T9+3.86e-02*T943+1.47e-02*T953)+1.89e+05/T923
| *exp(-19.431/T913-(T9/1.142)**2)*(1.0+0.021*T913+2.13*T923
| +0.320*T9+2.80*T943+1.07*T953)+8.45e+03/T954*exp(-7.64/T9)
r99 = f99*7.39e+09*T932*exp(-25.519/T9)
NE20(P,G)NA21
Q100 = 2.431
f100 = 9.55e+06*exp(-19.447/T913)/(T9**2*(1.0+0.0127/T923)**2)
| +2.05e+08/T923*exp(-19.447/T913)*(1.0+2.67
| *exp(-(T9/0.210)**(0.5)))+1.80e+01/T932*exp(-4.242/T9)+1.02e+01
| /T932*exp(-4.607/T9)+3.60e+04/T914*exp(-11.249/T9)
r100 = f100*4.63e+09*T932*exp(-28.216/T9)
NE20(P,A)F17
Q101 = -4.134
T9A = T9/(1.00+6.12e-02*T9+1.30e-02*T953/(1.0+6.12e-02*T9)
| **(2.0/3.0))
T9A13 = T9A**(1.0/3.0)
T9A56 = T9A**(5.0/6.0)
f101 = 3.25e+19*(5.31+0.544*T9-0.0523*T9**2)*T9A56/T932
| *exp(-43.176/T9A13-47.969/T9)
r101 = 3.25e+19*(5.31+0.544*T9-0.0523*T9**2)*T9A56/T932
| *exp(-43.176/T9A13)*5.37e-02
NE20(A,G)MG24
Q102 = 9.312
f102 = (4.11e+11/T923*exp(-46.766/T913-(T9/2.219)**2)*(1.00+0.009
| *T913+0.882*T923+0.055*T9+0.749*T943+0.119*T953)+5.27e+03
| /T932*exp(-15.869/T9)+6.51e+03*T912*exp(-16.223/T9)+rc102
| *4.21e+01/T932*exp(-9.115/T9)+3.20e+01/T923*exp(-9.383/T9))
| /(1.0+5.0*exp(-18.960/T9))
r102 = f102*6.01e+10*T932*exp(-108.059/T9)
NE21(P,G)NA22
Q103 = 6.738
f103 = 4.37e+08/T923*exp(-19.462/T913)+5.85e+00/T932*exp(-1.399/T9)
| +1.29e+04/T932*exp(-3.009/T9)+3.15e+05/T935*exp(-5.763/T9)
| +rc103*2.95e+08/T923*exp(-19.462/T913-(T9/0.058)**2)*(1.00
| +0.021*T913+13.29*T923+1.99*T9+124.1*T943+47.29*T953)
| +7.80e-01/T932*exp(-1.085/T9)
r103 = f103*1.06e+10*T932*exp(-78.194/T9)
NE21(A,G)MG25
Q104 = 9.882
T9A = T9/(1.0+0.0537*T9)
T9A13 = T9A**(1.0/3.0)
T9A56 = T9A**(5.0/6.0)
f104 = (4.94e+19*T9A56/T932*exp(-46.890/T9A13)+2.66e+07/T932
| *exp(-22.049/T9))*(1.52e-04*exp(-46.90/T913*(8.72e-03*T9
| -6.87e-04*T9**2+2.15e-05*T9**3))/(1.0+1.5*exp(-4.068/T9)
| +2.0*exp(-20.258/T9)))
r104 = f104*4.06e+10*T932*exp(-114.676/T9)
NE21(A,N)MG24
Q105 = 2.551
T9A = T9/(1.0+0.0537*T9)
T9A13 = T9A**(1.0/3.0)
T9A56 = T9A**(5.0/6.0)
f105 = (4.94e+19*T9A56/T932*exp(-46.890/T9A13)+2.66e+07/T932
| *exp(-22.049/T9))/(1.0+1.50*exp(-4.068/T9)+2.0
| *exp(-20.258/T9))
r105 = f105*1.29e+01*exp(-29.606/T9)
NE22(P,G)NA23
Q106 = 8.794
f106 = 1.15e+09/T923*exp(-19.475/T913)+9.77e-12/T932
| *exp(-0.348/T9)+8.96e+03/T932*exp(-4.840/T9)+6.52e+04/T932
| *exp(-5.319/T9)+7.97e+05/T912*exp(-7.418/T9)+rc106*1.63e-01
| /T932*exp(-1.775/T9)
r106 = f106*4.67e+09*T932*exp(-102.048/T9)
NE22(A,G)MG26
Q107 = 10.612
T9A = T9/(1.0+0.0548*T9)
T9A13 = T9A**(1.0/3.0)
T9A56 = T9A**(5.0/6.0)
GT9 = 1.0+5.0*exp(-14.791/T9)
FT9A = exp(-(0.197/T9A)**4.82)
FPT9A = exp(-(T9A/0.249)**2.31)
f107 = 4.16e+19*FPT9A/GT9*T9A56/T932*exp(-47.004/T9A13)
| +2.08e+16*FT9A/GT9*T9A56/T932*exp(-47.004/T9A13)
r107 = f107*6.15e+10*T932*exp(-123.151/T9)
NE22(A,N)MG25
Q108 = -0.481
T9A = T9/(1.0+0.0548*T9)
T9A13 = T9A**(1.0/3.0)
T9A56 = T9A**(5.0/6.0)
GT9 = 1.0+5.0*exp(-14.791/T9)
FT9A = exp(-(0.197/T9A)**4.82)
T9min108 = 0.025
f108 = 4.16e+19*FT9A/GT9*T9A56/T932*exp(-47.004/T9A13)
| +1.44e-04/GT9*exp(-5.577/T9)
r108 = 5.44e-01*4.16e+19/GT9*T9A56/T932
| *exp(-(0.197/T9A)**4.82 - 47.004/T9A13 + 5.577/T9)
| +1.44e-04/GT9*5.44e-01
NA21(P,G)MG22
Q109 = 5.497
f109 = 1.41e+05/T923*exp(-20.739/T913-(T9/0.366)**2)*(1.0+0.020
| *T913+4.741*T923+0.667*T9+16.380*T943+5.858*T953)
| +6.72e+02/T934*exp(-2.436/T9)
r109 = f109*7.44e+10*T932*exp(-63.790/T9)
NA22(N,P)NE22
Q110 = 3.624
f110 = 1.24e+08*exp(1.0-3.037e-02*T9+8.380e-03*T9**2-7.101e-04
| *T9**3)
r110 = f110*7.01*exp(-42.059/T9)
NA22(N,A)F19
Q111 = 1.949
f111 = 1.21e+06*exp(1.00+8.955e-01*T9-5.645e-02*T9**2+7.302e-04
| *T9**3)
r111 = f111*1.10*exp(-22.620/T9)
NA22(P,G)MG23
Q112 = 7.578
f112 = 9.63e-05*T932*exp(-0.517/T9)+2.51e+04*T9*exp(-2.013/T9)
r112 = f112*3.27e+10*T932*exp(-87.933/T9)
NA23(P,G)MG24
Q113 = 11.691
f113 = (2.93e+08/T923*exp(-20.766/T913-(T9/0.297)**2)*(1.0+0.020
| *T913+1.61*T923+0.226*T9+4.94*T943+1.76*T953)+9.34e+01
| /T932*exp(-2.789/T9)+1.89e+04/T932*exp(-3.434/T9)+5.10e+04
| *T915*exp(-5.510/T9))/(1.0+1.5*exp(-5.105/T9))
r113 = f113*7.49e+10*T932*exp(-135.665/T9)
NA23(P,N)MG23
Q114 = -4.841
T9A = T9/(1.0+0.141*T9)
T9A32 = T9A**(3.0/2.0)
f114 = 9.29e+08*(1.0-0.881*T9A32/T932)*exp(-56.173/T9)
r114 = 9.29e+08*(1.0-0.881*T9A32/T932)*9.98e-01
NA23(P,A)NE20
Q115 = 2.379
f115= 8.56e+09/T923*exp(-20.766/T913-(T9/0.131)**2)*(1.0+0.020
| *T913+8.21*T923+1.15*T9+44.36*T943+15.84*T953)+4.02e+00
| /T932*exp(-1.990/T9)+1.18e+04/T954*exp(-3.148/T9)+8.59e+05
| *T943*exp(-4.375/T9)+rc115*3.06e-12/T932*exp(-0.447/T9)
r115 = f115*1.25e+00*exp(-27.606/T9)
NA23(A,N)AL26T
Q116 = -2.968
GT9 = 1.0+exp(-4.612/T9-5.623e-04+7.460e-02*T9)
GPT9 = 1.0+exp(-3.573/T9-1.008+0.1357*T9)
f116 = 3.38e+06*exp(0.388*T9+9.08e-03*T9**2-2.07e-03*T9**3)
| *GPT9/GT9*8.36e-01*exp(-34.442/T9)
Note: r116=f141; See AL26T(N,A)NA23
NA23(A,N)AL26M
Q117 = -3.196
f117 = 4.18e+06+5.43e+07/T932*exp(-0.9653/T9)+6.97e+07
| *T927*exp(-1.494/T9)*7.60e-02*exp(-37.093/T9)
Note: r117=f142; See AL26M(N,A)NA23
NA23(A,N)AL26G
Q118 = -2.968
GT9 = 1.0+exp(-4.612/T9-5.623e-04+7.460e-02*T9)
GPT9 = 1.0+exp(-3.573/T9-1.008+0.1357*T9)
f118 = 3.38e+06*exp(0.388*T9+9.08e-03*T9**2-2.07e-03*T9**3)
| *GPT9/GT9-4.18e+06+5.43e+07/T932*exp(-0.9653/T9)+6.97e+07
| *T927*exp(-1.494/T9)*9.09e-02*exp(-2.651/T9)*8.36e-01
| *exp(-34.442/T9)
Note: r118=f143; See AL27G(N,A)NA23
MG24(P,G)AL25
Q119 = 2.271
f119 = (5.60e+08/T923*exp(-22.019/T913)*(1.0+0.019*T913-0.173
| *T923-0.023*T9)+1.48e+03/T932*exp(-2.484/T9)+4.00e+03
| *exp(-4.180/T9))/(1.0+5.0*exp(-15.882/T9))
r119 = f119*3.13e+09*T932*exp(-26.358/T9)
MG24(P,A)NA21
Q120 = -6.880
T9A = T9/(1.0+0.127*T9)
T9A13 = T9A**(1.0/3.0)
T9A56 = T9A**(5.0/6.0)
f120 = 1.81e+21*(4.43+3.31*T9-0.229*T9**2)*T9A56/T932
| *exp(-49.967/T9A13-79.843/T9)
r120 = 1.81e+21*(4.43+3.31*T9-0.229*T9**2)*T9A56/T932
| *exp(-49.967/T9A13)*7.71e-02
MG24(A,G)SI28
Q121 = 9.984
f121 = (4.78e+01/T932*exp(-13.506/T9)+2.38e+03/T932
| *exp(-15.218/T9)+2.47e+02*T932*exp(-15.147/T9)+rc121
| *1.72e-09/T932*exp(-5.028/T9)+1.25e-03/T932*exp(-7.929/T9)
| +2.43e+01/T9*exp(-11.523/T9))/(1.0+5.0*exp(-15.882/T9))
r121 = f121*6.27e+10*T932*exp(-115.862/T9)
MG25(P,G)AL26T
Q122 = 6.306
f122 = 3.57e+09/T923*exp(-22.031/T913-(T9/0.06)**2)*(1.0+0.019
| *T913+7.669*T923+1.015*T9+167.4*T943+56.35*T953)+3.07e-13
| /T932*exp(-0.435/T9)+1.94e-07/T932*exp(-0.673/T9)+3.15e-05
| /T9**(3.40)*exp(-1.342/T9-(T9/13.0)**2)+1.77e+04*T958
| *exp(-3.049/T9-(T9/13.0)**2)
r122 = f122*1.03e+10*T932*exp(-73.183/T9)
MG25(P,G)AL26M
Q123 = 6.078
f123 = 0.20*(3.57e+09/T923*exp(-22.031/T913-(T9/0.06)**2)*(1.0
| +0.019*T913+7.669*T923+1.015*T9+167.4*T943+56.35*T953)
| +3.07e-13/T932*exp(-0.435/T9)+1.94e-07/T932*exp(-0.673/T9)
| +3.15e-05/T9**(3.40)*exp(-1.342/T9-(T9/13.0)**2)+1.77e+04
| *T958*exp(-3.049/T9-(T9/13.0)**2))
r123 = f123*1.13e+11*T932*exp(-70.532/T9)
MG25(P,G)AL26G
Q124 = 6.306
f124 = 0.80*(3.57e+09/T923*exp(-22.031/T913-(T9/0.06)**2)*(1.0
| +0.019*T913+7.669*T923+1.015*T9+167.4*T943+56.35*T953)
| +3.07e-13/T932*exp(-0.435/T9)+1.94e-07/T932*exp(-0.673/T9)
| +3.15e-05/T9**(3.40)*exp(-1.342/T9-(T9/13.0)**2)+1.77e+04
| *T958*exp(-3.049/T9-(T9/13.0)**2))
r124 = f124*1.03e+10*T932*exp(-73.183/T9)
MG25(A,G)SI29
Q125 = 11.127
T9A = T9/(1.0+0.0630*T9)
T9A13 = T9A**(1.0/3.0)
T9A56 = T9A**(5.0/6.0)
GT9=1.0+10.0*exp(-13.180/T9)/3.0
f125 = 3.59e+20/GT9*T9A56/T932*exp(-53.410/T9A13)*5.87e-04
| *exp(-53.42/T913*(0.0156*T9-1.79e-03*T9**2+9.08e-05*T9**3))
r125 = f125*1.90e+11*T932*exp(-129.128/T9)
MG25(A,N)SI28
Q126 = 2.653
T9A = T9/(1.0+0.063*T9)
T9A13 = T9A**(1.0/3.0)
T9A56 = T9A**(5.0/6.0)
GT9=1.0+10.0*exp(-13.180/T9)/3.0
f126 = 3.59e+20/GT9*T9A56/T932*exp(-53.410/T9A13)
r126 = f126*2.00e+01*exp(-30.792/T9)
MG25(A,P)AL28
Q127 = -1.206
f127 = 3.23e+08/T923*exp(-23.271/T913-13.995/T9+6.46*T9-2.39*T9**2
| +0.506*T9**3-6.04e-02*T9**4+3.75e-03*T9**5-9.38e-05*T9**6)
r127 = 3.23e+08/T923*exp(-23.271/T913+6.46*T9-2.39*T9**2
| +0.506*T9**3-6.04e-02*T9**4+3.75e-03
| *T9**5-9.38e-05*T9**6)*2.86e+00
MG26(P,G)AL27
Q128 = 8.272
f128 = (7.39e+08/T923*exp(-22.042/T913-(T9/0.299)**2)*(1.0+0.019
| *T913+3.61*T923+0.478*T9+9.78*T943+3.29*T953)+1.32e-10
| /T932*exp(-0.603/T9)+2.90e-05/T932*exp(-1.056/T9)+6.45e-05
| /T932*exp(-1.230/T9)+5.64e-02/T932*exp(-1.694/T9)+2.86e+03
| /T932*exp(-3.265/T9)+7.99e+04/T932*exp(-3.781/T9)+4.23e+04
| *T912*exp(-3.661/T9))/(1.0+5.0*exp(-20.990/T9))
r128 = f128*3.14e+09*T932*exp(-95.990/T9)
MG26(P,N)AL26T
Q129 = -4.786
GT9 = 1.0+exp(-19.30/T9+0.6642+0.1386*T9)
GPT9 = 1.0+exp(-3.573/T9-1.008+0.1357*T9)
f129 = 3.09e+07*exp(0.0731*T9+0.0381*T9**2-3.22e-03*T9**3)
| *GPT9/GT9*11.0*exp(-55.544/T9)
Note: r129=f138; See AL26T(N,P)MG26
MG26(P,N)AL26M
Q130 = -5.015
rc130 = rc139
f130 = 1.84e+05/T932*exp(-0.043/T9)+2.28e+07/T932*exp(-0.348/T9)
| +6.54e+08/T938*exp(-0.826/T9)+rc130*1.60e+04/T932
| *exp(-58.195/T9)
Note: r130=f139; See AL26M(N,P)MG26
MG26(P,N)AL26G
Q131 = -4.786
GT9 = 1.0+exp(-19.30/T9+0.6642+0.1386*T9)
GPT9 = 1.0+exp(-3.573/T9-1.008+0.1357*T9)
rc131 = rc140
f131 = 3.09e+07*exp(0.0731*T9+0.0381*T9**2-3.22e-03*T9**3)
| *GPT9/GT9-1.84e+05/T932*exp(-0.043/T9)+2.28e+07/T932
| *exp(-0.348/T9)+6.54e+08/T938*exp(-0.826/T9)+rc131*1.60e+04
| /T932*9.09e-02*exp(-2.651/T9)*11.0*exp(-55.544/T9)
Note: r131=f140; See AL26G(N,P)MG26
MG26(A,G)SI30
Q132 = 10.644
T9A = T9/(1.0+0.0628*T9)
T9A13 = T9A**(1.0/3.0)
T9A56 = T9A**(5.0/6.0)
GT9 = 1.0+5.0*exp(-20.990/T9)
f132 = 2.93e+20/GT9*T9A56/T932*exp(-53.505/T9A13)*4.55e-02
| *exp(-53.51/T913*(0.0751*T9-0.0105*T9**2+5.57e-04*T9**3))
r132 = f132*6.38e+10*T932*exp(-123.520/T9)
MG26(A,N)SI29
Q133 = 0.035
T9A = T9/(1.0+0.0628*T9)
T9A13 = T9A**(1.0/3.0)
T9A56 = T9A**(5.0/6.0)
GT9 = 1.0+5.0*exp(-20.990/T9)
f133 = 2.93e+20/GT9*T9A56/T932*exp(-53.505/T9A13)
r133 = f133*1.68e+00*exp(-0.401/T9)
MG26(A,P)AL29
Q134 = -2.863
f134 = 9.72e+10/T923*exp(-23.280/T913-33.225/T9)+2.20e+07/T9**4
| *exp(-41.340/T9)+2.50e+05*T974*exp(-38.825/T9)
r134 = 9.72e+10/T923*exp(-23.280/T913)*5.60e-01
| +2.20e+07/T9**4*5.60e-01*exp(-8.115/T9)
| +2.50e+05*T974*5.60e-01*exp(-5.60/T9)
AL26T(G,P)MG25
Q135 = -6.306
f135 = 3.57e+09/T923*exp(-22.031/T913-(T9/0.06)**2)*(1.00+0.019
| *T913+7.669*T923+1.015*T9+167.4*T943+56.35*T953)
| +3.07e-013/T932*exp(-0.435/T9)+1.94e-07/T932*exp(-0.673/T9)
| +3.15e-05/T9**(3.40)*exp(-1.342/T9-(T9/13.0)**2)+1.77e+04
| *T958*exp(-3.049/T9-(T9/13.0)**2)*1.03e+10*T932
| *exp(-73.183/T9)
r135 = f122
AL26M(G,P)MG25
Q136 = -6.078
f136 = 0.20*(3.57e+09/T923*exp(-22.031/T913-(T9/0.06)**2)*(1.0
| +0.019*T913+7.669*T923+1.015*T9+167.4*T943+56.35*T953)
| +3.07e-13/T932*exp(-0.435/T9)+1.94e-07/T932*exp(-0.673/T9)
| +3.15e-05/T9**(3.40)*exp(-1.342/T9-(T9/13.0)**2)+1.77e+04
| *T958*exp(-3.049/T9-(T9/13.0)**2))*1.13e+11*T932
| *exp(-70.532/T9)
r136 = f123
AL26G(G,P)MG25
Q137 = -6.306
f137 = 0.80*(3.57e+09/T923*exp(-22.031/T913-(T9/0.06)**2)*(1.0
| +0.019*T913+7.669*T923+1.015*T9+167.4*T943+56.35
| *T953)+3.07e-13/T932*exp(-0.435/T9)+1.94e-07/T932
| *exp(-0.673/T9)+3.15e-05/T9**(3.40)*exp(-1.342/T9
| -(T9/13.0)**2)+1.77e+04*T958*exp(-3.049/T9-(T9/13.0)**2))
| *1.03e+10*T932*exp(-73.183/T9)
r137 = f124
AL26T(N,P)MG26
Q138 = 4.786
f138 = 3.09e+07*exp(0.0731*T9+0.0381*T9**2-3.22e-03*T9**3)
r138=f129
r129=f138
AL26M(N,P)MG26
Q139 = 5.015
f139 = 1.84e+05/T932*exp(-0.043/T9)+2.28e+07/T932*exp(-0.348/T9)
| +6.54e+08/T938*exp(-0.826/T9)+rc139*1.60e+04/T932
r139=f130
r130=f139
AL26G(N,P)MG26
Q140 = 4.786
GT9 = 1.0+exp(-19.30/T9+0.6642+0.1386*T9)
GPT9 = 1.0+exp(-3.573/T9-1.008+0.1357*T9)
rc140 = rc139
f140 = 3.09e+07*exp(0.0731*T9+0.0381*T9**2-3.22e-03*T9**3)
| *GPT9/GT9-1.84e+05/T932*exp(-0.043/T9)+2.28e+07/T932
| *exp(-0.348/T9)+6.54e+08/T938*exp(-0.826/T9)+rc140*1.60e+04
| /T932*9.09e-02*exp(-2.651/T9)
r140=f131
r131=f140
AL26T(N,A)NA23
Q141 = 2.968
f141 = 3.38e+06*exp(0.388*T9+9.08e-03*T9**2-2.07e-03*T9**3)
r141=f116
r116=f141
AL26M(N,A)NA23
Q142 = 3.196
f142 = 4.18e+06+5.43e+07/T932*exp(-0.9653/T9)+6.97e+07
| *T927*exp(-1.494/T9)
r142=f117
r117=f142
AL26G(N,A)NA23
Q143 = 2.968
GT9 = 1.0+exp(-4.612/T9-5.623e-04+7.460e-02*T9)
GPT9 = 1.0+exp(-3.573/T9-1.008+0.1357*T9)
f143 = 3.38e+06*exp(0.388*T9+9.08e-03*T9**2-2.07e-03*T9**3)
| *GPT9/GT9-4.18e+06+5.43e+07/T932*exp(-0.9653/T9)
| +6.97e+07*T927*exp(-1.494/T9)*9.09e-02*exp(-2.651/T9)
r143=f118
r118=f143
AL26T(P,G)SI27
Q144 = 7.464
f144 = 6.78e+13/T923*exp((-23.261/T913)*(1.0+2.004e-01*T9
| -1.538e-02*T9**2+5.723e-04*T9**3))+6.13e+02/T932
| *exp(-3.202/T9)+9.45e+03/T9*exp(-4.008/T9)
r144 = f144*3.46e+10*T932*exp(-86.621/T9)
AL26M(P,G)SI27
Q145 = 7.693
f145 = 1.36e+14/T923*exp((-23.261/T913)*(1.0+2.004e-01*T9
| -1.538e-02*T9**2+5.723e-04*T9**3))
r145 = f145*3.46e+10/11.0*T932*exp(-89.272/T9)
AL26G(P,G)SI27
Q146 = 7.464
f146 = 6.78e+13/T923*exp((-23.261/T913)*(1.0+2.004e-01*T9
| -1.538e-02*T9**2+5.723e-04*T9**3))+6.13e+02/T932
| *exp(-3.202/T9)+9.45e+03/T9*exp(-4.008/T9)-1.0/11.0
| *exp(-2.651/T9)*1.36e+14/T923*exp((-23.261/T913)*(1.0
| +2.004e-01*T9-1.538e-02*T9**2+5.723e-04*T9**3))
r146 = f146*3.46e+10*T932*exp(-86.621/T9)
AL27(P,G)SI28
Q147 = 11.585
f147 = (1.67e+08/T923*exp(-23.261/T913-(T9/0.155)**2)*(1.00+0.018
| *T913+5.81*T923+0.728*T9+27.31*T943+8.71*T953)+2.20e+00
| /T932*exp(-2.269/T9)+1.22e+01/T932*exp(-2.491/T9)+1.50e+04*T9
| *exp(-4.112/T9)+rc147*6.50e-10/T932*exp(-0.853/T9)+1.63e-10
| /T932*exp(-1.001/T9))/(1.0+exp(-9.792/T9)/3.0+2.0
| *exp(-11.773/T9)/3.0)
r147 = f147*1.13e+11*T932*exp(-134.434/T9)
AL27(P,A)MG24
Q148 = 1.600
f148 = (1.10e+08/T923*exp(-23.261/T913-(T9/0.157)**2)*(1.00+0.018
| *T913+12.85*T923+1.61*T9+89.87*T943+28.66*T953)+1.29e+02
| /T932*exp(-2.517/T9)+5.66e+03*T972*exp(-3.421/T9)+rc148
| *3.89e-08/T932*exp(-0.853/T9)+8.18e-09/T932*exp(-1.001/T9))
| /(1.00+exp(-9.792/T9)/3.0+2.0*exp(-11.773/T9)/3.0)
r148 = f148*1.81e+00*exp(-18.572/T9)
AL27(A,N)P30
Q149 = -2.636
f149 = 8.20e+04*exp(-30.588/T9)+5.21e+05*T974*exp(-33.554/T9)
r149 = 8.20e+04*6.75e+00+5.21e+05*T974*exp(-2.966/T9)*6.75e+00
SI27(G,P)AL26T
Q150 = -7.464
f150 = 6.78e+13/T923*exp((-23.261/T913)*(1.00+2.004e-01*T9
| -1.538e-02*T9**2+5.723e-04*T9**3))+6.13e+02/T932
| *exp(-3.202/T9)+9.45e+03/T9*exp(-4.008/T9)*3.46e+10
| *T932*exp(-86.621/T9)
r150 = f144
SI27(G,P)AL26M
Q151 = -7.693
f151 = 1.36e+14/T923*exp((-23.261/T913)*(1.00+2.004e-01*T9
| -1.538e-02*T9**2+5.723e-04*T9**3))*3.46e+10/11.0*T932
| *exp(-89.272/T9)
r151 = f145
SI27(G,P)AL26G
Q152 = -7.464
f152 = 6.78e+13/T923*exp((-23.261/T913)*(1.00+2.004e-01*T9
| -1.538e-02*T9**2+5.723e-04*T9**3))+6.13e+02/T932
| *exp(-3.202/T9)+9.45e+03/T9*exp(-4.008/T9)-1.0/11.0
| *exp(-2.651/T9)*1.36e+14/T923*exp((-23.261/T913)
| *(1.00+2.004e-01*T9-1.538e-02*T9**2+5.723 e-04*T9**3))
| *3.46e+10*T932*exp(-86.621/T9)
r152 = f146
SI27(P,G)P28
Q153 = 2.065
f153 = 1.64e+09/T923*exp(-24.439/T913)+2.00e-08/T932
| *exp(-0.928/T9)+1.95e-02/T932*exp(-1.857/T9)+3.70e+02
| /T947*exp(-3.817/T9)
r153 = f153*1.62e+10*T932*exp(-23.960/T9)
SI28(P,G)P29
Q154 = 2.747
f154 = 1.64e+08/T923*exp(-24.449/T913-(T9/2.91)**2)*(1.0+0.017
| *T913-4.11*T923-0.491*T9+5.22*T943+1.58*T953)+3.52e+02
| /T932*exp(-4.152/T9)+6.30e+05/T932*exp(-18.505/T9)+1.69e+03
| *exp(-14.518/T9)
r154 = f154*9.46e+09*T932*exp(-31.879/T9)
SI29(P,G)P30
Q155 = 5.601
f155 = 3.26e+09/T923*exp(-24.459/T913-(T9/0.256)**2)*(1.0+0.017
| *T913+4.27*T923+0.509*T9+15.40*T943+4.67*T953)
| +2.98e+03/T932*exp(-3.667/T9)+3.94e+04/T932*exp(-4.665/T9)
| +2.08e+04*T912*exp(-8.657/T9)
r155 = f155*1.26e+10*T932*exp(-65.002/T9)
SI30(P,G)P31
Q156 = 7.297
f156 = 4.25e+08/T923*exp(-24.468/T913-(T9/0.670)**2)*(1.0+0.017
| *T913+0.150*T923+0.018*T9+5.53*T943+1.68*T953)+1.86e+04
| /T932*exp(-5.601/T9)+3.15e+05/T932*exp(-6.961/T9)
| +2.75e+05/T912*exp(-10.062/T9)
r156 = f156*9.50e+09*T932*exp(-84.673/T9)
C12+C12 (MG24)
Q157 = 13.933
T9A=T9/(1.0+0.0396*T9)
T9A13=T9A**(1.0/3.0)
T9A56=T9A**(5.0/6.0)
f157 = 4.27e+26*T9A56/T932*exp(-84.165/T9A13-2.12e-03*T9**3)
C12+O16 (SI28)
Q158 = 16.755
T9A=T9/(1.0+0.055*T9)
T9A13=T9A**(1.0/3.0)
T9A56=T9A**(5.0/6.0)
T9A23=T9A**(2.0/3.0)
T9MAX158 = 0.5
f158 = 1.72e+31*T9A56/T932*exp(-106.594/T9A13)/(exp(-0.180*T9A
| **2)+1.06e-03*exp(2.562*T9A23))
O16+O16 (S32)
Q159 = 16.542
f159 = 7.10e+36/T923*exp(-135.93/T913-0.629*T923-0.445*T943
| +0.0103*T9**2)