Caughlan and Fowler 1988
Reaction Rates
Fortran Program
program rates
parameter(izmax=16,iamax=32)
real*8 temp(60),f(60,159),r(60,159)
integer iz1,ia1,ix1,iz2,ia2,ix2,sign
character *40, c(159)
real*8 ff(izmax,iamax,izmax,iamax)
real qq(izmax,iamax,izmax,iamax)
real*8 t9
| ,t912,t913,t914,t915,t916
| ,t923,t927
| ,t932,t934,t935,t938
| ,t943,t945,t947
| ,t952,t953,t954,t956,t957,t958
| ,t965
| ,t972,t973,t974,t976
| ,t995
| ,t9f56
| ,t9m1,t9m12,t9m13,t9m14,t9m15,t9m16
| ,t9m23,t9m27
| ,t9m32,t9m34,t9m35,t9m38
| ,t9m43,t9m45,t9m47
| ,t9m52,t9m53,t9m54,t9m56,t9m57,t9m58
| ,t9m65
| ,t9m72,t9m73,t9m74,t9m76
| ,t9m95
| ,t9a,t9a13,t9a23,t9a32,t9a56
| ,gt9,gpt9
| ,ft9a,fpt9a
real*8 onethrd,twothrd,threetwo,forthrd,fivthrd,fivsix
real*8 rc19, rc28, rc78, rc87, rc88, rc93,rc102,rc103,
| rc106, rc115,
| rc121, rc130, rc131, rc139, rc140, rc147, rc148
real*8 T9max,T9max2,T9max3,T9max13,T9max14,
|T9max158,T9min26,T9max41
data temp/
| 0.001,0.002,0.003,0.004,0.005,0.006,0.007,0.008,0.009,0.010,
| 0.011,0.012,0.013,0.014,0.015,0.016,0.018,0.020,0.025,0.030,
| 0.040,0.050,0.060,0.070,0.080,0.090,0.100,0.110,0.120,0.130,
| 0.140,0.150,0.160,0.180,0.200,0.250,0.300,0.350,0.400,0.450,
| 0.500,0.600,0.700,0.800,0.900,1.000,1.250,1.500,1.750,2.000,
| 2.500,3.000,3.500,4.000,5.000,6.000,7.000,8.000,9.000,10.000/
c open the output file
open(unit=20,file='num.dat',status='unknown')
write(20,*)'H1(E-,NU)N1'
write(20,*)'Q1 = -0.782'
write(20,*)'SEE AP.J.SUPPL. 48,279(1982) & AP.J.293,1(1985)'
write(20,*)
write(20,*)
write(20,*)
write(20,*)
write(20,*)
onethrd=1.0d0/3.0d0
twothrd=2.0d0/3.0d0
forthrd=4.0d0/3.0d0
fivthrd=5.0d0/3.0d0
fivsix=5.0d0/6.0d0
threetwo=3.0d0/2.0d0
do 88 j=1,60
t9=temp(j)
c print*,onethrd,twothrd,threetwo,forthrd,fivthrd,fivsix
t912 = t9**(1.0d0/2.0d0) !t9**1/2
t913 = t9**onethrd !t9**1/3
t914 = t9**(0.25d0) !t9**1/4
t915 = t9**(0.2d0) !t9**1/5
t916 = t9**(1.0d0/6.0d0) !t9**1/6
t923 = t9**twothrd !t9**2/3
t927 = t9**(2.0d0/7.0d0) !t9**2/7
t932 = t9**(1.5d0) !t9**3/2
t934 = t9**(0.75d0) !t9**3/4
t935 = t9**(0.6d0) !t9** 3/5
t938 = t9**(3.0d0/8.0d0) !t9**3/8
t943 = t9**forthrd !t9**4/3
t945 = t9**(0.8d0) !t9**4/5
t947 = t9**(4.0d0/7.0d0) !t9**4/7
t952 = t9**(2.5d0) !t9**5/2
t953 = t9**fivthrd !t9**5/3
t954 = t9**(1.25d0) !t9**5/4
t956 = t9**fivsix !t9**5/6
t958 = t9**(5.0d0/8.0d0) !t9**5/8
t965 = t9**(1.2d0) !t9**6/5
t972 = t9**(3.5d0) !t9**7/2
t973 = t9**(7.0d0/3.0d0) !t9**7/3
t974 = t9**(7.0d0/4.0d0) !t9**7/4
t976 = t9**(7.0d0/6.0d0) !t9**7/6
t995 = t9**(9.0d0/5.0d0) !t9**9/5
c t9f56 = t9f**fivsix !t9f**5/6
t9m1 = 1.0d0/t9 !t9**-1
t9m12 = t9**(-0.5d0) !t9**-1/2
t9m13 = t9**(-onethrd) !t9**-1/3
t9m14 = t9**(-0.25d0) !t9**-1/4
t9m15 = t9**(-0.2d0) !t9**-1/5
t9m16 = t9**(-1.0d0/6.0d0) !t9**-1/6
t9m23 = t9**(-twothrd) !t9**-2/3
t9m27 = t9**(-2.0d0/7.0d0) !t9**-2/7
t9m32 = t9**(-1.5d0) !t9**-3/2
t9m34 = t9**(-0.75d0) !t9**-3/4
t9m35 = t9**(-0.6d0) !t9**-3/5
t9m38 = t9**(-0.375d0) !t9**3/8
t9m43 = t9**(-forthrd) !t9**-4/3
t9m45 = t9**(-0.8d0) !t9**-4/5
t9m47 = t9**(-4.0d0/7.0d0) !t9**-4/7
t9m52 = t9**(-2.5d0) !t9**-5/2
t9m53 = t9**(-fivthrd) !t9**-5/3
t9m54 = t9**(-1.25d0) !t9**-5/4
t9m56 = t9**(-5.0d0/6.0d0) !t9**-5/6
t9m58 = t9**(-0.625d0) !t9**-5/8
t9m65 = t9**(-1.2d0) !t9**-6/5
t9m72 = t9**(-3.5d0) !t9**-7/2
t9m73 = t9**(-7.0d0/3.0d0) !t9**-7/3
t9m74 = t9**(-1.75d0) !t9**-7/4
t9m76 = t9**(-7.0d0/6.0d0) !t9**-7/6
t9m95 = t9**(-1.8d0) !t9**-9/5
c Rate coefficient of jth reaction (varies from 0 to 1; initially set to
c intermediate value of 0.1).
rc19=0.10d0 ! Rate coeffiecient for HE4(NN,G)HE6
rc28=0.10d0 ! Rate coeffiecient for HE4(2A,G)C12
rc78=0.10d0 ! Rate coeffiecient for N15(P,A)C12
rc87=0.10d0 ! Rate coeffiecient for O17(P,G)F18
rc88=rc87 ! Rate coeffiecient for O17(P,A)N14
rc93=0.10d0 ! Rate coeffiecient for O18(A,G)NE22
rc102=0.10d0 ! Rate coeffiecient for NE20(A,G)MG24
rc103=0.10d0 ! Rate coeffiecient for NE21(P,G)NA22
rc106=0.10d0 ! Rate coeffiecient for NE22(P,G)NA23
rc115=0.10d0 ! Rate coeffiecient for NA23(P,A)NE20
rc121=0.10d0 ! Rate coeffiecient for MG24(A,G)SI28
rc130=0.10d0 ! Rate coeffiecient for MG26(P,N)AL26M
rc131=0.10d0 ! Rate coeffiecient for MG26(P,N)AL26G
rc139=rc130 ! Rate coeffiecient for AL26M(N,P)MG26
rc140=rc131 ! Rate coeffiecient for AL26G(P,N)MG26
rc147=0.10d0 ! Rate coeffiecient for AL27(P,G)SI28
rc148=0.10d0 ! Rate coeffiecient for AL27(P,A)MG24
c(1)='H1(E-,NU)N1'
c f1
qq( 1, 1, 0, 1) = -0.782
c SEE AP.J.SUPPL. 48,279(1982) & AP.J. 293,1(1985)
ff( 1, 1, 0, 1) = 0.00
ff( 0, 1, 1, 1) = 0.00
f(j,1)=ff( 1, 1, 0, 1)
r(j,1)=ff( 0, 1, 1, 1)
c(2)= 'H1(P,E+NU)H2'
c f2
qq( 1, 1, 1, 2) = 1.442
T9max2 = 3
if (T9.le.T9max2) then
ff( 1, 1, 1, 2)
| = 4.01e-15/T923*exp(-3.380/T913)*(1.0+0.123*T913+1.09*T923
| +0.938*T9)
else if (T9.gt.T9max2) then
ff( 1, 1, 1, 2) = 1.1581136e-15
else
endif
ff( 1, 2, 1, 1)=0.0
f(j,2)=ff( 1, 1, 1, 2)
r(j,2)=ff( 1, 2, 1, 1)
c(3)= 'H1(E-P,NU)H2'
c f3
qq( 1, 1, 1, 2) = 1.442
T9max3 = 3
if (T9.le.T9max3) then
ff( 1, 1, 1, 2) = 1.36e-20/T976*exp(-3.380/T913)
| *(1.0-0.729*T913+9.82*T923)
else if (T9.gt.T9max3) then
ff( 1, 1, 1, 2) = 7.3554548e-21
else
endif
ff( 1, 2, 1, 1)=0.00
f(j,3)=ff( 1, 1, 1, 2)
r(j,3)=ff( 1, 2, 1, 1)
c(4)= 'H2(P,G)HE3'
c f4
qq( 1, 2, 2, 3) = 5.494
ff( 1, 2, 2, 3)
| = 2.24e+03/T923*exp(-3.720/T913)*(1.0+0.112*T913+3.38*T923
| +2.65*T9)
ff( 2, 3, 1, 2) = ff( 1, 2, 2, 3)*1.63e+10*T932*exp(-63.750/T9)
f(j,4)=ff( 1, 2, 2, 3)
r(j,4)=ff( 2, 3, 1, 2)
c(5)= 'H2(P,N)2H1'
c f5
qq( 1, 2, 0, 1) = -2.225
ff( 1, 2, 0, 1)
| = 3.35e+07*exp(-3.720/T913-25.815/T9)*(1.0+0.784*T913+0.346
| *T923+0.690*T9)
ff( 0, 1, 1, 2)
| = 3.35e+07*exp(-3.720/T913)*(1.0+0.784*T913+0.346
| *T923+0.690*T9)*4.24e-10/T932
f(j,5)=ff( 1, 2, 0, 1)
r(j,5)=ff( 0, 1, 1, 2)
c(6)= 'H2(D,G)HE4'
c f6
qq( 1, 2, 2, 4) = 23.847
ff( 1, 2, 2, 4)
| = 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)
ff( 2, 4, 1, 2) = ff( 1, 2, 2, 4)*4.53e+10*T932*exp(-276.729/T9)
f(j,6)=ff( 1, 2, 2, 4)
r(j,6)=ff( 2, 4, 1, 2)
c(7)= 'H2(D,N)HE3'
c f7
qq( 1, 2, 2, 3) = 3.269
ff( 1, 2, 2, 3)
| = 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)
ff( 2, 3, 1, 2) = ff( 1, 2, 2, 3)*1.73e+00*exp(-37.935/T9)
f(j,7)=ff( 1, 2, 2, 3)
r(j,7)=ff( 2, 3, 1, 2)
c(8)= 'H2(D,P)H3'
c f8
qq( 1, 2, 1, 3) = 4.033
ff( 1, 2, 1, 3)
| = 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)
ff( 1, 3, 1, 2)= ff( 1, 2, 1, 3)*1.73e+00*exp(-46.798/T9)
f(j,8)=ff( 1, 2, 1, 3)
r(j,8)=ff( 1, 3, 1, 2)
c(9)= 'H3(P,G)HE4'
c f9
qq( 1, 3, 2, 4) = 19.814
ff( 1, 3, 2, 4)
| = 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)
ff( 2, 4, 1, 3)= ff( 1, 3, 2, 4)*2.61e+10*T932*exp(-229.932/T9)
f(j,9)=ff( 1, 3, 2, 4)
r(j,9)=ff( 2, 4, 1, 3)
c(10)= 'H3(P,N)HE3'
c f10
qq( 1, 3, 2, 3) = -0.764
ff( 1, 3, 2, 3)
| = 7.07e+08*(1.0-0.150*T912+0.098*T9)*exp(-8.863/T9)
ff( 2, 3, 1, 3) = 7.07e+08*(1.0-0.150*T912+0.098*T9)*9.98e-01
f(j,10)=ff( 1, 3, 2, 3)
r(j,10)=ff( 2, 3, 1, 3)
c(11)= 'H3(D,N)HE4'
c f11
qq( 1, 3, 2, 4) = 17.589
ff( 1, 3, 2, 4)
| = 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)
ff( 2, 4, 1, 3) = ff( 1, 3, 2, 4)*5.54e+00*exp(-204.117/T9)
f(j,11)=ff( 1, 3, 2, 4)
r(j,11)=ff( 2, 4, 1, 3)
c(12)= 'H3(T,2N)HE4'
c f12
qq( 1, 3, 2, 4) = 11.332
ff( 1, 3, 2, 4)
| = 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)
ff( 2, 4, 1, 3) = ff( 1, 3, 2, 4)*3.38e-10/T932*exp(-131.504/T9)
f(j,12)=ff( 1, 3, 2, 4)
r(j,12)=ff( 2, 4, 1, 3)
c(13)= 'HE3(E-,NU)H3'
c f13
qq( 2, 3, 1, 3) = -0.019
T9max13 = 3
if (T9.le.T9max13) then
ff( 2, 3, 1, 3)
| = 7.71e-12*T932*(1.0+6.48*T9+7.48*T9**2+2.91*T9**3)
| *exp(-0.2158/T9)
else if (T9.gt.T9max13) then
ff( 2, 3, 1, 3) = 6.2010694e-09
else
endif
ff( 1, 3, 2, 3) = 0.00
f(j,13)=ff( 2, 3, 1, 3)
r(j,13)=ff( 1, 3, 2, 3)
c(14)= 'HE3(P,E+NU)HE4'
c f14
qq( 2, 3, 2, 4) = 19.795
T9max14 = 3
if (T9.le.T9max14) then
ff( 2, 3, 2, 4) = 8.78e-13/T923*exp(-6.141/T913)
else if (T9.gt.T9max14) then
ff( 2, 3, 2, 4) = 5.9733437e-15
else
endif
ff( 2, 4, 2, 3) = 0.00
qq( 2, 3, 2, 4) = 10.155 ! EXCLUSIVE OF NU-ENERGY
f(j,14)=ff( 2, 3, 2, 4)
r(j,14)=ff( 2, 4, 2, 3)
c(15)= 'HE3(D,P)HE4'
c f15
qq( 2, 3, 2, 4) = 18.353
ff( 2, 3, 2, 4)
| = 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)
ff( 2, 4, 2, 3) = ff( 2, 3, 2, 4)*5.55e+00*exp(-212.980/T9)
f(j,15)=ff( 2, 3, 2, 4)
r(j,15)=ff( 2, 4, 2, 3)
c(16)= 'HE3(T,D)HE4'
c f16
qq( 2, 3, 2, 4) = 14.320
T9A = T9/(1.0+0.128*T9)
T9A13 = T9A**(1.0/3.0)
T9A56 = T9A**(5.0/6.0)
ff( 2, 3, 2, 4) = 5.46e+09*T9A56/T932*exp(-7.733/T9A13)
ff( 2, 4, 2, 3) = ff( 2, 3, 2, 4)*1.6e+00*exp(-166.182/T9)
f(j,16)=ff( 2, 3, 2, 4)
r(j,16)=ff( 2, 4, 2, 3)
c(17)= 'HE3(T,NP)HE4'
c f17
qq( 2, 3, 2, 4) = 12.096
T9A = T9/(1.0+0.115*T9)
T9A13 = T9A**(1.0/3.0)
T9A56 = T9A**(5.0/6.0)
ff( 2, 3, 2, 4) = 7.71e+09*T9A56/T932*exp(-7.733/T9A13)
ff( 2, 4, 2, 3) = ff( 2, 3, 2, 4)*3.39e-10/T932*exp(-140.367/T9)
f(j,17)=ff( 2, 3, 2, 4)
r(j,17)=ff( 2, 4, 2, 3)
c(18)= 'HE3(HE3,2P)HE4'
c f18
qq( 2, 3, 2, 4) = 12.860
ff( 2, 3, 2, 4)
| = 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)
ff( 2, 4, 2, 3) = ff( 2, 3, 2, 4)*3.39e-10/T932*exp(-149.230/T9)
f(j,18)=ff( 2, 3, 2, 4)
r(j,18)=ff( 2, 4, 2, 3)
c(19)= 'HE4(NN,G)HE6'
c f19
qq( 2, 4, 2, 6) = 0.975
ff( 2, 4, 2, 6)
| = rc19*4.04e-11/T9**2*(1.0+0.138*T9)*exp(-9.585/T9)
ff( 2, 6, 2, 4) = ff( 2, 4, 2, 6)*1.08e+20*T9**3*exp(-11.319/T9)
f(j,19)=ff( 2, 4, 2, 6)
r(j,19)=ff( 2, 6, 2, 4)
c(20)= 'HE4(NP,G)LI6'
c f20
qq( 2, 4, 3, 6) = 3.700
ff( 2, 4, 3, 6)
| = 4.62e-06/T9**2*(1.0+0.075*T9)*exp(-19.353/T9)
ff( 3, 6, 2, 4) = ff( 2, 4, 3, 6)*7.22e+19*T9**3*exp(-42.933/T9)
f(j,20)=ff( 2, 4, 3, 6)
r(j,20)=ff( 3, 6, 2, 4)
c(21)= 'HE4(D,G)LI6'
c f21
qq( 2, 4, 3, 6) = 1.475
ff( 2, 4, 3, 6)
| = 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)
ff( 3, 6, 2, 4) = ff( 2, 4, 3, 6)*1.53e+10*T932*exp(-17.118/T9)
f(j,21)=ff( 2, 4, 3, 6)
r(j,21)=ff( 3, 6, 2, 4)
c(22)= 'HE4(T,G)LI7'
c f22
qq( 2, 4, 3, 7) = 2.468
ff( 2, 4, 3, 7)
| = 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)
ff( 3, 7, 2, 4) = ff( 2, 4, 3, 7)*1.11e+10*T932*exp(-28.640/T9)
f(j,22)=ff( 2, 4, 3, 7)
r(j,22)=ff( 3, 7, 2, 4)
c(23)= 'HE4(T,N)LI6'
c f23
qq( 2, 4, 3, 6) = -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)
ff( 2, 4, 3, 6)
| = 1.80e+08*exp(-55.494/T9)*(1.0-0.261*T9A32/T932)+2.72e+09
| /T932*exp(-57.884/T9)
ff( 3, 6, 2, 4)
| = 1.80e+08*(1.0-0.261*T9A32/T932) *9.35e-01+2.72e+09/T932
| *exp(-2.39/T9)*9.35e-01
f(j,23)=ff( 2, 4, 3, 6)
r(j,23)=ff( 3, 6, 2, 4)
c(24)= 'HE4(HE3,G)BE7'
c f24
qq( 2, 4, 4, 7) = 1.588
T9A = T9/(1.0+4.95e-02*T9)
T9A13 = T9A**(1.0/3.0)
T9A56 = T9A**(5.0/6.0)
ff( 2, 4, 4, 7) = 5.61e+06*T9A56/T932*exp(-12.826/T9A13)
ff( 4, 7, 2, 4) = ff( 2, 4, 4, 7)*1.11e+10*T932*exp(-18.423/T9)
f(j,24)=ff( 2, 4, 4, 7)
r(j,24)=ff( 4, 7, 2, 4)
c(25)= 'HE4(AN,G)BE9'
c f25
qq( 2, 4, 4, 9) = 1.574
ff( 2, 4, 4, 9)
| = 2.59e-06/(T9**2*(1.0+0.344*T9))*exp(-1.062/T9)
ff( 4, 9, 2, 4) = ff( 2, 4, 4, 9)*5.84e+19*T9**3*exp(-18.260/T9)
f(j,25)=ff( 2, 4, 4, 9)
r(j,25)=ff( 4, 9, 2, 4)
c(26)= 'HE4(A)BE8'
c f26
qq( 2, 4, 4, 8) = -0.092
T9min26 = 0.01
ff( 2, 4, 4, 8)
| = 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)
c if (T9.ge.T9min26) then
c
c ff( 2, 4, 4, 8)
c | = 7.40e+05/T932*1.40e+10*T932+4.164e+09/T923
c | *exp(-13.490/T913-(T9/0.098)**2+1.0663/T9)
c | *(1.0+0.031*T913+8.009*T923
c | +1.732*T9+49.883*T943+27.426*T953)*1.40e+10*T932
c
c else if (T9.lt.T9min26) then
c
c ff( 4, 8, 2, 4) = 0.2458047e+38
c
c else
c endif
if(t9.ge.0.01) then
ff( 4, 8, 2, 4) = ff( 2, 4, 4, 8)*1.40e+10*T932*exp(1.0663/T9)
else
ff( 4, 8, 2, 4) = 1.20e+34
endif
f(j,26)=ff( 2, 4, 4, 8)
r(j,26)=ff( 4, 8, 2, 4)
c(27)= 'BE8(A,G)C12'
c f27
qq( 4, 8, 6,12) = 7.367
ff( 4, 8, 6,12)
| = 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)
ff( 6,12, 4, 8) = ff( 4, 8, 6,12)*4.30e+10*T932*exp(-85.486/T9)
f(j,27)=ff( 4, 8, 6,12)
r(j,27)=ff( 6,12, 4, 8)
c(28)= 'HE4(2A,G)C12'
c f28
qq( 2, 4, 6,12) = 7.275
if (T9.le.0.08d0) then !Low Temperature Rate
ff( 2, 4, 6,12)
| = 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)
else if (T9.gt.0.08d0) then !High Temperature Rate
ff( 2, 4, 6,12) = 2.79e-08/T9**3*exp(-4.4027/T9)+rc28
| *1.35e-07/T932*exp(-24.811/T9)
else
endif
c Q = 14.437 ! IF Cl2(A,G)Ol6 ALWAYS FOLLOWS
ff( 6,12, 2, 4) = ff( 2, 4, 6,12)*2.00e+20*T9**3*exp(-84.420/T9)
f(j,28)=ff( 2, 4, 6,12)
r(j,28)=ff( 6,12, 2, 4)
c(29)= 'LI6(P,G)BE7'
c f29
qq( 3, 6, 4, 7) = 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)
ff( 3, 6, 4, 7) = 6.69e+05*T9A56/T932*exp(-8.413/T9A13)
ff( 4, 7, 3, 6) = ff( 3, 6, 4, 7)*1.19e+10*T932*exp(-65.054/T9)
f(j,29)=ff( 3, 6, 4, 7)
r(j,29)=ff( 4, 7, 3, 6)
c(30)= 'LI6(P,HE3)HE4'
c f30
qq( 3, 6, 2, 3) = 4.018
ff( 3, 6, 2, 3)
| = 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)
ff( 2, 3, 3, 6)= ff( 3, 6, 2, 3)*1.07e+00*exp(-46.631/T9)
f(j,30)=ff( 3, 6, 2, 3)
r(j,30)=ff( 2, 3, 3, 6)
c(31)= 'LI6(A,G)B10'
c f31
qq( 3, 6, 5,10) = 4.460
ff( 3, 6, 5,10)
| = 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)
ff( 5,10, 3, 6) = ff( 3, 6, 5,10)*1.58e+10*T932*exp(-51.753/T9)
f(j,31)=ff( 3, 6, 5,10)
r(j,31)=ff( 5,10, 3, 6)
c(32)= 'LI7(P,N)BE7'
c f32
qq( 3, 7, 4, 7) = -1.644
ff( 3, 7, 4, 7)
| = 5.15e+09*exp(-1.167*T913-19.081/T9)+7.84e+09/T932
| *exp(-22.832/T9)
ff( 4, 7, 3, 7) = 5.15e+09*exp(-1.167*T913)+7.84e+09/T932
| *exp(-3.751/T9)*9.98e-01
f(j,32)=ff( 3, 7, 4, 7)
r(j,32)=ff( 4, 7, 3, 7)
c(33)= 'LI7(P,G)BE8'
c f33
qq( 3, 7, 4, 8) = 17.254
ff( 3, 7, 4, 8)
| = 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)
ff( 4, 8, 3, 7)= ff( 3, 7, 4, 8)*6.55e+10*T932*exp(-200.225/T9)
f(j,33)=ff( 3, 7, 4, 8)
r(j,33)=ff( 4, 8, 3, 7)
c(34)= 'LI7(P,A)HE4'
c f34
qq( 3, 7, 2, 4) = 17.346
T9A = T9/(1.0+0.759*T9)
T9A13 = T9A**(1.0/3.0)
T9A56 = T9A**(5.0/6.0)
ff( 3, 7, 2, 4)
| = 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)
ff( 2, 4, 3, 7) = ff( 3, 7, 2, 4)*4.69e+00*exp(-201.291/T9)
c ff( 3, 7, 2, 4) APPLIES TO BOTH LI7(P,A)HE4 AND LI7(P,AG)HE4
f(j,34)=ff( 3, 7, 2, 4)
r(j,34)=ff( 2, 4, 3, 7)
c(35)= 'LI7(P,AG)HE4'
c f35
qq( 3, 7, 2, 4) = 17.346
T9A = T9/(1.0+0.759*T9)
T9A13 = T9A**(1.0/3.0)
T9A56 = T9A**(5.0/6.0)
c Note: f35=f33+f34
ff( 3, 7, 2, 4)
| = (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))
cc r35 = r34
ff( 2, 4, 3, 7) = ff( 3, 7, 2, 4)*4.69e+00*exp(-201.291/T9)
f(j,35)=ff( 3, 7, 2, 4)
r(j,35)=ff( 2, 4, 3, 7)
c(36)= 'LI7(D,N)2HE4'
c f36
qq( 3, 7, 2, 4) = 15.121
ff( 3, 7, 2, 4) = 2.92e+11/T923*exp(-10.259/T913)
ff( 2, 4, 3, 7) = ff( 3, 7, 2, 4)*9.95e-10/T932*exp(-175.476/T9)
f(j,36)=ff( 3, 7, 2, 4)
r(j,36)=ff( 2, 4, 3, 7)
c(37)= 'LI7(T,2N)2HE4'
c f37
qq( 3, 7, 2, 4) = 8.864
ff( 3, 7, 2, 4) = 8.81e+11/T923*exp(-11.333/T913)
ff( 2, 4, 3, 7) = ff( 3, 7, 2, 4)*1.22e-19/T9**3*exp(-102.864/T9)
f(j,37)=ff( 3, 7, 2, 4)
r(j,37)=ff( 2, 4, 3, 7)
c(38)= 'LI7(HE3,NP)2HE4'
c f38
qq( 3, 7, 4, 8) = 9.628
ff( 3, 7, 4, 8) = 1.11e+13/T923*exp(-17.989/T913)
ff( 4, 8, 3, 7) = ff( 3, 7, 4, 8)*6.09e-20/T9**3*exp(-111.727/T9)
f(j,38)=ff( 3, 7, 4, 8)
r(j,38)=ff( 4, 8, 3, 7)
c(39)= 'LI7(A,G)B11'
c f39
qq( 3, 7, 5,11) = 8.664
ff( 3, 7, 5,11)
| = 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)
ff( 5,11, 3, 7) = ff( 3, 7, 5,11)*4.02e+10*T932*exp(-100.538/T9)
f(j,39)=ff( 3, 7, 5,11)
r(j,39)=ff( 5,11, 3, 7)
c(40)= 'LI7(A,N)B10'
c f40
qq( 3, 7, 5,10) = -2.790
ff( 3, 7, 5,10) = 3.84e+08*exp(-32.382/T9)
ff( 5,10, 3, 7) = 3.84e+08*1.32e+00
f(j,40)=ff( 3, 7, 5,10)
r(j,40)=ff( 5,10, 3, 7)
c(41)= 'BE7(E-,NU+G)LI7'
c f41
qq( 4, 7, 3, 7) = 0.862
T9max41 = 3
if (T9.le.T9max41) then
ff( 4, 7, 3, 7)
| = 1.34e-10/T912*(1.0-0.537*T913+3.86*T923+0.0027/T9
| *exp(2.515e-03/T9))
else if (T9.gt.T9max41) then
ff( 4, 7, 3, 7) = 0.00
else
endif
ff( 3, 7, 4, 7) = 0.00
f(j,41)=ff( 4, 7, 3, 7)
r(j,41)=ff( 3, 7, 4, 7)
c Q = 0.049 EXCLUSIVE OF NU ENERGY
c RATE MUST NOT EXCEED l.51e-07/(RHO*(1.0+X)/2.0)
c F0R T9 LESS THAN 0.001.
c(42)= 'BE7(P,G)B8'
c f42
qq( 4, 7, 5, 8) = 0.137
ff( 4, 7, 5, 8)
| = 3.11e+05/T923*exp(-10.262/T913)+2.53e+03/T932
| *exp(-7.306/T9)
ff( 5, 8, 4, 7) = ff( 4, 7, 5, 8)*1.30e+10*T932*exp(-1.595/T9)
f(j,42)=ff( 4, 7, 5, 8)
r(j,42)=ff( 5, 8, 4, 7)
c(43)= 'BE7(D,P)2HE4'
c f43
qq( 4, 7, 4, 8) = 16.766
ff( 4, 7, 4, 8) = 1.07e+12/T923*exp(-12.428/T913)
ff( 4, 8, 4, 7) = ff( 4, 7, 4, 8)*9.97e-10/T932*exp(-194.557/T9)
f(j,43)=ff( 4, 7, 4, 8)
r(j,43)=ff( 4, 8, 4, 7)
c(44)= 'BE7(T,NP)2HE4'
c f44
qq( 4, 7, 4, 8) = 10.508
ff( 4, 7, 4, 8) = 2.91e+12/T923*exp(-13.729/T913)
ff( 4, 8, 4, 7) = ff( 4, 7, 4, 8)*6.09e-20/T9**3*exp(-121.944/T9)
f(j,44)=ff( 4, 7, 4, 8)
r(j,44)=ff( 4, 8, 4, 7)
c(45)= 'BE7(HE3,2P)2HE4'
c f45
qq( 4, 7, 4, 8) = 11.272
ff( 4, 7, 4, 8) = 6.11e+13/T923*exp(-21.793/T913)
ff( 4, 8, 4, 7) = ff( 4, 7, 4, 8)*1.22e-19/T9**3*exp(-130.807/T9)
f(j,45)=ff( 4, 7, 4, 8)
r(j,45)=ff( 4, 8, 4, 7)
c(46)= 'BE7(A,G)C11'
c f46
qq( 4, 7, 6,11) = 7.544
ff( 4, 7, 6,11)
| = 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)
ff( 6,11, 4, 7) = ff( 4, 7, 6,11)*4.02e+10*T932*exp(-87.539/T9)
f(j,46)=ff( 4, 7, 6,11)
r(j,46)=ff( 6,11, 4, 7)
c(47)= 'BE9(P,G)B10'
c f47
qq( 4, 9, 5,10) = 6.586
ff( 4, 9, 5,10)
| = 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)
ff( 5,10, 4, 9) = ff( 4, 9, 5,10)*9.73e+09*T932*exp(-76.427/T9)
f(j,47)=ff( 4, 9, 5,10)
r(j,47)=ff( 5,10, 4, 9)
c(48)= 'BE9(P,N)B9'
c f48
qq( 4, 9, 5, 9) = -1.850
ff( 4, 9, 5, 9)
| = 5.58e+07*(1.0+0.042*T912+0.985*T9)*exp(-21.473/T9)
| +1.02e+09/T932*exp(-26.725/T9)
ff( 5, 9, 4, 9)
| = 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
f(j,48)=ff( 4, 9, 5,9)
r(j,48)=ff( 5,9, 4, 9)
c(49)= 'BE9(P,D)2HE4'
c f49
qq( 4, 9, 4, 8) = 0.651
ff( 4, 9, 4, 8)
| = 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)
ff( 4, 8, 4, 9) = ff( 4, 9, 4, 8)*8.07e-11/T932*exp(-7.555/T9)
f(j,49)= ff( 4, 9, 4, 8)
r(j,49)=ff( 4, 8, 4, 9)
c(50)= 'BE9(P,A)LI6'
c f50
qq( 4, 9, 3, 6) = 2.126
ff( 4, 9, 3, 6)
| = 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)
ff( 3, 6, 4, 9) = ff( 4, 9, 3, 6)*6.18e-01*exp(-24.674/T9)
f(j,50)=ff( 4, 9, 3, 6)
r(j,50)=ff( 3, 6, 4, 9)
c(51)= 'BE9(A,N)C12'
c f51
qq( 4, 9, 6,12) = 5.701
ff( 4, 9, 6,12)
| = 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)
ff( 6,12, 4, 9)= ff( 4, 9, 6,12)*1.03e+01*exp(-66.160/T9)
f(j,51)= ff( 4, 9, 6,12)
r(j,51)=ff( 6,12, 4, 9)
c(52)= 'B10(P,G)C11'
c f52
qq( 5,10, 6,11) = 8.690
ff( 5,10, 6,11)
| = 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)
ff( 6,11, 5,10) = ff( 5,10, 6,11)*3.03e+10*T932*exp(-100.840/T9)
f(j,52)=ff( 5,10, 6,11)
r(j,52)=ff( 6,11, 5,10)
c(53)= 'B10(P,A)BE7'
c f53
qq( 5,10, 4, 7) = 1.146
ff( 5,10, 4, 7)
| = 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)
ff( 4, 7, 5,10) = ff( 5,10, 4, 7)*7.54e-01*exp(-13.301/T9)
f(j,53)=ff( 5,10, 4, 7)
r(j,53)=ff( 4, 7, 5,10)
c(54)= 'B10(A,N)N13'
c f54
qq( 5,10, 7,13) = 1.059
ff( 5,10, 7,13) = 1.20e+13/T923*exp(-27.989/T913-(T9/9.589)**2)
ff( 7,13, 5,10) = ff( 5,10, 7,13)*9.34e+00*exp(-12.287/T9)
f(j,54)=ff( 5,10, 7,13)
r(j,54)=ff( 7,13, 5,10)
c(55)= 'B11(P,G)C12'
c f55
qq( 5,11, 6,12) = 15.957
ff( 5,11, 6,12)
| = 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)
ff( 6,12, 5,11) = ff( 5,11, 6,12)*7.01e+10*T932*exp(-185.173/T9)
f(j,55)=ff( 5,11, 6,12)
r(j,55)=ff( 6,12, 5,11)
c(56)= 'B11(P,N)C11'
c f56
qq( 5,11, 6,11) = -2.764
ff( 5,11, 6,11)
| = 1.69e+08*(1.0-0.048*T912+0.010*T9)*exp(-32.080/T9)
ff( 6,11, 5,11) = 1.69e+08*(1.0-0.048*T912+0.010*T9)*9.98e-01
f(j,56)=ff( 5,11, 6,11)
r(j,56)=ff( 6,11, 5,11)
c(57)= 'B11(P,A)2HE4'
c f57
qq( 5,11, 4,8) = 8.682
ff( 5,11, 4,8)
|= 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)
ff( 4, 8, 5,11) = ff( 5,11, 4,8)*3.50e-10/T932*exp(-100.753/T9)
f(j,57)=ff( 5,11, 4,8)
r(j,57)=ff( 4, 8, 5,11)
c(58)= 'B11(A,N)N14'
c f58
qq( 5,11, 7,14) = 0.158
ff( 5,11, 7,14)
|= 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)
ff( 7,14, 5,11) = ff( 5,11, 7,14)*3.67e+00*exp(-1.835/T9)
f(j,58)=ff( 5,11, 7,14)
r(j,58)=ff( 7,14, 5,11)
c(59)= 'B11(A,P)C14'
c f59
qq( 5,11, 6,14) = 0.784
ff( 5,11, 6,14)
|= 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)
ff( 6,14, 5,11) = ff( 5,11, 6,14)*1.10e+01*exp(-9.098/T9)
f(j,59)=ff( 5,11, 6,14)
r(j,59)=ff( 6,14, 5,11)
c(60)= 'C11(P,G)N12'
c f60
qq( 6,11, 7,12) = 0.601
ff( 6,11, 7,12)
| = 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)
ff( 7,12, 6,11) = ff( 6,11, 7,12)*2.33e+10*T932*exp(-6.975/T9)
f(j,60)=ff( 6,11, 7,12)
r(j,60)=ff( 7,12, 6,11)
c(61)= 'C12(P,G)N13'
c f61
qq( 6,12, 7,13) = 1.944
ff( 6,12, 7,13)
| = 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)
ff( 7,13, 6,12) = ff( 6,12, 7,13)*8.84e+09*T932*exp(-22.553/T9)
f(j,61)=ff( 6,12, 7,13)
r(j,61)=ff( 7,13, 6,12)
c(62)= 'C12(A,G)O16'
c f62
qq( 6,12, 8,16) = 7.162
ff( 6,12, 8,16)
| = 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)
ff( 8,16, 6,12) = ff( 6,12, 8,16)*5.13e+10*T932*exp(-83.111/T9)
f(j,62)=ff( 6,12, 8,16)
r(j,62)=ff( 8,16, 6,12)
c(63)= 'C12(A,N)O15'
c f63
qq( 6,12, 8,15) = -8.502
ff( 6,12, 8,15)
| = 2.48e+07*(1.00+0.188*T912+0.015*T9)*exp(-98.661/T9)
ff( 8,15, 6,12) = 2.48e+07*(1.00+0.188*T912+0.015*T9)*1.41e+00
f(j,63)=ff( 6,12, 8,15)
r(j,63)=ff( 8,15, 6,12)
c(64)= 'C13(P,G)N14'
c f64
qq( 6,13, 7,14) = 7.551
ff( 6,13, 7,14)
| = 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)
ff( 7,14, 6,13) = ff( 6,13, 7,14)*1.19e+10*T932*exp(-87.621/T9)
f(j,64)=ff( 6,13, 7,14)
r(j,64)=ff( 7,14, 6,13)
c(65)= 'C13(P,N)N13'
c f65
qq( 6,13, 7,13) = -3.003
ff( 6,13, 7,13)
| = 1.88e+08*(1.0-0.167*T912+0.037*T9)*exp(-34.846/T9)
ff( 7,13, 6,13) = 1.88e+08*(1.0-0.167*T912+0.037*T9)*9.98e-01
f(j,65)=ff( 6,13, 7,13)
r(j,65)=ff( 7,13, 6,13)
c(66)= 'C13(A,N)O16'
c f66
qq( 6,13, 8,16) = 2.216
ff( 6,13, 8,16)
| = 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)
ff( 8,16, 6,13) = ff( 6,13, 8,16)*5.79e+00*exp(-25.711/T9)
f(j,66)=ff( 6,13, 8,16)
r(j,66)=ff( 8,16, 6,13)
c(67)= 'C14(P,G)N15'
c f67
qq( 6,14, 7,15) = 10.207
ff( 6,14, 7,15)
| = 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)
ff( 7,15, 6,14) = ff( 6,14, 7,15)*9.00e+09*T932*exp(-118.452/T9)
f(j,67)=ff( 6,14, 7,15)
r(j,67)=ff( 7,15, 6,14)
c(68)= 'C14(P,N)N14'
c f68
qq( 6,14, 7,14) = -0.626
ff( 6,14, 7,14)
| = 7.19e+05*(1.0+0.361*T912+0.502*T9)*exp(-7.263/T9)
| +3.34e+08/T912*exp(-12.246/T9)
ff( 7,14, 6,14) = 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
f(j,68)=ff( 6,14, 7,14)
r(j,68)=ff( 7,14, 6,14)
c(69)= 'C14(A,G)O18'
c f69
qq( 6,14, 8,18) = 6.227
ff( 6,14, 8,18)
| = 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)
ff( 8,18, 6,14) = ff( 6,14, 8,18)*5.42e+10*T932*exp(-72.262/T9)
f(j,69)=ff( 6,14, 8,18)
r(j,69)=ff( 8,18, 6,14)
c(70)= ' N13(P,G)O14'
c f70
qq( 7,13, 8,14) = 4.628
ff( 7,13, 8,14)
| = 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)
ff( 8,14, 7,13) = ff( 7,13, 8,14)*3.57e+10*T932*exp(-53.706/T9)
f(j,70)=ff( 7,13, 8,14)
r(j,70)=ff( 8,14, 7,13)
c(71)= 'N14(P,G)O15'
c f71
qq( 7,14, 8,15) = 7.297
ff( 7,14, 8,15)
| = 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)
ff( 8,15, 7,14) = ff( 7,14, 8,15)*2.70e+10*T932*exp(-84.678/T9)
f(j,71)=ff( 7,14, 8,15)
r(j,71)=ff( 8,15, 7,14)
c(72)= 'N14(P,N)O14'
c f72
qq( 7,14, 8,14) = -5.925
ff( 7,14, 8,14)
| = 6.74e+07*(1.0+0.658*T912+0.379*T9)*exp(-68.762/T9)
ff( 8,14, 7,14) = 6.74e+07*(1.0+0.658*T912+0.379*T9)*2.99e+00
f(j,72)=ff( 7,14, 8,14)
r(j,72)=ff( 8,14, 7,14)
c(73)= 'N14(P,A)C11'
c f73
qq( 7,14, 6,11) = -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)
ff( 7,14, 6,11)
| = 2.63e+16*T9A56/T932*exp(-31.883/T9A13-33.915/T9)
ff( 6,11, 7,14)
| = 2.63e+16*T9A56/T932*exp(-31.883/T9A13)*2.72e-01
f(j,73)=ff( 7,14, 6,11)
r(j,73)=ff( 6,11, 7,14)
c(74)= 'N14(A,G)F18'
c f74
qq( 7,14, 9,18) = 4.415
ff( 7,14, 9,18)
| = 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)
ff( 9,18, 7,14) = ff( 7,14, 9,18)*5.42e+10*T932*exp(-51.236/T9)
f(j,74)=ff( 7,14, 9,18)
r(j,74)=ff( 9,18, 7,14)
c(75)= 'N14(A,N)F17'
c f75
qq( 7,14, 9,17) = -4.735
ff( 7,14, 9,17)
| = 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)
ff( 9,17, 7,14)
| = 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
f(j,75)=ff( 7,14, 9,17)
r(j,75)=ff( 9,17, 7,14)
c(76)= 'N15(P,G)O16'
c f76
qq( 7,15, 8,16) = 12.128
ff( 7,15, 8,16)
| = 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)
ff( 8,16, 7,15) = ff( 7,15, 8,16)*3.62e+10*T932*exp(-140.734/T9)
f(j,76)=ff( 7,15, 8,16)
r(j,76)=ff( 8,16, 7,15)
c(77)= 'N15(P,N)O15'
c f77
qq( 7,15, 8,15) = -3.536
ff( 7,15, 8,15)
| = 3.51e+08*(1.0+0.452*T912-0.191*T9)*exp(-41.037/T9)
ff( 8,15, 7,15) = 3.51e+08*(1.0+0.452*T912-0.191*T9)*9.98e-01
f(j,77)=ff( 7,15, 8,15)
r(j,77)=ff( 8,15, 7,15)
c(78)= 'N15(P,A)C12'
c f78
qq( 7,15, 6,12) = 4.966
ff( 7,15, 6,12)
| = 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))
ff( 6,12, 7,15) = ff( 7,15, 6,12)*7.06e-01*exp(-57.623/T9)
f(j,78)=ff( 7,15, 6,12)
r(j,78)=ff( 6,12, 7,15)
c(79)= 'N15(A,G)F19'
c f79
qq( 7,15, 9,19) = 4.014
ff( 7,15, 9,19)
| = 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)
ff( 9,19, 7,15) = ff( 7,15, 9,19)*5.54e+10*T932*exp(-46.578/T9)
f(j,79)=ff( 7,15, 9,19)
r(j,79)=ff( 9,19, 7,15)
c(80)= ' N15(A,N)F18'
c f80
qq( 7,15, 9,18) = -6.418
ff( 7,15, 9,18)
| = 3.14e+08*(1.0-0.641*T912+0.108*T9)*exp(-74.479/T9)
ff( 9,18, 7,15) = 3.14e+08*(1.0-0.641*T912+0.108*T9)*2.00e+00
f(j,80)=ff( 7,15, 9,18)
r(j,80)=ff( 9,18, 7,15)
c(81)= 'O14(A,G)NE18'
c f81
qq( 8,14,10,18) = 5.112
ff( 8,14,10,18)
| = 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)
ff(10,18, 8,14) = ff( 8,14,10,18)*5.42e+10*T932*exp(-59.328/T9)
f(j,81)=ff( 8,14,10,18)
r(j,81)=ff(10,18, 8,14)
c(82)= 'O14(A,P)F17'
c f82
qq( 8,14, 9,17) = 1.191
ff( 8,14, 9,17)
| = 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)
ff( 9,17, 8,14) = ff( 8,14, 9,17)*4.93e-01*exp(-13.820/T9)
f(j,82)=ff( 8,14, 9,17)
r(j,82)=ff( 9,17, 8,14)
c(83)= 'O15(A,G)NE19'
c f83
qq( 8,15,10,19) = 3.529
ff( 8,15,10,19)
| = 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)
ff(10,19, 8,15) = ff( 8,15,10,19)*5.54e+10*T932*exp(-40.957/T9)
f(j,83)=ff( 8,15,10,19)
r(j,83)=ff(10,19, 8,15)
c(84)= 'O16(P,G)F17'
c f84
qq( 8,16, 9,17) = 0.600
ff( 8,16, 9,17)
| = 1.50e+08/(T923*(1.0+2.13*(1.0-exp(-0.728*T923))))
| *exp(-16.692/T913)
ff( 9,17, 8,16) = ff( 8,16, 9,17)*3.03e+09*T932*exp(-6.968/T9)
f(j,84)=ff( 8,16, 9,17)
r(j,84)=ff( 9,17, 8,16)
c(85)= 'O16(P,A)N13'
c f85
qq( 8,16, 7,13) = -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)
ff( 8,16, 7,13)
| = 1.88e+18*T9A56/T932*exp(-35.829/T9A13-60.557/T9)
ff( 7,13, 8,16)
| = 1.88e+18*T9A56/T932*exp(-35.829/T9A13)*1.72e-01
f(j,85)=ff( 8,16, 7,13)
r(j,85)=ff( 7,13, 8,16)
c(86)= 'O16(A,G)NE20'
c f86
qq( 8,16,10,20) = 4.734
ff( 8,16,10,20)
| = 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)
ff(10,20, 8,16) = ff( 8,16,10,20)*5.65e+10*T932*exp(-54.937/T9)
tf2=ff( 8,16,10,20)
f(j,86)=ff( 8,16,10,20)
r(j,86)=ff(10,20, 8,16)
c(87)= 'O17(P,G)F18'
c f87
qq( 8,17, 9,18) = 5.607
T9A = T9/(1.0+2.69*T9)
T9A13 = T9A**(1.0/3.0)
T9A56 = T9A**(5.0/6.0)
ff( 8,17, 9,18)
|= 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)
c tf=ff( 8,17, 9,18)
c write(*,*) 't9=',t9,'tf=',tf,'tf2=',tf2
ff( 9,18, 8,17) = ff( 8,17, 9,18)*3.66e+10*T932*exp(-65.061/T9)
f(j,87)=ff( 8,17, 9,18)
r(j,87)=ff( 9,18, 8,17)
c(88)= 'O17(P,A)N14'
c f88
qq( 8,17, 7,14) = 1.191
rc88=rc87
ff( 8,17, 7,14)
| = 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)
ff( 7,14, 8,17) = ff( 8,17, 7,14)*6.76e-01*exp(-13.825/T9)
f(j,88)=ff( 8,17, 7,14)
r(j,88)= ff( 7,14, 8,17)
c(89)= 'O17(A,G)NE21'
c f89
qq( 8,17,10,21) = 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)
ff( 8,17,10,21)
| = 1.73e+17*FPT9A/GT9*T9A56/T932*exp(-39.914/T9A13)+3.50e+15
| *FT9A/GT9*T9A56/T932*exp(-39.914/T9A13)
ff(10,21, 8,17) = ff( 8,17,10,21)*8.63e+10*T932*exp(-85.305/T9)
f(j,89)=ff( 8,17,10,21)
r(j,89)=ff(10,21, 8,17)
c(90)= 'O17(A,N)NE20'
c f90
qq( 8,17,10,20) = 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
ff( 8,17,10,20) = 1.03e+18/GT9*T9A56/T932*exp(-39.914/T9A13)
ff(10,20, 8,17) = ff( 8,17,10,20)*1.86e+01*exp(-6.852/T9)
f(j,90)=ff( 8,17,10,20)
r(j,90)=ff(10,20, 8,17)
c(91)= 'O18(P,G)F19'
c f91
qq( 8,16, 9,19) = 7.994
ff( 8,16, 9,19)
| = 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)
ff( 9,19, 8,16) = ff( 8,16, 9,19)*9.20e+09*T932*exp(-92.769/T9)
f(j,91)=ff( 8,16, 9,19)
r(j,91)=ff( 9,19, 8,16)
c(92)= 'O18(P,A)N15'
c f92
qq( 8,18, 7,15) = 3.980
ff( 8,18, 7,15)
| = 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)
ff( 7,15, 8,18) = ff( 8,18, 7,15)*1.66e-01*exp(-46.191/T9)
f(j,92)=ff( 8,18, 7,15)
r(j,92)=ff( 7,15, 8,18)
c(93)= 'O18(A,G)NE22'
c f93
qq( 8,18,10,22) = 9.669
ff( 8,18,10,22)
| = 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)
ff(10,22, 8,18) = ff( 8,18,10,22)*5.85e+10*T932*exp(-112.208/T9)
f(j,93)=ff( 8,18,10,22)
r(j,93)=ff(10,22, 8,18)
c(94)= 'O18(A,N)NE21'
c f94
qq( 8,18,10,21) = -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)
ff( 8,18,10,21)
| = 7.22e+17*FT9A/GT9*T9A56/T932*exp(-40.056/T9A13)+150.31
| /GT9*exp(-8.045/T9)
ff(10,21, 8,18) = 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
f(j,94)=ff( 8,18,10,21)
r(j,94)=ff(10,21, 8,18)
c(95)= 'F19(P,G)NE20'
c f95
qq( 9,19,10,20) = 12.848
ff( 9,19,10,20)
| = (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))
ff(10,20, 9,19) = ff( 9,19,10,20)*3.70e+10*T932*exp(-149.093/T9)
f(j,95)= ff( 9,19,10,20)
r(j,95)=ff(10,20, 9,19)
c(96)= 'F19(P,N)NE19'
c f96
qq( 9,19,10,19) = -4.021
ff( 9,19,10,19)
| = 1.27e+08*(1.0-0.147*T912+0.069*T9)*exp(-46.659/T9)
ff(10,19, 9,19) = 1.27e+08*(1.0-0.147*T912+0.069*T9)*9.98e-01
f(j,96)=ff( 9,19,10,19)
r(j,96)=ff(10,19, 9,19)
c(97)= 'F19(P,A)O16'
c f97
qq( 9,19, 8,16) = 8.114
ff( 9,19, 8,16)
| = (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))
ff( 8,16, 9,19) = ff( 9,19, 8,16)*6.54e-01*exp(-94.156/T9)
f(j,97)=ff( 9,19, 8,16)
r(j,97)=ff( 8,16, 9,19)
c(98)= 'F19(A,P)NE22'
c f98
qq( 9,19,10,22) = 1.675
ff( 9,19,10,22)
| = 4.50e+18/T923*exp(-43.467/T913-(T9/0.637)**2)+7.98e+04*T932
| *exp(-12.760/T9)
ff(10,22, 9,19) = ff( 9,19,10,22)*6.36e+00*exp(-19.439/T9)
f(j,98)=ff( 9,19,10,22)
r(j,98)=ff(10,22, 9,19)
c(99)= 'NE19(P,G)NA20'
c f99
qq(10,19,11,20) = 2.199
ff(10,19,11,20)
| = 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)
ff(11,20,10,19) = ff(10,19,11,20)*7.39e+09*T932*exp(-25.519/T9)
f(j,99)=ff(10,19,11,20)
r(j,99)=ff(11,20,10,19)
c(100)= 'NE20(P,G)NA21'
c f100
qq(10,20,11,21) = 2.431
ff(10,20,11,21)
| = 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)
ff(11,21,10,20) = ff(10,20,11,21)*4.63e+09*T932*exp(-28.216/T9)
f(j,100)=ff(10,20,11,21)
r(j,100)=ff(11,21,10,20)
c(101)= 'NE20(P,A)F17'
c f101
qq(10,20, 9,17) = -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)
ff(10,20, 9,17)
| = 3.25e+19*(5.31+0.544*T9-0.0523*T9**2)*T9A56/T932
| *exp(-43.176/T9A13-47.969/T9)
ff( 9,17,10,20)
| = 3.25e+19*(5.31+0.544*T9-0.0523*T9**2)*T9A56/T932
| *exp(-43.176/T9A13)*5.37e-02
f(j,101)=ff(10,20,9,17)
r(j,101)=ff(9,17,10,20)
c(102)= 'NE20(A,G)MG24'
c f102
qq(10,20,12,24) = 9.312
ff(10,20,12,24)
| = (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)+rc102*3.20e+01/
| T923*exp(-9.383/T9))
| /(1.00+5.0*exp(-18.960/T9))
ff(12,24,10,20) = ff(10,20,12,24)*6.01e+10*T932*exp(-108.059/T9)
f(j,102)=ff(10,20,12,24)
r(j,102)=ff(12,24,10,20)
c(103)= 'NE21(P,G)NA22'
c f103
qq(10,21,11,22) = 6.738
ff(10,21,11,22)
| = 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)+rc103*7.80e-01/T932*exp(-1.085/T9)
ff(11,22,10,21) = ff(10,21,11,22)*1.06e+10*T932*exp(-78.194/T9)
f(j,103)=ff(10,21,11,22)
r(j,103)=ff(11,22,10,21)
c(104)= 'NE21(A,G)MG25'
c f104
qq(10,21,12,25)= 9.882
T9A = T9/(1.0+0.0537*T9)
T9A13 = T9A**(1.0/3.0)
T9A56 = T9A**(5.0/6.0)
ff(10,21,12,25)
| = (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)))
ff(12,25,10,21) = ff(10,21,12,25)*4.06e+10*T932*exp(-114.676/T9)
f(j,104)=ff(10,21,12,25)
r(j,104)=ff(12,25,10,21)
c(105)= 'NE21(A,N)MG24'
c f105
qq(10,21,12,24) = 2.551
T9A = T9/(1.0+0.0537*T9)
T9A13 = T9A**(1.0/3.0)
T9A56 = T9A**(5.0/6.0)
ff(10,21,12,24)
| = (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))
ff(12,24,10,21) = ff(10,21,12,24)*1.29e+01*exp(-29.606/T9)
f(j,105)=ff(10,21,12,24)
r(j,105)=ff(12,24,10,21)
c(106)= 'NE22(P,G)NA23'
c f106
qq(10,22,11,23) = 8.794
ff(10,22,11,23)
| = 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)
ff(11,23,10,22) = ff(10,22,11,23)*4.67e+09*T932*exp(-102.048/T9)
f(j,106)=ff(10,22,11,23)
r(j,106)=ff(11,23,10,22)
c(107)= 'NE22(A,G)MG26'
c f107
qq(10,22,12,26) = 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)
ff(10,22,12,26)
| = 4.16e+19*FPT9A/GT9*T9A56/T932*exp(-47.004/T9A13)
| +2.08e+16*FT9A/GT9*T9A56/T932*exp(-47.004/T9A13)
ff(12,26,10,22) = ff(10,22,12,26)*6.15e+10*T932*exp(-123.151/T9)
f(j,107)= ff(10,22,12,26)
r(j,107)= ff(12,26,10,22)
c(108)= 'NE22(A,N)MG25'
c f108
qq(10,22,12,25) = -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)
ff(10,22,12,25)
| = 4.16e+19*FT9A/GT9*T9A56/T932*exp(-47.004/T9A13)
| +1.44e-04/GT9*exp(-5.577/T9)
ff(12,25,10,22) = 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
f(j,108)= ff(10,22,12,25)
r(j,108)= ff(12,25,10,22)
c(109)= ' NA21(P,G)MG22'
c f109
qq(11,21,12,22) = 5.497
ff(11,21,12,22)
| = 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)
ff(12,22,11,21) = ff(11,21,12,22)*7.44e+10*T932*exp(-63.790/T9)
f(j,109)=ff(11,21,12,22)
r(j,109)=ff(12,22,11,21)
c(110)= 'NA22(N,P)NE22'
c f110
qq(11,22,10,22) = 3.624
ff(11,22,10,22)
| = 1.24e+08*exp(1.0-3.037e-02*T9+8.380e-03*T9**2-7.101e-04
| *T9**3)
ff(10,22,11,22) = ff(11,22,10,22)*7.01*exp(-42.059/T9)
f(j,110)=ff(11,22,10,22)
r(j,110)=ff(10,22,11,22)
c(111)= 'NA22(N,A)F19'
c f111
qq(11,22, 9,19) = 1.949
ff(11,22, 9,19)
| = 1.21e+06*exp(1.00+8.955e-01*T9-5.645e-02*T9**2+7.302e-04
| *T9**3)
ff( 9,19,11,22) = ff(11,22, 9,19)*1.10*exp(-22.620/T9)
f(j,111)=ff(11,22, 9,19)
r(j,111)=ff( 9,19,11,22)
c(112)= 'NA22(P,G)MG23'
c f112
qq(11,22,12,23) = 7.578
ff(11,22,12,23)
| = 9.63e-05*T932*exp(-0.517/T9)+2.51e+04*T9*exp(-2.013/T9)
ff(12,23,11,22) = ff(11,22,12,23)*3.27e+10*T932*exp(-87.933/T9)
f(j,112)=ff(11,22,12,23)
r(j,112)=ff(12,23,11,22)
c(113)= 'NA23(P,G)MG24'
c f113
qq(11,23,12,24) = 11.691
ff(11,23,12,24)
| = (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))
ff(12,24,11,23) = ff(11,23,12,24)*7.49e+10*T932*exp(-135.665/T9)
f(j,113)=ff(11,23,12,24)
r(j,113)=ff(12,24,11,23)
c(114)= 'NA23(P,N)MG23'
c f114
qq(11,23,12,23) = -4.841
T9A = T9/(1.0+0.141*T9)
T9A32 = T9A**(3.0/2.0)
ff(11,23,12,23) = 9.29e+08*(1.0-0.881*T9A32/T932)*exp(-56.173/T9)
ff(12,23,11,23) = 9.29e+08*(1.0-0.881*T9A32/T932)*9.98e-01
f(j,114)=ff(11,23,12,23)
r(j,114)=ff(12,23,11,23)
c(115)= 'NA23(P,A)NE20'
c f115
qq(11,23,10,20) = 2.379
ff(11,23,10,20)
| = 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)
ff(10,20,11,23) = ff(11,23,10,20)*1.25e+00*exp(-27.606/T9)
f(j,115)=ff(11,23,10,20)
r(j,115)=ff(10,20,11,23)
c(116)= 'NA23(A,N)AL26T'
c f116
qq(11,23,13,26) = -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)
ff(11,23,13,26)
| = 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)
c Note: r116=f141 See AL26T(N,A)NA23
ff(13,26,11,23)
| = 3.38e+06*exp(0.388*T9+9.08e-03*T9**2-2.07e-03*T9**3)
f(j,116)=ff(11,23,13,26)
r(j,116)=ff(13,26,11,23)
c(117)= 'NA23(A,N)AL26M'
c f117
qq(11,23,13,26) = -3.196
ff(11,23,13,26)
| = (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)
c Note: r117=f142; See AL26M(N,A)NA23
ff(13,26,11,23)
| = 4.18e+06+5.43e+07/T932*exp(-0.9653/T9)+6.97e+07
| *T927*exp(-1.494/T9)
f(j,117)=ff(11,23,13,26)
r(j,117)=ff(13,26,11,23)
c(118)= 'NA23(A,N)AL26G'
c f118
qq(11,23,13,26) = -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)
ff(11,23,13,26)
| = ((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)
c Note: r118=f143; See AL27G(N,A)NA23
ff(13,26,11,23)
| = (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)
f(j,118)=ff(11,23,13,26)
r(j,118)=ff(13,26,11,23)
c(119)= 'MG24(P,G)AL25'
c f119
qq(12,24,13,25) = 2.271
ff(12,24,13,25)
| = (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))
ff(13,25,12,24) = ff(12,24,13,25)*3.13e+09*T932*exp(-26.358/T9)
f(j,119)=ff(12,24,13,25)
r(j,119)=ff(13,25,12,24)
c(120)= 'MG24(P,A)NA21'
c f120
qq(12,24,11,21) = -6.880
T9A = T9/(1.0+0.127*T9)
T9A13 = T9A**(1.0/3.0)
T9A56 = T9A**(5.0/6.0)
ff(12,24,11,21)
| = 1.81e+21*(4.43+3.31*T9-0.229*T9**2)*T9A56/T932
| *exp(-49.967/T9A13-79.843/T9)
ff(11,21,12,24)
| = 1.81e+21*(4.43+3.31*T9-0.229*T9**2)*T9A56/T932
| *exp(-49.967/T9A13)*7.71e-02
f(j,120)=ff(12,24,11,21)
r(j,120)=ff(11,21,12,24)
c(121)= 'MG24(A,G)SI28'
c f121
qq(12,24,14,28) = 9.984
ff(12,24,14,28)
| = (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)+rc121*1.25e-03/T932*
| exp(-7.929/T9)
| +rc121*2.43e+01/T9*exp(-11.523/T9))/(1.0+5.0*
| exp(-15.882/T9))
ff(14,28,12,24) = ff(12,24,14,28)*6.27e+10*T932*exp(-115.862/T9)
f(j,121)=ff(12,24,14,28)
r(j,121)=ff(14,28,12,24)
c(122)= 'MG25(P,G)AL26T'
c f122
qq(12,25,13,26) = 6.306
ff(12,25,13,26)
| = 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)
ff(13,26,12,25) = ff(12,25,13,26)*1.03e+10*T932*exp(-73.183/T9)
f(j,122)= ff(12,25,13,26)
r(j,122)=ff(13,26,12,25)
c(123)= 'MG25(P,G)AL26M'
c f123
qq(12,25,13,26) = 6.078
ff(12,25,13,26)
| = 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))
ff(13,26,12,25) = ff(12,25,13,26)*1.13e+11*T932*exp(-70.532/T9)
f(j,123)= ff(12,25,13,26)
r(j,123)=ff(13,26,12,25)
c(124)= 'MG25(P,G)AL26G'
c f124
qq(12,25,13,26) = 6.306
ff(12,25,13,26)
| = 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))
ff(13,26,12,25) = ff(12,25,13,26)*1.03e+10*T932*exp(-73.183/T9)
f(j,124)= ff(12,25,13,26)
r(j,124)=ff(13,26,12,25)
c(125)= 'MG25(A,G)SI29'
c f125
qq(12,25,14,29) = 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
ff(12,25,14,29)
| = 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))
ff(14,29,12,25) = ff(12,25,14,29)*1.90e+11*T932*exp(-129.128/T9)
f(j,125)=ff(12,25,14,29)
r(j,125)= ff(14,29,12,25)
c(126)= 'MG25(A,N)SI28'
c f126
qq(12,25,14,28) = 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
ff(12,25,14,28) = 3.59e+20/GT9*T9A56/T932*exp(-53.410/T9A13)
ff(14,28,12,25) = ff(12,25,14,28)*2.00e+01*exp(-30.792/T9)
f(j,126)=ff(12,25,14,28)
r(j,126)= ff(14,28,12,25)
c(127)= 'MG25(A,P)AL28'
c f127
qq(12,25,13,28) = -1.206
ff(12,25,13,28)
| = 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)
ff(13,28,12,25)
| = 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
f(j,127)=ff(12,25,13,28)
r(j,127)= ff(13,28,12,25)
c(128)= 'MG26(P,G)AL27'
c f128
qq(12,26,13,27) = 8.272
ff(12,26,13,27)
| = (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))
ff(13,27,12,26) = ff(12,26,13,27)*3.14e+09*T932*exp(-95.990/T9)
f(j,128)=ff(12,26,13,27)
r(j,128)=ff(13,27,12,26)
c(129)= ' MG26(P,N)AL26T'
c f129
qq(12,26,13,26) = -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)
ff(12,26,13,26)
| = 3.09e+07*exp(0.0731*T9+0.0381*T9**2-3.22e-03*T9**3)
| *GPT9/GT9*11.0*exp(-55.544/T9)
c Note: r129=f138; See AL26T(N,P)MG26
ff(13,26,12,26)
| = 3.09e+07*exp(0.0731*T9+0.0381*T9**2-3.22e-03*T9**3)
f(j,129)=ff(12,26,13,26)
r(j,129)=ff(13,26,12,26)
c(130)= 'MG26(P,N)AL26M'
c f130
qq(12,26,13,26) = -5.015
rc130 = rc139
ff(12,26,13,26)
| = (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)
c Note: r130=f139; See AL26M(N,P)MG26
ff(13,26,12,26)
| = 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
f(j,130)=ff(12,26,13,26)
r(j,130)=ff(13,26,12,26)
c(131)= 'MG26(P,N)AL26G'
c f131
qq(12,26,13,26) = -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
ff(12,26,13,26)
| = ((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)
c Note: r131=f140; See AL26G(N,P)MG26
ff(13,26,12,26)
| = 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)
f(j,131)=ff(12,26,13,26)
r(j,131)=ff(13,26,12,26)
c(132)= 'MG26(A,G)SI30'
c f132
qq(12,26,14,30) = 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)
ff(12,26,14,30)
| = 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))
ff(14,30,12,26) = ff(12,26,14,30)*6.38e+10*T932*exp(-123.520/T9)
f(j,132)=ff(12,26,14,30)
r(j,132)=ff(14,30,12,26)
c(133)= 'MG26(A,N)SI29'
c f133
qq(12,26,14,29) = 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)
ff(12,26,14,29) = 2.93e+20/GT9*T9A56/T932*exp(-53.505/T9A13)
ff(14,29,12,26) = ff(12,26,14,29)*1.68e+00*exp(-0.401/T9)
f(j,133)=ff(12,26,14,29)
r(j,133)=ff(14,29,12,26)
c(134)= 'MG26(A,P)AL29'
c f134
qq(12,26,13,29) = -2.863
ff(12,26,13,29)
| = 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)
ff(13,29,12,26) = 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.6/T9)
f(j,134)=ff(12,26,13,29)
r(j,134)=ff(13,29,12,26)
c(135)= ' AL26T(G,P)MG25'
c f135
qq(13,26,12,25) = -6.306
ff(13,26,12,25)
| = (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)**2)+1.77e+04
| *T958*exp(-3.049/T9-(T9/13.0)**2))*1.03e+10*T932
| *exp(-73.183/T9)
c r135 = f122
f(j,135)=ff(13,26,12,25)
r(j,135)=f(j,122)
c(136)= 'AL26M(G,P)MG25'
c f136
qq(13,26,12,25) = -6.078
ff(13,26,12,25)
| = 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)
c r136 = f123
f(j,136)=ff(13,26,12,25)
r(j,136)=f(j,123)
c(137)= 'AL26G(G,P)MG25'
c f137
qq(13,26,12,25) = -6.306
ff(13,26,12,25)
| = 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)
cc r137 = f124
f(j,137)=ff(13,26,12,25)
r(j,137)=f(j,124)
c(138)= 'AL26T(N,P)MG26'
c f138
qq(13,26,12,26) = 4.786
ff(13,26,12,26)
| = 3.09e+07*exp(0.0731*T9+0.0381*T9**2-3.22e-03*T9**3)
c r138 = f129
c r129 = f138
f(j,138)=ff(13,26,12,26)
r(j,138)=f(j,129)
c(139)= ' AL26M(N,P)MG26'
c f139
qq(13,26,12,26) = 5.015
ff(13,26,12,26)
| = 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
c r139 = f130
c r130 = f139
f(j,139)=ff(13,26,12,26)
r(j,139)=f(j,130)
c(140)= 'AL26G(N,P)MG26'
c f140
qq(12,26,13,26) = 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
ff(12,26,13,26)
| = (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)
c r140 = f131
c r131 = f140
f(j,140)=ff(12,26,13,26)
r(j,140)=f(j,131)
c(141)= ' AL26T(N,A)NA23'
c f141
qq(13,26,11,23) = 2.968
ff(13,26,11,23)
| = 3.38e+06*exp(0.388*T9+9.08e-03*T9**2-2.07e-03*T9**3)
c r141=f116
c r116=f141
f(j,141)=ff(13,26,11,23)
r(j,141)=f(j,116)
c(142)= 'AL26M(N,A)NA23'
c f142
qq(13,26,11,23) = 3.196
ff(13,26,11,23)
| = 4.18e+06+5.43e+07/T932*exp(-0.9653/T9)+6.97e+07
| *T927*exp(-1.494/T9)
c r142=f117
f(j,142)=ff(13,26,11,23)
r(j,142)=f(j,117)
c(143)= ' AL26G(N,A)NA23'
c f143
qq(13,26,11,23) = 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)
ff(13,26,11,23)
| = 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)
c r143=f118
c r118=f143
c f(143)=ff(13,26,11,23)
f(j,143)=r(j,118)
r(j,143)=f(j,118)
c(144)= ' AL26T(P,G)SI27'
c f144
qq(13,26,14,27) = 7.464
ff(13,26,14,27)
| = 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)
ff(14,27,13,26) = ff(13,26,14,27)*3.46e+10*T932*exp(-86.621/T9)
f(j,144)=ff(13,26,14,27)
r(j,144)=ff(14,27,13,26)
c(145)= ' AL26M(P,G)SI27'
c f145
qq(13,26,14,27) = 7.693
ff(13,26,14,27)
| = 1.36e+14/T923*exp((-23.261/T913)*(1.0+2.004e-01*T9
| -1.538e-02*T9**2+5.723e-04*T9**3))
ff(14,27,13,26) = ff(13,26,14,27)*3.46e+10/11.0*T932
| *exp(-89.272/T9)
f(j,145)=ff(13,26,14,27)
r(j,145)=ff(14,27,13,26)
c(146)= ' AL26G(P,G)SI27'
c f146
qq(13,26,14,27) = 7.464
ff(13,26,14,27)
| = f(j,144)-1.0/11.0*exp(-2.651/T9)*f(j,145)
ff(14,27,13,26) = ff(13,26,14,27)*3.46e+10*T932*exp(-86.621/T9)
f(j,146)=ff(13,26,14,27)
r(j,146)=ff(14,27,13,26)
c(147)= 'AL27(P,G)SI28'
c f147
qq(13,27,14,28) = 11.585
ff(13,27,14,28)
| = (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)
| +rc147*1.63e-10/T932*exp(-1.001/T9))
| /(1.0+exp(-9.792/T9)/3.0+2.0*exp(-11.773/T9)/3.0)
ff(14,28,13,27) = ff(13,27,14,28)*1.13e+11*T932*exp(-134.434/T9)
f(j,147)=ff(13,27,14,28)
r(j,147)=ff(14,28,13,27)
c(148)= 'AL27(P,A)MG24'
c f148
qq(13,27,12,24) = 1.600
ff(13,27,12,24)
| = (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)+rc148*8.18e-09/T932*exp(-1.001/T9))
| /(1.00+exp(-9.792/T9)/3.0+2.0*exp(-11.773/T9)/3.0)
ff(12,24,13,27) = ff(13,27,12,24)*1.81e+00*exp(-18.572/T9)
f(j,148)=ff(13,27,12,24)
r(j,148)=ff(12,24,13,27)
c(149)= 'AL27(A,N)P30'
c f149
qq(13,27,15,30) = -2.636
ff(13,27,15,30)
| = 8.20e+04*exp(-30.588/T9)+5.21e+05*T974*exp(-33.554/T9)
ff(15,30,23,27)
| = 8.20e+04*6.75e+00+5.21e+05*T974*exp(-2.966/T9)*6.75e+00
f(j,149)=ff(13,27,15,30)
r(j,149)=ff(15,30,23,27)
c(150)='SI27(G,P)AL26T'
c f150
qq(14,27,13,26) = -7.464
ff(14,27,13,26)
| = (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)
f(j,150)=ff(14,27,13,26)
r(j,150) = f(j,144)
c(151)= 'SI27(G,P)AL26M'
c f151
qq(14,27,13,26) = -7.693
ff(14,27,13,26)
| = 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)
f(j,151)=ff(14,27,13,26)
c r151 = f145
r(j,151)=f(j,145)
c(152)= 'SI27(G,P)AL26G'
c f152
qq(14,27,13,26) = -7.464
ff(14,27,13,26)
| = (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)
f(j,152)=r(j,146)
r(j,152)=f(j,146)
c r152 = f146
c(153)= ' SI27(P,G)P28'
c f153
qq(14,27,15,28) = 2.065
ff(14,27,15,28)
| = 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)
ff(15,28,14,27) = ff(14,27,15,28)*1.62e+10*T932*exp(-23.960/T9)
f(j,153)=ff(14,27,15,28)
r(j,153)=ff(15,28,14,27)
c(154)= 'SI28(P,G)P29'
c f154
qq(14,28,15,29) = 2.747
ff(14,28,15,29)
| = 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)
ff(15,29,14,28) = ff(14,28,15,29)*9.46e+09*T932*exp(-31.879/T9)
f(j,154)=ff(14,28,15,29)
r(j,154)=ff(15,29,14,28)
c(155)= 'SI29(P,G)P30'
c f155
qq(14,29,15,30) = 5.601
ff(14,29,15,30)
| = 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)
ff(15,30,14,29) = ff(14,29,15,30)*1.26e+10*T932*exp(-65.002/T9)
f(j,155)=ff(14,29,15,30)
r(j,155)=ff(15,30,14,29)
c(156)= 'SI30(P,G)P31'
c f156
qq (14,30,15,31)= 7.297
ff(14,30,15,31)
| = 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)
ff(15,31,14,30) = ff(14,30,15,31)*9.50e+09*T932*exp(-84.673/T9)
f(j,156)=ff(14,30,15,31)
r(j,156)=ff(15,31,14,30)
c(157)= 'C12+C12 (MG24)'
c f157
qq( 6,12,12,24) = 13.933
T9A=T9/(1.0+0.0396*T9)
T9A13=T9A**(1.0/3.0)
T9A56=T9A**(5.0/6.0)
ff( 6,12,12,24)
| = 4.27e+26*T9A56/T932*exp(-84.165/T9A13-2.12e-03*T9**3)
ff(12,24, 6,12) = 0.00
f(j,157)=ff( 6,12,12,24)
r(j,157)=ff(12,24, 6,12)
c(158)= 'C12+O16 (SI28)'
c f158
qq( 6,12,14,28) = 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)
if (T9.ge.0.5) then
c! or ff( 8,16,14,28)
ff( 6,12,14,28)
| = 1.72e+31*T9A56/T932*exp(-106.594/T9A13)/(exp(-0.180*T9A
| **2)+1.06e-03*exp(2.562*T9A23))
else if (T9.lt.0.5) then
ff( 6,12,14,28) = 2.6288035e-29
else
endif
ff(14,28, 6,12) = 0.00
T9MAX158 = 0.5
f(j,158)=ff( 6,12,14,28)
r(j,158)=ff(14,28, 6,12)
c(159)='O16+O16 (S32)'
c f159
qq( 8,16,16,32) = 16.542
ff( 8,16,16,32)
| = 7.10e+36/T923*exp(-135.93/T913-0.629*T923-0.445*T943
| +0.0103*T9**2)
ff(16,32, 8,16) = 0.00
f(j,159)=ff( 8,16,16,32)
r(j,159)=ff(16,32, 8,16)
88 continue
c output dat file
sign=0
do 99 i=2,159
do 188 j=1,60
t9=temp(j)
if((i.eq.27).or.(i.eq.28)) then
if(f(j,i).lt.1.0e-99) f(j,i)=0.0e+0
else
if(f(j,i).lt.1.0e-50) f(j,i)=0.0e+0
endif
if(sign.eq.0) then
write(20,299)c(i)
299 format(5x,A40)
write(20,*)' t9',' reaction rate'
sign=1
endif
if((i.ge.68).and.(i.le.97)) then
if(t9.ge.0.002) then
write(20,77)t9,f(j,i)
endif
elseif((i.ge.98).and.(i.le.109)) then
if(t9.ge.0.003) then
write(20,77)t9,f(j,i)
endif
elseif(((i.ge.119).and.(i.le.128)).or.((i.ge.132).and.
|(i.le.137)).or.((i.ge.144).and.(i.le.148)).or.((i.ge.150).
|and.(i.le.152))) then
if(t9.ge.0.004) then
write(20,77)t9,f(j,i)
endif
elseif((i.eq.149).or.((i.ge.153).and.(i.le.156))) then
if(t9.ge.0.006) then
write(20,77)t9,f(j,i)
endif
elseif((i.eq.2).or.(i.eq.3).or.(i.eq.13).or.(i.eq.14)
|.or.(i.eq.41)) then
if(t9.le.3.000) then
write(20,77)t9,f(j,i)
endif
elseif((i.eq.157).or.(i.eq.159)) then
if(t9.ge.0.110) then
write(20,77)t9,f(j,i)
endif
elseif(i.eq.158) then
if(t9.ge.0.500) then
write(20,77)t9,f(j,i)
endif
else
write(20,77)t9,f(j,i)
endif
77 format(1x,F6.3,8x,1PE8.2)
188 continue
sign=0
write(20,*)
write(20,*)
99 continue
close(unit=20) ! close the .dat file
stop
end
Last modified March 30, 1998.
Any questions or comments should be directed to Michael Smith
at msmith@mail.phy.ornl.gov