Caughlan and Fowler 1988
Reaction Rates Temperature Derivatives
Fortran Program
program derivatives
c these "analytic" derivative expressions are from Mathematica.
real*8 temp(60),dfdt(60,159),drdt(60,159),E
integer sign
character *40, c(159)
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/
E=dexp(1.0d0)
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
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 open the output file
open(unit=21,file='derivatives.dat',status='unknown')
write(21,*)'H1(E-,NU)N1'
write(21,*)'Q1 = -0.782'
write(21,*)'SEE AP.J.SUPPL. 48,279(1982) & AP.J.293,1(1985)'
write(21,*)
write(21,*)
do 88 j=1,60
t9=temp(j)
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
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(1)= 'H1(E-,NU)N1'
c SEE AP.J.SUPPL. 48,279(1982) & AP.J. 293,1(1985)
c(2)= 'H1(P,E+NU)H2'
dfdt(j,2)= 1.2538e-15*(3.6034 - 1.689*T9**(1.d0/3) +
| 3.7966*T9**(2.d0/3) + 3.38*T9 + 1.*T9**(4.d0/3))*
| exp(-3.38*T9**(-1.d0/3))*T9**(-2)
c(3)= 'H1(E-P,NU)H2'
dfdt(j,3)= 6.6776e-20*(2.2946e-1 - 4.0489e-1*T9**(1.d0/3) +
| 2.3771*T9**(2.d0/3) - 1.*T9)*exp(-3.38*T9**(-1.d0/3))*
| T9**(-5.d0/2)
c(4)= 'H2(P,G)HE3'
dfdt(j,4)= 1.9787e3*(1.4038 - 5.9749e-1*T9**(1.d0/3) +
| 4.7025*T9**(2.d0/3) + 3.72*T9 + 1.*T9**(4.d0/3))*
| exp(-3.72*T9**(-1.d0/3))*T9**(-2)
drdt(j,4)= 3.6512e13*(2.65 + 3.7333e-2*T9**(-2.d0/3) +
| 2.2533*T9**(-1.d0/3))*
| exp(-6.375e1*(1/T9) - 3.72*T9**(-1.d0/3))*T9**(5.d0/6) +
| 3.0427e13*(1. + 1.12e-1*T9**(1.d0/3) + 3.38*T9**(2.d0/3) +
| 2.65*T9)*exp(-6.375e1*(1/T9) - 3.72*T9**(-1.d0/3))*
| T9**(-1.d0/6) + 3.6512e13*
| (6.375e1*T9**(-2) + 1.24*T9**(-4.d0/3))*
| (1. + 1.12e-1*T9**(1.d0/3) + 3.38*T9**(2.d0/3) + 2.65*T9)*
| exp(-6.375e1*(1/T9) - 3.72*T9**(-1.d0/3))*T9**(5.d0/6)
c(5)= 'H2(P,N)2H1'
dfdt(j,5)= 2.3115e7*(3.7413e1 + 2.9332e1*T9**(1.d0/3) +
| 1.4742e1*T9**(2.d0/3) + 2.7224e1*T9 +
| 1.0005*T9**(4.d0/3) + 1.5743*T9**(5.d0/3) + 1.*T9**2)*
| exp(-2.5815e1*(1/T9) - 3.72*T9**(-1.d0/3))*T9**(-2)
drdt(j,5)= 1.7613e-2*(1. + 7.84e-1*T9**(1.d0/3) +
| 3.46e-1*T9**(2.d0/3) +
| 6.9e-1*T9)*exp(-3.72*T9**(-1.d0/3))*T9**(-17.d0/6) -
| 2.1306e-2*(1. + 7.84e-1*T9**(1.d0/3) +
| 3.46e-1*T9**(2.d0/3) + 6.9e-1*T9)*
| exp(-3.72*T9**(-1.d0/3))*T9**(-5.d0/2) +
| 1.4204e-2*(6.9e-1 + 2.6133e-1*T9**(-2.d0/3) +
| 2.3067e-1*T9**(-1.d0/3))*exp(-3.72*T9**(-1.d0/3))*
| T9**(-3.d0/2)
c(6)= 'H2(D,G)HE4'
dfdt(j,6)= 8.954*(7.6721 - 2.8517*T9**(1.d0/3)-
| 1.734*T9**(2.d0/3)-
| 1.0664*T9 + 5.6279e-1*T9**(4.d0/3) +
| 1.8013*T9**(5.d0/3) + 1.*T9**2)*exp(-4.258*T9**(-1.d0/3))*
| T9**(-2)
drdt(j,6)= 2.1925e12*(-1.39e-1 + 3.2667e-2*T9**(-2.d0/3) -
| 1.3533e-1*T9**(-1.d0/3) + 1.4133e-1*T9**(1.d0/3) +
| 3.0833e-1*T9**(2.d0/3))*
| exp(-2.7673e2*(1/T9) - 4.258*T9**(-1.d0/3))*T9**(5.d0/6) +
| 1.8271e12*(1. + 9.8e-2*T9**(1.d0/3) -
| 2.03e-1*T9**(2.d0/3) - 1.39e-1*T9 +
| 1.06e-1*T9**(4.d0/3) + 1.85e-1*T9**(5.d0/3))*
| exp(-2.7673e2*(1/T9) - 4.258*T9**(-1.d0/3))*T9**(-1.d0/6)
| + 2.1925e12*(2.7673e2*T9**(-2) + 1.4193*T9**(-4.d0/3))*
| (1. + 9.8e-2*T9**(1.d0/3) - 2.03e-1*T9**(2.d0/3) -
| 1.39e-1*T9 + 1.06e-1*T9**(4.d0/3) + 1.85e-1*T9**(5.d0/3))
| *exp(-2.7673e2*(1/T9) - 4.258*T9**(-1.d0/3))*T9**(5.d0/6)
c(7)= 'H2(D,N)HE3'
dfdt(j,7)= 4.3146e8*(1.2764 - 4.7444e-1*T9**(1.d0/3) +
| 5.0415e-1*T9**(2.d0/3) + 3.6632e-1*T9 +
| 9.0036e-1*T9**(4.d0/3) + 1.8018*T9**(5.d0/3) + 1.*T9**2)*
| exp(-4.258*T9**(-1.d0/3))*T9**(-2)
drdt(j,7)= 6.7124e8*(2.87e-1 + 3.2667e-2*T9**(-2.d0/3) +
| 2.7867e-1*T9**(-1.d0/3) + 8.5067e-1*T9**(1.d0/3) +
| 1.8533*T9**(2.d0/3))*
| exp(-3.7935e1*(1/T9) - 4.258*T9**(-1.d0/3))*T9**(-2.d0/3)
| - 4.4749e8*(1. + 9.8e-2*T9**(1.d0/3) +
| 4.18e-1*T9**(2.d0/3) + 2.87e-1*T9 +
| 6.38e-1*T9**(4.d0/3) + 1.112*T9**(5.d0/3))*
| exp(-3.7935e1*(1/T9) - 4.258*T9**(-1.d0/3))*T9**(-5.d0/3)
| + 6.7124e8*(3.7935e1*T9**(-2) + 1.4193*T9**(-4.d0/3))*
| (1. + 9.8e-2*T9**(1.d0/3) + 4.18e-1*T9**(2.d0/3) +
| 2.87e-1*T9 + 6.38e-1*T9**(4.d0/3) + 1.112*T9**(5.d0/3))*
| exp(-3.7935e1*(1/T9) - 4.258*T9**(-1.d0/3))*T9**(-2.d0/3)
c(8)= 'H2(D,P)H3'
dfdt(j,8)= 3.907e8*(1.5004 - 5.5769e-1*T9**(1.d0/3) +
| 3.1334e-2*T9**(2.d0/3) + 4.5161e-2*T9 +
| 8.253e-1*T9**(4.d0/3) + 1.802*T9**(5.d0/3) + 1.*T9**2)*
| exp(-4.258*T9**(-1.d0/3))*T9**(-2)
drdt(j,8)= 7.1449e8*(3.01e-2 + 3.2667e-2*T9**(-2.d0/3) +
| 2.9267e-2*T9**(-1.d0/3) + 7.24e-1*T9**(1.d0/3) +
| 1.5767*T9**(2.d0/3))*
| exp(-4.6798e1*(1/T9) - 4.258*T9**(-1.d0/3))*T9**(-2.d0/3)
| - 4.7633e8*(1. + 9.8e-2*T9**(1.d0/3) +
| 4.39e-2*T9**(2.d0/3) + 3.01e-2*T9 +
| 5.43e-1*T9**(4.d0/3) + 9.46e-1*T9**(5.d0/3))*
| exp(-4.6798e1*(1/T9) - 4.258*T9**(-1.d0/3))*T9**(-5.d0/3)
| + 7.1449e8*(4.6798e1*T9**(-2) + 1.4193*T9**(-4.d0/3))*
| (1. + 9.8e-2*T9**(1.d0/3) + 4.39e-2*T9**(2.d0/3) +
| 3.01e-2*T9 + 5.43e-1*T9**(4.d0/3) + 9.46e-1*T9**(5.d0/3))
| *exp(-4.6798e1*(1/T9) - 4.258*T9**(-1.d0/3))*T9**(-2.d0/3)
c(9)= 'H3(P,G)HE4'
dfdt(j,9)= 2.332e4*(1.2167 - 4.9753e-1*T9**(1.d0/3) +
| 2.01*T9**(2.d0/3) + 1.533*T9 + 1.0666*T9**(4.d0/3) +
| 1.6362*T9**(5.d0/3) + 1.*T9**2)*exp(-3.869*T9**(-1.d0/3))*
| T9**(-2)
drdt(j,9)= 5.742e14*(1.26 + 3.6e-2*T9**(-2.d0/3) +
| 1.12*T9**(-1.d0/3) +
| 7.3467e-1*T9**(1.d0/3) + 1.7667*T9**(2.d0/3))*
| exp(-2.2993e2*(1/T9) - 3.869*T9**(-1.d0/3))*T9**(5.d0/6) +
| 4.785e14*(1. + 1.08e-1*T9**(1.d0/3) + 1.68*T9**(2.d0/3) +
| 1.26*T9 + 5.51e-1*T9**(4.d0/3) + 1.06*T9**(5.d0/3))*
| exp(-2.2993e2*(1/T9) - 3.869*T9**(-1.d0/3))*T9**(-1.d0/6)
| + 5.742e14*(2.2993e2*T9**(-2) + 1.2897*T9**(-4.d0/3))*
| (1. + 1.08e-1*T9**(1.d0/3) + 1.68*T9**(2.d0/3) + 1.26*T9 +
| 5.51e-1*T9**(4.d0/3) + 1.06*T9**(5.d0/3))*
| exp(-2.2993e2*(1/T9) - 3.869*T9**(-1.d0/3))*T9**(5.d0/6)
c(10)= 'H3(P,N)HE3'
dfdt(j,10)= 6.9286e7*(9.0439e1 - 1.3566e1*Sqrt(T9)
| + 8.863*T9 -
| 7.6531e-1*T9**(3.d0/2) + 1.*T9**2)*exp(-8.863*(1/T9))*
| T9**(-2)
drdt(j,10)= 7.0559e8*(9.8e-2 - 7.5e-2*(1/Sqrt(T9)))
c(11)= 'H3(D,N)HE4'
dfdt(j,11)= 4.5658e8*exp(-5.23e-1*(1/T9))*T9**(-8.d0/3) -
| 5.82e8*exp(-5.23e-1*(1/T9))*T9**(-5.d0/3) +
| 2.3245e12*(1.0673e-3 + 4.1763e-2*T9**(1.d0/3) +
| 4.0371e-2*T9**(2.d0/3) + 4.8817e-1*T9 + 1.*T9**(4.d0/3))*
| exp(-4.524*T9**(-1.d0/3) - 6.9444e1*T9**2)*T9**(-4.d0/3) -
| 5.3933e10*(1. + 9.2e-2*T9**(1.d0/3) + 1.8*T9**(2.d0/3) +
| 1.16*T9 + 1.052e1*T9**(4.d0/3) + 1.724e1*T9**(5.d0/3))*
| exp(-4.524*T9**(-1.d0/3) - 6.9444e1*T9**2)*T9**(-5.d0/3) +
| 1.9371e14*(5.8005e-2 + 5.3364e-3*T9**(1.d0/3) +
| 1.0441e-1*T9**(2.d0/3) + 6.7285e-2*T9 +
| 6.1021e-1*T9**(4.d0/3) + 1.*T9**(5.d0/3))*
| (1.0858e-2 - 1.*T9**(7.d0/3))*
| exp(-4.524*T9**(-1.d0/3) - 6.9444e1*T9**2)*T9**(-2)
drdt(j,11)= 1.1308e3*(8.73e8*exp(-5.23e-1*(1/T9))
| *T9**(-2.d0/3) +
| 8.09e10*(1. + 9.2e-2*T9**(1.d0/3) + 1.8*T9**(2.d0/3) +
| 1.16*T9 + 1.052e1*T9**(4.d0/3) + 1.724e1*T9**(5.d0/3))
| *exp(-4.524*T9**(-1.d0/3) - 6.9444e1*T9**2)*
| T9**(-2.d0/3))*exp(-2.0412e2*(1/T9))*T9**(-2) +
| 5.54*(4.5658e8*exp(-5.23e-1*(1/T9))*T9**(-8.d0/3) -
| 5.82e8*exp(-5.23e-1*(1/T9))*T9**(-5.d0/3) +
| 8.09e10*(1.16 + 3.0667e-2*T9**(-2.d0/3) +
| 1.2*T9**(-1.d0/3) + 1.4027e1*T9**(1.d0/3) +
| 2.8733e1*T9**(2.d0/3))*
| exp(-4.524*T9**(-1.d0/3) - 6.9444e1*T9**2)*T9**(-2.d0/3)
| - 5.3933e10*(1. + 9.2e-2*T9**(1.d0/3) +
| 1.8*T9**(2.d0/3) + 1.16*T9 + 1.052e1*T9**(4.d0/3) +
| 1.724e1*T9**(5.d0/3))*
| exp(-4.524*T9**(-1.d0/3) - 6.9444e1*T9**2)*T9**(-5.d0/3)
| + 8.09e10*(1.508*T9**(-4.d0/3) - 1.3889e2*T9)*
| (1. + 9.2e-2*T9**(1.d0/3) + 1.8*T9**(2.d0/3) +
| 1.16*T9 + 1.052e1*T9**(4.d0/3) + 1.724e1*T9**(5.d0/3))
| *exp(-4.524*T9**(-1.d0/3) - 6.9444e1*T9**2)*
| T9**(-2.d0/3))*exp(-2.0412e2*(1/T9))
c(12)= 'H3(T,2N)HE4'
dfdt(j,12)= 3.7575e8*(7.2178 - 2.3422*T9**(1.d0/3) -
| 3.4115*T9**(2.d0/3) - 1.9632*T9 +
| 6.6527e-1*T9**(4.d0/3) + 2.0625*T9**(5.d0/3) + 1.*T9**2)*
| exp(-4.872*T9**(-1.d0/3))*T9**(-2)
drdt(j,12)= 5.6446e-1*(-2.72e-1 + 2.8667e-2*T9**(-2.d0/3) -
| 3.0333e-1*T9**(-1.d0/3) + 1.9733e-1*T9**(1.d0/3) +
| 3.75e-1*T9**(2.d0/3))*
| exp(-1.315e2*(1/T9) - 4.872*T9**(-1.d0/3))*T9**(-13.d0/6)
| - 1.223*(1. + 8.6e-2*T9**(1.d0/3) - 4.55e-1*T9**(2.d0/3) -
| 2.72e-1*T9 + 1.48e-1*T9**(4.d0/3) + 2.25e-1*T9**(5.d0/3))
| *exp(-1.315e2*(1/T9) - 4.872*T9**(-1.d0/3))*T9**(-19.d0/6)
| + 5.6446e-1*(1.315e2*T9**(-2) + 1.624*T9**(-4.d0/3))*
| (1. + 8.6e-2*T9**(1.d0/3) - 4.55e-1*T9**(2.d0/3) -
| 2.72e-1*T9 + 1.48e-1*T9**(4.d0/3) + 2.25e-1*T9**(5.d0/3))
| *exp(-1.315e2*(1/T9) - 4.872*T9**(-1.d0/3))*T9**(-13.d0/6)
c(13)= 'HE3(E-,NU)H3'
dfdt(j,13)= 1.0096e-10*(1.648e-2 + 2.2134e-1*T9 +
| 1.3604*T9**2 +
| 2.0472*T9**3 + 1.*T9**4)*exp(-2.158e-1*(1/T9))*
| (1/Sqrt(T9))
c(14)= 'HE3(P,E+NU)HE4'
dfdt(j,14)= 5.8533e-13*(3.0705 - 1.*T9**(1.d0/3))*
| exp(-6.141*T9**(-1.d0/3))*T9**(-2)
c(15)= 'HE3(D,P)HE4'
dfdt(j,15)= 7.4992e8*exp(-1.72*(1/T9))*T9**(-5.d0/2) -
| 2.18e8*exp(-1.72*(1/T9))*T9**(-3.d0/2) +
| 2.2659e11*(5.e-3 + 2.4483e-2*T9**(1.d0/3) +
| 1.4948e-2*T9**(2.d0/3) + 7.7586e-1*T9 + 1.*T9**(4.d0/3))*
| exp(-7.181*T9**(-1.d0/3) - 1.0078e1*T9**2)*T9**(-4.d0/3) -
| 3.9067e10*(1. + 5.8e-2*T9**(1.d0/3) +
| 1.42e-1*T9**(2.d0/3) + 5.78e-2*T9 + 2.25*T9**(4.d0/3) +
| 2.32*T9**(5.d0/3))*exp
| (-7.181*T9**(-1.d0/3) - 1.0078e1*T9**2)*T9**(-5.d0/3) +
| 2.7403e12*(4.3103e-1 + 2.5e-2*T9**(1.d0/3) +
| 6.1207e-2*T9**(2.d0/3) + 2.4914e-2*T9 +
| 9.6983e-1*T9**(4.d0/3) + 1.*T9**(5.d0/3))*
| (1.1876e-1 - 1.*T9**(7.d0/3))*
| exp(-7.181*T9**(-1.d0/3) - 1.0078e1*T9**2)*T9**(-2)
drdt(j,15)= 1.182e3*(4.36e8*exp(-1.72*(1/T9))
| *(1/Sqrt(T9))+
| 5.86e10*(1. + 5.8e-2*T9**(1.d0/3) +
| 1.42e-1*T9**(2.d0/3) + 5.78e-2*T9 +
| 2.25*T9**(4.d0/3) + 2.32*T9**(5.d0/3))*
| exp(-7.181*T9**(-1.d0/3) - 1.0078e1*T9**2)*T9**(-2.d0/3)
| )*exp(-2.1298e2*(1/T9))*T9**(-2) +
| 5.55*(7.4992e8*exp(-1.72*(1/T9))*T9**(-5.d0/2) -
| 2.18e8*exp(-1.72*(1/T9))*T9**(-3.d0/2) +
| 5.86e10*(5.78e-2 + 1.9333e-2*T9**(-2.d0/3) +
| 9.4667e-2*T9**(-1.d0/3) + 3.*T9**(1.d0/3) +
| 3.8667*T9**(2.d0/3))*
| exp(-7.181*T9**(-1.d0/3) - 1.0078e1*T9**2)*T9**(-2.d0/3)
| - 3.9067e10*(1. + 5.8e-2*T9**(1.d0/3) +
| 1.42e-1*T9**(2.d0/3) + 5.78e-2*T9 +
| 2.25*T9**(4.d0/3) + 2.32*T9**(5.d0/3))*
| exp(-7.181*T9**(-1.d0/3) - 1.0078e1*T9**2)*T9**(-5.d0/3)
| + 5.86e10*(2.3937*T9**(-4.d0/3) - 2.0156e1*T9)*
| (1. + 5.8e-2*T9**(1.d0/3) + 1.42e-1*T9**(2.d0/3) +
| 5.78e-2*T9 + 2.25*T9**(4.d0/3) + 2.32*T9**(5.d0/3))*
| exp(-7.181*T9**(-1.d0/3) - 1.0078e1*T9**2)*T9**(-2.d0/3)
| )*exp(-2.1298e2*(1/T9))
|
c(16)= 'HE3(T,D)HE4'
dfdt(j,16)= 8.19e9*(1.0489e2*T9 + 3.3908e1*T9*
| (T9*(1/(1. + 1.28e-1*T9)))**(1.d0/3) -
| 6.1035e1*(T9*(1/(1. + 1.28e-1*T9)))**(4.d0/3) -
| 1.5625e1*T9*(T9*(1/(1. + 1.28e-1*T9)))**(4.d0/3) -
| 1.*T9**2*(T9*(1/(1. + 1.28e-1*T9)))**(4.d0/3))*
| exp(-7.733*(T9*(1/(1. + 1.28e-1*T9)))**(-1.d0/3) )*
| T9**(-5.d0/2)*(T9*(1/(1. + 1.28e-1*T9)))**(-5.e-1)*
| (7.8125 + 1.*T9)**(-2)
drdt(j,16)= -1.3104e10*exp
| (-1.6618e2*(1/T9) -
| 7.733*(T9*(1/(1. + 1.28e-1*T9)))**(-1.d0/3) )*
| T9**(-5.d0/2)*(T9*(1/(1. + 1.28e-1*T9)))**(2.5d0/3) +
| 7.28e9*(1/(1. + 1.28e-1*T9) -
| 1.28e-1*T9*(1. + 1.28e-1*T9)**(-2))*
| exp(-1.6618e2*(1/T9) -
| 7.733*(T9*(1/(1. + 1.28e-1*T9)))**(-1.d0/3) )*
| T9**(-3.d0/2)*(T9*(1/(1. + 1.28e-1*T9)))**(-.5d0/3) +
| 8.736e9*(1.6618e2*T9**(-2) +
| 2.5777*(1/(1. + 1.28e-1*T9) -
| 1.28e-1*T9*(1. + 1.28e-1*T9)**(-2))*
| (T9*(1/(1. + 1.28e-1*T9)))**(-4.d0/3))*
| exp(-1.6618e2*(1/T9) -
| 7.733*(T9*(1/(1. + 1.28e-1*T9)))**(-1.d0/3) )*
| T9**(-3.d0/2)*(T9*(1/(1. + 1.28e-1*T9)))**(2.5d0/3)
c(17)= 'HE3(T,NP)HE4'
dfdt(j,17)= 1.1565e10*(1.2994e2*T9 +
| 4.2008e1*T9*(T9*(1/(1. + 1.15e-1*T9)))**(1.d0/3) -
| 7.5614e1*(T9*(1/(1. + 1.15e-1*T9)))**(4.d0/3) -
| 1.7391e1*T9*(T9*(1/(1. + 1.15e-1*T9)))**(4.d0/3) -
| 1.*T9**2*(T9*(1/(1. + 1.15e-1*T9)))**(4.d0/3))*
| exp(-7.733*(T9*(1/(1. + 1.15e-1*T9)))**(-1.d0/3) )*
| T9**(-5.d0/2)*(T9*(1/(1. + 1.15e-1*T9)))**(-5.e-1)*
| (8.6957 + 1.*T9)**(-2)
drdt(j,17)= -7.8411*exp(-1.4037e2*(1/T9) -
| 7.733*(T9*(1/(1. + 1.15e-1*T9)))**(-1.d0/3) )*
| T9**(-4)*(T9*(1/(1. + 1.15e-1*T9)))**(2.5d0/3) +
| 2.1781*(1/(1. + 1.15e-1*T9) -
| 1.15e-1*T9*(1. + 1.15e-1*T9)**(-2))*
| exp(-1.4037e2*(1/T9) -
| 7.733*(T9*(1/(1. + 1.15e-1*T9)))**(-1.d0/3) )*
| T9**(-3)*(T9*(1/(1. + 1.15e-1*T9)))**(-.5d0/3) +
| 2.6137*(1.4037e2*T9**(-2) +
| 2.5777*(1/(1. + 1.15e-1*T9) -
| 1.15e-1*T9*(1. + 1.15e-1*T9)**(-2))*
| (T9*(1/(1. + 1.15e-1*T9)))**(-4.d0/3))*
| exp(-1.4037e2*(1/T9) -
| 7.733*(T9*(1/(1. + 1.15e-1*T9)))**(-1.d0/3) )*
| T9**(-3)*(T9*(1/(1. + 1.15e-1*T9)))**(2.5d0/3)
c(18)= 'HE3(HE3,2P)HE4'
dfdt(j,18)= 1.2865e10*(1.9211e1 - 2.4767*T9**(1.d0/3) -
| 1.0081e1*T9**(2.d0/3) - 2.3822*T9 + 6.5875*T9**(4.d0/3) +
| 5.1969*T9**(5.d0/3) + 1.*T9**2)*
| exp(-1.2276e1*T9**(-1.d0/3))*T9**(-2)
drdt(j,18)= 2.0476e1*(-1.24e-1 + 1.1333e-2*T9**(-2.d0/3) -
| 3.48e-1*T9**(-1.d0/3) + 4.7067e-1*T9**(1.d0/3) +
| 3.55e-1*T9**(2.d0/3))*
| exp(-1.4923e2*(1/T9) - 1.2276e1*T9**(-1.d0/3))*
| T9**(-13.d0/6) - 4.4364e1*
| (1. + 3.4e-2*T9**(1.d0/3) - 5.22e-1*T9**(2.d0/3) -
| 1.24e-1*T9 + 3.53e-1*T9**(4.d0/3) + 2.13e-1*T9**(5.d0/3))
| *exp(-1.4923e2*(1/T9) - 1.2276e1*T9**(-1.d0/3))*
| T9**(-19.d0/6) + 2.0476e1*
| (1.4923e2*T9**(-2) + 4.092*T9**(-4.d0/3))*
| (1. + 3.4e-2*T9**(1.d0/3) - 5.22e-1*T9**(2.d0/3) -
| 1.24e-1*T9 + 3.53e-1*T9**(4.d0/3) + 2.13e-1*T9**(5.d0/3))
| *exp(-1.4923e2*(1/T9) - 1.2276e1*T9**(-1.d0/3))*
| T9**(-13.d0/6)
c(19)= 'HE4(NN,G)HE6'
dfdt(j,19)= 5.5752e-12*(6.9457e1*rc19 - 4.9078*rc19*T9 -
| 1.*rc19*T9**2)*exp(-9.585*(1/T9))*T9**(-4)
drdt(j,19)= 4.3632e9*rc19*(1. + 1.38e-1*T9)*
| exp(-2.0904e1*(1/T9)) +
| 9.1208e10*rc19*(1. + 1.38e-1*T9)*
| exp(-2.0904e1*(1/T9))*(1/T9) +
| 6.0212e8*rc19*T9*exp(-2.0904e1*(1/T9))
c(20)= 'HE4(NP,G)LI6'
dfdt(j,20)= 3.465e-7*(2.5804e2 - 7.3137*T9 - 1.*T9**2)*
| exp(-1.9353e1*(1/T9))*T9**(-4)
drdt(j,20)= 3.3356e14*(1. + 7.5e-2*T9)*
| exp(-6.2286e1*(1/T9)) +
| 2.0776e16*(1. + 7.5e-2*T9)*exp(-6.2286e1*(1/T9))*
| (1/T9) + 2.5017e13*T9*exp(-6.2286e1*(1/T9))
c(21)= 'HE4(D,G)LI6'
dfdt(j,21)= -1.7578e1*exp(-7.423*T9**(-1.d0/3)) +
| 7.0349e2*exp(-8.228*(1/T9))*T9**(-7.d0/2) -
| 1.2825e2*exp(-8.228*(1/T9))*T9**(-5.d0/2) +
| 7.4477e1*exp(-7.423*T9**(-1.d0/3))*T9**(-2) -
| 1.5896e1*exp(-7.423*T9**(-1.d0/3))*T9**(-5.d0/3) -
| 3.6178e2*exp(-7.423*T9**(-1.d0/3))*T9**(-4.d0/3) +
| 6.5913e2*exp(-7.423*T9**(-1.d0/3))*(1/T9) +
| 4.5226e1*exp(-7.423*T9**(-1.d0/3))*T9**(-2.d0/3) -
| 5.5234e1*exp(-7.423*T9**(-1.d0/3))*T9**(-1.d0/3)
drdt(j,21)= 1.53e10*(7.0349e2*exp(-8.228*(1/T9))*
| T9**(-7.d0/2) -
| 1.2825e2*exp(-8.228*(1/T9))*T9**(-5.d0/2) +
| 3.01e1*(8.85 + 1.8667e-2*T9**(-2.d0/3) -
| 3.2333*T9**(-1.d0/3) - 7.8e-1*T9**(1.d0/3) -
| 9.7333e-1*T9**(2.d0/3))*exp(-7.423*T9**(-1.d0/3))*
| T9**(-2.d0/3) + 7.4477e1*
| (1. + 5.6e-2*T9**(1.d0/3) - 4.85*T9**(2.d0/3) +
| 8.85*T9 - 5.85e-1*T9**(4.d0/3) - 5.84e-1*T9**(5.d0/3))
| *exp(-7.423*T9**(-1.d0/3))*T9**(-2) -
| 2.0067e1*(1. + 5.6e-2*T9**(1.d0/3) - 4.85*T9**(2.d0/3) +
| 8.85*T9 - 5.85e-1*T9**(4.d0/3) - 5.84e-1*T9**(5.d0/3))
| *exp(-7.423*T9**(-1.d0/3))*T9**(-5.d0/3))*
| exp(-1.7118e1*(1/T9))*T9**(3.d0/2) +
| 2.6191e11*(8.55e1*exp(-8.228*(1/T9))*T9**(-3.d0/2) +
| 3.01e1*(1. + 5.6e-2*T9**(1.d0/3) - 4.85*T9**(2.d0/3) +
| 8.85*T9 - 5.85e-1*T9**(4.d0/3) - 5.84e-1*T9**(5.d0/3))
| *exp(-7.423*T9**(-1.d0/3))*T9**(-2.d0/3))*
| exp(-1.7118e1*(1/T9))*(1/Sqrt(T9)) +
| 2.295e10*(8.55e1*exp(-8.228*(1/T9))*T9**(-3.d0/2) +
| 3.01e1*(1. + 5.6e-2*T9**(1.d0/3) - 4.85*T9**(2.d0/3) +
| 8.85*T9 - 5.85e-1*T9**(4.d0/3) - 5.84e-1*T9**(5.d0/3))
| *exp(-7.423*T9**(-1.d0/3))*T9**(-2.d0/3))*
| exp(-1.7118e1*(1/T9))*Sqrt(T9)
c(22)= 'HE4(T,G)LI7'
dfdt(j,22)= 1.1618e5*(2.01e1 - 3.93*T9**(1.d0/3)
| - 9.1339*T9**(2.d0/3) -
| 3.3164*T9 + 2.4839*T9**(4.d0/3) + 3.4098*T9**(5.d0/3) +
| 1.*T9**2)*exp(-8.08*T9**(-1.d0/3))*T9**(-2)
drdt(j,22)= 9.6237e15*(-1.65e-1 + 1.7333e-2*T9**(-2.d0/3) -
| 2.9867e-1*T9**(-1.d0/3) + 1.92e-1*T9**(1.d0/3) +
| 2.2333e-1*T9**(2.d0/3))*
| exp(-2.864e1*(1/T9) - 8.08*T9**(-1.d0/3))*T9**(5.d0/6) +
| 8.0198e15*(1. + 5.2e-2*T9**(1.d0/3) -
| 4.48e-1*T9**(2.d0/3) - 1.65e-1*T9 +
| 1.44e-1*T9**(4.d0/3) + 1.34e-1*T9**(5.d0/3))*
| exp(-2.864e1*(1/T9) - 8.08*T9**(-1.d0/3))*T9**(-1.d0/6) +
| 9.6237e15*(2.864e1*T9**(-2) + 2.6933*T9**(-4.d0/3))*
| (1. + 5.2e-2*T9**(1.d0/3) - 4.48e-1*T9**(2.d0/3) -
| 1.65e-1*T9 + 1.44e-1*T9**(4.d0/3) + 1.34e-1*T9**(5.d0/3))
| *exp(-2.864e1*(1/T9) - 8.08*T9**(-1.d0/3))*T9**(5.d0/6)
c(23)= 'HE4(T,N)LI6'
dfdt(j,23)= 1.5744e11*exp(-5.7884e1*(1/T9))*T9**(-7.d0/2) -
| 4.08e9*exp(-5.7884e1*(1/T9))*T9**(-5.d0/2) +
| 1.8e8*(-1.3583e-18*(1.1917e14 - 1.*T9)*T9**(-3.d0/2)*
| (2.0333e-2 + 1.*T9)**(-2)*
| (T9*(1/(1. + 4.918e1*T9)))**5.e-1 +
| 3.915e-1*T9**(-5.d0/2)*(T9*(1/(1. + 4.918e1*T9)))**1.5
| )*exp(-5.5494e1*(1/T9)) +
| 9.9889e9*(1. - 2.61e-1*T9**(-3.d0/2)*
| (T9*(1/(1. + 4.918e1*T9)))**1.5)*
| exp(-5.5494e1*(1/T9))*T9**(-2)
drdt(j,23)= 6.0782e9*exp(-2.39*(1/T9))*T9**(-7.d0/2) -
| 3.8148e9*exp(-2.39*(1/T9))*T9**(-5.d0/2) +
| 1.683e8*(3.915e-1*T9**(-5.d0/2)*
| (T9*(1/(1. + 4.918e1*T9)))**1.5 -
| 3.915e-1*(-4.918e1*T9*(1. + 4.918e1*T9)**(-2) +
| 1/(1. + 4.918e1*T9))*T9**(-3.d0/2)*
| (T9*(1/(1. + 4.918e1*T9)))**5.e-1)
c(24)= 'HE4(HE3,G)BE7'
dfdt(j,24)= 8.415e6*(1.1632e3*T9 +
| 2.2673e2*T9*(T9*(1/(1. + 4.95e-2*T9)))**(1.d0/3) -
| 4.0812e2*(T9*(1/(1. + 4.95e-2*T9)))**(4.d0/3) -
| 4.0404e1*T9*(T9*(1/(1. + 4.95e-2*T9)))**(4.d0/3) -
| 1.*T9**2*(T9*(1/(1. + 4.95e-2*T9)))**(4.d0/3))*
| exp(-1.2826e1*(T9*(1/(1. + 4.95e-2*T9)))**(-1.d0/3) )*
| T9**(-5.d0/2)*(T9*(1/(1. + 4.95e-2*T9)))**(-5.e-1)*
| (2.0202e1 + 1.*T9)**(-2)
drdt(j,24)= 5.1892e16*(1/(1. + 4.95e-2*T9) -
| 4.95e-2*T9*(1. + 4.95e-2*T9)**(-2))*
| exp(-1.8423e1*(1/T9) -
| 1.2826e1*(T9*(1/(1. + 4.95e-2*T9)))**(-1.d0/3) )*
| (T9*(1/(1. + 4.95e-2*T9)))**(-.5d0/3) +
| 6.2271e16*(1.8423e1*T9**(-2) +
| 4.2753*(1/(1. + 4.95e-2*T9) -
| 4.95e-2*T9*(1. + 4.95e-2*T9)**(-2))*
| (T9*(1/(1. + 4.95e-2*T9)))**(-4.d0/3))*
| exp(-1.8423e1*(1/T9) -
| 1.2826e1*(T9*(1/(1. + 4.95e-2*T9)))**(-1.d0/3) )*
| (T9*(1/(1. + 4.95e-2*T9)))**(2.5d0/3)
c(25)= 'HE4(AN,G)BE9'
dfdt(j,25)= 2.2587e-5*(1.0291 - 1.584*T9 - 1.*T9**2)*
| exp(-1.062*(1/T9))*T9**(-4)*(2.907 + 1.*T9)**(-2)
drdt(j,25)= 1.5126e14*exp(-1.9322e1*(1/T9))*(1/(1. +
| 3.44e-1*T9)) +
| 2.9226e15*exp(-1.9322e1*(1/T9))*(1/(1. + 3.44e-1*T9))*
| (1/T9) - 5.2032e13*T9*exp(-1.9322e1*(1/T9))*
| (1. + 3.44e-1*T9)**(-2)
c(26)= 'HE4(A)BE8'
dfdt(j,26)=
|7.89062e5*exp(-1.0663*(1/T9))*T9**(-7.d0/2) -
| 1.11e6*exp(-1.0663*(1/T9))*T9**(-5.d0/2) +
| 1.90336e11*(2.26063e-4 + 1.16809e-1*T9**(1.d0/3) +
| 3.78911e-2*T9**(2.d0/3) + 1.45506*T9 + 1.*T9**(4.d0/3))*
| exp(-1.349e1*T9**(-1.d0/3) - 1.04123e2*T9**2)*T9**(-4.d0/3)
| - 2.776e9*(1. + 3.1e-2*T9**(1.d0/3) +
| 8.009*T9**(2.d0/3) + 1.732*T9 + 4.9883e1*T9**(4.d0/3) +
| 2.7426e1*T9**(5.d0/3))*
| exp(-1.349e1*T9**(-1.d0/3) - 1.04123e2*T9**2)*T9**(-5.d0/3)
| + 2.37821e13*(3.64618e-2 + 1.13031e-3*T9**(1.d0/3) +
| 2.92022e-1*T9**(2.d0/3) + 6.31518e-2*T9 +
| 1.81882*T9**(4.d0/3) + 1.*T9**(5.d0/3))*
| (2.1593e-2 - 1.*T9**(7.d0/3))*
| exp(-1.349e1*T9**(-1.d0/3) - 1.04123e2*T9**2)*T9**(-2)
if (t9.le.0.005) then
drdt(j,26)=0.0
else
drdt(j,26)=
| -1.49282e10*(7.4e5*exp(-1.0663*(1/T9))*T9**(-3.d0/2) +
| 4.164e9*(1. + 3.1e-2*T9**(1.d0/3) + 8.009*T9**(2.d0/3) +
| 1.732*T9 + 4.9883e1*T9**(4.d0/3) +
| 2.7426e1*T9**(5.d0/3))*
| exp(-1.349e1*T9**(-1.d0/3) - 1.041233e2*T9**2)*
| T9**(-2.d0/3))*exp(1.0663*(1/T9))*(1/Sqrt(T9)) +
| 2.1e10*(7.4e5*exp(-1.0663*(1/T9))*T9**(-3.d0/2) +
| 4.164e9*(1. + 3.1e-2*T9**(1.d0/3) + 8.009*T9**(2.d0/3) +
| 1.732*T9 + 4.9883e1*T9**(4.d0/3) +
| 2.7426e1*T9**(5.d0/3))*
| exp(-1.349e1*T9**(-1.d0/3) - 1.041233e2*T9**2)*
| T9**(-2.d0/3))*exp(1.0663*(1/T9))*Sqrt(T9) +
| 1.4e10*(7.89062e5*exp(-1.0663*(1/T9))*T9**(-7.d0/2) -
| 1.11e6*exp(-1.0663*(1/T9))*T9**(-5.d0/2) +
| 4.164e9*(1.732 + 1.033333e-2*T9**(-2.d0/3) +
| 5.339333*T9**(-1.d0/3) + 6.651067e1*T9**(1.d0/3) +
| 4.571e1*T9**(2.d0/3))*
| exp(-1.349e1*T9**(-1.d0/3) - 1.041233e2*T9**2)*
| T9**(-2.d0/3) - 2.776e9*
| (1. + 3.1e-2*T9**(1.d0/3) + 8.009*T9**(2.d0/3) +
| 1.732*T9 + 4.9883e1*T9**(4.d0/3) +
| 2.7426e1*T9**(5.d0/3))*
| exp(-1.349e1*T9**(-1.d0/3) - 1.041233e2*T9**2)*
| T9**(-5.d0/3) + 4.164e9*
| (4.496667*T9**(-4.d0/3) - 2.082466e2*T9)*
| (1. + 3.1e-2*T9**(1.d0/3) + 8.009*T9**(2.d0/3) +
| 1.732*T9 + 4.9883e1*T9**(4.d0/3) +
| 2.7426e1*T9**(5.d0/3))*
| exp(-1.349e1*T9**(-1.d0/3) - 1.041233e2*T9**2)*
| T9**(-2.d0/3))*exp(1.0663*(1/T9))*T9**(3.d0/2)
endif
c(27)= 'BE8(A,G)C12'
dfdt(j,27)= 4.3373e2*exp(-3.3364*(1/T9))*T9**(-7.d0/2) -
| 1.95e2*exp(-3.3364*(1/T9))*T9**(-5.d0/2) +
| 2.5242e8*(5.9662e-4 + 3.4796e-1*T9**(1.d0/3) +
| 6.4634e-2*T9**(2.d0/3) + 2.5424*T9 + 1.*T9**(4.d0/3))*
| exp(-2.357e1*T9**(-1.d0/3) - 1.8108e1*T9**2)*T9**(-4.d0/3)
| - 1.6733e7*(1. + 1.8e-2*T9**(1.d0/3) +
| 5.249*T9**(2.d0/3) + 6.5e-1*T9 + 1.9176e1*T9**(4.d0/3) +
| 6.034*T9**(5.d0/3))*exp
| (-2.357e1*T9**(-1.d0/3) - 1.8108e1*T9**2)*T9**(-5.d0/3) +
| 5.485e9*(1.6573e-1 + 2.9831e-3*T9**(1.d0/3) +
| 8.699e-1*T9**(2.d0/3) + 1.0772e-1*T9 +
| 3.178*T9**(4.d0/3) + 1.*T9**(5.d0/3))*
| (2.1694e-1 - 1.*T9**(7.d0/3))*
| exp(-2.357e1*T9**(-1.d0/3) - 1.8108e1*T9**2)*T9**(-2)
drdt(j,27)= 3.6759e12*(1.3e2*exp(-3.3364*(1/T9))
| *T9**(-3.d0/2) +
| 2.51e7*(1. + 1.8e-2*T9**(1.d0/3) + 5.249*T9**(2.d0/3) +
| 6.5e-1*T9 + 1.9176e1*T9**(4.d0/3) + 6.034*T9**(5.d0/3)
| )*exp(-2.357e1*T9**(-1.d0/3) - 1.8108e1*T9**2)*
| T9**(-2.d0/3))*exp(-8.5486e1*(1/T9))*(1/Sqrt(T9)) +
| 6.45e10*(1.3e2*exp(-3.3364*(1/T9))*T9**(-3.d0/2) +
| 2.51e7*(1. + 1.8e-2*T9**(1.d0/3) + 5.249*T9**(2.d0/3) +
| 6.5e-1*T9 + 1.9176e1*T9**(4.d0/3) + 6.034*T9**(5.d0/3)
| )*exp(-2.357e1*T9**(-1.d0/3) - 1.8108e1*T9**2)*
| T9**(-2.d0/3))*exp(-8.5486e1*(1/T9))*Sqrt(T9) +
| 4.3e10*(4.3373e2*exp(-3.3364*(1/T9))*T9**(-7.d0/2) -
| 1.95e2*exp(-3.3364*(1/T9))*T9**(-5.d0/2) +
| 2.51e7*(6.5e-1 + 6.e-3*T9**(-2.d0/3) +
| 3.4993*T9**(-1.d0/3) + 2.5568e1*T9**(1.d0/3) +
| 1.0057e1*T9**(2.d0/3))*
| exp(-2.357e1*T9**(-1.d0/3) - 1.8108e1*T9**2)*
| T9**(-2.d0/3) - 1.6733e7*
| (1. + 1.8e-2*T9**(1.d0/3) + 5.249*T9**(2.d0/3) +
| 6.5e-1*T9 + 1.9176e1*T9**(4.d0/3) + 6.034*T9**(5.d0/3)
| )*exp(-2.357e1*T9**(-1.d0/3) - 1.8108e1*T9**2)*
| T9**(-5.d0/3) + 2.51e7*
| (7.8567*T9**(-4.d0/3) - 3.6215e1*T9)*
| (1. + 1.8e-2*T9**(1.d0/3) + 5.249*T9**(2.d0/3) +
| 6.5e-1*T9 + 1.9176e1*T9**(4.d0/3) + 6.034*T9**(5.d0/3)
| )*exp(-2.357e1*T9**(-1.d0/3) - 1.8108e1*T9**2)*
| T9**(-2.d0/3))*exp(-8.5486e1*(1/T9))*T9**(3.d0/2)
c(28)= 'HE4(2A,G)C12'
if (T9.le.0.08d0) then
dfdt(j,28)= 3.3495e-6*rc28*exp(-2.4811e1*(1/T9))
| *T9**(-7.d0/2) -
| 2.025e-7*rc28*exp(-2.4811e1*(1/T9))*T9**(-5.d0/2) +
| 2.9e-16*(1.5459e-5*exp(-5.9221e-6*(1/T9)**3.263)*
| (1/(1. + 4.*exp(-6.0563e14*T9**9.227)))*
| (1/T9)**2.263*T9**(-2) +
| 4.4705e15*(1. + 4.*exp(-5.9221e-6*(1/T9)**3.263))*
| exp(-6.0563e14*T9**9.227)*
| (1. + 4.*exp(-6.0563e14*T9**9.227))**(-2)*
| T9**8.227)*(1.3e2*exp(-3.3364*(1/T9))*
| T9**(-3.d0/2) + 2.51e7*
| (1. + 1.8e-2*T9**(1.d0/3) + 5.249*T9**(2.d0/3) +
| 6.5e-1*T9 + 1.9176e1*T9**(4.d0/3) + 6.034*T9**(5.d0/3)
| )*exp(-2.357e1*T9**(-1.d0/3) - 1.8108e1*T9**2)*
| T9**(-2.d0/3))*(7.4e5*exp(-1.0663*(1/T9))*
| T9**(-3.d0/2) + 4.16e9*
| (1. + 3.1e-2*T9**(1.d0/3) + 8.009*T9**(2.d0/3) +
| 1.732*T9 + 4.9883e1*T9**(4.d0/3) +
| 2.7426e1*T9**(5.d0/3))*
| exp(-1.349e1*T9**(-1.d0/3) - 1.0412e2*T9**2)*
| T9**(-2.d0/3)) + 2.9e-16*
| (1.e-2 + 2.e-1*(1. +
| 4.*exp(-5.9221e-6*(1/T9)**3.263))*
| (1/(1. + 4.*exp(-6.0563e14*T9**9.227))))*
| (4.3373e2*exp(-3.3364*(1/T9))*T9**(-7.d0/2) -
| 1.95e2*exp(-3.3364*(1/T9))*T9**(-5.d0/2) +
| 2.51e7*(6.5e-1 + 6.e-3*T9**(-2.d0/3) +
| 3.4993*T9**(-1.d0/3) + 2.5568e1*T9**(1.d0/3) +
| 1.0057e1*T9**(2.d0/3))*
| exp(-2.357e1*T9**(-1.d0/3) - 1.8108e1*T9**2)*
| T9**(-2.d0/3) - 1.6733e7*
| (1. + 1.8e-2*T9**(1.d0/3) + 5.249*T9**(2.d0/3) +
| 6.5e-1*T9 + 1.9176e1*T9**(4.d0/3) + 6.034*T9**(5.d0/3)
| )*exp(-2.357e1*T9**(-1.d0/3) - 1.8108e1*T9**2)*
| T9**(-5.d0/3) + 2.51e7*
| (7.8567*T9**(-4.d0/3) - 3.6215e1*T9)*
| (1. + 1.8e-2*T9**(1.d0/3) + 5.249*T9**(2.d0/3) +
| 6.5e-1*T9 + 1.9176e1*T9**(4.d0/3) + 6.034*T9**(5.d0/3)
| )*exp(-2.357e1*T9**(-1.d0/3) - 1.8108e1*T9**2)*
| T9**(-2.d0/3))*(7.4e5*exp(-1.0663*(1/T9))*
| T9**(-3.d0/2) + 4.16e9*
| (1. + 3.1e-2*T9**(1.d0/3) + 8.009*T9**(2.d0/3) +
| 1.732*T9 + 4.9883e1*T9**(4.d0/3) +
| 2.7426e1*T9**(5.d0/3))*
| exp(-1.349e1*T9**(-1.d0/3) - 1.0412e2*T9**2)*
| T9**(-2.d0/3)) + 2.9e-16*
| (1.e-2 + 2.e-1*(1. +
| 4.*exp(-5.9221e-6*(1/T9)**3.263))*
| (1/(1. + 4.*exp(-6.0563e14*T9**9.227))))*
| (1.3e2*exp(-3.3364*(1/T9))*T9**(-3.d0/2) +
| 2.51e7*(1. + 1.8e-2*T9**(1.d0/3) + 5.249*T9**(2.d0/3) +
| 6.5e-1*T9 + 1.9176e1*T9**(4.d0/3) + 6.034*T9**(5.d0/3)
| )*exp(-2.357e1*T9**(-1.d0/3) - 1.8108e1*T9**2)*
| T9**(-2.d0/3))*(7.8906e5*exp(-1.0663*(1/T9))*
| T9**(-7.d0/2) - 1.11e6*exp(-1.0663*(1/T9))*
| T9**(-5.d0/2) + 4.16e9*
| (1.732 + 1.0333e-2*T9**(-2.d0/3) +
| 5.3393*T9**(-1.d0/3) + 6.6511e1*T9**(1.d0/3) +
| 4.571e1*T9**(2.d0/3))*
| exp(-1.349e1*T9**(-1.d0/3) - 1.0412e2*T9**2)*
| T9**(-2.d0/3) - 2.7733e9*
| (1. + 3.1e-2*T9**(1.d0/3) + 8.009*T9**(2.d0/3) +
| 1.732*T9 + 4.9883e1*T9**(4.d0/3) +
| 2.7426e1*T9**(5.d0/3))*
| exp(-1.349e1*T9**(-1.d0/3) - 1.0412e2*T9**2)*
| T9**(-5.d0/3) + 4.16e9*
| (4.4967*T9**(-4.d0/3) - 2.0825e2*T9)*
| (1. + 3.1e-2*T9**(1.d0/3) + 8.009*T9**(2.d0/3) +
| 1.732*T9 + 4.9883e1*T9**(4.d0/3) +
| 2.7426e1*T9**(5.d0/3))*
| exp(-1.349e1*T9**(-1.d0/3) - 1.0412e2*T9**2)*
| T9**(-2.d0/3))
drdt(j,28)= 1.6884e22*T9*(1.35e-7*rc28*
| exp(-2.4811e1*(1/T9))*
| T9**(-3.d0/2) + 2.9e-16*
| (1.e-2 + 2.e-1*(1. +
| 4.*exp(-5.9221e-6*(1/T9)**3.263))*
| (1/(1. + 4.*exp(-6.0563e14*T9**9.227))))*
| (1.3e2*exp(-3.3364*(1/T9))*T9**(-3.d0/2) +
| 2.51e7*(1. + 1.8e-2*T9**(1.d0/3) +
| 5.249*T9**(2.d0/3) + 6.5e-1*T9 +
| 1.9176e1*T9**(4.d0/3) + 6.034*T9**(5.d0/3))*
| exp(-2.357e1*T9**(-1.d0/3) - 1.8108e1*T9**2)*
| T9**(-2.d0/3))*(7.4e5*exp(-1.0663*(1/T9))*
| T9**(-3.d0/2) + 4.16e9*
| (1. + 3.1e-2*T9**(1.d0/3) + 8.009*T9**(2.d0/3) +
| 1.732*T9 + 4.9883e1*T9**(4.d0/3) +
| 2.7426e1*T9**(5.d0/3))*
| exp(-1.349e1*T9**(-1.d0/3) - 1.0412e2*T9**2)*
| T9**(-2.d0/3)))*exp(-8.442e1*(1/T9)) +
| 6.e20*(1.35e-7*rc28*exp(-2.4811e1*(1/T9))*
| T9**(-3.d0/2) + 2.9e-16*
| (1.e-2 + 2.e-1*(1. +
| 4.*exp(-5.9221e-6*(1/T9)**3.263))*
| (1/(1. + 4.*exp(-6.0563e14*T9**9.227))))*
| (1.3e2*exp(-3.3364*(1/T9))*T9**(-3.d0/2) +
| 2.51e7*(1. + 1.8e-2*T9**(1.d0/3) +
| 5.249*T9**(2.d0/3) + 6.5e-1*T9 +
| 1.9176e1*T9**(4.d0/3) + 6.034*T9**(5.d0/3))*
| exp(-2.357e1*T9**(-1.d0/3) - 1.8108e1*T9**2)*
| T9**(-2.d0/3))*(7.4e5*exp(-1.0663*(1/T9))*
| T9**(-3.d0/2) + 4.16e9*
| (1. + 3.1e-2*T9**(1.d0/3) + 8.009*T9**(2.d0/3) +
| 1.732*T9 + 4.9883e1*T9**(4.d0/3) +
| 2.7426e1*T9**(5.d0/3))*
| exp(-1.349e1*T9**(-1.d0/3) - 1.0412e2*T9**2)*
| T9**(-2.d0/3)))*exp(-8.442e1*(1/T9))*T9**2 +
| 2.e20*(3.3495e-6*rc28*exp(-2.4811e1*(1/T9))*
| T9**(-7.d0/2) - 2.025e-7*rc28*exp(-2.4811e1*(1/T9))*
| T9**(-5.d0/2) + 2.9e-16*
| (1.5459e-5*exp(-5.9221e-6*(1/T9)**3.263)*
| (1/(1. + 4.*exp(-6.0563e14*T9**9.227)))*
| (1/T9)**2.263*T9**(-2) +
| 4.4705e15*(1. +
| 4.*exp(-5.9221e-6*(1/T9)**3.263))*
| exp(-6.0563e14*T9**9.227)*
| (1. + 4.*exp(-6.0563e14*T9**9.227))**(-2)*
| T9**8.227)*(1.3e2*exp(-3.3364*(1/T9))*
| T9**(-3.d0/2) + 2.51e7*
| (1. + 1.8e-2*T9**(1.d0/3) + 5.249*T9**(2.d0/3) +
| 6.5e-1*T9 + 1.9176e1*T9**(4.d0/3) +
| 6.034*T9**(5.d0/3))*
| exp(-2.357e1*T9**(-1.d0/3) - 1.8108e1*T9**2)*
| T9**(-2.d0/3))*(7.4e5*exp(-1.0663*(1/T9))*
| T9**(-3.d0/2) + 4.16e9*
| (1. + 3.1e-2*T9**(1.d0/3) + 8.009*T9**(2.d0/3) +
| 1.732*T9 + 4.9883e1*T9**(4.d0/3) +
| 2.7426e1*T9**(5.d0/3))*
| exp(-1.349e1*T9**(-1.d0/3) - 1.0412e2*T9**2)*
| T9**(-2.d0/3)) +
| 2.9e-16*(1.e-2 + 2.e-1*
| (1. + 4.*exp(-5.9221e-6*(1/T9)**3.263))*
| (1/(1. + 4.*exp(-6.0563e14*T9**9.227))))*
| (4.3373e2*exp(-3.3364*(1/T9))*T9**(-7.d0/2) -
| 1.95e2*exp(-3.3364*(1/T9))*T9**(-5.d0/2) +
| 2.51e7*(6.5e-1 + 6.e-3*T9**(-2.d0/3) +
| 3.4993*T9**(-1.d0/3) + 2.5568e1*T9**(1.d0/3) +
| 1.0057e1*T9**(2.d0/3))*
| exp(-2.357e1*T9**(-1.d0/3) - 1.8108e1*T9**2)*
| T9**(-2.d0/3) - 1.6733e7*
| (1. + 1.8e-2*T9**(1.d0/3) + 5.249*T9**(2.d0/3) +
| 6.5e-1*T9 + 1.9176e1*T9**(4.d0/3) +
| 6.034*T9**(5.d0/3))*
| exp(-2.357e1*T9**(-1.d0/3) - 1.8108e1*T9**2)*
| T9**(-5.d0/3) + 2.51e7*
| (7.8567*T9**(-4.d0/3) - 3.6215e1*T9)*
| (1. + 1.8e-2*T9**(1.d0/3) + 5.249*T9**(2.d0/3) +
| 6.5e-1*T9 + 1.9176e1*T9**(4.d0/3) +
| 6.034*T9**(5.d0/3))*
| exp(-2.357e1*T9**(-1.d0/3) - 1.8108e1*T9**2)*
| T9**(-2.d0/3))*(7.4e5*exp(-1.0663*(1/T9))*
| T9**(-3.d0/2) + 4.16e9*
| (1. + 3.1e-2*T9**(1.d0/3) + 8.009*T9**(2.d0/3) +
| 1.732*T9 + 4.9883e1*T9**(4.d0/3) +
| 2.7426e1*T9**(5.d0/3))*
| exp(-1.349e1*T9**(-1.d0/3) - 1.0412e2*T9**2)*
| T9**(-2.d0/3)) +
| 2.9e-16*(1.e-2 + 2.e-1*
| (1. + 4.*exp(-5.9221e-6*(1/T9)**3.263))*
| (1/(1. + 4.*exp(-6.0563e14*T9**9.227))))*
| (1.3e2*exp(-3.3364*(1/T9))*T9**(-3.d0/2) +
| 2.51e7*(1. + 1.8e-2*T9**(1.d0/3) +
| 5.249*T9**(2.d0/3) + 6.5e-1*T9 +
| 1.9176e1*T9**(4.d0/3) + 6.034*T9**(5.d0/3))*
| exp(-2.357e1*T9**(-1.d0/3) - 1.8108e1*T9**2)*
| T9**(-2.d0/3))*(7.8906e5*exp(-1.0663*(1/T9))*
| T9**(-7.d0/2) - 1.11e6*exp(-1.0663*(1/T9))*
| T9**(-5.d0/2) + 4.16e9*
| (1.732 + 1.0333e-2*T9**(-2.d0/3) +
| 5.3393*T9**(-1.d0/3) + 6.6511e1*T9**(1.d0/3) +
| 4.571e1*T9**(2.d0/3))*
| exp(-1.349e1*T9**(-1.d0/3) - 1.0412e2*T9**2)*
| T9**(-2.d0/3) - 2.7733e9*
| (1. + 3.1e-2*T9**(1.d0/3) + 8.009*T9**(2.d0/3) +
| 1.732*T9 + 4.9883e1*T9**(4.d0/3) +
| 2.7426e1*T9**(5.d0/3))*
| exp(-1.349e1*T9**(-1.d0/3) - 1.0412e2*T9**2)*
| T9**(-5.d0/3) + 4.16e9*
| (4.4967*T9**(-4.d0/3) - 2.0825e2*T9)*
| (1. + 3.1e-2*T9**(1.d0/3) + 8.009*T9**(2.d0/3) +
| 1.732*T9 + 4.9883e1*T9**(4.d0/3) +
| 2.7426e1*T9**(5.d0/3))*
| exp(-1.349e1*T9**(-1.d0/3) - 1.0412e2*T9**2)*
| T9**(-2.d0/3)))*exp(-8.442e1*(1/T9))*T9**3
else if (T9.gt.0.08d0) then
dfdt(j,28)=
| 8.37e-8*(1.46757*exp(2.04083e1*(1/T9)) -
| 1.*T9*exp(2.04083e1*(1/T9)) +
| 4.00177e1*rc28*T9**(3.d0/2) - 2.41935*rc28*T9**(5.d0/2))*
| exp(-2.4811e1*(1/T9))*T9**(-5)
drdt(j,28)= 1.6884e22*(2.79e-8*exp(-4.4027*(1/T9))
| *T9**(-3) +
| 1.35e-7*rc28*exp(-2.4811e1*(1/T9))*T9**(-3.d0/2))*T9*
| exp(-8.442e1*(1/T9)) +
| 6.e20*(2.79e-8*exp(-4.4027*(1/T9))*T9**(-3) +
| 1.35e-7*rc28*exp(-2.4811e1*(1/T9))*T9**(-3.d0/2))*
| exp(-8.442e1*(1/T9))*T9**2 +
| 2.e20*(1.2284e-7*exp(-4.4027*(1/T9))*T9**(-5) -
| 8.37e-8*exp(-4.4027*(1/T9))*T9**(-4) +
| 3.3495e-6*rc28*exp(-2.4811e1*(1/T9))*T9**(-7.d0/2) -
| 2.025e-7*rc28*exp(-2.4811e1*(1/T9))*T9**(-5.d0/2))*
| exp(-8.442e1*(1/T9))*T9**3
endif
c(29)= 'LI6(P,G)BE7'
dfdt(j,29)= -1.0035e6*exp
| (-8.413*(T9*(1/
| (1. - 9.69e-2*T9 +
| 2.84e-2*(1. - 9.69e-2*T9)**(-2.d0/3)*
| T9**(5.d0/3))))**(-1.d0/3) )*T9**(-5.d0/2)*
| (T9*(1/(1. - 9.69e-2*T9 +
| 2.84e-2*(1. - 9.69e-2*T9)**(-2.d0/3)*
| T9**(5.d0/3))))**(2.5d0/3) +
| 1.8761e6*(-(T9*(-9.69e-2 +
| 4.7333e-2*(1. - 9.69e-2*T9)**(-2.d0/3)*
| T9**(2.d0/3) +
| 1.8346e-3*(1. - 9.69e-2*T9)**(-5.d0/3)*T9**(5.d0/3)
| )*(1. - 9.69e-2*T9 +
| 2.84e-2*(1. - 9.69e-2*T9)**(-2.d0/3)*
| T9**(5.d0/3))**(-2)) +
| 1/(1. - 9.69e-2*T9 +
| 2.84e-2*(1. - 9.69e-2*T9)**(-2.d0/3)*T9**(5.d0/3))
| )*exp(-8.413*(T9*(1/
| (1. - 9.69e-2*T9 +
| 2.84e-2*(1. - 9.69e-2*T9)**(-2.d0/3)*
| T9**(5.d0/3))))**(-1.d0/3) )*T9**(-3.d0/2)*
| (T9*(1/(1. - 9.69e-2*T9 +
| 2.84e-2*(1. - 9.69e-2*T9)**(-2.d0/3)*
| T9**(5.d0/3))))**(-5.e-1) +
| 5.575e5*(-(T9*(-9.69e-2 +
| 4.7333e-2*(1. - 9.69e-2*T9)**(-2.d0/3)*
| T9**(2.d0/3) +
| 1.8346e-3*(1. - 9.69e-2*T9)**(-5.d0/3)*T9**(5.d0/3)
| )*(1. - 9.69e-2*T9 +
| 2.84e-2*(1. - 9.69e-2*T9)**(-2.d0/3)*
| T9**(5.d0/3))**(-2)) +
| 1/(1. - 9.69e-2*T9 +
| 2.84e-2*(1. - 9.69e-2*T9)**(-2.d0/3)*T9**(5.d0/3))
| )*exp(-8.413*(T9*(1/
| (1. - 9.69e-2*T9 +
| 2.84e-2*(1. - 9.69e-2*T9)**(-2.d0/3)*
| T9**(5.d0/3))))**(-1.d0/3))*T9**(-3.d0/2)*
| (T9*(1/(1. - 9.69e-2*T9 +
| 2.84e-2*(1. - 9.69e-2*T9)**(-2.d0/3)*
| T9**(5.d0/3))))**(-.5d0/3)
drdt(j,29)= 6.6342e15*(-(T9*(-9.69e-2 +
| 4.7333e-2*(1. - 9.69e-2*T9)**(-2.d0/3)*
| T9**(2.d0/3) +
| 1.8346e-3*(1. - 9.69e-2*T9)**(-5.d0/3)*T9**(5.d0/3)
| )*(1. - 9.69e-2*T9 +
| 2.84e-2*(1. - 9.69e-2*T9)**(-2.d0/3)*
| T9**(5.d0/3))**(-2)) +
| 1/(1. - 9.69e-2*T9 +
| 2.84e-2*(1. - 9.69e-2*T9)**(-2.d0/3)*T9**(5.d0/3))
| )*exp(-6.5054e1*(1/T9) -
| 8.413*(T9*(1/
| (1. - 9.69e-2*T9 +
| 2.84e-2*(1. - 9.69e-2*T9)**(-2.d0/3)*
| T9**(5.d0/3))))**(-1.d0/3))*
| (T9*(1/(1. - 9.69e-2*T9 +
| 2.84e-2*(1. - 9.69e-2*T9)**(-2.d0/3)*
| T9**(5.d0/3))))**(-.5d0/3) +
| 7.9611e15*(6.5054e1*T9**(-2) +
| 2.8043*(-(T9*(-9.69e-2 +
| 4.7333e-2*(1. - 9.69e-2*T9)**(-2.d0/3)*
| T9**(2.d0/3) +
| 1.8346e-3*(1. - 9.69e-2*T9)**(-5.d0/3)*
| T9**(5.d0/3))*
| (1. - 9.69e-2*T9 +
| 2.84e-2*(1. - 9.69e-2*T9)**(-2.d0/3)*
| T9**(5.d0/3))**(-2)) +
| 1/(1. - 9.69e-2*T9 +
| 2.84e-2*(1. - 9.69e-2*T9)**(-2.d0/3)*
| T9**(5.d0/3)))*
| (T9*(1/
| (1. - 9.69e-2*T9 +
| 2.84e-2*(1. - 9.69e-2*T9)**(-2.d0/3)*
| T9**(5.d0/3))))**(-4.d0/3))*
| exp(-6.5054e1*(1/T9) -
| 8.413*(T9*(1/
| (1. - 9.69e-2*T9 +
| 2.84e-2*(1. - 9.69e-2*T9)**(-2.d0/3)*
| T9**(5.d0/3))))**(-1.d0/3))*
| (T9*(1/(1. - 9.69e-2*T9 +
| 2.84e-2*(1. - 9.69e-2*T9)**(-2.d0/3)*
| T9**(5.d0/3))))**(2.5d0/3)
c(30)= 'LI6(P,HE3)HE4'
dfdt(j,30)= 2.3625e11*exp(-1.7763e1*(1/T9))*T9**(-7.d0/2) +
| 2.8148e10*exp(-2.182e1*(1/T9))*T9**(-3) -
| 1.995e10*exp(-1.7763e1*(1/T9))*T9**(-5.d0/2) -
| 1.29e9*exp(-2.182e1*(1/T9))*T9**(-2) +
| 3.1083e8*(2. - 4.88*T9**(1.d0/3) - 2.52*T9**(2.d0/3) +
| 9.6e-1*T9 + 1.*T9**(4.d0/3))*
| exp(-8.413*T9**(-1.d0/3) - 3.3058e-2*T9**2)*T9**(-4.d0/3)
| + 1.2433e8*(-2.e2 - 1.e1*T9**(1.d0/3) +
| 1.22e1*T9**(2.d0/3) + 4.2*T9 - 1.2*T9**(4.d0/3) -
| 1.*T9**(5.d0/3))*exp
| (-8.413*T9**(-1.d0/3) - 3.3058e-2*T9**2)*T9**(-5.d0/3) +
| 1.2331e7*(2.e2 + 1.e1*T9**(1.d0/3) - 1.22e1*T9**(2.d0/3) -
| 4.2*T9 + 1.2*T9**(4.d0/3) + 1.*T9**(5.d0/3))*
| (4.2416e1 - 1.*T9**(7.d0/3))*
| exp(-8.413*T9**(-1.d0/3) - 3.3058e-2*T9**2)*T9**(-2)
drdt(j,30)= 4.9895e1*(1.33e10*exp(-1.7763e1*(1/T9))
| *T9**(-3.d0/2) +
| 1.29e9*exp(-2.182e1*(1/T9))*(1/T9) +
| 3.73e10*(1. + 5.e-2*T9**(1.d0/3) - 6.1e-2*T9**(2.d0/3) -
| 2.1e-2*T9 + 6.e-3*T9**(4.d0/3) + 5.e-3*T9**(5.d0/3))*
| exp(-8.413*T9**(-1.d0/3) - 3.3058e-2*T9**2)*
| T9**(-2.d0/3))*exp(-4.6631e1*(1/T9))*T9**(-2) +
| 1.07*(2.3625e11*exp(-1.7763e1*(1/T9))*T9**(-7.d0/2) +
| 2.8148e10*exp(-2.182e1*(1/T9))*T9**(-3) -
| 1.995e10*exp(-1.7763e1*(1/T9))*T9**(-5.d0/2) -
| 1.29e9*exp(-2.182e1*(1/T9))*T9**(-2) +
| 3.73e10*(-2.1e-2 + 1.6667e-2*T9**(-2.d0/3) -
| 4.0667e-2*T9**(-1.d0/3) + 8.e-3*T9**(1.d0/3) +
| 8.3333e-3*T9**(2.d0/3))*
| exp(-8.413*T9**(-1.d0/3) - 3.3058e-2*T9**2)*
| T9**(-2.d0/3) - 2.4867e10*
| (1. + 5.e-2*T9**(1.d0/3) - 6.1e-2*T9**(2.d0/3) -
| 2.1e-2*T9 + 6.e-3*T9**(4.d0/3) + 5.e-3*T9**(5.d0/3))*
| exp(-8.413*T9**(-1.d0/3) - 3.3058e-2*T9**2)*
| T9**(-5.d0/3) + 3.73e10*
| (2.8043*T9**(-4.d0/3) - 6.6116e-2*T9)*
| (1. + 5.e-2*T9**(1.d0/3) - 6.1e-2*T9**(2.d0/3) -
| 2.1e-2*T9 + 6.e-3*T9**(4.d0/3) + 5.e-3*T9**(5.d0/3))*
| exp(-8.413*T9**(-1.d0/3) - 3.3058e-2*T9**2)*
| T9**(-2.d0/3))*exp(-4.6631e1*(1/T9))
|
c(31)= 'LI6(A,G)B10'
dfdt(j,31)= 6.6544e3*exp(-3.484*(1/T9))*T9**(-7.d0/2) +
| 7.3417e4*exp(-7.269*(1/T9))*T9**(-3) -
| 2.865e3*exp(-3.484*(1/T9))*T9**(-5.d0/2) -
| 1.01e4*exp(-7.269*(1/T9))*T9**(-2) +
| 5.8802e6*(5.0633e-3 + 7.0886e-1*T9**(1.d0/3) +
| 1.6502e-1*T9**(2.d0/3) + 2.0253*T9 + 1.*T9**(4.d0/3))*
| exp(-1.879e1*T9**(-1.d0/3) - 5.6874e-1*T9**2)*T9**(-4.d0/3)
| - 2.7067e6*(1. + 2.2e-2*T9**(1.d0/3) +
| 1.54*T9**(2.d0/3) + 2.39e-1*T9 + 2.2*T9**(4.d0/3) +
| 8.69e-1*T9**(5.d0/3))*
| exp(-1.879e1*T9**(-1.d0/3) - 5.6874e-1*T9**2)*T9**(-5.d0/3)
| + 4.0132e6*(1.1507 + 2.5316e-2*T9**(1.d0/3) +
| 1.7722*T9**(2.d0/3) + 2.7503e-1*T9 +
| 2.5316*T9**(4.d0/3) + 1.*T9**(5.d0/3))*
| (5.5063 - 1.*T9**(7.d0/3))*
| exp(-1.879e1*T9**(-1.d0/3) - 5.6874e-1*T9**2)*T9**(-2)
drdt(j,31)= 8.177e11*(1.91e3*exp(-3.484*(1/T9))
| *T9**(-3.d0/2) +
| 1.01e4*exp(-7.269*(1/T9))*(1/T9) +
| 4.06e6*(1. + 2.2e-2*T9**(1.d0/3) + 1.54*T9**(2.d0/3) +
| 2.39e-1*T9 + 2.2*T9**(4.d0/3) + 8.69e-1*T9**(5.d0/3))*
| exp(-1.879e1*T9**(-1.d0/3) - 5.6874e-1*T9**2)*
| T9**(-2.d0/3))*exp(-5.1753e1*(1/T9))*(1/Sqrt(T9)) +
| 2.37e10*(1.91e3*exp(-3.484*(1/T9))*T9**(-3.d0/2) +
| 1.01e4*exp(-7.269*(1/T9))*(1/T9) +
| 4.06e6*(1. + 2.2e-2*T9**(1.d0/3) + 1.54*T9**(2.d0/3) +
| 2.39e-1*T9 + 2.2*T9**(4.d0/3) + 8.69e-1*T9**(5.d0/3))*
| exp(-1.879e1*T9**(-1.d0/3) - 5.6874e-1*T9**2)*
| T9**(-2.d0/3))*exp(-5.1753e1*(1/T9))*Sqrt(T9) +
| 1.58e10*(6.6544e3*exp(-3.484*(1/T9))*T9**(-7.d0/2) +
| 7.3417e4*exp(-7.269*(1/T9))*T9**(-3) -
| 2.865e3*exp(-3.484*(1/T9))*T9**(-5.d0/2) -
| 1.01e4*exp(-7.269*(1/T9))*T9**(-2) +
| 4.06e6*(2.39e-1 + 7.3333e-3*T9**(-2.d0/3) +
| 1.0267*T9**(-1.d0/3) + 2.9333*T9**(1.d0/3) +
| 1.4483*T9**(2.d0/3))*
| exp(-1.879e1*T9**(-1.d0/3) - 5.6874e-1*T9**2)*
| T9**(-2.d0/3) - 2.7067e6*
| (1. + 2.2e-2*T9**(1.d0/3) + 1.54*T9**(2.d0/3) +
| 2.39e-1*T9 + 2.2*T9**(4.d0/3) + 8.69e-1*T9**(5.d0/3))*
| exp(-1.879e1*T9**(-1.d0/3) - 5.6874e-1*T9**2)*
| T9**(-5.d0/3) + 4.06e6*
| (6.2633*T9**(-4.d0/3) - 1.1375*T9)*
| (1. + 2.2e-2*T9**(1.d0/3) + 1.54*T9**(2.d0/3) +
| 2.39e-1*T9 + 2.2*T9**(4.d0/3) + 8.69e-1*T9**(5.d0/3))*
| exp(-1.879e1*T9**(-1.d0/3) - 5.6874e-1*T9**2)*
| T9**(-2.d0/3))*exp(-5.1753e1*(1/T9))*T9**(3.d0/2)
c(32)= 'LI7(P,N)BE7'
dfdt(j,32)= 1.79e11*exp(-2.2832e1*(1/T9))*T9**(-7.d0/2) -
| 1.176e10*exp(-2.2832e1*(1/T9))*T9**(-5.d0/2) +
| 2.0033e9*(4.9051e1 - 1.*T9**(4.d0/3))*
| exp(-1.9081e1*(1/T9) - 1.167*T9**(1.d0/3))*T9**(-2)
drdt(j,32)= 2.9349e10*exp(-3.751*(1/T9))*T9**(-7.d0/2) -
| 1.1736e10*exp(-3.751*(1/T9))*T9**(-5.d0/2) -
| 2.0033e9*exp(-1.167*T9**(1.d0/3))*T9**(-2.d0/3)
c(33)= 'LI7(P,G)BE8'
dfdt(j,33)= 6.9409e6*exp(-4.478*(1/T9))*T9**(-7.d0/2) -
| 2.325e6*exp(-4.478*(1/T9))*T9**(-5.d0/2) +
| 8.008e5*(3.1818e-3 + 3.2442e-1*T9**(1.d0/3) +
| 1.6753e-1*T9**(2.d0/3) + 9.1377e-1*T9 + 1.*T9**(4.d0/3))*
| exp(-8.472*T9**(-1.d0/3) - 3.4765e-1*T9**2)*T9**(-4.d0/3)
| - 1.04e5*(1. + 4.9e-2*T9**(1.d0/3) + 2.498*T9**(2.d0/3) +
| 8.6e-1*T9 + 3.518*T9**(4.d0/3) + 3.08*T9**(5.d0/3))*
| exp(-8.472*T9**(-1.d0/3) - 3.4765e-1*T9**2)*T9**(-5.d0/3)
| + 3.3408e5*(3.2468e-1 + 1.5909e-2*T9**(1.d0/3) +
| 8.1104e-1*T9**(2.d0/3) + 2.7922e-1*T9 +
| 1.1422*T9**(4.d0/3) + 1.*T9**(5.d0/3))*
| (4.0615 - 1.*T9**(7.d0/3))*
| exp(-8.472*T9**(-1.d0/3) - 3.4765e-1*T9**2)*T9**(-2)
drdt(j,33)= 1.3115e13*(1.55e6*exp(-4.478*(1/T9))
| *T9**(-3.d0/2) +
| 1.56e5*(1. + 4.9e-2*T9**(1.d0/3) + 2.498*T9**(2.d0/3) +
| 8.6e-1*T9 + 3.518*T9**(4.d0/3) + 3.08*T9**(5.d0/3))*
| exp(-8.472*T9**(-1.d0/3) - 3.4765e-1*T9**2)*
| T9**(-2.d0/3))*exp(-2.0022e2*(1/T9))*(1/Sqrt(T9)) +
| 9.825e10*(1.55e6*exp(-4.478*(1/T9))*T9**(-3.d0/2) +
| 1.56e5*(1. + 4.9e-2*T9**(1.d0/3) + 2.498*T9**(2.d0/3) +
| 8.6e-1*T9 + 3.518*T9**(4.d0/3) + 3.08*T9**(5.d0/3))*
| exp(-8.472*T9**(-1.d0/3) - 3.4765e-1*T9**2)*
| T9**(-2.d0/3))*exp(-2.0022e2*(1/T9))*Sqrt(T9) +
| 6.55e10*(6.9409e6*exp(-4.478*(1/T9))*T9**(-7.d0/2) -
| 2.325e6*exp(-4.478*(1/T9))*T9**(-5.d0/2) +
| 1.56e5*(8.6e-1 + 1.6333e-2*T9**(-2.d0/3) +
| 1.6653*T9**(-1.d0/3) + 4.6907*T9**(1.d0/3) +
| 5.1333*T9**(2.d0/3))*
| exp(-8.472*T9**(-1.d0/3) - 3.4765e-1*T9**2)*
| T9**(-2.d0/3) - 1.04e5*
| (1. + 4.9e-2*T9**(1.d0/3) + 2.498*T9**(2.d0/3) +
| 8.6e-1*T9 + 3.518*T9**(4.d0/3) + 3.08*T9**(5.d0/3))*
| exp(-8.472*T9**(-1.d0/3) - 3.4765e-1*T9**2)*
| T9**(-5.d0/3) + 1.56e5*
| (2.824*T9**(-4.d0/3) - 6.9531e-1*T9)*
| (1. + 4.9e-2*T9**(1.d0/3) + 2.498*T9**(2.d0/3) +
| 8.6e-1*T9 + 3.518*T9**(4.d0/3) + 3.08*T9**(5.d0/3))*
| exp(-8.472*T9**(-1.d0/3) - 3.4765e-1*T9**2)*
| T9**(-2.d0/3))*exp(-2.0022e2*(1/T9))*T9**(3.d0/2)
c(34)= 'LI7(P,A)HE4'
dfdt(j,34)= 3.2269e11*exp(-3.0442e1*(1/T9))*T9**(-7.d0/2) -
| 1.59e10*exp(-3.0442e1*(1/T9))*T9**(-5.d0/2) +
| 3.0951e9*exp(-8.472*T9**(-1.d0/3))*T9**(-2) -
| 7.3067e8*exp(-8.472*T9**(-1.d0/3))*T9**(-5.d0/3) +
| 7.245e8*exp
| (-8.472*(T9*(1/(1. + 7.59e-1*T9)))**(-1.d0/3))*
| T9**(-5.d0/2)*(T9*(1/(1. + 7.59e-1*T9)))**(2.5d0/3) -
| 3.0287e-7*(7.8176e15 + 1.*T9)*
| exp(-8.472*(T9*(1/(1. + 7.59e-1*T9)))**(-1.d0/3))*
| T9**(-3.d0/2)*(T9*(1/(1. + 7.59e-1*T9)))**(-5.e-1)*
| (1.3175 + 1.*T9)**(-2) -
| 8.9373e-8*(7.8176e15 + 1.*T9)*
| exp(-8.472*(T9*(1/(1. + 7.59e-1*T9)))**(-1.d0/3))*
| T9**(-3.d0/2)*(T9*(1/(1. + 7.59e-1*T9)))**(-.5d0/3)*
| (1.3175 + 1.*T9)**(-2)
drdt(j,34)= 9.4405e2*(1.06e10*exp(-3.0442e1*(1/T9))
| *T9**(-3.d0/2) +
| 1.096e9*exp(-8.472*T9**(-1.d0/3))*T9**(-2.d0/3) -
| 4.83e8*exp
| (-8.472*(T9*(1/(1. + 7.59e-1*T9)))**(-1.d0/3))*
| T9**(-3.d0/2)*(T9*(1/(1. + 7.59e-1*T9)))**(2.5d0/3))*
| exp(-2.0129e2*(1/T9))*T9**(-2) +
| 4.69*(3.2269e11*exp(-3.0442e1*(1/T9))*T9**(-7.d0/2) -
| 1.59e10*exp(-3.0442e1*(1/T9))*T9**(-5.d0/2) +
| 3.0951e9*exp(-8.472*T9**(-1.d0/3))*T9**(-2) -
| 7.3067e8*exp(-8.472*T9**(-1.d0/3))*T9**(-5.d0/3) +
| 7.245e8*exp
| (-8.472*(T9*(1/(1. + 7.59e-1*T9)))**(-1.d0/3))*
| T9**(-5.d0/2)*(T9*(1/(1. + 7.59e-1*T9)))**(2.5d0/3)
| - 1.364e9*(1/(1. + 7.59e-1*T9) -
| 7.59e-1*T9*(1. + 7.59e-1*T9)**(-2))*
| exp(-8.472*(T9*(1/(1. + 7.59e-1*T9)))**
| (-1.d0/3) )*T9**(-3.d0/2)*
| (T9*(1/(1. + 7.59e-1*T9)))**(-5.e-1) -
| 4.025e8*(1/(1. + 7.59e-1*T9) -
| 7.59e-1*T9*(1. + 7.59e-1*T9)**(-2))*
| exp(-8.472*(T9*(1/(1. + 7.59e-1*T9)))**
| (-1.d0/3) )*T9**(-3.d0/2)*
| (T9*(1/(1. + 7.59e-1*T9)))**(-.5d0/3))*
| exp(-2.0129e2*(1/T9))
c(35)= 'LI7(P,AG)HE4'
dfdt(j,35)= 3.2269e11*exp(-3.0442e1*(1/T9))*T9**(-7.d0/2) +
| 6.9409e6*exp(-4.478*(1/T9))*T9**(-7.d0/2) -
| 1.59e10*exp(-3.0442e1*(1/T9))*T9**(-5.d0/2) -
| 2.325e6*exp(-4.478*(1/T9))*T9**(-5.d0/2) +
| 3.0951e9*exp(-8.472*T9**(-1.d0/3))*T9**(-2) -
| 7.3067e8*exp(-8.472*T9**(-1.d0/3))*T9**(-5.d0/3) +
| 7.245e8*exp
| (-8.472*(T9*(1/(1. + 7.59e-1*T9)))**(-1.d0/3))*
| T9**(-5.d0/2)*(T9*(1/(1. + 7.59e-1*T9)))**(2.5d0/3) -
| 3.0287e-7*(7.8176e15 + 1.*T9)*
| exp(-8.472*(T9*(1/(1. + 7.59e-1*T9)))**(-1.d0/3))*
| T9**(-3.d0/2)*(T9*(1/(1. + 7.59e-1*T9)))**(-5.e-1)*
| (1.3175 + 1.*T9)**(-2) -
| 8.9373e-8*(7.8176e15 + 1.*T9)*
| exp(-8.472*(T9*(1/(1. + 7.59e-1*T9)))**(-1.d0/3))*
| T9**(-3.d0/2)*(T9*(1/(1. + 7.59e-1*T9)))**(-.5d0/3)*
| (1.3175 + 1.*T9)**(-2) +
| 8.008e5*(3.1818e-3 + 3.2442e-1*T9**(1.d0/3) +
| 1.6753e-1*T9**(2.d0/3) + 9.1377e-1*T9 + 1.*T9**(4.d0/3))*
| exp(-8.472*T9**(-1.d0/3) - 3.4765e-1*T9**2)*T9**(-4.d0/3)
| - 1.04e5*(1. + 4.9e-2*T9**(1.d0/3) + 2.498*T9**(2.d0/3) +
| 8.6e-1*T9 + 3.518*T9**(4.d0/3) + 3.08*T9**(5.d0/3))*
| exp(-8.472*T9**(-1.d0/3) - 3.4765e-1*T9**2)*T9**(-5.d0/3)
| + 3.3408e5*(3.2468e-1 + 1.5909e-2*T9**(1.d0/3) +
| 8.1104e-1*T9**(2.d0/3) + 2.7922e-1*T9 +
| 1.1422*T9**(4.d0/3) + 1.*T9**(5.d0/3))*
| (4.0615 - 1.*T9**(7.d0/3))*
| exp(-8.472*T9**(-1.d0/3) - 3.4765e-1*T9**2)*T9**(-2)
drdt(j,35)= 9.4405e2*(1.06e10*exp(-3.0442e1*(1/T9))*
| T9**(-3.d0/2) +
| 1.55e6*exp(-4.478*(1/T9))*T9**(-3.d0/2) +
| 1.096e9*exp(-8.472*T9**(-1.d0/3))*T9**(-2.d0/3) -
| 4.83e8*exp
| (-8.472*(T9*(1/(1. + 7.59e-1*T9)))**(-1.d0/3))*
| T9**(-3.d0/2)*(T9*(1/(1. + 7.59e-1*T9)))**(2.5d0/3)
| + 1.56e5*(1. + 4.9e-2*T9**(1.d0/3) +
| 2.498*T9**(2.d0/3) + 8.6e-1*T9 + 3.518*T9**(4.d0/3) +
| 3.08*T9**(5.d0/3))*
| exp(-8.472*T9**(-1.d0/3) - 3.4765e-1*T9**2)*
| T9**(-2.d0/3))*exp(-2.0129e2*(1/T9))*T9**(-2) +
| 4.69*(3.2269e11*exp(-3.0442e1*(1/T9))*T9**(-7.d0/2) +
| 6.9409e6*exp(-4.478*(1/T9))*T9**(-7.d0/2) -
| 1.59e10*exp(-3.0442e1*(1/T9))*T9**(-5.d0/2) -
| 2.325e6*exp(-4.478*(1/T9))*T9**(-5.d0/2) +
| 3.0951e9*exp(-8.472*T9**(-1.d0/3))*T9**(-2) -
| 7.3067e8*exp(-8.472*T9**(-1.d0/3))*T9**(-5.d0/3) +
| 1.56e5*(8.6e-1 + 1.6333e-2*T9**(-2.d0/3) +
| 1.6653*T9**(-1.d0/3) + 4.6907*T9**(1.d0/3) +
| 5.1333*T9**(2.d0/3))*
| exp(-8.472*T9**(-1.d0/3) - 3.4765e-1*T9**2)*
| T9**(-2.d0/3) + 7.245e8*
| exp(-8.472*(T9*(1/(1. + 7.59e-1*T9)))**
| (-1.d0/3) )*T9**(-5.d0/2)*
| (T9*(1/(1. + 7.59e-1*T9)))**(2.5d0/3) -
| 1.364e9*(1/(1. + 7.59e-1*T9) -
| 7.59e-1*T9*(1. + 7.59e-1*T9)**(-2))*
| exp(-8.472*(T9*(1/(1. + 7.59e-1*T9)))**
| (-1.d0/3) )*T9**(-3.d0/2)*
| (T9*(1/(1. + 7.59e-1*T9)))**(-5.e-1) -
| 4.025e8*(1/(1. + 7.59e-1*T9) -
| 7.59e-1*T9*(1. + 7.59e-1*T9)**(-2))*
| exp(-8.472*(T9*(1/(1. + 7.59e-1*T9)))**
| (-1.d0/3) )*T9**(-3.d0/2)*
| (T9*(1/(1. + 7.59e-1*T9)))**(-.5d0/3) -
| 1.04e5*(1. + 4.9e-2*T9**(1.d0/3) + 2.498*T9**(2.d0/3) +
| 8.6e-1*T9 + 3.518*T9**(4.d0/3) + 3.08*T9**(5.d0/3))*
| exp(-8.472*T9**(-1.d0/3) - 3.4765e-1*T9**2)*
| T9**(-5.d0/3) + 1.56e5*
| (2.824*T9**(-4.d0/3) - 6.9531e-1*T9)*
| (1. + 4.9e-2*T9**(1.d0/3) + 2.498*T9**(2.d0/3) +
| 8.6e-1*T9 + 3.518*T9**(4.d0/3) + 3.08*T9**(5.d0/3))*
| exp(-8.472*T9**(-1.d0/3) - 3.4765e-1*T9**2)*
| T9**(-2.d0/3))*exp(-2.0129e2*(1/T9))
c(36)= 'LI7(D,N)2HE4'
dfdt(j,36)= 1.9467e11*(5.1295 - 1.*T9**(1.d0/3))*
| exp(-1.0259e1*T9**(-1.d0/3))*T9**(-2)
drdt(j,36)= -6.295e2*exp(-1.7548e2*(1/T9) - 1.0259e1*
| T9**(-1.d0/3))*
| T9**(-19.d0/6) + 2.9054e2*
| (1.7548e2*T9**(-2) + 3.4197*T9**(-4.d0/3))*
| exp(-1.7548e2*(1/T9) - 1.0259e1*T9**(-1.d0/3))*
| T9**(-13.d0/6)
c(37)= 'LI7(T,2N)2HE4'
dfdt(j,37)= 5.8733e11*(5.6665 - 1.*T9**(1.d0/3))*
| exp(-1.1333e1*T9**(-1.d0/3))*T9**(-2)
drdt(j,37)= -3.941e-7*exp(-1.0286e2*(1/T9) -
| 1.1333e1*T9**(-1.d0/3))*
| T9**(-14.d0/3) + 1.0748e-7*
| (1.0286e2*T9**(-2) + 3.7777*T9**(-4.d0/3))*
| exp(-1.0286e2*(1/T9) - 1.1333e1*T9**(-1.d0/3))*
| T9**(-11.d0/3)
c(38)= 'LI7(HE3,NP)2HE4'
dfdt(j,38)= 7.4e12*(8.9945 - 1.*T9**(1.d0/3))*
| exp(-1.7989e1*T9**(-1.d0/3))*T9**(-2)
drdt(j,38)= -2.4786e-6*exp(-1.1173e2*(1/T9) -
| 1.7989e1*T9**(-1.d0/3))*
| T9**(-14.d0/3) + 6.7599e-7*
| (1.1173e2*T9**(-2) + 5.9963*T9**(-4.d0/3))*
| exp(-1.1173e2*(1/T9) - 1.7989e1*T9**(-1.d0/3))*
| T9**(-11.d0/3)
c(39)= 'LI7(A,G)B11'
dfdt(j,39)= 9.9134e2*exp(-2.977*(1/T9))*T9**(-7.d0/2) +
| 2.5531e5*exp(-6.227*(1/T9))*T9**(-3) -
| 4.995e2*exp(-2.977*(1/T9))*T9**(-5.d0/2) -
| 4.1e4*exp(-6.227*(1/T9))*T9**(-2) +
| 2.0235e7*(1.2865e-2 + 9.0643e-1*T9**(1.d0/3) +
| 2.0702e-1*T9**(2.d0/3) + 2.0678*T9 + 1.*T9**(4.d0/3))*
| exp(-1.9161e1*T9**(-1.d0/3) - 5.6825e-2*T9**2)*
| T9**(-4.d0/3) - 2.3667e7*
| (1. + 2.2e-2*T9**(1.d0/3) + 7.75e-1*T9**(2.d0/3) +
| 1.18e-1*T9 + 8.84e-1*T9**(4.d0/3) + 3.42e-1*T9**(5.d0/3))
| *exp(-1.9161e1*T9**(-1.d0/3) - 5.6825e-2*T9**2)*
| T9**(-5.d0/3) + 1.3798e6*
| (2.924 + 6.4327e-2*T9**(1.d0/3) + 2.2661*T9**(2.d0/3) +
| 3.4503e-1*T9 + 2.5848*T9**(4.d0/3) + 1.*T9**(5.d0/3))*
| (5.6199e1 - 1.*T9**(7.d0/3))*
| exp(-1.9161e1*T9**(-1.d0/3) - 5.6825e-2*T9**2)*T9**(-2)
drdt(j,39)= 4.0416e12*(3.33e2*exp(-2.977*(1/T9))*
| T9**(-3.d0/2) +
| 4.1e4*exp(-6.227*(1/T9))*(1/T9) +
| 3.55e7*(1. + 2.2e-2*T9**(1.d0/3) +
| 7.75e-1*T9**(2.d0/3) + 1.18e-1*T9 +
| 8.84e-1*T9**(4.d0/3) + 3.42e-1*T9**(5.d0/3))*
| exp(-1.9161e1*T9**(-1.d0/3) - 5.6825e-2*T9**2)*
| T9**(-2.d0/3))*exp(-1.0054e2*(1/T9))*(1/Sqrt(T9)) +
| 6.03e10*(3.33e2*exp(-2.977*(1/T9))*T9**(-3.d0/2) +
| 4.1e4*exp(-6.227*(1/T9))*(1/T9) +
| 3.55e7*(1. + 2.2e-2*T9**(1.d0/3) +
| 7.75e-1*T9**(2.d0/3) + 1.18e-1*T9 +
| 8.84e-1*T9**(4.d0/3) + 3.42e-1*T9**(5.d0/3))*
| exp(-1.9161e1*T9**(-1.d0/3) - 5.6825e-2*T9**2)*
| T9**(-2.d0/3))*exp(-1.0054e2*(1/T9))*Sqrt(T9) +
| 4.02e10*(9.9134e2*exp(-2.977*(1/T9))*T9**(-7.d0/2) +
| 2.5531e5*exp(-6.227*(1/T9))*T9**(-3) -
| 4.995e2*exp(-2.977*(1/T9))*T9**(-5.d0/2) -
| 4.1e4*exp(-6.227*(1/T9))*T9**(-2) +
| 3.55e7*(1.18e-1 + 7.3333e-3*T9**(-2.d0/3) +
| 5.1667e-1*T9**(-1.d0/3) + 1.1787*T9**(1.d0/3) +
| 5.7e-1*T9**(2.d0/3))*
| exp(-1.9161e1*T9**(-1.d0/3) - 5.6825e-2*T9**2)*
| T9**(-2.d0/3) - 2.3667e7*
| (1. + 2.2e-2*T9**(1.d0/3) + 7.75e-1*T9**(2.d0/3) +
| 1.18e-1*T9 + 8.84e-1*T9**(4.d0/3) +
| 3.42e-1*T9**(5.d0/3))*
| exp(-1.9161e1*T9**(-1.d0/3) - 5.6825e-2*T9**2)*
| T9**(-5.d0/3) + 3.55e7*
| (6.387*T9**(-4.d0/3) - 1.1365e-1*T9)*
| (1. + 2.2e-2*T9**(1.d0/3) + 7.75e-1*T9**(2.d0/3) +
| 1.18e-1*T9 + 8.84e-1*T9**(4.d0/3) +
| 3.42e-1*T9**(5.d0/3))*
| exp(-1.9161e1*T9**(-1.d0/3) - 5.6825e-2*T9**2)*
| T9**(-2.d0/3))*exp(-1.0054e2*(1/T9))*T9**(3.d0/2)
c(40)= 'LI7(A,N)B10'
dfdt(j,40)= 1.2435e10*exp(-3.2382e1*(1/T9))*T9**(-2)
drdt(j,40)= 0
c(41)= 'BE7(E-,NU+G)LI7'
dfdt(j,41)= 5.427e-13*(-1.6767e-3*exp(2.515e-3*(1/T9)) -
| 1.*T9*exp(2.515e-3*(1/T9)) - 1.2346e2*T9**2 +
| 2.2099e1*T9**(7.d0/3) + 1.5885e2*T9**(8.d0/3))
| *T9**(-7.d0/2)
c(42)= 'BE7(P,G)B8'
dfdt(j,42)= 1.8484e4*exp(-7.306*(1/T9))*T9**(-7.d0/2) -
| 3.795e3*exp(-7.306*(1/T9))*T9**(-5.d0/2) +
| 1.0638e6*exp(-1.0262e1*T9**(-1.d0/3))*T9**(-2) -
| 2.0733e5*exp(-1.0262e1*T9**(-1.d0/3))*T9**(-5.d0/3)
drdt(j,42)= 2.0735e10*(2.53e3*exp(-7.306*(1/T9))*
| T9**(-3.d0/2) +
| 3.11e5*exp(-1.0262e1*T9**(-1.d0/3))*T9**(-2.d0/3))*
| exp(-1.595*(1/T9))*(1/Sqrt(T9)) +
| 1.95e10*(2.53e3*exp(-7.306*(1/T9))*T9**(-3.d0/2) +
| 3.11e5*exp(-1.0262e1*T9**(-1.d0/3))*T9**(-2.d0/3))*
| exp(-1.595*(1/T9))*Sqrt(T9) +
| 1.3e10*(1.8484e4*exp(-7.306*(1/T9))*T9**(-7.d0/2) -
| 3.795e3*exp(-7.306*(1/T9))*T9**(-5.d0/2) +
| 1.0638e6*exp(-1.0262e1*T9**(-1.d0/3))*T9**(-2) -
| 2.0733e5*exp(-1.0262e1*T9**(-1.d0/3))*T9**(-5.d0/3))*
| exp(-1.595*(1/T9))*T9**(3.d0/2)
c(43)= 'BE7(D,P)2HE4'
dfdt(j,43)= 7.1333e11*(6.214 - 1.*T9**(1.d0/3))*
| exp(-1.2428e1*T9**(-1.d0/3))*T9**(-2)
drdt(j,43)= -2.3114e3*exp(-1.9456e2*(1/T9) -
| 1.2428e1*T9**(-1.d0/3))*
| T9**(-19.d0/6) + 1.0668e3*
| (1.9456e2*T9**(-2) + 4.1427*T9**(-4.d0/3))*
| exp(-1.9456e2*(1/T9) - 1.2428e1*T9**(-1.d0/3))*
| T9**(-13.d0/6)
c(44)= 'BE7(T,NP)2HE4'
dfdt(j,44)= 1.94e12*(6.8645 - 1.*T9**(1.d0/3))*
| exp(-1.3729e1*T9**(-1.d0/3))*T9**(-2)
drdt(j,44)= -6.498e-7*exp(-1.2194e2*(1/T9) -
| 1.3729e1*T9**(-1.d0/3))*
| T9**(-14.d0/3) + 1.7722e-7*
| (1.2194e2*T9**(-2) + 4.5763*T9**(-4.d0/3))*
| exp(-1.2194e2*(1/T9) - 1.3729e1*T9**(-1.d0/3))*
| T9**(-11.d0/3)
c(45)= 'BE7(HE3,2P)2HE4'
dfdt(j,45)= 4.0733e13*(1.0896e1 - 1.*T9**(1.d0/3))*
| exp(-2.1793e1*T9**(-1.d0/3))*T9**(-2)
drdt(j,45)= -2.7332e-5*exp(-1.3081e2*(1/T9) -
| 2.1793e1*T9**(-1.d0/3))*
| T9**(-14.d0/3) + 7.4542e-6*
| (1.3081e2*T9**(-2) + 7.2643*T9**(-4.d0/3))*
| exp(-1.3081e2*(1/T9) - 2.1793e1*T9**(-1.d0/3))*
| T9**(-11.d0/3)
|
c(46)= 'BE7(A,G)C11'
dfdt(j,46)= 8.1375e4*exp(-6.51*(1/T9))*T9**(-7.d0/2) +
| 1.295e6*exp(-1.0039e1*(1/T9))*T9**(-13.d0/4) -
| 1.875e4*exp(-6.51*(1/T9))*T9**(-5.d0/2) -
| 1.6125e5*exp(-1.0039e1*(1/T9))*T9**(-9.d0/4) +
| 1.3379e7*(3.7895e-2 + 2.0547*T9**(1.d0/3) +
| 3.8526e-1*T9**(2.d0/3) + 2.4926*T9 + 1.*T9**(4.d0/3))*
| exp(-2.3212e1*T9**(-1.d0/3) - 4.3969e-2*T9**2)*
| T9**(-4.d0/3) - 5.6333e7*
| (1. + 1.8e-2*T9**(1.d0/3) + 4.88e-1*T9**(2.d0/3) +
| 6.1e-2*T9 + 2.96e-1*T9**(4.d0/3) + 9.5e-2*T9**(5.d0/3))*
| exp(-2.3212e1*T9**(-1.d0/3) - 4.3969e-2*T9**2)*
| T9**(-5.d0/3) + 7.0592e5*
| (1.0526e1 + 1.8947e-1*T9**(1.d0/3) + 5.1368*T9**(2.d0/3) +
| 6.4211e-1*T9 + 3.1158*T9**(4.d0/3) + 1.*T9**(5.d0/3))*
| (8.7986e1 - 1.*T9**(7.d0/3))*
| exp(-2.3212e1*T9**(-1.d0/3) - 4.3969e-2*T9**2)*T9**(-2)
drdt(j,46)= 3.5191e12*(1.25e4*exp(-6.51*(1/T9))
| *T9**(-3.d0/2) +
| 1.29e5*exp(-1.0039e1*(1/T9))*T9**(-5.d0/4) +
| 8.45e7*(1. + 1.8e-2*T9**(1.d0/3) +
| 4.88e-1*T9**(2.d0/3) + 6.1e-2*T9 +
| 2.96e-1*T9**(4.d0/3) + 9.5e-2*T9**(5.d0/3))*
| exp(-2.3212e1*T9**(-1.d0/3) - 4.3969e-2*T9**2)*
| T9**(-2.d0/3))*exp(-8.7539e1*(1/T9))*(1/Sqrt(T9)) +
| 6.03e10*(1.25e4*exp(-6.51*(1/T9))*T9**(-3.d0/2) +
| 1.29e5*exp(-1.0039e1*(1/T9))*T9**(-5.d0/4) +
| 8.45e7*(1. + 1.8e-2*T9**(1.d0/3) +
| 4.88e-1*T9**(2.d0/3) + 6.1e-2*T9 +
| 2.96e-1*T9**(4.d0/3) + 9.5e-2*T9**(5.d0/3))*
| exp(-2.3212e1*T9**(-1.d0/3) - 4.3969e-2*T9**2)*
| T9**(-2.d0/3))*exp(-8.7539e1*(1/T9))*Sqrt(T9) +
| 4.02e10*(8.1375e4*exp(-6.51*(1/T9))*T9**(-7.d0/2) +
| 1.295e6*exp(-1.0039e1*(1/T9))*T9**(-13.d0/4) -
| 1.875e4*exp(-6.51*(1/T9))*T9**(-5.d0/2) -
| 1.6125e5*exp(-1.0039e1*(1/T9))*T9**(-9.d0/4) +
| 8.45e7*(6.1e-2 + 6.e-3*T9**(-2.d0/3) +
| 3.2533e-1*T9**(-1.d0/3) + 3.9467e-1*T9**(1.d0/3) +
| 1.5833e-1*T9**(2.d0/3))*
| exp(-2.3212e1*T9**(-1.d0/3) - 4.3969e-2*T9**2)*
| T9**(-2.d0/3) - 5.6333e7*
| (1. + 1.8e-2*T9**(1.d0/3) + 4.88e-1*T9**(2.d0/3) +
| 6.1e-2*T9 + 2.96e-1*T9**(4.d0/3) + 9.5e-2*T9**(5.d0/3)
| )*exp(-2.3212e1*T9**(-1.d0/3) - 4.3969e-2*T9**2)*
| T9**(-5.d0/3) + 8.45e7*
| (7.7373*T9**(-4.d0/3) - 8.7938e-2*T9)*
| (1. + 1.8e-2*T9**(1.d0/3) + 4.88e-1*T9**(2.d0/3) +
| 6.1e-2*T9 + 2.96e-1*T9**(4.d0/3) + 9.5e-2*T9**(5.d0/3)
| )*exp(-2.3212e1*T9**(-1.d0/3) - 4.3969e-2*T9**2)*
| T9**(-2.d0/3))*exp(-8.7539e1*(1/T9))*T9**(3.d0/2)
c(47)= 'BE9(P,G)B10'
dfdt(j,47)= 2.8886e7*exp(-1.062e1*(1/T9))*T9**(-7.d0/2) +
| 3.321e5*exp(-3.445*(1/T9))*T9**(-7.d0/2) -
| 4.08e6*exp(-1.062e1*(1/T9))*T9**(-5.d0/2) -
| 1.446e5*exp(-3.445*(1/T9))*T9**(-5.d0/2) +
| 3.4137e7*(5.1948e-3 + 3.9481e-1*T9**(1.d0/3) +
| 1.6675e-1*T9**(2.d0/3) + 1.1169*T9 + 1.*T9**(4.d0/3))*
| exp(-1.0359e1*T9**(-1.d0/3) - 1.3972*T9**2)*T9**(-4.d0/3)
| - 8.8667e6*(1. + 4.e-2*T9**(1.d0/3) + 1.52*T9**(2.d0/3) +
| 4.28e-1*T9 + 2.15*T9**(4.d0/3) + 1.54*T9**(5.d0/3))*
| exp(-1.0359e1*T9**(-1.d0/3) - 1.3972*T9**2)*T9**(-5.d0/3)
| + 5.7235e7*(6.4935e-1 + 2.5974e-2*T9**(1.d0/3) +
| 9.8701e-1*T9**(2.d0/3) + 2.7792e-1*T9 +
| 1.3961*T9**(4.d0/3) + 1.*T9**(5.d0/3))*
| (1.2357 - 1.*T9**(7.d0/3))*
| exp(-1.0359e1*T9**(-1.d0/3) - 1.3972*T9**2)*T9**(-2)
drdt(j,47)= 7.4363e11*(2.72e6*exp(-1.062e1*(1/T9))*
| T9**(-3.d0/2) +
| 9.64e4*exp(-3.445*(1/T9))*T9**(-3.d0/2) +
| 1.33e7*(1. + 4.e-2*T9**(1.d0/3) + 1.52*T9**(2.d0/3) +
| 4.28e-1*T9 + 2.15*T9**(4.d0/3) + 1.54*T9**(5.d0/3))*
| exp(-1.0359e1*T9**(-1.d0/3) - 1.3972*T9**2)*
| T9**(-2.d0/3))*exp(-7.6427e1*(1/T9))*(1/Sqrt(T9)) +
| 1.4595e10*(2.72e6*exp(-1.062e1*(1/T9))*T9**(-3.d0/2) +
| 9.64e4*exp(-3.445*(1/T9))*T9**(-3.d0/2) +
| 1.33e7*(1. + 4.e-2*T9**(1.d0/3) + 1.52*T9**(2.d0/3) +
| 4.28e-1*T9 + 2.15*T9**(4.d0/3) + 1.54*T9**(5.d0/3))*
| exp(-1.0359e1*T9**(-1.d0/3) - 1.3972*T9**2)*
| T9**(-2.d0/3))*exp(-7.6427e1*(1/T9))*Sqrt(T9) +
| 9.73e9*(2.8886e7*exp(-1.062e1*(1/T9))*T9**(-7.d0/2) +
| 3.321e5*exp(-3.445*(1/T9))*T9**(-7.d0/2) -
| 4.08e6*exp(-1.062e1*(1/T9))*T9**(-5.d0/2) -
| 1.446e5*exp(-3.445*(1/T9))*T9**(-5.d0/2) +
| 1.33e7*(4.28e-1 + 1.3333e-2*T9**(-2.d0/3) +
| 1.0133*T9**(-1.d0/3) + 2.8667*T9**(1.d0/3) +
| 2.5667*T9**(2.d0/3))*
| exp(-1.0359e1*T9**(-1.d0/3) - 1.3972*T9**2)*
| T9**(-2.d0/3) - 8.8667e6*
| (1. + 4.e-2*T9**(1.d0/3) + 1.52*T9**(2.d0/3) +
| 4.28e-1*T9 + 2.15*T9**(4.d0/3) + 1.54*T9**(5.d0/3))*
| exp(-1.0359e1*T9**(-1.d0/3) - 1.3972*T9**2)*
| T9**(-5.d0/3) + 1.33e7*
| (3.453*T9**(-4.d0/3) - 2.7944*T9)*
| (1. + 4.e-2*T9**(1.d0/3) + 1.52*T9**(2.d0/3) +
| 4.28e-1*T9 + 2.15*T9**(4.d0/3) + 1.54*T9**(5.d0/3))*
| exp(-1.0359e1*T9**(-1.d0/3) - 1.3972*T9**2)*
| T9**(-2.d0/3))*exp(-7.6427e1*(1/T9))*T9**(3.d0/2)
c(48)= 'BE9(P,N)B9'
dfdt(j,48)= 5.58e7*(9.85e-1 + 2.1e-2*(1/Sqrt(T9)))*
| exp(-2.1473e1*(1/T9)) +
| 2.7259e10*exp(-2.6725e1*(1/T9))*T9**(-7.d0/2) -
| 1.53e9*exp(-2.6725e1*(1/T9))*T9**(-5.d0/2) +
| 1.1982e9*(1. + 4.2e-2*Sqrt(T9) + 9.85e-1*T9)*
| exp(-2.1473e1*(1/T9))*T9**(-2)
drdt(j,48)= 5.5688e7*(9.85e-1 + 2.1e-2*(1/Sqrt(T9))) +
| 5.3463e9*exp(-5.252*(1/T9))*T9**(-7.d0/2) -
| 1.5269e9*exp(-5.252*(1/T9))*T9**(-5.d0/2)
c(49)= 'BE9(P,D)2HE4'
dfdt(j,49)= 1.7636e9*exp(-3.046*(1/T9))*T9**(-3) +
| 4.93e9*exp(-5.8*(1/T9))*T9**(-11.d0/4) -
| 5.79e8*exp(-3.046*(1/T9))*T9**(-2) -
| 6.375e8*exp(-5.8*(1/T9))*T9**(-7.d0/4) +
| 8.0883e11*(3.4783e-3 + 1.8957e-1*T9**(1.d0/3) +
| 8.0087e-2*T9**(2.d0/3) + 1.1165*T9 + 1.*T9**(4.d0/3))*
| exp(-1.0359e1*T9**(-1.d0/3) - 3.6982*T9**2)*T9**(-4.d0/3)
| - 1.4067e11*(1. + 4.e-2*T9**(1.d0/3) + 1.09*T9**(2.d0/3) +
| 3.07e-1*T9 + 3.21*T9**(4.d0/3) + 2.3*T9**(5.d0/3))*
| exp(-1.0359e1*T9**(-1.d0/3) - 3.6982*T9**2)*T9**(-5.d0/3)
| + 3.5895e12*(4.3478e-1 + 1.7391e-2*T9**(1.d0/3) +
| 4.7391e-1*T9**(2.d0/3) + 1.3348e-1*T9 +
| 1.3957*T9**(4.d0/3) + 1.*T9**(5.d0/3))*
| (4.6685e-1 - 1.*T9**(7.d0/3))*
| exp(-1.0359e1*T9**(-1.d0/3) - 3.6982*T9**2)*T9**(-2)
drdt(j,49)= 6.0969e-10*(5.79e8*exp(-3.046*(1/T9))
| *(1/T9) +
| 8.5e8*exp(-5.8*(1/T9))*T9**(-3.d0/4) +
| 2.11e11*(1. + 4.e-2*T9**(1.d0/3) + 1.09*T9**(2.d0/3) +
| 3.07e-1*T9 + 3.21*T9**(4.d0/3) + 2.3*T9**(5.d0/3))*
| exp(-1.0359e1*T9**(-1.d0/3) - 3.6982*T9**2)*
| T9**(-2.d0/3))*exp(-7.555*(1/T9))*T9**(-7.d0/2) -
| 1.2105e-10*(5.79e8*exp(-3.046*(1/T9))*(1/T9) +
| 8.5e8*exp(-5.8*(1/T9))*T9**(-3.d0/4) +
| 2.11e11*(1. + 4.e-2*T9**(1.d0/3) + 1.09*T9**(2.d0/3) +
| 3.07e-1*T9 + 3.21*T9**(4.d0/3) + 2.3*T9**(5.d0/3))*
| exp(-1.0359e1*T9**(-1.d0/3) - 3.6982*T9**2)*
| T9**(-2.d0/3))*exp(-7.555*(1/T9))*T9**(-5.d0/2) +
| 8.07e-11*(1.7636e9*exp(-3.046*(1/T9))*T9**(-3) +
| 4.93e9*exp(-5.8*(1/T9))*T9**(-11.d0/4) -
| 5.79e8*exp(-3.046*(1/T9))*T9**(-2) -
| 6.375e8*exp(-5.8*(1/T9))*T9**(-7.d0/4) +
| 2.11e11*(3.07e-1 + 1.3333e-2*T9**(-2.d0/3) +
| 7.2667e-1*T9**(-1.d0/3) + 4.28*T9**(1.d0/3) +
| 3.8333*T9**(2.d0/3))*
| exp(-1.0359e1*T9**(-1.d0/3) - 3.6982*T9**2)*
| T9**(-2.d0/3) - 1.4067e11*
| (1. + 4.e-2*T9**(1.d0/3) + 1.09*T9**(2.d0/3) +
| 3.07e-1*T9 + 3.21*T9**(4.d0/3) + 2.3*T9**(5.d0/3))*
| exp(-1.0359e1*T9**(-1.d0/3) - 3.6982*T9**2)*
| T9**(-5.d0/3) + 2.11e11*
| (3.453*T9**(-4.d0/3) - 7.3964*T9)*
| (1. + 4.e-2*T9**(1.d0/3) + 1.09*T9**(2.d0/3) +
| 3.07e-1*T9 + 3.21*T9**(4.d0/3) + 2.3*T9**(5.d0/3))*
| exp(-1.0359e1*T9**(-1.d0/3) - 3.6982*T9**2)*
| T9**(-2.d0/3))*exp(-7.555*(1/T9))*T9**(-3.d0/2)
c(50)= 'BE9(P,A)LI6'
dfdt(j,50)= 1.3737e9*exp(-3.046*(1/T9))*T9**(-3) +
| 3.4572e9*exp(-5.16*(1/T9))*T9**(-11.d0/4) -
| 4.51e8*exp(-3.046*(1/T9))*T9**(-2) -
| 5.025e8*exp(-5.16*(1/T9))*T9**(-7.d0/4) +
| 8.0883e11*(3.4783e-3 + 1.8957e-1*T9**(1.d0/3) +
| 8.0087e-2*T9**(2.d0/3) + 1.1165*T9 + 1.*T9**(4.d0/3))*
| exp(-1.0359e1*T9**(-1.d0/3) - 3.6982*T9**2)*T9**(-4.d0/3)
| - 1.4067e11*(1. + 4.e-2*T9**(1.d0/3) + 1.09*T9**(2.d0/3) +
| 3.07e-1*T9 + 3.21*T9**(4.d0/3) + 2.3*T9**(5.d0/3))*
| exp(-1.0359e1*T9**(-1.d0/3) - 3.6982*T9**2)*T9**(-5.d0/3)
| + 3.5895e12*(4.3478e-1 + 1.7391e-2*T9**(1.d0/3) +
| 4.7391e-1*T9**(2.d0/3) + 1.3348e-1*T9 +
| 1.3957*T9**(4.d0/3) + 1.*T9**(5.d0/3))*
| (4.6685e-1 - 1.*T9**(7.d0/3))*
| exp(-1.0359e1*T9**(-1.d0/3) - 3.6982*T9**2)*T9**(-2)
drdt(j,50)= 1.5249e1*(4.51e8*exp(-3.046*(1/T9))*(1/T9) +
| 6.7e8*exp(-5.16*(1/T9))*T9**(-3.d0/4) +
| 2.11e11*(1. + 4.e-2*T9**(1.d0/3) + 1.09*T9**(2.d0/3) +
| 3.07e-1*T9 + 3.21*T9**(4.d0/3) + 2.3*T9**(5.d0/3))*
| exp(-1.0359e1*T9**(-1.d0/3) - 3.6982*T9**2)*
| T9**(-2.d0/3))*exp(-2.4674e1*(1/T9))*T9**(-2) +
| 6.18e-1*(1.3737e9*exp(-3.046*(1/T9))*T9**(-3) +
| 3.4572e9*exp(-5.16*(1/T9))*T9**(-11.d0/4) -
| 4.51e8*exp(-3.046*(1/T9))*T9**(-2) -
| 5.025e8*exp(-5.16*(1/T9))*T9**(-7.d0/4) +
| 2.11e11*(3.07e-1 + 1.3333e-2*T9**(-2.d0/3) +
| 7.2667e-1*T9**(-1.d0/3) + 4.28*T9**(1.d0/3) +
| 3.8333*T9**(2.d0/3))*
| exp(-1.0359e1*T9**(-1.d0/3) - 3.6982*T9**2)*
| T9**(-2.d0/3) - 1.4067e11*
| (1. + 4.e-2*T9**(1.d0/3) + 1.09*T9**(2.d0/3) +
| 3.07e-1*T9 + 3.21*T9**(4.d0/3) + 2.3*T9**(5.d0/3))*
| exp(-1.0359e1*T9**(-1.d0/3) - 3.6982*T9**2)*
| T9**(-5.d0/3) + 2.11e11*
| (3.453*T9**(-4.d0/3) - 7.3964*T9)*
| (1. + 4.e-2*T9**(1.d0/3) + 1.09*T9**(2.d0/3) +
| 3.07e-1*T9 + 3.21*T9**(4.d0/3) + 2.3*T9**(5.d0/3))*
| exp(-1.0359e1*T9**(-1.d0/3) - 3.6982*T9**2)*
| T9**(-2.d0/3))*exp(-2.4674e1*(1/T9))
c(51)= 'BE9(A,N)C12'
dfdt(j,51)= 5.2655e5*exp(-4.179*(1/T9))*T9**(-7.d0/2) +
| 4.1632e-1*exp(-1.834*(1/T9))*T9**(-7.d0/2) +
| 8.6906e-5*exp(-1.184*(1/T9))*T9**(-7.d0/2) -
| 1.89e5*exp(-4.179*(1/T9))*T9**(-5.d0/2) -
| 3.405e-1*exp(-1.834*(1/T9))*T9**(-5.d0/2) -
| 1.101e-4*exp(-1.184*(1/T9))*T9**(-5.d0/2) +
| 3.0557e9*exp(-1.2732e1*(1/T9))*T9**(-2) +
| 1.7826e15*(1.4687e-4 + 1.4808e-1*T9**(1.d0/3) +
| 2.7214e-2*T9**(2.d0/3) + 2.5749*T9 + 1.*T9**(4.d0/3))*
| exp(-2.387e1*T9**(-1.d0/3) - 4.1649e2*T9**2)*T9**(-4.d0/3)
| - 3.08e13*(1. + 1.7e-2*T9**(1.d0/3) + 8.57*T9**(2.d0/3) +
| 1.05*T9 + 7.451e1*T9**(4.d0/3) + 2.315e1*T9**(5.d0/3))*
| exp(-2.387e1*T9**(-1.d0/3) - 4.1649e2*T9**2)*T9**(-5.d0/3)
| + 8.909e17*(4.3197e-2 + 7.3434e-4*T9**(1.d0/3) +
| 3.7019e-1*T9**(2.d0/3) + 4.5356e-2*T9 +
| 3.2186*T9**(4.d0/3) + 1.*T9**(5.d0/3))*
| (9.552e-3 - 1.*T9**(7.d0/3))*
| exp(-2.387e1*T9**(-1.d0/3) - 4.1649e2*T9**2)*T9**(-2)
drdt(j,51)= 6.8145e2*(2.4e8*exp(-1.2732e1*(1/T9)) +
| 1.26e5*exp(-4.179*(1/T9))*T9**(-3.d0/2) +
| 2.27e-1*exp(-1.834*(1/T9))*T9**(-3.d0/2) +
| 7.34e-5*exp(-1.184*(1/T9))*T9**(-3.d0/2) +
| 4.62e13*(1. + 1.7e-2*T9**(1.d0/3) + 8.57*T9**(2.d0/3) +
| 1.05*T9 + 7.451e1*T9**(4.d0/3) + 2.315e1*T9**(5.d0/3))
| *exp(-2.387e1*T9**(-1.d0/3) - 4.1649e2*T9**2)*
| T9**(-2.d0/3))*exp(-6.616e1*(1/T9))*T9**(-2) +
| 1.03e1*(5.2655e5*exp(-4.179*(1/T9))*T9**(-7.d0/2) +
| 4.1632e-1*exp(-1.834*(1/T9))*T9**(-7.d0/2) +
| 8.6906e-5*exp(-1.184*(1/T9))*T9**(-7.d0/2) -
| 1.89e5*exp(-4.179*(1/T9))*T9**(-5.d0/2) -
| 3.405e-1*exp(-1.834*(1/T9))*T9**(-5.d0/2) -
| 1.101e-4*exp(-1.184*(1/T9))*T9**(-5.d0/2) +
| 3.0557e9*exp(-1.2732e1*(1/T9))*T9**(-2) +
| 4.62e13*(1.05 + 5.6667e-3*T9**(-2.d0/3) +
| 5.7133*T9**(-1.d0/3) + 9.9347e1*T9**(1.d0/3) +
| 3.8583e1*T9**(2.d0/3))*
| exp(-2.387e1*T9**(-1.d0/3) - 4.1649e2*T9**2)*
| T9**(-2.d0/3) - 3.08e13*
| (1. + 1.7e-2*T9**(1.d0/3) + 8.57*T9**(2.d0/3) +
| 1.05*T9 + 7.451e1*T9**(4.d0/3) + 2.315e1*T9**(5.d0/3))
| *exp(-2.387e1*T9**(-1.d0/3) - 4.1649e2*T9**2)*
| T9**(-5.d0/3) + 4.62e13*
| (7.9567*T9**(-4.d0/3) - 8.3299e2*T9)*
| (1. + 1.7e-2*T9**(1.d0/3) + 8.57*T9**(2.d0/3) +
| 1.05*T9 + 7.451e1*T9**(4.d0/3) + 2.315e1*T9**(5.d0/3))
| *exp(-2.387e1*T9**(-1.d0/3) - 4.1649e2*T9**2)*
| T9**(-2.d0/3))*exp(-6.616e1*(1/T9))
c(52)= 'B10(P,G)C11'
dfdt(j,52)= 1.8427e5*exp(-1.6164e1*(1/T9))*T9**(-7.d0/2) +
| 2.3239e6*exp(-1.2041e1*(1/T9))*T9**(-7.d0/2) -
| 1.71e4*exp(-1.6164e1*(1/T9))*T9**(-5.d0/2) -
| 2.895e5*exp(-1.2041e1*(1/T9))*T9**(-5.d0/2) +
| 1.3292e5*(4.0462e-2 + 9.8497e-1*T9**(1.d0/3) +
| 3.5723e-1*T9**(2.d0/3) + 1.2994*T9 + 1.*T9**(4.d0/3))*
| exp(-1.2062e1*T9**(-1.d0/3) - 5.1606e-2*T9**2)*
| T9**(-4.d0/3) - 3.0733e5*
| (1. + 3.5e-2*T9**(1.d0/3) + 4.26e-1*T9**(2.d0/3) +
| 1.03e-1*T9 + 2.81e-1*T9**(4.d0/3) + 1.73e-1*T9**(5.d0/3))
| *exp(-1.2062e1*T9**(-1.d0/3) - 5.1606e-2*T9**2)*
| T9**(-5.d0/3) + 8.2315e3*
| (5.7803 + 2.0231e-1*T9**(1.d0/3) + 2.4624*T9**(2.d0/3) +
| 5.9538e-1*T9 + 1.6243*T9**(4.d0/3) + 1.*T9**(5.d0/3))*
| (3.8955e1 - 1.*T9**(7.d0/3))*
| exp(-1.2062e1*T9**(-1.d0/3) - 5.1606e-2*T9**2)*T9**(-2)
drdt(j,52)= 3.0555e12*(1.14e4*exp(-1.6164e1*(1/T9))
| *T9**(-3.d0/2) +
| 1.93e5*exp(-1.2041e1*(1/T9))*T9**(-3.d0/2) +
| 4.61e5*(1. + 3.5e-2*T9**(1.d0/3) +
| 4.26e-1*T9**(2.d0/3) + 1.03e-1*T9 +
| 2.81e-1*T9**(4.d0/3) + 1.73e-1*T9**(5.d0/3))*
| exp(-1.2062e1*T9**(-1.d0/3) - 5.1606e-2*T9**2)*
| T9**(-2.d0/3))*exp(-1.0084e2*(1/T9))*(1/Sqrt(T9)) +
| 4.545e10*(1.14e4*exp(-1.6164e1*(1/T9))*T9**(-3.d0/2) +
| 1.93e5*exp(-1.2041e1*(1/T9))*T9**(-3.d0/2) +
| 4.61e5*(1. + 3.5e-2*T9**(1.d0/3) +
| 4.26e-1*T9**(2.d0/3) + 1.03e-1*T9 +
| 2.81e-1*T9**(4.d0/3) + 1.73e-1*T9**(5.d0/3))*
| exp(-1.2062e1*T9**(-1.d0/3) - 5.1606e-2*T9**2)*
| T9**(-2.d0/3))*exp(-1.0084e2*(1/T9))*Sqrt(T9) +
| 3.03e10*(1.8427e5*exp(-1.6164e1*(1/T9))*T9**(-7.d0/2) +
| 2.3239e6*exp(-1.2041e1*(1/T9))*T9**(-7.d0/2) -
| 1.71e4*exp(-1.6164e1*(1/T9))*T9**(-5.d0/2) -
| 2.895e5*exp(-1.2041e1*(1/T9))*T9**(-5.d0/2) +
| 4.61e5*(1.03e-1 + 1.1667e-2*T9**(-2.d0/3) +
| 2.84e-1*T9**(-1.d0/3) + 3.7467e-1*T9**(1.d0/3) +
| 2.8833e-1*T9**(2.d0/3))*
| exp(-1.2062e1*T9**(-1.d0/3) - 5.1606e-2*T9**2)*
| T9**(-2.d0/3) - 3.0733e5*
| (1. + 3.5e-2*T9**(1.d0/3) + 4.26e-1*T9**(2.d0/3) +
| 1.03e-1*T9 + 2.81e-1*T9**(4.d0/3) +
| 1.73e-1*T9**(5.d0/3))*
| exp(-1.2062e1*T9**(-1.d0/3) - 5.1606e-2*T9**2)*
| T9**(-5.d0/3) + 4.61e5*
| (4.0207*T9**(-4.d0/3) - 1.0321e-1*T9)*
| (1. + 3.5e-2*T9**(1.d0/3) + 4.26e-1*T9**(2.d0/3) +
| 1.03e-1*T9 + 2.81e-1*T9**(4.d0/3) +
| 1.73e-1*T9**(5.d0/3))*
| exp(-1.2062e1*T9**(-1.d0/3) - 5.1606e-2*T9**2)*
| T9**(-2.d0/3))*exp(-1.0084e2*(1/T9))*T9**(3.d0/2)
c(53)= 'B10(P,A)BE7'
dfdt(j,53)= 3.1753e10*exp(-1.226e1*(1/T9))*T9**(-3) -
| 2.59e9*exp(-1.226e1*(1/T9))*T9**(-2) +
| 3.864e10*(3.8043e-2 - 1.0826*T9**(1.d0/3) -
| 3.9457e-1*T9**(2.d0/3) + 1.3043*T9 + 1.*T9**(4.d0/3))*
| exp(-1.2062e1*T9**(-1.d0/3) - 5.1606e-2*T9**2)*
| T9**(-4.d0/3) + 1.5456e10*
| (-5.4348 - 1.9022e-1*T9**(1.d0/3) + 2.7065*T9**(2.d0/3) +
| 6.5761e-1*T9 - 1.6304*T9**(4.d0/3) - 1.*T9**(5.d0/3))*
| exp(-1.2062e1*T9**(-1.d0/3) - 5.1606e-2*T9**2)*
| T9**(-5.d0/3) + 2.3929e9*
| (5.4348 + 1.9022e-1*T9**(1.d0/3) - 2.7065*T9**(2.d0/3) -
| 6.5761e-1*T9 + 1.6304*T9**(4.d0/3) + 1.*T9**(5.d0/3))*
| (3.8955e1 - 1.*T9**(7.d0/3))*
| exp(-1.2062e1*T9**(-1.d0/3) - 5.1606e-2*T9**2)*T9**(-2)
drdt(j,53)= 1.0029e1*(2.59e9*exp(-1.226e1*(1/T9))*(1/T9) +
| 1.26e11*(1. + 3.5e-2*T9**(1.d0/3) -
| 4.98e-1*T9**(2.d0/3) - 1.21e-1*T9 +
| 3.e-1*T9**(4.d0/3) + 1.84e-1*T9**(5.d0/3))*
| exp(-1.2062e1*T9**(-1.d0/3) - 5.1606e-2*T9**2)*
| T9**(-2.d0/3))*exp(-1.3301e1*(1/T9))*T9**(-2) +
| 7.54e-1*(3.1753e10*exp(-1.226e1*(1/T9))*T9**(-3) -
| 2.59e9*exp(-1.226e1*(1/T9))*T9**(-2) +
| 1.26e11*(-1.21e-1 + 1.1667e-2*T9**(-2.d0/3) -
| 3.32e-1*T9**(-1.d0/3) + 4.e-1*T9**(1.d0/3) +
| 3.0667e-1*T9**(2.d0/3))*
| exp(-1.2062e1*T9**(-1.d0/3) - 5.1606e-2*T9**2)*
| T9**(-2.d0/3) - 8.4e10*
| (1. + 3.5e-2*T9**(1.d0/3) - 4.98e-1*T9**(2.d0/3) -
| 1.21e-1*T9 + 3.e-1*T9**(4.d0/3) + 1.84e-1*T9**(5.d0/3)
| )*exp(-1.2062e1*T9**(-1.d0/3) - 5.1606e-2*T9**2)*
| T9**(-5.d0/3) + 1.26e11*
| (4.0207*T9**(-4.d0/3) - 1.0321e-1*T9)*
| (1. + 3.5e-2*T9**(1.d0/3) - 4.98e-1*T9**(2.d0/3) -
| 1.21e-1*T9 + 3.e-1*T9**(4.d0/3) + 1.84e-1*T9**(5.d0/3)
| )*exp(-1.2062e1*T9**(-1.d0/3) - 5.1606e-2*T9**2)*
| T9**(-2.d0/3))*exp(-1.3301e1*(1/T9))
c(54)= 'B10(A,N)N13'
dfdt(j,54)= 2.6101e11*(4.2893e2 - 3.065e1*T9**(1.d0/3) -
| 1.*T9**(7.d0/3))*
| exp(-2.7989e1*T9**(-1.d0/3) - 1.0876e-2*T9**2)*T9**(-2)
drdt(j,54)= -7.472e13*exp
| (-1.2287e1*(1/T9) - 2.7989e1*T9**(-1.d0/3) -
| 1.0876e-2*T9**2)*T9**(-5.d0/3) +
| 1.1208e14*(1.2287e1*T9**(-2) + 9.3297*T9**(-4.d0/3) -
| 2.1751e-2*T9)*exp
| (-1.2287e1*(1/T9) - 2.7989e1*T9**(-1.d0/3) -
| 1.0876e-2*T9**2)*T9**(-2.d0/3)
c(55)= 'B11(P,G)C12'
dfdt(j,55)= 1.3673e4*exp(-1.733*(1/T9))*T9**(-7.d0/2) -
| 1.1835e4*exp(-1.733*(1/T9))*T9**(-5.d0/2) +
| 5.4373e5*exp(-5.617*(1/T9))*T9**(-11.d0/5) -
| 1.936e4*exp(-5.617*(1/T9))*T9**(-6.d0/5) +
| 4.6739e8*(1.1532e-3 + 1.9769e-1*T9**(1.d0/3) +
| 7.1466e-2*T9**(2.d0/3) + 1.3061*T9 + 1.*T9**(4.d0/3))*
| exp(-1.2095e1*T9**(-1.d0/3) - 1.7507e1*T9**2)*T9**(-4.d0/3)
| - 3.08e7*(1. + 3.5e-2*T9**(1.d0/3) + 3.*T9**(2.d0/3) +
| 7.23e-1*T9 + 9.91*T9**(4.d0/3) + 6.07*T9**(5.d0/3))*
| exp(-1.2095e1*T9**(-1.d0/3) - 1.7507e1*T9**2)*T9**(-5.d0/3)
| + 9.8189e9*(1.6474e-1 + 5.7661e-3*T9**(1.d0/3) +
| 4.9423e-1*T9**(2.d0/3) + 1.1911e-1*T9 +
| 1.6326*T9**(4.d0/3) + 1.*T9**(5.d0/3))*
| (1.1515e-1 - 1.*T9**(7.d0/3))*
| exp(-1.2095e1*T9**(-1.d0/3) - 1.7507e1*T9**2)*T9**(-2)
drdt(j,55)= 1.2981e13*(7.89e3*exp(-1.733*(1/T9))*
| T9**(-3.d0/2) +
| 9.68e4*exp(-5.617*(1/T9))*T9**(-1.d0/5) +
| 4.62e7*(1. + 3.5e-2*T9**(1.d0/3) + 3.*T9**(2.d0/3) +
| 7.23e-1*T9 + 9.91*T9**(4.d0/3) + 6.07*T9**(5.d0/3))*
| exp(-1.2095e1*T9**(-1.d0/3) - 1.7507e1*T9**2)*
| T9**(-2.d0/3))*exp(-1.8517e2*(1/T9))*(1/Sqrt(T9)) +
| 1.0515e11*(7.89e3*exp(-1.733*(1/T9))*T9**(-3.d0/2) +
| 9.68e4*exp(-5.617*(1/T9))*T9**(-1.d0/5) +
| 4.62e7*(1. + 3.5e-2*T9**(1.d0/3) + 3.*T9**(2.d0/3) +
| 7.23e-1*T9 + 9.91*T9**(4.d0/3) + 6.07*T9**(5.d0/3))*
| exp(-1.2095e1*T9**(-1.d0/3) - 1.7507e1*T9**2)*
| T9**(-2.d0/3))*exp(-1.8517e2*(1/T9))*Sqrt(T9) +
| 7.01e10*(1.3673e4*exp(-1.733*(1/T9))*T9**(-7.d0/2) -
| 1.1835e4*exp(-1.733*(1/T9))*T9**(-5.d0/2) +
| 5.4373e5*exp(-5.617*(1/T9))*T9**(-11.d0/5) -
| 1.936e4*exp(-5.617*(1/T9))*T9**(-6.d0/5) +
| 4.62e7*(7.23e-1 + 1.1667e-2*T9**(-2.d0/3) +
| 2.*T9**(-1.d0/3) + 1.3213e1*T9**(1.d0/3) +
| 1.0117e1*T9**(2.d0/3))*
| exp(-1.2095e1*T9**(-1.d0/3) - 1.7507e1*T9**2)*
| T9**(-2.d0/3) - 3.08e7*
| (1. + 3.5e-2*T9**(1.d0/3) + 3.*T9**(2.d0/3) +
| 7.23e-1*T9 + 9.91*T9**(4.d0/3) + 6.07*T9**(5.d0/3))*
| exp(-1.2095e1*T9**(-1.d0/3) - 1.7507e1*T9**2)*
| T9**(-5.d0/3) + 4.62e7*
| (4.0317*T9**(-4.d0/3) - 3.5013e1*T9)*
| (1. + 3.5e-2*T9**(1.d0/3) + 3.*T9**(2.d0/3) +
| 7.23e-1*T9 + 9.91*T9**(4.d0/3) + 6.07*T9**(5.d0/3))*
| exp(-1.2095e1*T9**(-1.d0/3) - 1.7507e1*T9**2)*
| T9**(-2.d0/3))*exp(-1.8517e2*(1/T9))*T9**(3.d0/2)
c(56)= 'B11(P,N)C11'
dfdt(j,56)= 1.69e6*(3.208e3 - 1.5398e2*Sqrt(T9) +
| 3.208e1*T9 -
| 2.4*T9**(3.d0/2) + 1.*T9**2)*exp(-3.208e1*(1/T9))*
| T9**(-2)
drdt(j,56)= 1.6866e8*(1.e-2 - 2.4e-2*(1/Sqrt(T9)))
c(57)= 'B11(P,A)2HE4'
dfdt(j,57)= 4.2143e10*exp(-6.262*(1/T9))*T9**(-7.d0/2) +
| 6.988e6*exp(-1.734*(1/T9))*T9**(-7.d0/2) +
| 5.4918e10*exp(-1.4154e1*(1/T9))*T9**(-3) -
| 1.0095e10*exp(-6.262*(1/T9))*T9**(-5.d0/2) -
| 6.045e6*exp(-1.734*(1/T9))*T9**(-5.d0/2) -
| 3.88e9*exp(-1.4154e1*(1/T9))*T9**(-2) +
| 4.2533e11*(5.8621e-2 + 4.8276e-1*T9**(1.d0/3) +
| 1.7586e-1*T9**(2.d0/3) + 1.3103*T9 + 1.*T9**(4.d0/3))*
| exp(-1.2095e1*T9**(-1.d0/3) - 3.7e-1*T9**2)*T9**(-4.d0/3)
| - 1.4667e12*(1. + 3.4e-2*T9**(1.d0/3) +
| 1.4e-1*T9**(2.d0/3) + 3.4e-2*T9 + 1.9e-1*T9**(4.d0/3) +
| 1.16e-1*T9**(5.d0/3))*
| exp(-1.2095e1*T9**(-1.d0/3) - 3.7e-1*T9**2)*T9**(-5.d0/3)
| + 1.8885e11*(8.6207 + 2.931e-1*T9**(1.d0/3) +
| 1.2069*T9**(2.d0/3) + 2.931e-1*T9 +
| 1.6379*T9**(4.d0/3) + 1.*T9**(5.d0/3))*
| (5.4483 - 1.*T9**(7.d0/3))*
| exp(-1.2095e1*T9**(-1.d0/3) - 3.7e-1*T9**2)*T9**(-2)
drdt(j,57)= 3.5264e-8*(6.73e9*exp(-6.262*(1/T9))*
| T9**(-3.d0/2) +
| 4.03e6*exp(-1.734*(1/T9))*T9**(-3.d0/2) +
| 3.88e9*exp(-1.4154e1*(1/T9))*(1/T9) +
| 2.2e12*(1. + 3.4e-2*T9**(1.d0/3) + 1.4e-1*T9**(2.d0/3) +
| 3.4e-2*T9 + 1.9e-1*T9**(4.d0/3) + 1.16e-1*T9**(5.d0/3)
| )*exp(-1.2095e1*T9**(-1.d0/3) - 3.7e-1*T9**2)*
| T9**(-2.d0/3))*exp(-1.0075e2*(1/T9))*T9**(-7.d0/2) -
| 5.25e-10*(6.73e9*exp(-6.262*(1/T9))*T9**(-3.d0/2) +
| 4.03e6*exp(-1.734*(1/T9))*T9**(-3.d0/2) +
| 3.88e9*exp(-1.4154e1*(1/T9))*(1/T9) +
| 2.2e12*(1. + 3.4e-2*T9**(1.d0/3) + 1.4e-1*T9**(2.d0/3) +
| 3.4e-2*T9 + 1.9e-1*T9**(4.d0/3) + 1.16e-1*T9**(5.d0/3)
| )*exp(-1.2095e1*T9**(-1.d0/3) - 3.7e-1*T9**2)*
| T9**(-2.d0/3))*exp(-1.0075e2*(1/T9))*T9**(-5.d0/2) +
| 3.5e-10*(4.2143e10*exp(-6.262*(1/T9))*T9**(-7.d0/2) +
| 6.988e6*exp(-1.734*(1/T9))*T9**(-7.d0/2) +
| 5.4918e10*exp(-1.4154e1*(1/T9))*T9**(-3) -
| 1.0095e10*exp(-6.262*(1/T9))*T9**(-5.d0/2) -
| 6.045e6*exp(-1.734*(1/T9))*T9**(-5.d0/2) -
| 3.88e9*exp(-1.4154e1*(1/T9))*T9**(-2) +
| 2.2e12*(3.4e-2 + 1.1333e-2*T9**(-2.d0/3) +
| 9.3333e-2*T9**(-1.d0/3) + 2.5333e-1*T9**(1.d0/3) +
| 1.9333e-1*T9**(2.d0/3))*
| exp(-1.2095e1*T9**(-1.d0/3) - 3.7e-1*T9**2)*
| T9**(-2.d0/3) - 1.4667e12*
| (1. + 3.4e-2*T9**(1.d0/3) + 1.4e-1*T9**(2.d0/3) +
| 3.4e-2*T9 + 1.9e-1*T9**(4.d0/3) + 1.16e-1*T9**(5.d0/3)
| )*exp(-1.2095e1*T9**(-1.d0/3) - 3.7e-1*T9**2)*
| T9**(-5.d0/3) + 2.2e12*
| (4.0317*T9**(-4.d0/3) - 7.3999e-1*T9)*
| (1. + 3.4e-2*T9**(1.d0/3) + 1.4e-1*T9**(2.d0/3) +
| 3.4e-2*T9 + 1.9e-1*T9**(4.d0/3) + 1.16e-1*T9**(5.d0/3)
| )*exp(-1.2095e1*T9**(-1.d0/3) - 3.7e-1*T9**2)*
| T9**(-2.d0/3))*exp(-1.0075e2*(1/T9))*T9**(-3.d0/2)
c(58)= 'B11(A,N)N14'
dfdt(j,58)= 8.8544e3*exp(-5.178*(1/T9))*T9**(-7.d0/2) +
| 5.0603*exp(-2.827*(1/T9))*T9**(-7.d0/2) -
| 2.565e3*exp(-5.178*(1/T9))*T9**(-5.d0/2) -
| 2.685*exp(-2.827*(1/T9))*T9**(-5.d0/2) +
| 3.8596e7*exp(-8.596*(1/T9))*T9**(-7.d0/5) +
| 2.694e6*exp(-8.596*(1/T9))*T9**(-2.d0/5) +
| 1.2146e14*(2.8692e-4 + 3.1044e-1*T9**(1.d0/3) +
| 4.8087e-2*T9**(2.d0/3) + 3.0454*T9 + 1.*T9**(4.d0/3))*
| exp(-2.8234e1*T9**(-1.d0/3) - 5.102e1*T9**2)*T9**(-4.d0/3)
| - 4.6467e12*(1. + 1.5e-2*T9**(1.d0/3) +
| 8.115*T9**(2.d0/3) + 8.38e-1*T9 +
| 3.9804e1*T9**(4.d0/3) + 1.0456e1*T9**(5.d0/3))*
| exp(-2.8234e1*T9**(-1.d0/3) - 5.102e1*T9**2)*T9**(-5.d0/3)
| + 7.4366e15*(9.5639e-2 + 1.4346e-3*T9**(1.d0/3) +
| 7.7611e-1*T9**(2.d0/3) + 8.0145e-2*T9 +
| 3.8068*T9**(4.d0/3) + 1.*T9**(5.d0/3))*
| (9.2231e-2 - 1.*T9**(7.d0/3))*
| exp(-2.8234e1*T9**(-1.d0/3) - 5.102e1*T9**2)*T9**(-2)
drdt(j,58)= 6.7344*(1.71e3*exp(-5.178*(1/T9))
| *T9**(-3.d0/2) +
| 1.79*exp(-2.827*(1/T9))*T9**(-3.d0/2) +
| 4.49e6*exp(-8.596*(1/T9))*T9**(3.d0/5) +
| 6.97e12*(1. + 1.5e-2*T9**(1.d0/3) + 8.115*T9**(2.d0/3) +
| 8.38e-1*T9 + 3.9804e1*T9**(4.d0/3) +
| 1.0456e1*T9**(5.d0/3))*
| exp(-2.8234e1*T9**(-1.d0/3) - 5.102e1*T9**2)*
| T9**(-2.d0/3))*exp(-1.835*(1/T9))*T9**(-2) +
| 3.67*(8.8544e3*exp(-5.178*(1/T9))*T9**(-7.d0/2) +
| 5.0603*exp(-2.827*(1/T9))*T9**(-7.d0/2) -
| 2.565e3*exp(-5.178*(1/T9))*T9**(-5.d0/2) -
| 2.685*exp(-2.827*(1/T9))*T9**(-5.d0/2) +
| 3.8596e7*exp(-8.596*(1/T9))*T9**(-7.d0/5) +
| 6.97e12*(8.38e-1 + 5.e-3*T9**(-2.d0/3) +
| 5.41*T9**(-1.d0/3) + 5.3072e1*T9**(1.d0/3) +
| 1.7427e1*T9**(2.d0/3))*
| exp(-2.8234e1*T9**(-1.d0/3) - 5.102e1*T9**2)*
| T9**(-2.d0/3) + 2.694e6*exp(-8.596*(1/T9))*
| T9**(-2.d0/5) - 4.6467e12*
| (1. + 1.5e-2*T9**(1.d0/3) + 8.115*T9**(2.d0/3) +
| 8.38e-1*T9 + 3.9804e1*T9**(4.d0/3) +
| 1.0456e1*T9**(5.d0/3))*
| exp(-2.8234e1*T9**(-1.d0/3) - 5.102e1*T9**2)*
| T9**(-5.d0/3) + 6.97e12*
| (9.4113*T9**(-4.d0/3) - 1.0204e2*T9)*
| (1. + 1.5e-2*T9**(1.d0/3) + 8.115*T9**(2.d0/3) +
| 8.38e-1*T9 + 3.9804e1*T9**(4.d0/3) +
| 1.0456e1*T9**(5.d0/3))*
| exp(-2.8234e1*T9**(-1.d0/3) - 5.102e1*T9**2)*
| T9**(-2.d0/3))*exp(-1.835*(1/T9))
c(59)= 'B11(A,P)C14'
dfdt(j,59)= 1.7398e3*exp(-5.178*(1/T9))*T9**(-7.d0/2) +
| 1.5379e-2*exp(-2.827*(1/T9))*T9**(-7.d0/2) -
| 5.04e2*exp(-5.178*(1/T9))*T9**(-5.d0/2) -
| 8.16e-3*exp(-2.827*(1/T9))*T9**(-5.d0/2) +
| 6.1802e7*exp(-1.1617e1*(1/T9))*T9**(-19.d0/8) -
| 1.995e6*exp(-1.1617e1*(1/T9))*T9**(-11.d0/8) +
| 3.7357e12*(7.1874e-4 + 5.3426e-1*T9**(1.d0/3) +
| 8.2798e-2*T9**(2.d0/3) + 3.0451*T9 + 1.*T9**(4.d0/3))*
| exp(-2.8234e1*T9**(-1.d0/3) - 8.305*T9**2)*T9**(-4.d0/3) -
| 3.58e11*(1. + 1.5e-2*T9**(1.d0/3) + 5.575*T9**(2.d0/3) +
| 5.76e-1*T9 + 1.5888e1*T9**(4.d0/3) + 4.174*T9**(5.d0/3))*
| exp(-2.8234e1*T9**(-1.d0/3) - 8.305*T9**2)*T9**(-5.d0/3) +
| 3.723e13*(2.3958e-1 + 3.5937e-3*T9**(1.d0/3) +
| 1.3356*T9**(2.d0/3) + 1.38e-1*T9 + 3.8064*T9**(4.d0/3) +
| 1.*T9**(5.d0/3))*(5.666e-1 - 1.*T9**(7.d0/3))*
| exp(-2.8234e1*T9**(-1.d0/3) - 8.305*T9**2)*T9**(-2)
drdt(j,59)= 1.0008e2*(3.36e2*exp(-5.178*(1/T9))*
| T9**(-3.d0/2) +
| 5.44e-3*exp(-2.827*(1/T9))*T9**(-3.d0/2) +
| 5.32e6*exp(-1.1617e1*(1/T9))*T9**(-3.d0/8) +
| 5.37e11*(1. + 1.5e-2*T9**(1.d0/3) + 5.575*T9**(2.d0/3) +
| 5.76e-1*T9 + 1.5888e1*T9**(4.d0/3) +
| 4.174*T9**(5.d0/3))*
| exp(-2.8234e1*T9**(-1.d0/3) - 8.305*T9**2)*T9**(-2.d0/3)
| )*exp(-9.098*(1/T9))*T9**(-2) +
| 1.1e1*(1.7398e3*exp(-5.178*(1/T9))*T9**(-7.d0/2) +
| 1.5379e-2*exp(-2.827*(1/T9))*T9**(-7.d0/2) -
| 5.04e2*exp(-5.178*(1/T9))*T9**(-5.d0/2) -
| 8.16e-3*exp(-2.827*(1/T9))*T9**(-5.d0/2) +
| 6.1802e7*exp(-1.1617e1*(1/T9))*T9**(-19.d0/8) -
| 1.995e6*exp(-1.1617e1*(1/T9))*T9**(-11.d0/8) +
| 5.37e11*(5.76e-1 + 5.e-3*T9**(-2.d0/3) +
| 3.7167*T9**(-1.d0/3) + 2.1184e1*T9**(1.d0/3) +
| 6.9567*T9**(2.d0/3))*
| exp(-2.8234e1*T9**(-1.d0/3) - 8.305*T9**2)*T9**(-2.d0/3)
| - 3.58e11*(1. + 1.5e-2*T9**(1.d0/3) +
| 5.575*T9**(2.d0/3) + 5.76e-1*T9 +
| 1.5888e1*T9**(4.d0/3) + 4.174*T9**(5.d0/3))*
| exp(-2.8234e1*T9**(-1.d0/3) - 8.305*T9**2)*T9**(-5.d0/3)
| + 5.37e11*(9.4113*T9**(-4.d0/3) - 1.661e1*T9)*
| (1. + 1.5e-2*T9**(1.d0/3) + 5.575*T9**(2.d0/3) +
| 5.76e-1*T9 + 1.5888e1*T9**(4.d0/3) +
| 4.174*T9**(5.d0/3))*
| exp(-2.8234e1*T9**(-1.d0/3) - 8.305*T9**2)*T9**(-2.d0/3)
| )*exp(-9.098*(1/T9))
c(60)= 'C11(P,G)N12'
dfdt(j,60)= 6.2066e4*exp(-7.021*(1/T9))*T9**(-7.d0/2) -
| 1.326e4*exp(-7.021*(1/T9))*T9**(-5.d0/2) +
| 1.7667e5*(2.48e-3 + 4.976e-1*T9**(1.d0/3) +
| 1.596e-1*T9**(2.d0/3) + 1.4752*T9 + 1.*T9**(4.d0/3))*
| exp(-1.3658e1*T9**(-1.d0/3) - 3.7777e-1*T9**2)*
| T9**(-4.d0/3) - 2.8267e4*
| (1. + 3.1e-2*T9**(1.d0/3) + 3.11*T9**(2.d0/3) +
| 6.65e-1*T9 + 4.61*T9**(4.d0/3) + 2.5*T9**(5.d0/3))*
| exp(-1.3658e1*T9**(-1.d0/3) - 3.7777e-1*T9**2)*
| T9**(-5.d0/3) + 8.0087e4*
| (4.e-1 + 1.24e-2*T9**(1.d0/3) + 1.244*T9**(2.d0/3) +
| 2.66e-1*T9 + 1.844*T9**(4.d0/3) + 1.*T9**(5.d0/3))*
| (6.0257 - 1.*T9**(7.d0/3))*
| exp(-1.3658e1*T9**(-1.d0/3) - 3.7777e-1*T9**2)*T9**(-2)
drdt(j,60)= 1.6252e11*(8.84e3*exp(-7.021*(1/T9))*
| T9**(-3.d0/2) +
| 4.24e4*(1. + 3.1e-2*T9**(1.d0/3) + 3.11*T9**(2.d0/3) +
| 6.65e-1*T9 + 4.61*T9**(4.d0/3) + 2.5*T9**(5.d0/3))*
| exp(-1.3658e1*T9**(-1.d0/3) - 3.7777e-1*T9**2)*
| T9**(-2.d0/3))*exp(-6.975*(1/T9))*(1/Sqrt(T9)) +
| 3.495e10*(8.84e3*exp(-7.021*(1/T9))*T9**(-3.d0/2) +
| 4.24e4*(1. + 3.1e-2*T9**(1.d0/3) + 3.11*T9**(2.d0/3) +
| 6.65e-1*T9 + 4.61*T9**(4.d0/3) + 2.5*T9**(5.d0/3))*
| exp(-1.3658e1*T9**(-1.d0/3) - 3.7777e-1*T9**2)*
| T9**(-2.d0/3))*exp(-6.975*(1/T9))*Sqrt(T9) +
| 2.33e10*(6.2066e4*exp(-7.021*(1/T9))*T9**(-7.d0/2) -
| 1.326e4*exp(-7.021*(1/T9))*T9**(-5.d0/2) +
| 4.24e4*(6.65e-1 + 1.0333e-2*T9**(-2.d0/3) +
| 2.0733*T9**(-1.d0/3) + 6.1467*T9**(1.d0/3) +
| 4.1667*T9**(2.d0/3))*
| exp(-1.3658e1*T9**(-1.d0/3) - 3.7777e-1*T9**2)*
| T9**(-2.d0/3) - 2.8267e4*
| (1. + 3.1e-2*T9**(1.d0/3) + 3.11*T9**(2.d0/3) +
| 6.65e-1*T9 + 4.61*T9**(4.d0/3) + 2.5*T9**(5.d0/3))*
| exp(-1.3658e1*T9**(-1.d0/3) - 3.7777e-1*T9**2)*
| T9**(-5.d0/3) + 4.24e4*
| (4.5527*T9**(-4.d0/3) - 7.5554e-1*T9)*
| (1. + 3.1e-2*T9**(1.d0/3) + 3.11*T9**(2.d0/3) +
| 6.65e-1*T9 + 4.61*T9**(4.d0/3) + 2.5*T9**(5.d0/3))*
| exp(-1.3658e1*T9**(-1.d0/3) - 3.7777e-1*T9**2)*
| T9**(-2.d0/3))*exp(-6.975*(1/T9))*T9**(3.d0/2)
|
c(61)= 'C12(P,G)N13'
dfdt(j,61)= 3.9085e6*exp(-1.8179e1*(1/T9))*T9**(-7.d0/2) +
| 5.319e5*exp(-4.925*(1/T9))*T9**(-7.d0/2) -
| 3.225e5*exp(-1.8179e1*(1/T9))*T9**(-5.d0/2) -
| 1.62e5*exp(-4.925*(1/T9))*T9**(-5.d0/2) +
| 3.808e7*(5.3571e-3 + 4.25e-1*T9**(1.d0/3) +
| 1.3607e-1*T9**(2.d0/3) + 1.4714*T9 + 1.*T9**(4.d0/3))*
| exp(-1.369e1*T9**(-1.d0/3) - 4.4444e-1*T9**2)*T9**(-4.d0/3)
| - 1.36e7*(1. + 3.e-2*T9**(1.d0/3) + 1.19*T9**(2.d0/3) +
| 2.54e-1*T9 + 2.06*T9**(4.d0/3) + 1.12*T9**(5.d0/3))*
| exp(-1.369e1*T9**(-1.d0/3) - 4.4444e-1*T9**2)*T9**(-5.d0/3)
| + 2.0309e7*(8.9286e-1 + 2.6786e-2*T9**(1.d0/3) +
| 1.0625*T9**(2.d0/3) + 2.2679e-1*T9 +
| 1.8393*T9**(4.d0/3) + 1.*T9**(5.d0/3))*
| (5.1337 - 1.*T9**(7.d0/3))*
| exp(-1.369e1*T9**(-1.d0/3) - 4.4444e-1*T9**2)*T9**(-2)
drdt(j,61)= 1.9937e11*(2.15e5*exp(-1.8179e1*(1/T9))*
| T9**(-3.d0/2) +
| 1.08e5*exp(-4.925*(1/T9))*T9**(-3.d0/2) +
| 2.04e7*(1. + 3.e-2*T9**(1.d0/3) + 1.19*T9**(2.d0/3) +
| 2.54e-1*T9 + 2.06*T9**(4.d0/3) + 1.12*T9**(5.d0/3))*
| exp(-1.369e1*T9**(-1.d0/3) - 4.4444e-1*T9**2)*
| T9**(-2.d0/3))*exp(-2.2553e1*(1/T9))*(1/Sqrt(T9)) +
| 1.326e10*(2.15e5*exp(-1.8179e1*(1/T9))*T9**(-3.d0/2) +
| 1.08e5*exp(-4.925*(1/T9))*T9**(-3.d0/2) +
| 2.04e7*(1. + 3.e-2*T9**(1.d0/3) + 1.19*T9**(2.d0/3) +
| 2.54e-1*T9 + 2.06*T9**(4.d0/3) + 1.12*T9**(5.d0/3))*
| exp(-1.369e1*T9**(-1.d0/3) - 4.4444e-1*T9**2)*
| T9**(-2.d0/3))*exp(-2.2553e1*(1/T9))*Sqrt(T9) +
| 8.84e9*(3.9085e6*exp(-1.8179e1*(1/T9))*T9**(-7.d0/2) +
| 5.319e5*exp(-4.925*(1/T9))*T9**(-7.d0/2) -
| 3.225e5*exp(-1.8179e1*(1/T9))*T9**(-5.d0/2) -
| 1.62e5*exp(-4.925*(1/T9))*T9**(-5.d0/2) +
| 2.04e7*(2.54e-1 + 1.e-2*T9**(-2.d0/3) +
| 7.9333e-1*T9**(-1.d0/3) + 2.7467*T9**(1.d0/3) +
| 1.8667*T9**(2.d0/3))*
| exp(-1.369e1*T9**(-1.d0/3) - 4.4444e-1*T9**2)*
| T9**(-2.d0/3) - 1.36e7*
| (1. + 3.e-2*T9**(1.d0/3) + 1.19*T9**(2.d0/3) +
| 2.54e-1*T9 + 2.06*T9**(4.d0/3) + 1.12*T9**(5.d0/3))*
| exp(-1.369e1*T9**(-1.d0/3) - 4.4444e-1*T9**2)*
| T9**(-5.d0/3) + 2.04e7*
| (4.5633*T9**(-4.d0/3) - 8.8889e-1*T9)*
| (1. + 3.e-2*T9**(1.d0/3) + 1.19*T9**(2.d0/3) +
| 2.54e-1*T9 + 2.06*T9**(4.d0/3) + 1.12*T9**(5.d0/3))*
| exp(-1.369e1*T9**(-1.d0/3) - 4.4444e-1*T9**2)*
| T9**(-2.d0/3))*exp(-2.2553e1*(1/T9))*T9**(3.d0/2)
c(62)= 'C12(A,G)O16'
dfdt(j,62)= 3.4374e4*exp(-2.7499e1*(1/T9))*T9**(-7.d0/2) -
| 2.08e8*exp(-3.212e1*T9**(-1.d0/3) - 8.182e-2*T9**2)*
| (1. + 4.89e-2*T9**(-2.d0/3))**(-2)*T9**(-3) -
| 3.52e8*exp(-3.212e1*T9**(-1.d0/3))*
| (1. + 2.654e-1*T9**(-2.d0/3))**(-2)*T9**(-3) -
| 1.875e3*exp(-2.7499e1*(1/T9))*T9**(-5.d0/2) +
| 1.8844e9*exp(-3.212e1*T9**(-1.d0/3))*
| (2.654e-1 + 1.*T9**(2.d0/3))**(-2)*T9**(-2) +
| 6.7808e6*exp(-3.212e1*T9**(-1.d0/3) - 8.182e-2*T9**2)*
| (4.89e-2 + 1.*T9**(2.d0/3))**(-3)*T9**(-5.d0/3) +
| 6.2281e7*exp(-3.212e1*T9**(-1.d0/3))*
| (2.654e-1 + 1.*T9**(2.d0/3))**(-3)*T9**(-5.d0/3) +
| 2.2224e-1*exp(-1.5541e1*(1/T9))*T9**3 +
| 7.15e-2*exp(-1.5541e1*(1/T9))*T9**4 +
| 1.7018e7*(6.5429e1 - 1.*T9**(7.d0/3))*
| exp(-3.212e1*T9**(-1.d0/3) - 8.182e-2*T9**2)*
| (4.89e-2 + 1.*T9**(2.d0/3))**(-2)*T9**(-2)
drdt(j,62)= 5.13e10*(6.7808e6*exp
| (-3.212e1*T9**(-1.d0/3) - 8.182e-2*T9**2)*
| (1. + 4.89e-2*T9**(-2.d0/3))**(-3)*T9**(-11.d0/3) +
| 6.2281e7*exp(-3.212e1*T9**(-1.d0/3))*
| (1. + 2.654e-1*T9**(-2.d0/3))**(-3)*T9**(-11.d0/3) +
| 3.4374e4*exp(-2.7499e1*(1/T9))*T9**(-7.d0/2) +
| 1.8844e9*exp(-3.212e1*T9**(-1.d0/3))*
| (1. + 2.654e-1*T9**(-2.d0/3))**(-2)*T9**(-10.d0/3) -
| 2.08e8*exp(-3.212e1*T9**(-1.d0/3) - 8.182e-2*T9**2)*
| (1. + 4.89e-2*T9**(-2.d0/3))**(-2)*T9**(-3) -
| 3.52e8*exp(-3.212e1*T9**(-1.d0/3))*
| (1. + 2.654e-1*T9**(-2.d0/3))**(-2)*T9**(-3) -
| 1.875e3*exp(-2.7499e1*(1/T9))*T9**(-5.d0/2) +
| 1.04e8*(1.0707e1*T9**(-4.d0/3) - 1.6364e-1*T9)*
| exp(-3.212e1*T9**(-1.d0/3) - 8.182e-2*T9**2)*
| (1. + 4.89e-2*T9**(-2.d0/3))**(-2)*T9**(-2) +
| 2.2224e-1*exp(-1.5541e1*(1/T9))*T9**3 +
| 7.15e-2*exp(-1.5541e1*(1/T9))*T9**4)*
| exp(-8.3111e1*(1/T9))*T9**(3.d0/2) +
| 4.2636e12*(1.04e8*exp
| (-3.212e1*T9**(-1.d0/3) - 8.182e-2*T9**2)*
| (1. + 4.89e-2*T9**(-2.d0/3))**(-2)*T9**(-2) +
| 1.76e8*exp(-3.212e1*T9**(-1.d0/3))*
| (1. + 2.654e-1*T9**(-2.d0/3))**(-2)*T9**(-2) +
| 1.25e3*exp(-2.7499e1*(1/T9))*T9**(-3.d0/2) +
| 1.43e-2*exp(-1.5541e1*(1/T9))*T9**5)*
| exp(-8.3111e1*(1/T9))*(1/Sqrt(T9)) +
| 7.695e10*(1.04e8*exp
| (-3.212e1*T9**(-1.d0/3) - 8.182e-2*T9**2)*
| (1. + 4.89e-2*T9**(-2.d0/3))**(-2)*T9**(-2) +
| 1.76e8*exp(-3.212e1*T9**(-1.d0/3))*
| (1. + 2.654e-1*T9**(-2.d0/3))**(-2)*T9**(-2) +
| 1.25e3*exp(-2.7499e1*(1/T9))*T9**(-3.d0/2) +
| 1.43e-2*exp(-1.5541e1*(1/T9))*T9**5)*
| exp(-8.3111e1*(1/T9))*Sqrt(T9)
c(63)= 'C12(A,N)O15'
dfdt(j,63)= 3.72e5*(6.5774e3 + 1.2366e3*Sqrt(T9) +
| 9.8661e1*T9 +
| 6.2667*T9**(3.d0/2) + 1.*T9**2)*exp(-9.8661e1*(1/T9))*
| T9**(-2)
drdt(j,63)= 3.4968e7*(1.5e-2 + 9.4e-2*(1/Sqrt(T9)))
c(64)= 'C13(P,G)N14'
dfdt(j,64)= 6.8982e6*exp(-5.701*(1/T9))*T9**(-16.d0/5) -
| 1.452e6*exp(-5.701*(1/T9))*T9**(-11.d0/5) +
| 1.0053e8*(7.9681e-3 + 5.089e-1*T9**(1.d0/3) +
| 1.6255e-1*T9**(2.d0/3) + 1.4768*T9 + 1.*T9**(4.d0/3))*
| exp(-1.3717e1*T9**(-1.d0/3) - 2.5e-1*T9**2)*T9**(-4.d0/3)
| - 5.34e7*(1. + 3.e-2*T9**(1.d0/3) + 9.58e-1*T9**(2.d0/3) +
| 2.04e-1*T9 + 1.39*T9**(4.d0/3) + 7.53e-1*T9**(5.d0/3))*
| exp(-1.3717e1*T9**(-1.d0/3) - 2.5e-1*T9**2)*T9**(-5.d0/3)
| + 3.0158e7*(1.328 + 3.9841e-2*T9**(1.d0/3) +
| 1.2722*T9**(2.d0/3) + 2.7092e-1*T9 +
| 1.8459*T9**(4.d0/3) + 1.*T9**(5.d0/3))*
| (9.1447 - 1.*T9**(7.d0/3))*
| exp(-1.3717e1*T9**(-1.d0/3) - 2.5e-1*T9**2)*T9**(-2)
drdt(j,64)= 1.0427e12*(1.21e6*exp(-5.701*(1/T9))*
| T9**(-6.d0/5) +
| 8.01e7*(1. + 3.e-2*T9**(1.d0/3) + 9.58e-1*T9**(2.d0/3) +
| 2.04e-1*T9 + 1.39*T9**(4.d0/3) + 7.53e-1*T9**(5.d0/3))
| *exp(-1.3717e1*T9**(-1.d0/3) - 2.5e-1*T9**2)*
| T9**(-2.d0/3))*exp(-8.7621e1*(1/T9))*(1/Sqrt(T9)) +
| 1.785e10*(1.21e6*exp(-5.701*(1/T9))*T9**(-6.d0/5) +
| 8.01e7*(1. + 3.e-2*T9**(1.d0/3) + 9.58e-1*T9**(2.d0/3) +
| 2.04e-1*T9 + 1.39*T9**(4.d0/3) + 7.53e-1*T9**(5.d0/3))
| *exp(-1.3717e1*T9**(-1.d0/3) - 2.5e-1*T9**2)*
| T9**(-2.d0/3))*exp(-8.7621e1*(1/T9))*Sqrt(T9) +
| 1.19e10*(6.8982e6*exp(-5.701*(1/T9))*T9**(-16.d0/5) -
| 1.452e6*exp(-5.701*(1/T9))*T9**(-11.d0/5) +
| 8.01e7*(2.04e-1 + 1.e-2*T9**(-2.d0/3) +
| 6.3867e-1*T9**(-1.d0/3) + 1.8533*T9**(1.d0/3) +
| 1.255*T9**(2.d0/3))*
| exp(-1.3717e1*T9**(-1.d0/3) - 2.5e-1*T9**2)*
| T9**(-2.d0/3) - 5.34e7*
| (1. + 3.e-2*T9**(1.d0/3) + 9.58e-1*T9**(2.d0/3) +
| 2.04e-1*T9 + 1.39*T9**(4.d0/3) + 7.53e-1*T9**(5.d0/3))
| *exp(-1.3717e1*T9**(-1.d0/3) - 2.5e-1*T9**2)*
| T9**(-5.d0/3) + 8.01e7*
| (4.5723*T9**(-4.d0/3) - 5.e-1*T9)*
| (1. + 3.e-2*T9**(1.d0/3) + 9.58e-1*T9**(2.d0/3) +
| 2.04e-1*T9 + 1.39*T9**(4.d0/3) + 7.53e-1*T9**(5.d0/3))
| *exp(-1.3717e1*T9**(-1.d0/3) - 2.5e-1*T9**2)*
| T9**(-2.d0/3))*exp(-8.7621e1*(1/T9))*T9**(3.d0/2)
c(65)= 'C13(P,N)N13'
dfdt(j,65)= 6.956e6*(9.4178e2 - 1.5728e2*Sqrt(T9) +
| 3.4846e1*T9 -
| 2.2568*T9**(3.d0/2) + 1.*T9**2)*exp(-3.4846e1*(1/T9))*
| T9**(-2)
drdt(j,65)= 1.8762e8*(3.7e-2 - 8.35e-2*(1/Sqrt(T9)))
c(66)= 'C13(A,N)O16'
dfdt(j,66)= 8.5506e10*exp(-2.9283e1*(1/T9))*T9**(-7.d0/2) +
| 4.0818e10*exp(-2.0409e1*(1/T9))*T9**(-7.d0/2) +
| 1.6741e7*exp(-1.1873e1*(1/T9))*T9**(-7.d0/2) +
| 3.5805e6*exp(-9.373*(1/T9))*T9**(-7.d0/2) -
| 4.38e9*exp(-2.9283e1*(1/T9))*T9**(-5.d0/2) -
| 3.e9*exp(-2.0409e1*(1/T9))*T9**(-5.d0/2) -
| 2.115e6*exp(-1.1873e1*(1/T9))*T9**(-5.d0/2) -
| 5.73e5*exp(-9.373*(1/T9))*T9**(-5.d0/2) -
| 4.5133e15*(1. + 1.3e-2*T9**(1.d0/3) + 2.04*T9**(2.d0/3) +
| 1.84e-1*T9)*exp
| (-3.2329e1*T9**(-1.d0/3) - 6.0655e-1*T9**2)*T9**(-5.d0/3)
| + 1.2457e15*(2.3551e-2 + 7.3913*T9**(1.d0/3) +
| 1.*T9**(2.d0/3))*exp
| (-3.2329e1*T9**(-1.d0/3) - 6.0655e-1*T9**2)*T9**(-4.d0/3)
| + 1.5111e15*(5.4348 + 7.0652e-2*T9**(1.d0/3) +
| 1.1087e1*T9**(2.d0/3) + 1.*T9)*
| (8.8832 - 1.*T9**(7.d0/3))*
| exp(-3.2329e1*T9**(-1.d0/3) - 6.0655e-1*T9**2)*T9**(-2)
drdt(j,66)= 5.79*(8.5506e10*exp(-2.9283e1*(1/T9))*
| T9**(-7.d0/2) +
| 4.0818e10*exp(-2.0409e1*(1/T9))*T9**(-7.d0/2) +
| 1.6741e7*exp(-1.1873e1*(1/T9))*T9**(-7.d0/2) +
| 3.5805e6*exp(-9.373*(1/T9))*T9**(-7.d0/2) -
| 4.38e9*exp(-2.9283e1*(1/T9))*T9**(-5.d0/2) -
| 3.e9*exp(-2.0409e1*(1/T9))*T9**(-5.d0/2) -
| 2.115e6*exp(-1.1873e1*(1/T9))*T9**(-5.d0/2) -
| 5.73e5*exp(-9.373*(1/T9))*T9**(-5.d0/2) -
| 4.5133e15*(1. + 1.3e-2*T9**(1.d0/3) +
| 2.04*T9**(2.d0/3) + 1.84e-1*T9)*
| exp(-3.2329e1*T9**(-1.d0/3) - 6.0655e-1*T9**2)*
| T9**(-5.d0/3) + 6.77e15*
| (1.84e-1 + 4.3333e-3*T9**(-2.d0/3) +
| 1.36*T9**(-1.d0/3))*
| exp(-3.2329e1*T9**(-1.d0/3) - 6.0655e-1*T9**2)*
| T9**(-2.d0/3) + 6.77e15*
| (1.0776e1*T9**(-4.d0/3) - 1.2131*T9)*
| (1. + 1.3e-2*T9**(1.d0/3) + 2.04*T9**(2.d0/3) +
| 1.84e-1*T9)*exp
| (-3.2329e1*T9**(-1.d0/3) - 6.0655e-1*T9**2)*
| T9**(-2.d0/3))*exp(-2.5711e1*(1/T9)) +
| 1.4887e2*(2.92e9*exp(-2.9283e1*(1/T9))*T9**(-3.d0/2) +
| 2.e9*exp(-2.0409e1*(1/T9))*T9**(-3.d0/2) +
| 1.41e6*exp(-1.1873e1*(1/T9))*T9**(-3.d0/2) +
| 3.82e5*exp(-9.373*(1/T9))*T9**(-3.d0/2) +
| 6.77e15*(1. + 1.3e-2*T9**(1.d0/3) + 2.04*T9**(2.d0/3) +
| 1.84e-1*T9)*exp
| (-3.2329e1*T9**(-1.d0/3) - 6.0655e-1*T9**2)*
| T9**(-2.d0/3))*exp(-2.5711e1*(1/T9))*T9**(-2)
c(67)= 'C14(P,G)N15'
dfdt(j,67)= 2.0427e4*exp(-3.811*(1/T9))*T9**(-7.d0/2) -
| 8.04e3*exp(-3.811*(1/T9))*T9**(-5.d0/2) +
| 4.6537e5*exp(-4.739*(1/T9))*T9**(-7.d0/3) -
| 3.2733e4*exp(-4.739*(1/T9))*T9**(-4.d0/3) +
| 1.3033e6*(5.2174e-2 + 1.7496*T9**(1.d0/3) +
| 5.5826e-1*T9**(2.d0/3) + 1.4817*T9 + 1.*T9**(4.d0/3))*
| exp(-1.3741e1*T9**(-1.d0/3) - 3.0553e-2*T9**2)*
| T9**(-4.d0/3) - 4.5333e6*
| (1. + 3.e-2*T9**(1.d0/3) + 5.03e-1*T9**(2.d0/3) +
| 1.07e-1*T9 + 2.13e-1*T9**(4.d0/3) + 1.15e-1*T9**(5.d0/3))
| *exp(-1.3741e1*T9**(-1.d0/3) - 3.0553e-2*T9**2)*
| T9**(-5.d0/3) + 4.7785e4*
| (8.6957 + 2.6087e-1*T9**(1.d0/3) + 4.3739*T9**(2.d0/3) +
| 9.3043e-1*T9 + 1.8522*T9**(4.d0/3) + 1.*T9**(5.d0/3))*
| (7.4957e1 - 1.*T9**(7.d0/3))*
| exp(-1.3741e1*T9**(-1.d0/3) - 3.0553e-2*T9**2)*T9**(-2)
drdt(j,67)= 1.0661e12*(5.36e3*exp(-3.811*(1/T9))*
| T9**(-3.d0/2) +
| 9.82e4*exp(-4.739*(1/T9))*T9**(-1.d0/3) +
| 6.8e6*(1. + 3.e-2*T9**(1.d0/3) + 5.03e-1*T9**(2.d0/3) +
| 1.07e-1*T9 + 2.13e-1*T9**(4.d0/3) +
| 1.15e-1*T9**(5.d0/3))*
| exp(-1.3741e1*T9**(-1.d0/3) - 3.0553e-2*T9**2)*
| T9**(-2.d0/3))*exp(-1.1845e2*(1/T9))*(1/Sqrt(T9)) +
| 1.35e10*(5.36e3*exp(-3.811*(1/T9))*T9**(-3.d0/2) +
| 9.82e4*exp(-4.739*(1/T9))*T9**(-1.d0/3) +
| 6.8e6*(1. + 3.e-2*T9**(1.d0/3) + 5.03e-1*T9**(2.d0/3) +
| 1.07e-1*T9 + 2.13e-1*T9**(4.d0/3) +
| 1.15e-1*T9**(5.d0/3))*
| exp(-1.3741e1*T9**(-1.d0/3) - 3.0553e-2*T9**2)*
| T9**(-2.d0/3))*exp(-1.1845e2*(1/T9))*Sqrt(T9) +
| 9.e9*(2.0427e4*exp(-3.811*(1/T9))*T9**(-7.d0/2) -
| 8.04e3*exp(-3.811*(1/T9))*T9**(-5.d0/2) +
| 4.6537e5*exp(-4.739*(1/T9))*T9**(-7.d0/3) -
| 3.2733e4*exp(-4.739*(1/T9))*T9**(-4.d0/3) +
| 6.8e6*(1.07e-1 + 1.e-2*T9**(-2.d0/3) +
| 3.3533e-1*T9**(-1.d0/3) + 2.84e-1*T9**(1.d0/3) +
| 1.9167e-1*T9**(2.d0/3))*
| exp(-1.3741e1*T9**(-1.d0/3) - 3.0553e-2*T9**2)*
| T9**(-2.d0/3) - 4.5333e6*
| (1. + 3.e-2*T9**(1.d0/3) + 5.03e-1*T9**(2.d0/3) +
| 1.07e-1*T9 + 2.13e-1*T9**(4.d0/3) +
| 1.15e-1*T9**(5.d0/3))*
| exp(-1.3741e1*T9**(-1.d0/3) - 3.0553e-2*T9**2)*
| T9**(-5.d0/3) + 6.8e6*
| (4.5803*T9**(-4.d0/3) - 6.1106e-2*T9)*
| (1. + 3.e-2*T9**(1.d0/3) + 5.03e-1*T9**(2.d0/3) +
| 1.07e-1*T9 + 2.13e-1*T9**(4.d0/3) +
| 1.15e-1*T9**(5.d0/3))*
| exp(-1.3741e1*T9**(-1.d0/3) - 3.0553e-2*T9**2)*
| T9**(-2.d0/3))*exp(-1.1845e2*(1/T9))*T9**(3.d0/2)
c(68)= 'C14(P,N)N14'
dfdt(j,68)= 7.19e5*(5.02e-1 + 1.805e-1*(1/Sqrt(T9)))*
| exp(-7.263*(1/T9)) +
| 4.0902e9*exp(-1.2246e1*(1/T9))*T9**(-5.d0/2) +
| 5.2221e6*(1. + 3.61e-1*Sqrt(T9) + 5.02e-1*T9)*
| exp(-7.263*(1/T9))*T9**(-2) -
| 1.67e8*exp(-1.2246e1*(1/T9))*T9**(-3.d0/2)
drdt(j,68)= 2.3943e5*(5.02e-1 + 1.805e-1*(1/Sqrt(T9))) +
| 5.5422e8*exp(-4.983*(1/T9))*T9**(-5.d0/2) -
| 5.5611e7*exp(-4.983*(1/T9))*T9**(-3.d0/2)
c(69)= 'C14(A,G)O18'
dfdt(j,69)= 1.9026e-7*exp(-2.048*(1/T9))*T9**(-7.d0/2) +
| 3.6577e9*exp(-3.2513e1*T9**(-1.d0/3))*T9**(-10.d0/3) -
| 6.75e8*exp(-3.2513e1*T9**(-1.d0/3))*T9**(-3) +
| 2.7357e4*exp(-9.876*(1/T9))*T9**(-14.d0/5) -
| 1.3935e-7*exp(-2.048*(1/T9))*T9**(-5.d0/2) -
| 2.216e3*exp(-9.876*(1/T9))*T9**(-9.d0/5) +
| 1.8642e8*(3.4973e-2 - 4.7486*T9**(1.d0/3) -
| 6.3852e-1*T9**(2.d0/3) + 3.5082*T9 + 1.*T9**(4.d0/3))*
| exp(-3.2513e1*T9**(-1.d0/3) - 1.4112e-1*T9**2)*
| T9**(-4.d0/3) + 7.4566e7*
| (-1.3661e1 - 1.7486e-1*T9**(1.d0/3) +
| 1.1872e1*T9**(2.d0/3) + 1.0642*T9 -
| 4.3852*T9**(4.d0/3) - 1.*T9**(5.d0/3))*
| exp(-3.2513e1*T9**(-1.d0/3) - 1.4112e-1*T9**2)*
| T9**(-5.d0/3) + 3.1568e7*
| (1.3661e1 + 1.7486e-1*T9**(1.d0/3) -
| 1.1872e1*T9**(2.d0/3) - 1.0642*T9 +
| 4.3852*T9**(4.d0/3) + 1.*T9**(5.d0/3))*
| (3.8399e1 - 1.*T9**(7.d0/3))*
| exp(-3.2513e1*T9**(-1.d0/3) - 1.4112e-1*T9**2)*T9**(-2)
drdt(j,69)= 3.9166e12*(3.375e8*exp(-3.2513e1*T9**(-1.d0/3))
| *T9**(-2) +
| 9.29e-8*exp(-2.048*(1/T9))*T9**(-3.d0/2) +
| 2.77e3*exp(-9.876*(1/T9))*T9**(-4.d0/5) +
| 1.528e9*(1. + 1.28e-2*T9**(1.d0/3) -
| 8.69e-1*T9**(2.d0/3) - 7.79e-2*T9 +
| 3.21e-1*T9**(4.d0/3) + 7.32e-2*T9**(5.d0/3))*
| exp(-3.2513e1*T9**(-1.d0/3) - 1.4112e-1*T9**2)*
| T9**(-2.d0/3))*exp(-7.2262e1*(1/T9))*(1/Sqrt(T9)) +
| 8.13e10*(3.375e8*exp(-3.2513e1*T9**(-1.d0/3))*T9**(-2) +
| 9.29e-8*exp(-2.048*(1/T9))*T9**(-3.d0/2) +
| 2.77e3*exp(-9.876*(1/T9))*T9**(-4.d0/5) +
| 1.528e9*(1. + 1.28e-2*T9**(1.d0/3) -
| 8.69e-1*T9**(2.d0/3) - 7.79e-2*T9 +
| 3.21e-1*T9**(4.d0/3) + 7.32e-2*T9**(5.d0/3))*
| exp(-3.2513e1*T9**(-1.d0/3) - 1.4112e-1*T9**2)*
| T9**(-2.d0/3))*exp(-7.2262e1*(1/T9))*Sqrt(T9) +
| 5.42e10*(1.9026e-7*exp(-2.048*(1/T9))*T9**(-7.d0/2) +
| 3.6577e9*exp(-3.2513e1*T9**(-1.d0/3))*T9**(-10.d0/3) -
| 6.75e8*exp(-3.2513e1*T9**(-1.d0/3))*T9**(-3) +
| 2.7357e4*exp(-9.876*(1/T9))*T9**(-14.d0/5) -
| 1.3935e-7*exp(-2.048*(1/T9))*T9**(-5.d0/2) -
| 2.216e3*exp(-9.876*(1/T9))*T9**(-9.d0/5) +
| 1.528e9*(-7.79e-2 + 4.2667e-3*T9**(-2.d0/3) -
| 5.7933e-1*T9**(-1.d0/3) + 4.28e-1*T9**(1.d0/3) +
| 1.22e-1*T9**(2.d0/3))*
| exp(-3.2513e1*T9**(-1.d0/3) - 1.4112e-1*T9**2)*
| T9**(-2.d0/3) - 1.0187e9*
| (1. + 1.28e-2*T9**(1.d0/3) - 8.69e-1*T9**(2.d0/3) -
| 7.79e-2*T9 + 3.21e-1*T9**(4.d0/3) +
| 7.32e-2*T9**(5.d0/3))*
| exp(-3.2513e1*T9**(-1.d0/3) - 1.4112e-1*T9**2)*
| T9**(-5.d0/3) + 1.528e9*
| (1.0838e1*T9**(-4.d0/3) - 2.8224e-1*T9)*
| (1. + 1.28e-2*T9**(1.d0/3) - 8.69e-1*T9**(2.d0/3) -
| 7.79e-2*T9 + 3.21e-1*T9**(4.d0/3) +
| 7.32e-2*T9**(5.d0/3))*
| exp(-3.2513e1*T9**(-1.d0/3) - 1.4112e-1*T9**2)*
| T9**(-2.d0/3))*exp(-7.2262e1*(1/T9))*T9**(3.d0/2)
c(70)= ' N13(P,G)O14'
dfdt(j,70)= 1.5426e6*exp(-6.348*(1/T9))*T9**(-7.d0/2) -
| 3.645e5*exp(-6.348*(1/T9))*T9**(-5.d0/2) +
| 1.6429e8*(2.2131e-3 - 1.3164e-1*T9**(1.d0/3) -
| 3.7869e-2*T9**(2.d0/3) + 1.6393*T9 + 1.*T9**(4.d0/3))*
| exp(-1.5202e1*T9**(-1.d0/3) - 7.0498e-1*T9**2)*
| T9**(-4.d0/3) + 6.5717e7*
| (-4.0984e-1 - 1.1066e-2*T9**(1.d0/3) +
| 3.291e-1*T9**(2.d0/3) + 6.3115e-2*T9 -
| 2.0492*T9**(4.d0/3) - 1.*T9**(5.d0/3))*
| exp(-1.5202e1*T9**(-1.d0/3) - 7.0498e-1*T9**2)*
| T9**(-5.d0/3) + 1.3899e8*
| (4.0984e-1 + 1.1066e-2*T9**(1.d0/3) -
| 3.291e-1*T9**(2.d0/3) - 6.3115e-2*T9 +
| 2.0492*T9**(4.d0/3) + 1.*T9**(5.d0/3))*
| (3.594 - 1.*T9**(7.d0/3))*
| exp(-1.5202e1*T9**(-1.d0/3) - 7.0498e-1*T9**2)*T9**(-2)
drdt(j,70)= 1.9173e12*(2.43e5*exp(-6.348*(1/T9))*
| T9**(-3.d0/2) +
| 4.04e7*(1. + 2.7e-2*T9**(1.d0/3) -
| 8.03e-1*T9**(2.d0/3) - 1.54e-1*T9 + 5.*T9**(4.d0/3) +
| 2.44*T9**(5.d0/3))*
| exp(-1.5202e1*T9**(-1.d0/3) - 7.0498e-1*T9**2)*
| T9**(-2.d0/3))*exp(-5.3706e1*(1/T9))*(1/Sqrt(T9)) +
| 5.355e10*(2.43e5*exp(-6.348*(1/T9))*T9**(-3.d0/2) +
| 4.04e7*(1. + 2.7e-2*T9**(1.d0/3) -
| 8.03e-1*T9**(2.d0/3) - 1.54e-1*T9 + 5.*T9**(4.d0/3) +
| 2.44*T9**(5.d0/3))*
| exp(-1.5202e1*T9**(-1.d0/3) - 7.0498e-1*T9**2)*
| T9**(-2.d0/3))*exp(-5.3706e1*(1/T9))*Sqrt(T9) +
| 3.57e10*(1.5426e6*exp(-6.348*(1/T9))*T9**(-7.d0/2) -
| 3.645e5*exp(-6.348*(1/T9))*T9**(-5.d0/2) +
| 4.04e7*(-1.54e-1 + 9.e-3*T9**(-2.d0/3) -
| 5.3533e-1*T9**(-1.d0/3) + 6.6667*T9**(1.d0/3) +
| 4.0667*T9**(2.d0/3))*
| exp(-1.5202e1*T9**(-1.d0/3) - 7.0498e-1*T9**2)*
| T9**(-2.d0/3) - 2.6933e7*
| (1. + 2.7e-2*T9**(1.d0/3) - 8.03e-1*T9**(2.d0/3) -
| 1.54e-1*T9 + 5.*T9**(4.d0/3) + 2.44*T9**(5.d0/3))*
| exp(-1.5202e1*T9**(-1.d0/3) - 7.0498e-1*T9**2)*
| T9**(-5.d0/3) + 4.04e7*(5.0673*T9**(-4.d0/3) - 1.41*T9)*
| (1. + 2.7e-2*T9**(1.d0/3) - 8.03e-1*T9**(2.d0/3) -
| 1.54e-1*T9 + 5.*T9**(4.d0/3) + 2.44*T9**(5.d0/3))*
| exp(-1.5202e1*T9**(-1.d0/3) - 7.0498e-1*T9**2)*
| T9**(-2.d0/3))*exp(-5.3706e1*(1/T9))*T9**(3.d0/2)
c(71)= 'N14(P,G)O15'
dfdt(j,71)= 7.1361e3*exp(-3.011*(1/T9))*T9**(-7.d0/2) -
| 3.555e3*exp(-3.011*(1/T9))*T9**(-5.d0/2) +
| 2.7441e5*exp(-1.253e1*(1/T9))*T9**(-2) +
| 1.0372e7*(4.252e-2 - 2.4504*T9**(1.d0/3) -
| 7.0394e-1*T9**(2.d0/3) + 1.6441*T9 + 1.*T9**(4.d0/3))*
| exp(-1.5228e1*T9**(-1.d0/3) - 9.2162e-2*T9**2)*
| T9**(-4.d0/3) + 4.1487e6*
| (-7.874 - 2.126e-1*T9**(1.d0/3) + 6.126*T9**(2.d0/3) +
| 1.1732*T9 - 2.0551*T9**(4.d0/3) - 1.*T9**(5.d0/3))*
| exp(-1.5228e1*T9**(-1.d0/3) - 9.2162e-2*T9**2)*
| T9**(-5.d0/3) + 1.1471e6*
| (7.874 + 2.126e-1*T9**(1.d0/3) - 6.126*T9**(2.d0/3) -
| 1.1732*T9 + 2.0551*T9**(4.d0/3) + 1.*T9**(5.d0/3))*
| (2.7538e1 - 1.*T9**(7.d0/3))*
| exp(-1.5228e1*T9**(-1.d0/3) - 9.2162e-2*T9**2)*T9**(-2)
drdt(j,71)= 2.2863e12*(2.19e4*exp(-1.253e1*(1/T9)) +
| 2.37e3*exp(-3.011*(1/T9))*T9**(-3.d0/2) +
| 4.9e7*(1. + 2.7e-2*T9**(1.d0/3) - 7.78e-1*T9**(2.d0/3) -
| 1.49e-1*T9 + 2.61e-1*T9**(4.d0/3) +
| 1.27e-1*T9**(5.d0/3))*
| exp(-1.5228e1*T9**(-1.d0/3) - 9.2162e-2*T9**2)*
| T9**(-2.d0/3))*exp(-8.4678e1*(1/T9))*(1/Sqrt(T9)) +
| 4.05e10*(2.19e4*exp(-1.253e1*(1/T9)) +
| 2.37e3*exp(-3.011*(1/T9))*T9**(-3.d0/2) +
| 4.9e7*(1. + 2.7e-2*T9**(1.d0/3) - 7.78e-1*T9**(2.d0/3) -
| 1.49e-1*T9 + 2.61e-1*T9**(4.d0/3) +
| 1.27e-1*T9**(5.d0/3))*
| exp(-1.5228e1*T9**(-1.d0/3) - 9.2162e-2*T9**2)*
| T9**(-2.d0/3))*exp(-8.4678e1*(1/T9))*Sqrt(T9) +
| 2.7e10*(7.1361e3*exp(-3.011*(1/T9))*T9**(-7.d0/2) -
| 3.555e3*exp(-3.011*(1/T9))*T9**(-5.d0/2) +
| 2.7441e5*exp(-1.253e1*(1/T9))*T9**(-2) +
| 4.9e7*(-1.49e-1 + 9.e-3*T9**(-2.d0/3) -
| 5.1867e-1*T9**(-1.d0/3) + 3.48e-1*T9**(1.d0/3) +
| 2.1167e-1*T9**(2.d0/3))*
| exp(-1.5228e1*T9**(-1.d0/3) - 9.2162e-2*T9**2)*
| T9**(-2.d0/3) - 3.2667e7*
| (1. + 2.7e-2*T9**(1.d0/3) - 7.78e-1*T9**(2.d0/3) -
| 1.49e-1*T9 + 2.61e-1*T9**(4.d0/3) +
| 1.27e-1*T9**(5.d0/3))*
| exp(-1.5228e1*T9**(-1.d0/3) - 9.2162e-2*T9**2)*
| T9**(-5.d0/3) + 4.9e7*
| (5.076*T9**(-4.d0/3) - 1.8432e-1*T9)*
| (1. + 2.7e-2*T9**(1.d0/3) - 7.78e-1*T9**(2.d0/3) -
| 1.49e-1*T9 + 2.61e-1*T9**(4.d0/3) +
| 1.27e-1*T9**(5.d0/3))*
| exp(-1.5228e1*T9**(-1.d0/3) - 9.2162e-2*T9**2)*
| T9**(-2.d0/3))*exp(-8.4678e1*(1/T9))*T9**(3.d0/2)
c(72)= 'N14(P,N)O14'
dfdt(j,72)= 2.5545e7*(1.8143e2 + 1.1938e2*Sqrt(T9) +
| 6.8762e1*T9 +
| 8.6807e-1*T9**(3.d0/2) + 1.*T9**2)*
| exp(-6.8762e1*(1/T9))*T9**(-2)
drdt(j,72)= 2.0153e8*(3.79e-1 + 3.29e-1*(1/Sqrt(T9)))
c(73)= 'N14(P,A)C11'
dfdt(j,73)= -3.945e16*exp
| (-3.3915e1*(1/T9) -
| 3.1883e1*(T9*(1/
| (1. + 4.78e-2*T9 +
| 7.56e-3*(1. + 4.78e-2*T9)**(-2.d0/3)*
| T9**(5.d0/3))))**(-1.d0/3))*T9**(-5.d0/2)*
| (T9*(1/(1. + 4.78e-2*T9 +
| 7.56e-3*(1. + 4.78e-2*T9)**(-2.d0/3)*
| T9**(5.d0/3))))**(2.5d0/3) +
| 2.1917e16*(-(T9*(4.78e-2 +
| 1.26e-2*(1. + 4.78e-2*T9)**(-2.d0/3)*
| T9**(2.d0/3) -
| 2.4091e-4*(1. + 4.78e-2*T9)**(-5.d0/3)*T9**(5.d0/3)
| )*(1. + 4.78e-2*T9 +
| 7.56e-3*(1. + 4.78e-2*T9)**(-2.d0/3)*
| T9**(5.d0/3))**(-2)) +
| 1/(1. + 4.78e-2*T9 +
| 7.56e-3*(1. + 4.78e-2*T9)**(-2.d0/3)*T9**(5.d0/3))
| )*exp(-3.3915e1*(1/T9) -
| 3.1883e1*(T9*(1/
| (1. + 4.78e-2*T9 +
| 7.56e-3*(1. + 4.78e-2*T9)**(-2.d0/3)*
| T9**(5.d0/3))))**(-1.d0/3))*T9**(-3.d0/2)*
| (T9*(1/(1. + 4.78e-2*T9 +
| 7.56e-3*(1. + 4.78e-2*T9)**(-2.d0/3)*
| T9**(5.d0/3))))**(-.5d0/3) +
| 2.63e16*(3.3915e1*T9**(-2) +
| 1.0628e1*(-(T9*(4.78e-2 +
| 1.26e-2*(1. + 4.78e-2*T9)**(-2.d0/3)*
| T9**(2.d0/3) -
| 2.4091e-4*(1. + 4.78e-2*T9)**(-5.d0/3)*
| T9**(5.d0/3))*
| (1. + 4.78e-2*T9 +
| 7.56e-3*(1. + 4.78e-2*T9)**(-2.d0/3)*
| T9**(5.d0/3))**(-2)) +
| 1/(1. + 4.78e-2*T9 +
| 7.56e-3*(1. + 4.78e-2*T9)**(-2.d0/3)*
| T9**(5.d0/3)))*
| (T9*(1/
| (1. + 4.78e-2*T9 +
| 7.56e-3*(1. + 4.78e-2*T9)**(-2.d0/3)*
| T9**(5.d0/3))))**(-4.d0/3))*
| exp(-3.3915e1*(1/T9) -
| 3.1883e1*(T9*(1/
| (1. + 4.78e-2*T9 +
| 7.56e-3*(1. + 4.78e-2*T9)**(-2.d0/3)*
| T9**(5.d0/3))))**(-1.d0/3))*T9**(-3.d0/2)*
| (T9*(1/(1. + 4.78e-2*T9 +
| 7.56e-3*(1. + 4.78e-2*T9)**(-2.d0/3)*
| T9**(5.d0/3))))**(2.5d0/3)
drdt(j,73)= -1.073e16*exp
| (-3.1883e1*(T9*(1/
| (1. + 4.78e-2*T9 +
| 7.56e-3*(1. + 4.78e-2*T9)**(-2.d0/3)*
| T9**(5.d0/3))))**(-1.d0/3))*T9**(-5.d0/2)*
| (T9*(1/(1. + 4.78e-2*T9 +
| 7.56e-3*(1. + 4.78e-2*T9)**(-2.d0/3)*
| T9**(5.d0/3))))**(2.5d0/3) +
| 7.6026e16*(-(T9*(4.78e-2 +
| 1.26e-2*(1. + 4.78e-2*T9)**(-2.d0/3)*
| T9**(2.d0/3) -
| 2.4091e-4*(1. + 4.78e-2*T9)**(-5.d0/3)*T9**(5.d0/3)
| )*(1. + 4.78e-2*T9 +
| 7.56e-3*(1. + 4.78e-2*T9)**(-2.d0/3)*
| T9**(5.d0/3))**(-2)) +
| 1/(1. + 4.78e-2*T9 +
| 7.56e-3*(1. + 4.78e-2*T9)**(-2.d0/3)*T9**(5.d0/3))
| )*exp(-3.1883e1*(T9*
| (1/
| (1. + 4.78e-2*T9 +
| 7.56e-3*(1. + 4.78e-2*T9)**(-2.d0/3)*
| T9**(5.d0/3))))**(-1.d0/3))*T9**(-3.d0/2)*
| (T9*(1/(1. + 4.78e-2*T9 +
| 7.56e-3*(1. + 4.78e-2*T9)**(-2.d0/3)*
| T9**(5.d0/3))))**(-5.e-1) +
| 5.9613e15*(-(T9*(4.78e-2 +
| 1.26e-2*(1. + 4.78e-2*T9)**(-2.d0/3)*
| T9**(2.d0/3) -
| 2.4091e-4*(1. + 4.78e-2*T9)**(-5.d0/3)*T9**(5.d0/3)
| )*(1. + 4.78e-2*T9 +
| 7.56e-3*(1. + 4.78e-2*T9)**(-2.d0/3)*
| T9**(5.d0/3))**(-2)) +
| 1/(1. + 4.78e-2*T9 +
| 7.56e-3*(1. + 4.78e-2*T9)**(-2.d0/3)*T9**(5.d0/3))
| )*exp(-3.1883e1*(T9*
| (1/
| (1. + 4.78e-2*T9 +
| 7.56e-3*(1. + 4.78e-2*T9)**(-2.d0/3)*
| T9**(5.d0/3))))**(-1.d0/3))*T9**(-3.d0/2)*
| (T9*(1/(1. + 4.78e-2*T9 +
| 7.56e-3*(1. + 4.78e-2*T9)**(-2.d0/3)*
| T9**(5.d0/3))))**(-.5d0/3)
c(74)= 'N14(A,G)F18'
dfdt(j,74)= 1.026e1*exp(-5.054*(1/T9))*T9**(-7.d0/2) +
| 6.6033e-10*exp(-2.798*(1/T9))*T9**(-7.d0/2) +
| 1.4157e5*exp(-1.231e1*(1/T9))*T9**(-8.d0/3) -
| 3.045*exp(-5.054*(1/T9))*T9**(-5.d0/2) -
| 3.54e-10*exp(-2.798*(1/T9))*T9**(-5.d0/2) -
| 7.6667e3*exp(-1.231e1*(1/T9))*T9**(-5.d0/3) +
| 5.2645e9*(5.9113e-3 + 1.4286*T9**(1.d0/3) +
| 1.7291e-1*T9**(2.d0/3) + 3.8818*T9 + 1.*T9**(4.d0/3))*
| exp(-3.6031e1*T9**(-1.d0/3) - 1.2884*T9**2)*T9**(-4.d0/3)
| - 5.1867e9*(1. + 1.2e-2*T9**(1.d0/3) + 1.45*T9**(2.d0/3) +
| 1.17e-1*T9 + 1.97*T9**(4.d0/3) + 4.06e-1*T9**(5.d0/3))*
| exp(-3.6031e1*T9**(-1.d0/3) - 1.2884*T9**2)*T9**(-5.d0/3)
| + 8.1392e9*(2.4631 + 2.9557e-2*T9**(1.d0/3) +
| 3.5714*T9**(2.d0/3) + 2.8818e-1*T9 +
| 4.8522*T9**(4.d0/3) + 1.*T9**(5.d0/3))*
| (4.661 - 1.*T9**(7.d0/3))*
| exp(-3.6031e1*T9**(-1.d0/3) - 1.2884*T9**2)*T9**(-2)
drdt(j,74)= 2.777e12*(2.03*exp(-5.054*(1/T9))
| *T9**(-3.d0/2) +
| 2.36e-10*exp(-2.798*(1/T9))*T9**(-3.d0/2) +
| 1.15e4*exp(-1.231e1*(1/T9))*T9**(-2.d0/3) +
| 7.78e9*(1. + 1.2e-2*T9**(1.d0/3) + 1.45*T9**(2.d0/3) +
| 1.17e-1*T9 + 1.97*T9**(4.d0/3) + 4.06e-1*T9**(5.d0/3))
| *exp(-3.6031e1*T9**(-1.d0/3) - 1.2884*T9**2)*
| T9**(-2.d0/3))*exp(-5.1236e1*(1/T9))*(1/Sqrt(T9)) +
| 8.13e10*(2.03*exp(-5.054*(1/T9))*T9**(-3.d0/2) +
| 2.36e-10*exp(-2.798*(1/T9))*T9**(-3.d0/2) +
| 1.15e4*exp(-1.231e1*(1/T9))*T9**(-2.d0/3) +
| 7.78e9*(1. + 1.2e-2*T9**(1.d0/3) + 1.45*T9**(2.d0/3) +
| 1.17e-1*T9 + 1.97*T9**(4.d0/3) + 4.06e-1*T9**(5.d0/3))
| *exp(-3.6031e1*T9**(-1.d0/3) - 1.2884*T9**2)*
| T9**(-2.d0/3))*exp(-5.1236e1*(1/T9))*Sqrt(T9) +
| 5.42e10*(1.026e1*exp(-5.054*(1/T9))*T9**(-7.d0/2) +
| 6.6033e-10*exp(-2.798*(1/T9))*T9**(-7.d0/2) +
| 1.4157e5*exp(-1.231e1*(1/T9))*T9**(-8.d0/3) -
| 3.045*exp(-5.054*(1/T9))*T9**(-5.d0/2) -
| 3.54e-10*exp(-2.798*(1/T9))*T9**(-5.d0/2) -
| 7.6667e3*exp(-1.231e1*(1/T9))*T9**(-5.d0/3) +
| 7.78e9*(1.17e-1 + 4.e-3*T9**(-2.d0/3) +
| 9.6667e-1*T9**(-1.d0/3) + 2.6267*T9**(1.d0/3) +
| 6.7667e-1*T9**(2.d0/3))*
| exp(-3.6031e1*T9**(-1.d0/3) - 1.2884*T9**2)*
| T9**(-2.d0/3) - 5.1867e9*
| (1. + 1.2e-2*T9**(1.d0/3) + 1.45*T9**(2.d0/3) +
| 1.17e-1*T9 + 1.97*T9**(4.d0/3) + 4.06e-1*T9**(5.d0/3))
| *exp(-3.6031e1*T9**(-1.d0/3) - 1.2884*T9**2)*
| T9**(-5.d0/3) + 7.78e9*
| (1.201e1*T9**(-4.d0/3) - 2.5768*T9)*
| (1. + 1.2e-2*T9**(1.d0/3) + 1.45*T9**(2.d0/3) +
| 1.17e-1*T9 + 1.97*T9**(4.d0/3) + 4.06e-1*T9**(5.d0/3))
| *exp(-3.6031e1*T9**(-1.d0/3) - 1.2884*T9**2)*
| T9**(-2.d0/3))*exp(-5.1236e1*(1/T9))*T9**(3.d0/2)
c(75)= 'N14(A,N)F17'
dfdt(j,75)= 5.24e9*(3.65e-1 - 5.75e-1*(1/Sqrt(T9)))*
| exp(-5.4942e1*(1/T9) - 1.2773e-1*T9**2) +
| 5.24e9*(5.4942e1*T9**(-2) - 2.5547e-1*T9)*
| (1. - 1.15*Sqrt(T9) + 3.65e-1*T9)*
| exp(-5.4942e1*(1/T9) - 1.2773e-1*T9**2) +
| 2.3192e12*exp(-7.0708e1*(1/T9))*T9**(-7.d0/2) -
| 4.92e10*exp(-7.0708e1*(1/T9))*T9**(-5.d0/2)
drdt(j,75)= 7.7552e9*(3.65e-1 - 5.75e-1*(1/Sqrt(T9)))*
| exp(-1.2773e-1*T9**2) +
| 7.6534e11*exp(-1.5766e1*(1/T9))*T9**(-7.d0/2) -
| 7.2816e10*exp(-1.5766e1*(1/T9))*T9**(-5.d0/2) -
| 1.9812e9*(1. - 1.15*Sqrt(T9) + 3.65e-1*T9)*T9*
| exp(-1.2773e-1*T9**2)
|
c(76)= 'N15(P,G)O16'
dfdt(j,76)= 4.1982e7*exp(-1.1048e1*(1/T9))*T9**(-7.d0/2) +
| 6.9508e4*exp(-4.665*(1/T9))*T9**(-7.d0/2) +
| 3.6941e4*exp(-3.328*(1/T9))*T9**(-7.d0/2) -
| 5.7e6*exp(-1.1048e1*(1/T9))*T9**(-5.d0/2) -
| 2.235e4*exp(-4.665*(1/T9))*T9**(-5.d0/2) -
| 1.665e4*exp(-3.328*(1/T9))*T9**(-5.d0/2) +
| 5.4116e9*(1.6265e-3 + 2.6386e-2*T9**(1.d0/3) +
| 7.5904e-3*T9**(2.d0/3) + 1.6458*T9 + 1.*T9**(4.d0/3))*
| exp(-1.5251e1*T9**(-1.d0/3) - 4.9383*T9**2)*T9**(-4.d0/3)
| - 6.52e8*(1. + 2.7e-2*T9**(1.d0/3) +
| 2.19e-1*T9**(2.d0/3) + 4.2e-2*T9 + 6.83*T9**(4.d0/3) +
| 3.32*T9**(5.d0/3))*exp
| (-1.5251e1*T9**(-1.d0/3) - 4.9383*T9**2)*T9**(-5.d0/3) +
| 3.2069e10*(3.012e-1 + 8.1325e-3*T9**(1.d0/3) +
| 6.5964e-2*T9**(2.d0/3) + 1.2651e-2*T9 +
| 2.0572*T9**(4.d0/3) + 1.*T9**(5.d0/3))*
| (5.1472e-1 - 1.*T9**(7.d0/3))*
| exp(-1.5251e1*T9**(-1.d0/3) - 4.9383*T9**2)*T9**(-2)
drdt(j,76)= 5.0946e12*(3.8e6*exp(-1.1048e1*(1/T9))*
| T9**(-3.d0/2) +
| 1.49e4*exp(-4.665*(1/T9))*T9**(-3.d0/2) +
| 1.11e4*exp(-3.328*(1/T9))*T9**(-3.d0/2) +
| 9.78e8*(1. + 2.7e-2*T9**(1.d0/3) +
| 2.19e-1*T9**(2.d0/3) + 4.2e-2*T9 +
| 6.83*T9**(4.d0/3) + 3.32*T9**(5.d0/3))*
| exp(-1.5251e1*T9**(-1.d0/3) - 4.9383*T9**2)*
| T9**(-2.d0/3))*exp(-1.4073e2*(1/T9))*(1/Sqrt(T9)) +
| 5.43e10*(3.8e6*exp(-1.1048e1*(1/T9))*T9**(-3.d0/2) +
| 1.49e4*exp(-4.665*(1/T9))*T9**(-3.d0/2) +
| 1.11e4*exp(-3.328*(1/T9))*T9**(-3.d0/2) +
| 9.78e8*(1. + 2.7e-2*T9**(1.d0/3) +
| 2.19e-1*T9**(2.d0/3) + 4.2e-2*T9 +
| 6.83*T9**(4.d0/3) + 3.32*T9**(5.d0/3))*
| exp(-1.5251e1*T9**(-1.d0/3) - 4.9383*T9**2)*
| T9**(-2.d0/3))*exp(-1.4073e2*(1/T9))*Sqrt(T9) +
| 3.62e10*(4.1982e7*exp(-1.1048e1*(1/T9))*T9**(-7.d0/2) +
| 6.9508e4*exp(-4.665*(1/T9))*T9**(-7.d0/2) +
| 3.6941e4*exp(-3.328*(1/T9))*T9**(-7.d0/2) -
| 5.7e6*exp(-1.1048e1*(1/T9))*T9**(-5.d0/2) -
| 2.235e4*exp(-4.665*(1/T9))*T9**(-5.d0/2) -
| 1.665e4*exp(-3.328*(1/T9))*T9**(-5.d0/2) +
| 9.78e8*(4.2e-2 + 9.e-3*T9**(-2.d0/3) +
| 1.46e-1*T9**(-1.d0/3) + 9.1067*T9**(1.d0/3) +
| 5.5333*T9**(2.d0/3))*
| exp(-1.5251e1*T9**(-1.d0/3) - 4.9383*T9**2)*
| T9**(-2.d0/3) - 6.52e8*
| (1. + 2.7e-2*T9**(1.d0/3) + 2.19e-1*T9**(2.d0/3) +
| 4.2e-2*T9 + 6.83*T9**(4.d0/3) + 3.32*T9**(5.d0/3))*
| exp(-1.5251e1*T9**(-1.d0/3) - 4.9383*T9**2)*
| T9**(-5.d0/3) + 9.78e8*
| (5.0837*T9**(-4.d0/3) - 9.8765*T9)*
| (1. + 2.7e-2*T9**(1.d0/3) + 2.19e-1*T9**(2.d0/3) +
| 4.2e-2*T9 + 6.83*T9**(4.d0/3) + 3.32*T9**(5.d0/3))*
| exp(-1.5251e1*T9**(-1.d0/3) - 4.9383*T9**2)*
| T9**(-2.d0/3))*exp(-1.4073e2*(1/T9))*T9**(3.d0/2)
c(77)= 'N15(P,N)O15'
dfdt(j,77)=