diff --git a/net/test/test_output b/net/test/test_output index 3f16af2df..8a8fa1607 100644 --- a/net/test/test_output +++ b/net/test/test_output @@ -9,7 +9,7 @@ eps_nuc 0.839654151E+13 - d_lneps_dlnT 9.328757835 + d_lneps_dlnT 9.328757799 d_lneps_dlnRho 1.009292766 @@ -34,7 +34,7 @@ eps_nuc 0.302124143E+11 - d_lneps_dlnT 11.356884345 + d_lneps_dlnT 11.356884263 d_lneps_dlnRho 1.015667220 @@ -59,7 +59,7 @@ eps_nuc 0.835361769E+13 - d_lneps_dlnT 9.326560302 + d_lneps_dlnT 9.326560266 d_lneps_dlnRho 1.009292884 @@ -84,7 +84,7 @@ eps_nuc 0.536017372E+09 - d_lneps_dlnT 11.559205087 + d_lneps_dlnT 11.559204995 d_lneps_dlnRho 1.013842798 @@ -109,7 +109,7 @@ eps_nuc 0.150660462E+25 - d_lneps_dlnT -0.253424018 + d_lneps_dlnT -0.253423995 d_lneps_dlnRho 1.008418067 @@ -151,7 +151,7 @@ eps_nuc 0.150660462E+25 - d_lneps_dlnT -0.253424018 + d_lneps_dlnT -0.253423995 d_lneps_dlnRho 1.008418067 @@ -194,22 +194,22 @@ test_one_zone_burn_small_net number of species 21 large final abundances 1.0000000000000000D-02 - c12 1 9.8714342254068388D-01 - he4 2 1.2802186836379657D-02 + c12 1 9.8714342283429213D-01 + he4 2 1.2802186549072116D-02 - xsum 1.0000000000006330D+00 + xsum 1.0000000000005884D+00 test_one_zone_burn_const_P number of species 21 large final abundances 1.0000000000000000D-02 - o16 1 7.1798810044269901D-01 - si28 2 1.6142410416068431D-01 - s32 3 6.0439378324441878D-02 - ar36 4 2.9692711889885499D-02 - ca40 5 1.6560528697491497D-02 + o16 1 7.1798810044342132D-01 + si28 2 1.6142410415591638D-01 + s32 3 6.0439378327693444D-02 + ar36 4 2.9692711893865516D-02 + ca40 5 1.6560528702223146D-02 - xsum 1.0000000000000000D+00 + xsum 9.9999999999999978D-01 diff --git a/rates/private/rates_support.f90 b/rates/private/rates_support.f90 index 9fc001279..46365c7fb 100644 --- a/rates/private/rates_support.f90 +++ b/rates/private/rates_support.f90 @@ -167,17 +167,14 @@ subroutine get_rates_from_table(r1, r2) do i = r1,r2 - rate_raw(i) = & - (rattab_f(1,k,i) + dt*(rattab_f(2,k,i) + & - dt*(rattab_f(3,k,i) + dt*rattab_f(4,k,i))) & - ) * dtab(i) + rate_raw(i) = exp10(rattab_f(1,k,i) + dt*(rattab_f(2,k,i) + & + dt*(rattab_f(3,k,i) + dt*rattab_f(4,k,i))) ) * dtab(i) - rate_raw_dRho(i) = rate_raw(i) * ddtab(i) / dtab(i) + ! Derivative with respect to density (dRho) + rate_raw_dRho(i) = rate_raw(i) * (ddtab(i) / dtab(i)) - rate_raw_dT(i) = & - (rattab_f(2,k,i) + 2*dt*(rattab_f(3,k,i) + & - 1.5d0*dt*rattab_f(4,k,i)) & - ) * dtab(i) / (btemp * ln10) + ! Derivative with respect to temperature (dT) + rate_raw_dT(i) = rate_raw(i) * (rattab_f(2,k,i) + 2*dt*(rattab_f(3,k,i) + 1.5d0*dt*rattab_f(4,k,i))) / (btemp) end do @@ -201,7 +198,7 @@ subroutine do_make_rate_tables( & integer, intent(out) :: ierr integer :: i, j, operr, num_to_add_to_cache,thread_num - real(dp) :: logT, btemp + real(dp) :: logT, btemp, rate_logR real(dp), pointer :: work1(:)=>null(), f1(:)=>null(), rattab_f(:,:,:)=>null() integer, pointer :: reaction_id(:) =>null() real(dp), allocatable, target :: work(:,:) @@ -211,7 +208,8 @@ subroutine do_make_rate_tables( & include 'formats' ierr = 0 - + rate_logR = 0d0 + rattab_f(1:4,1:nrattab,1:num_reactions) => & rattab_f1(1:4*nrattab*num_reactions) @@ -282,6 +280,20 @@ subroutine do_make_rate_tables( & write(*, '(a,i4,2x,a)') 'missing raw rate for ', & j, trim(reaction_Name(reaction_id(j))) a_okay = .false. + else + ! **New Addition: Store the logarithm of the rate** + if (rattab(j, i) > 0.0_dp) then ! Only take log of positive values + rate_logR = log10(rattab(j, i)) + else if (rattab(j, i) == 0.0_dp) then + rate_logR = -323d0 ! set to a tiny value that won't underflow when exponentiated + else + ! warning for rates set to -1 + ! like rni56ec_to_co56,rco56ec_to_fe56 + !write(*, '(a,i4,2x,a)') 'Warning: non-positive value for rate in ', & + ! j, trim(reaction_Name(reaction_id(j))) + call set_nan(rate_logR) + end if + rattab(j, i) = rate_logR ! Store value in rattab end if end do if (.not. a_okay) all_okay = .false. diff --git a/rates/public/rates_def.f90 b/rates/public/rates_def.f90 index 32f43a535..7f5e2c82a 100644 --- a/rates/public/rates_def.f90 +++ b/rates/public/rates_def.f90 @@ -1150,7 +1150,7 @@ subroutine do_start_rates_def_init(ierr) ierr = 0 call set_rattab_range(5.30102999566398d0, 10.301029995664d0) - reaclib_min_T9 = 1d-2 + reaclib_min_T9 = 1d-4 ! need <= 2d-3 for pre-ms li7 burning ! pre-ms deuterium burning needs much lower (4d-4) ! but that seems to cause problems during advanced burning. diff --git a/rates/test/test_output b/rates/test/test_output index 7aaf02fb6..770a17cc3 100644 --- a/rates/test/test_output +++ b/rates/test/test_output @@ -381,8 +381,8 @@ r_he3_ag_be7 1.0717812086522761D-47 r34_pp2 1.0717812086522761D-47 r34_pp3 1.0717812086522761D-47 - r_al27_pa_mg24 0.0000000000000000D+00 - r_ar36_ag_ca40 0.0000000000000000D+00 + r_al27_pa_mg24 1.1726448908365033-100 + r_ar36_ag_ca40 3.4668300436425714-284 r_ar36_ga_s32 0.0000000000000000D+00 r_b8_gp_be7 0.0000000000000000D+00 r_be7_pg_b8 6.9283961564215479D-38 @@ -390,63 +390,63 @@ r_c12_ap_n15 0.0000000000000000D+00 r_c12_pg_n13 6.9450591676016374D-51 r_c12_to_he4_he4_he4 0.0000000000000000D+00 - r_c13_an_o16 0.0000000000000000D+00 + r_c13_an_o16 1.8900211955401780-134 r_c13_pg_n14 2.4960422321922773D-50 - r_ca40_ag_ti44 0.0000000000000000D+00 + r_ca40_ag_ti44 3.8947751206152086-307 r_ca40_ga_ar36 0.0000000000000000D+00 r_cr48_ag_fe52 0.0000000000000000D+00 r_cr48_ga_ti44 0.0000000000000000D+00 - r_f17_ap_ne20 0.0000000000000000D+00 + r_f17_ap_ne20 2.7619647416652612-169 r_f17_gp_o16 0.0000000000000000D+00 r_f17_pa_o14 0.0000000000000000D+00 r_f18_gp_o17 0.0000000000000000D+00 - r_f18_pa_o15 0.0000000000000000D+00 - r_f17_pg_ne18 0.0000000000000000D+00 - r_f18_pg_ne19 0.0000000000000000D+00 - r_f19_ap_ne22 0.0000000000000000D+00 + r_f18_pa_o15 3.3473994880853350D-68 + r_f17_pg_ne18 4.6046359502054981D-72 + r_f18_pg_ne19 9.7254531143126770D-71 + r_f19_ap_ne22 2.2909817203896192-170 r_f19_gp_o18 0.0000000000000000D+00 r_f19_pa_o16 5.5147721045188877D-66 - r_f19_pg_ne20 0.0000000000000000D+00 + r_f19_pg_ne20 1.5919282201534221D-69 r_fe52_ag_ni56 0.0000000000000000D+00 r_fe52_ga_cr48 0.0000000000000000D+00 - r_h2_be7_to_h1_he4_he4 0.0000000000000000D+00 - r_h2_h2_to_he4 0.0000000000000000D+00 - r_h2_he3_to_h1_he4 0.0000000000000000D+00 - r_he3_be7_to_h1_h1_he4_he4 0.0000000000000000D+00 + r_h2_be7_to_h1_he4_he4 1.1357306727406292D-40 + r_h2_h2_to_he4 1.3471260714309430D-15 + r_h2_he3_to_h1_he4 3.5491543661924754D-19 + r_he3_be7_to_h1_h1_he4_he4 1.3811200901115312D-79 r_he3_he3_to_h1_h1_he4 2.6854658991383834D-41 r_h1_h1_he4_to_he3_he3 0.0000000000000000D+00 - r_he4_he4_he4_to_c12 0.0000000000000000D+00 + r_he4_he4_he4_to_c12 7.1309748215131773-157 r_li7_pa_he4 1.1481972998876799D-26 r_he4_ap_li7 0.0000000000000000D+00 r_mg24_ag_si28 0.0000000000000000D+00 r_mg24_ap_al27 0.0000000000000000D+00 r_mg24_ga_ne20 0.0000000000000000D+00 - r_n13_ap_o16 0.0000000000000000D+00 + r_n13_ap_o16 8.8438391452107862-137 r_n13_gp_c12 0.0000000000000000D+00 r_n13_pg_o14 3.7274186108029641D-57 - r_n14_ag_f18 0.0000000000000000D+00 + r_n14_ag_f18 8.2898440628435526-147 r_n14_ap_o17 0.0000000000000000D+00 r_n14_gp_c13 0.0000000000000000D+00 r_n14_pg_o15 3.4330959194227231D-57 - r_n15_ag_f19 0.0000000000000000D+00 + r_n15_ag_f19 4.7001337826755437-145 r_n15_ap_o18 0.0000000000000000D+00 - r_n15_pa_c12 0.0000000000000000D+00 - r_n15_pg_o16 0.0000000000000000D+00 - r_na23_pa_ne20 0.0000000000000000D+00 + r_n15_pa_c12 5.7564809766760218D-53 + r_n15_pg_o16 3.3386726887925942D-56 + r_na23_pa_ne20 7.4956336357959249D-77 r_ne18_gp_f17 0.0000000000000000D+00 r_ne19_ga_o15 0.0000000000000000D+00 r_ne19_gp_f18 0.0000000000000000D+00 - r_ne20_ag_mg24 0.0000000000000000D+00 + r_ne20_ag_mg24 1.0392014393758480-188 r_ne20_ap_na23 0.0000000000000000D+00 r_ne20_ga_o16 0.0000000000000000D+00 r_ne20_gp_f19 0.0000000000000000D+00 r_ne22_ag_mg26 0.0000000000000000D+00 - r_ne22_pg_na23 0.0000000000000000D+00 + r_ne22_pg_na23 4.2586804983614625D-74 r_ni56_ga_fe52 0.0000000000000000D+00 - r_o14_ag_ne18 0.0000000000000000D+00 - r_o14_ap_f17 0.0000000000000000D+00 + r_o14_ag_ne18 2.8524007639314185-158 + r_o14_ap_f17 8.5404735293280779-153 r_o14_gp_n13 0.0000000000000000D+00 - r_o15_ag_ne19 0.0000000000000000D+00 + r_o15_ag_ne19 3.4129338044899019-159 r_o15_ap_f18 0.0000000000000000D+00 r_o15_gp_n14 0.0000000000000000D+00 r_o16_ag_ne20 5.6578133965133730-161 @@ -454,19 +454,19 @@ r_o16_ga_c12 0.0000000000000000D+00 r_o16_gp_n15 0.0000000000000000D+00 r_o16_pg_f17 6.5714436667798447D-63 - r_o17_pa_n14 0.0000000000000000D+00 - r_o17_pg_f18 0.0000000000000000D+00 + r_o17_pa_n14 1.1841359437155215D-63 + r_o17_pg_f18 7.0786869127070877D-63 r_o18_ag_ne22 0.0000000000000000D+00 r_o18_pa_n15 1.2074711916489890D-59 - r_o18_pg_f19 0.0000000000000000D+00 - r_s32_ag_ar36 0.0000000000000000D+00 + r_o18_pg_f19 1.0370895228440540D-62 + r_s32_ag_ar36 3.8006708631673725-261 r_s32_ga_si28 0.0000000000000000D+00 - r_si28_ag_s32 0.0000000000000000D+00 + r_si28_ag_s32 4.4055940886229377-236 r_si28_ga_mg24 0.0000000000000000D+00 r_ti44_ag_cr48 0.0000000000000000D+00 r_ti44_ga_ca40 0.0000000000000000D+00 - ral27pa_aux 0.0000000000000000D+00 - ral27pg_aux 0.0000000000000000D+00 + ral27pa_aux 1.1726448908365033-100 + ral27pg_aux 2.0944916149633688D-90 rar36ap_aux 0.0000000000000000D+00 rar36ap_to_ca40 0.0000000000000000D+00 rar36gp_aux 0.0000000000000000D+00 @@ -481,16 +481,16 @@ rca40ap_to_ti44 0.0000000000000000D+00 rca40gp_aux 0.0000000000000000D+00 rca40gp_to_ar36 0.0000000000000000D+00 - rcl35pa_aux 0.0000000000000000D+00 - rcl35pg_aux 0.0000000000000000D+00 + rcl35pa_aux 4.4128898333781778-119 + rcl35pg_aux 4.5012891879051611-105 rco55gprot_aux 0.0000000000000000D+00 - rco55pg_aux 0.0000000000000000D+00 - rco55protg_aux 0.0000000000000000D+00 + rco55pg_aux 5.3902477026630388-148 + rco55protg_aux 5.3902477026630388-148 rcr48ap_aux 0.0000000000000000D+00 rcr48ap_to_fe52 0.0000000000000000D+00 rcr48gp_aux 0.0000000000000000D+00 rcr48gp_to_ti44 0.0000000000000000D+00 - rf19pg_aux 0.0000000000000000D+00 + rf19pg_aux 1.5919282201534221D-69 rfe52ap_aux 0.0000000000000000D+00 rfe52ap_to_ni56 0.0000000000000000D+00 rfe52aprot_aux 0.0000000000000000D+00 @@ -499,21 +499,21 @@ rfe52gp_aux 0.0000000000000000D+00 rfe52gp_to_cr48 0.0000000000000000D+00 rfe52neut_to_fe54 -1.0000000000000000D+00 - rfe52ng_aux 0.0000000000000000D+00 + rfe52ng_aux 7.6808904520025365D+06 rfe53gn_aux 0.0000000000000000D+00 - rfe53ng_aux 0.0000000000000000D+00 + rfe53ng_aux 1.2526517097161818D+07 rfe54a_to_ni56 0.0000000000000000D+00 rfe54an_aux 0.0000000000000000D+00 rfe54an_to_ni56 0.0000000000000000D+00 rfe54aprot_to_fe56 0.0000000000000000D+00 rfe54g_to_fe52 0.0000000000000000D+00 - rfe54ng_aux 0.0000000000000000D+00 + rfe54ng_aux 1.6116706715887450D+18 rfe54ng_to_fe56 -1.0000000000000000D+00 rfe54prot_to_fe52 -1.0000000000000000D+00 rfe54prot_to_ni56 -1.0000000000000000D+00 - rfe54protg_aux 0.0000000000000000D+00 + rfe54protg_aux 2.6041048156863966-144 rfe55gn_aux 0.0000000000000000D+00 - rfe55ng_aux 0.0000000000000000D+00 + rfe55ng_aux 7.8340729757598912D-18 rfe56ec_fake_to_mn56 -1.0000000000000000D+00 rfe56ec_fake_to_mn57 -1.0000000000000000D+00 rfe56ec_fake_to_cr56 -1.0000000000000000D+00 @@ -538,28 +538,28 @@ rhe4_breakup -1.0000000000000000D+00 rhe4_rebuild -1.0000000000000000D+00 rhe4g_neut_aux 0.0000000000000000D+00 - rk39pa_aux 0.0000000000000000D+00 - rk39pg_aux 0.0000000000000000D+00 + rk39pa_aux 4.0639163851271452-120 + rk39pg_aux 4.6186923666554649-113 rmg24ap_aux 0.0000000000000000D+00 rmg24ap_to_si28 0.0000000000000000D+00 rmg24gp_aux 0.0000000000000000D+00 rmg24gp_to_ne20 0.0000000000000000D+00 - rmn51pg_aux 0.0000000000000000D+00 + rmn51pg_aux 4.4198070517200279-140 rn14_to_c12 3.4330959194227231D-57 rn14_to_n15 3.4330959194227231D-57 rn14_to_o16 3.4330959194227231D-57 - rn14ag_lite 0.0000000000000000D+00 + rn14ag_lite 8.2898440628435526-147 rn14gc12 0.0000000000000000D+00 rn14pg_aux 3.4330959194227231D-57 - rn15pa_aux 0.0000000000000000D+00 - rn15pg_aux 0.0000000000000000D+00 - rna23pa_aux 0.0000000000000000D+00 - rna23pg_aux 0.0000000000000000D+00 - rne18ag_to_mg24 0.0000000000000000D+00 - rne18ap_to_mg22 0.0000000000000000D+00 - rne18ap_to_mg24 0.0000000000000000D+00 - rne19pg_to_mg22 0.0000000000000000D+00 - rne19pg_to_mg24 0.0000000000000000D+00 + rn15pa_aux 5.7564809766760218D-53 + rn15pg_aux 3.3386726887925942D-56 + rna23pa_aux 7.4956336357959249D-77 + rna23pg_aux 4.2510781286576797D-80 + rne18ag_to_mg24 2.7408283875886730-186 + rne18ap_to_mg22 1.9720554076528380-252 + rne18ap_to_mg24 1.9720554076528380-252 + rne19pg_to_mg22 3.6009620782629716D-76 + rne19pg_to_mg24 3.6009620782629716D-76 rne20ap_aux 0.0000000000000000D+00 rne20ap_to_mg24 0.0000000000000000D+00 rne20gp_aux 0.0000000000000000D+00 @@ -577,16 +577,16 @@ rni56gprot_to_fe52 -1.0000000000000000D+00 rni56gprot_to_fe54 -1.0000000000000000D+00 rni56ng_to_fe54 -1.0000000000000000D+00 - rni57na_aux 0.0000000000000000D+00 + rni57na_aux 2.1505573077156749D+06 ro16_to_n14 6.5714436667798447D-63 ro16_to_o17 6.5714436667798447D-63 ro16ap_aux 0.0000000000000000D+00 ro16ap_to_ne20 0.0000000000000000D+00 ro16gp_aux 0.0000000000000000D+00 ro16gp_to_c12 0.0000000000000000D+00 - ro17_to_o18 0.0000000000000000D+00 - rp31pa_aux 0.0000000000000000D+00 - rp31pg_aux 0.0000000000000000D+00 + ro17_to_o18 7.0786869127070877D-63 + rp31pa_aux 2.4985373343461576-114 + rp31pg_aux 4.4163367090472991-100 rpep_to_he3 9.2307370018631361D-32 rpp_to_he3 8.4882903687436408D-28 rprot_neutg_aux 4.6167849112068943D+04 @@ -595,8 +595,8 @@ rs32ap_to_ar36 0.0000000000000000D+00 rs32gp_aux 0.0000000000000000D+00 rs32gp_to_si28 0.0000000000000000D+00 - rsc43pa_aux 0.0000000000000000D+00 - rsc43pg_aux 0.0000000000000000D+00 + rsc43pa_aux 2.5781102957576896-123 + rsc43pg_aux 1.5690365347252879-122 rsi28ap_aux 0.0000000000000000D+00 rsi28ap_to_s32 0.0000000000000000D+00 rsi28gp_aux 0.0000000000000000D+00 @@ -606,54 +606,54 @@ rti44gp_aux 0.0000000000000000D+00 rti44gp_to_ca40 0.0000000000000000D+00 rv47pa_aux 0.0000000000000000D+00 - rv47pg_aux 0.0000000000000000D+00 + rv47pg_aux 4.1763717093215707D-57 r_h1_h1_wk_h2 8.4882903687436408D-28 r_h1_h1_ec_h2 9.2307370018631361D-32 - rn14ag_to_ne22 0.0000000000000000D+00 + rn14ag_to_ne22 8.2898440628435526-147 rf19pa_aux 5.5147721045188877D-66 r_b8_wk_he4_he4 9.0019114358434449D-01 rmn51pa_aux 0.0000000000000000D+00 rfe54gn_aux 0.0000000000000000D+00 rco55pa_aux 0.0000000000000000D+00 rco55prota_aux 0.0000000000000000D+00 - rn14ag_to_o18 0.0000000000000000D+00 + rn14ag_to_o18 8.2898440628435526-147 r_h1_he3_wk_he4 1.8770264056112323D-37 r_be7_wk_li7 1.4566132867780899D-07 - r_al27_pg_si28 0.0000000000000000D+00 + r_al27_pg_si28 2.0944916149633688D-90 r_si28_gp_al27 0.0000000000000000D+00 r_si28_ap_p31 0.0000000000000000D+00 - r_p31_pa_si28 0.0000000000000000D+00 - r_p31_pg_s32 0.0000000000000000D+00 + r_p31_pa_si28 2.4985373343461576-114 + r_p31_pg_s32 4.4163367090472991-100 r_s32_gp_p31 0.0000000000000000D+00 r_s32_ap_cl35 0.0000000000000000D+00 - r_cl35_pa_s32 0.0000000000000000D+00 - r_cl35_pg_ar36 0.0000000000000000D+00 + r_cl35_pa_s32 4.4128898333781778-119 + r_cl35_pg_ar36 4.5012891879051611-105 r_ar36_gp_cl35 0.0000000000000000D+00 r_ar36_ap_k39 0.0000000000000000D+00 - r_k39_pa_ar36 0.0000000000000000D+00 - r_k39_pg_ca40 0.0000000000000000D+00 + r_k39_pa_ar36 4.0639163851271452-120 + r_k39_pg_ca40 4.6186923666554649-113 r_ca40_gp_k39 0.0000000000000000D+00 r_ca40_ap_sc43 0.0000000000000000D+00 - r_sc43_pa_ca40 0.0000000000000000D+00 - r_sc43_pg_ti44 0.0000000000000000D+00 + r_sc43_pa_ca40 2.5781102957576896-123 + r_sc43_pg_ti44 1.5690365347252879-122 r_ti44_gp_sc43 0.0000000000000000D+00 r_ti44_ap_v47 0.0000000000000000D+00 r_v47_pa_ti44 0.0000000000000000D+00 - r_v47_pg_cr48 0.0000000000000000D+00 + r_v47_pg_cr48 4.1763717093215707D-57 r_cr48_gp_v47 0.0000000000000000D+00 r_cr48_ap_mn51 0.0000000000000000D+00 r_mn51_pa_cr48 0.0000000000000000D+00 - r_mn51_pg_fe52 0.0000000000000000D+00 + r_mn51_pg_fe52 4.4198070517200279-140 r_fe52_gp_mn51 0.0000000000000000D+00 r_fe52_ap_co55 0.0000000000000000D+00 r_co55_pa_fe52 0.0000000000000000D+00 - r_co55_pg_ni56 0.0000000000000000D+00 + r_co55_pg_ni56 5.3902477026630388-148 r_ni56_gp_co55 0.0000000000000000D+00 - r_fe52_ng_fe53 0.0000000000000000D+00 + r_fe52_ng_fe53 7.6808904520025365D+06 r_fe53_gn_fe52 0.0000000000000000D+00 - r_fe53_ng_fe54 0.0000000000000000D+00 + r_fe53_ng_fe54 1.2526517097161818D+07 r_fe54_gn_fe53 0.0000000000000000D+00 - r_fe54_pg_co55 0.0000000000000000D+00 + r_fe54_pg_co55 2.6041048156863966-144 r_co55_gp_fe54 0.0000000000000000D+00 r_he3_ng_he4 1.2611100000000000D+01 r_he4_gn_he3 0.0000000000000000D+00 @@ -661,13 +661,13 @@ r_h2_gn_h1 0.0000000000000000D+00 r_h2_pg_he3 1.2697222985127328D-11 r_he3_gp_h2 0.0000000000000000D+00 - r_fe54_ng_fe55 0.0000000000000000D+00 + r_fe54_ng_fe55 1.6116706715887450D+18 r_fe55_gn_fe54 0.0000000000000000D+00 - r_fe55_ng_fe56 0.0000000000000000D+00 + r_fe55_ng_fe56 7.8340729757598912D-18 r_fe56_gn_fe55 0.0000000000000000D+00 r_fe54_ap_co57 0.0000000000000000D+00 - r_co57_pa_fe54 0.0000000000000000D+00 - r_fe56_pg_co57 0.0000000000000000D+00 + r_co57_pa_fe54 6.2375165937633259-164 + r_fe56_pg_co57 2.0634641787669126-144 r_co57_gp_fe56 0.0000000000000000D+00 r_c12_c12_to_h1_na23 0.0000000000000000D+00 r_he4_ne20_to_c12_c12 0.0000000000000000D+00 @@ -681,8 +681,8 @@ r1212 1.8116989486805094-221 r1216 0.0000000000000000D+00 - r1216_to_mg24 0.0000000000000000D+00 - r1216_to_si28 0.0000000000000000D+00 + r1216_to_mg24 7.5518111961788331-261 + r1216_to_si28 5.7975894636903584-270 r1616 0.0000000000000000D+00 r1616a 0.0000000000000000D+00 r1616g 0.0000000000000000D+00 @@ -692,72 +692,72 @@ r_he3_ag_be7 2.8231403139027404D-30 r34_pp2 2.8231403139027404D-30 r34_pp3 2.8231403139027404D-30 - r_al27_pa_mg24 0.0000000000000000D+00 - r_ar36_ag_ca40 0.0000000000000000D+00 + r_al27_pa_mg24 1.9918779499772236D-64 + r_ar36_ag_ca40 3.6638259180956996-186 r_ar36_ga_s32 0.0000000000000000D+00 r_b8_gp_be7 1.4668059007716089-236 r_be7_pg_b8 5.1757312531762798D-24 - r_c12_ag_o16 0.0000000000000000D+00 + r_c12_ag_o16 9.7297725055896462-173 r_c12_ap_n15 0.0000000000000000D+00 r_c12_pg_n13 2.9425431722017017D-32 r_c12_to_he4_he4_he4 0.0000000000000000D+00 - r_c13_an_o16 0.0000000000000000D+00 + r_c13_an_o16 2.5450969261817348D-81 r_c13_pg_n14 1.1415893957710228D-31 - r_ca40_ag_ti44 0.0000000000000000D+00 + r_ca40_ag_ti44 1.3118106505932359-201 r_ca40_ga_ar36 0.0000000000000000D+00 - r_cr48_ag_fe52 0.0000000000000000D+00 + r_cr48_ag_fe52 4.8045041641220315-228 r_cr48_ga_ti44 0.0000000000000000D+00 - r_f17_ap_ne20 0.0000000000000000D+00 + r_f17_ap_ne20 9.2024141691416590-110 r_f17_gp_o16 0.0000000000000000D+00 r_f17_pa_o14 0.0000000000000000D+00 r_f18_gp_o17 0.0000000000000000D+00 - r_f18_pa_o15 0.0000000000000000D+00 - r_f17_pg_ne18 0.0000000000000000D+00 - r_f18_pg_ne19 0.0000000000000000D+00 - r_f19_ap_ne22 0.0000000000000000D+00 + r_f18_pa_o15 1.9715159138530976D-40 + r_f17_pg_ne18 3.6367066583240212D-43 + r_f18_pg_ne19 3.4388148063104511D-43 + r_f19_ap_ne22 2.6116557728083981-110 r_f19_gp_o18 0.0000000000000000D+00 r_f19_pa_o16 3.0412796782211088D-41 - r_f19_pg_ne20 0.0000000000000000D+00 - r_fe52_ag_ni56 0.0000000000000000D+00 + r_f19_pg_ne20 8.3225081090292684D-45 + r_fe52_ag_ni56 5.3901176951702360-242 r_fe52_ga_cr48 0.0000000000000000D+00 - r_h2_be7_to_h1_he4_he4 0.0000000000000000D+00 - r_h2_h2_to_he4 0.0000000000000000D+00 - r_h2_he3_to_h1_he4 0.0000000000000000D+00 - r_he3_be7_to_h1_h1_he4_he4 0.0000000000000000D+00 + r_h2_be7_to_h1_he4_he4 8.3930806127232442D-24 + r_h2_h2_to_he4 4.9326118728026079D-10 + r_h2_he3_to_h1_he4 1.4721607322617134D-09 + r_he3_be7_to_h1_h1_he4_he4 9.3597424900485328D-50 r_he3_he3_to_h1_h1_he4 1.2261267385307417D-24 r_h1_h1_he4_to_he3_he3 0.0000000000000000D+00 - r_he4_he4_he4_to_c12 0.0000000000000000D+00 + r_he4_he4_he4_to_c12 9.3309154251394056-106 r_li7_pa_he4 2.8536849590579373D-15 r_he4_ap_li7 0.0000000000000000D+00 r_mg24_ag_si28 0.0000000000000000D+00 r_mg24_ap_al27 0.0000000000000000D+00 r_mg24_ga_ne20 0.0000000000000000D+00 - r_n13_ap_o16 0.0000000000000000D+00 + r_n13_ap_o16 1.5618583622449647D-87 r_n13_gp_c12 0.0000000000000000D+00 r_n13_pg_o14 1.9381819262283661D-36 - r_n14_ag_f18 0.0000000000000000D+00 + r_n14_ag_f18 5.7595359837054342D-97 r_n14_ap_o17 0.0000000000000000D+00 r_n14_gp_c13 0.0000000000000000D+00 r_n14_pg_o15 1.9150164715337093D-36 - r_n15_ag_f19 0.0000000000000000D+00 + r_n15_ag_f19 3.0839883777838008D-95 r_n15_ap_o18 0.0000000000000000D+00 - r_n15_pa_c12 0.0000000000000000D+00 - r_n15_pg_o16 0.0000000000000000D+00 - r_na23_pa_ne20 0.0000000000000000D+00 + r_n15_pa_c12 3.7461382926599108D-32 + r_n15_pg_o16 1.9741989978881749D-35 + r_na23_pa_ne20 2.8956408617284671D-49 r_ne18_gp_f17 0.0000000000000000D+00 r_ne19_ga_o15 0.0000000000000000D+00 r_ne19_gp_f18 0.0000000000000000D+00 - r_ne20_ag_mg24 0.0000000000000000D+00 + r_ne20_ag_mg24 6.6602783188653664-125 r_ne20_ap_na23 0.0000000000000000D+00 r_ne20_ga_o16 0.0000000000000000D+00 r_ne20_gp_f19 0.0000000000000000D+00 - r_ne22_ag_mg26 0.0000000000000000D+00 - r_ne22_pg_na23 0.0000000000000000D+00 + r_ne22_ag_mg26 3.7722492592751900-162 + r_ne22_pg_na23 1.5522617227535292D-47 r_ni56_ga_fe52 0.0000000000000000D+00 - r_o14_ag_ne18 0.0000000000000000D+00 - r_o14_ap_f17 0.0000000000000000D+00 + r_o14_ag_ne18 4.2343952152106083-104 + r_o14_ap_f17 6.2223403374710201D-99 r_o14_gp_n13 0.0000000000000000D+00 - r_o15_ag_ne19 0.0000000000000000D+00 + r_o15_ag_ne19 9.5687159836124304-105 r_o15_ap_f18 0.0000000000000000D+00 r_o15_gp_n14 0.0000000000000000D+00 r_o16_ag_ne20 2.8337440971737175-106 @@ -765,19 +765,19 @@ r_o16_ga_c12 0.0000000000000000D+00 r_o16_gp_n15 0.0000000000000000D+00 r_o16_pg_f17 3.2893588509057685D-40 - r_o17_pa_n14 0.0000000000000000D+00 - r_o17_pg_f18 0.0000000000000000D+00 - r_o18_ag_ne22 0.0000000000000000D+00 + r_o17_pa_n14 1.3887592875354731D-40 + r_o17_pg_f18 2.0334669089529022D-40 + r_o18_ag_ne22 3.8529665567318468-125 r_o18_pa_n15 8.1448127072519293D-37 - r_o18_pg_f19 0.0000000000000000D+00 - r_s32_ag_ar36 0.0000000000000000D+00 + r_o18_pg_f19 6.7833491710710671D-40 + r_s32_ag_ar36 6.9053347096806561-171 r_s32_ga_si28 0.0000000000000000D+00 - r_si28_ag_s32 0.0000000000000000D+00 + r_si28_ag_s32 5.4756966972878739-154 r_si28_ga_mg24 0.0000000000000000D+00 - r_ti44_ag_cr48 0.0000000000000000D+00 + r_ti44_ag_cr48 3.0969058786136191-213 r_ti44_ga_ca40 0.0000000000000000D+00 - ral27pa_aux 0.0000000000000000D+00 - ral27pg_aux 0.0000000000000000D+00 + ral27pa_aux 1.9918779499772236D-64 + ral27pg_aux 1.3424791010808044D-58 rar36ap_aux 0.0000000000000000D+00 rar36ap_to_ca40 0.0000000000000000D+00 rar36gp_aux 0.0000000000000000D+00 @@ -792,39 +792,39 @@ rca40ap_to_ti44 0.0000000000000000D+00 rca40gp_aux 0.0000000000000000D+00 rca40gp_to_ar36 0.0000000000000000D+00 - rcl35pa_aux 0.0000000000000000D+00 - rcl35pg_aux 0.0000000000000000D+00 + rcl35pa_aux 1.2491308643655965D-76 + rcl35pg_aux 4.3247045064840110D-67 rco55gprot_aux 0.0000000000000000D+00 - rco55pg_aux 0.0000000000000000D+00 - rco55protg_aux 0.0000000000000000D+00 - rcr48ap_aux 0.0000000000000000D+00 - rcr48ap_to_fe52 0.0000000000000000D+00 + rco55pg_aux 1.3620833659952599D-95 + rco55protg_aux 1.3620833659952599D-95 + rcr48ap_aux 6.1613814457736093-230 + rcr48ap_to_fe52 6.1613814457736093-230 rcr48gp_aux 0.0000000000000000D+00 rcr48gp_to_ti44 0.0000000000000000D+00 - rf19pg_aux 0.0000000000000000D+00 - rfe52ap_aux 0.0000000000000000D+00 - rfe52ap_to_ni56 0.0000000000000000D+00 - rfe52aprot_aux 0.0000000000000000D+00 - rfe52aprot_to_fe54 0.0000000000000000D+00 - rfe52aprot_to_ni56 0.0000000000000000D+00 + rf19pg_aux 8.3225081090292684D-45 + rfe52ap_aux 2.4948190097338029-243 + rfe52ap_to_ni56 2.4948190097338029-243 + rfe52aprot_aux 2.4948190097338029-243 + rfe52aprot_to_fe54 2.4948190097338029-243 + rfe52aprot_to_ni56 2.4948190097338029-243 rfe52gp_aux 0.0000000000000000D+00 rfe52gp_to_cr48 0.0000000000000000D+00 rfe52neut_to_fe54 -1.0000000000000000D+00 - rfe52ng_aux 0.0000000000000000D+00 + rfe52ng_aux 7.5610631459147874D+06 rfe53gn_aux 0.0000000000000000D+00 - rfe53ng_aux 0.0000000000000000D+00 - rfe54a_to_ni56 0.0000000000000000D+00 + rfe53ng_aux 1.1905625028580947D+07 + rfe54a_to_ni56 3.8728400667314115-244 rfe54an_aux 0.0000000000000000D+00 rfe54an_to_ni56 0.0000000000000000D+00 rfe54aprot_to_fe56 0.0000000000000000D+00 rfe54g_to_fe52 0.0000000000000000D+00 - rfe54ng_aux 0.0000000000000000D+00 + rfe54ng_aux 1.8827867900769722D+08 rfe54ng_to_fe56 -1.0000000000000000D+00 rfe54prot_to_fe52 -1.0000000000000000D+00 rfe54prot_to_ni56 -1.0000000000000000D+00 - rfe54protg_aux 0.0000000000000000D+00 + rfe54protg_aux 3.3614445565734901D-93 rfe55gn_aux 0.0000000000000000D+00 - rfe55ng_aux 0.0000000000000000D+00 + rfe55ng_aux 1.1450537918056973D+03 rfe56ec_fake_to_mn56 -1.0000000000000000D+00 rfe56ec_fake_to_mn57 -1.0000000000000000D+00 rfe56ec_fake_to_cr56 -1.0000000000000000D+00 @@ -849,28 +849,28 @@ rhe4_breakup -1.0000000000000000D+00 rhe4_rebuild -1.0000000000000000D+00 rhe4g_neut_aux 0.0000000000000000D+00 - rk39pa_aux 0.0000000000000000D+00 - rk39pg_aux 0.0000000000000000D+00 + rk39pa_aux 1.1033282777003295D-78 + rk39pg_aux 7.2279010928964974D-72 rmg24ap_aux 0.0000000000000000D+00 rmg24ap_to_si28 0.0000000000000000D+00 rmg24gp_aux 0.0000000000000000D+00 rmg24gp_to_ne20 0.0000000000000000D+00 - rmn51pg_aux 0.0000000000000000D+00 + rmn51pg_aux 2.5610235199081346D-90 rn14_to_c12 1.9150164715337093D-36 rn14_to_n15 1.9150164715337093D-36 rn14_to_o16 1.9150164715337093D-36 - rn14ag_lite 0.0000000000000000D+00 + rn14ag_lite 5.7595359837054342D-97 rn14gc12 0.0000000000000000D+00 rn14pg_aux 1.9150164715337093D-36 - rn15pa_aux 0.0000000000000000D+00 - rn15pg_aux 0.0000000000000000D+00 - rna23pa_aux 0.0000000000000000D+00 - rna23pg_aux 0.0000000000000000D+00 - rne18ag_to_mg24 0.0000000000000000D+00 - rne18ap_to_mg22 0.0000000000000000D+00 - rne18ap_to_mg24 0.0000000000000000D+00 - rne19pg_to_mg22 0.0000000000000000D+00 - rne19pg_to_mg24 0.0000000000000000D+00 + rn15pa_aux 3.7461382926599108D-32 + rn15pg_aux 1.9741989978881749D-35 + rna23pa_aux 2.8956408617284671D-49 + rna23pg_aux 7.9640617295751365D-52 + rne18ag_to_mg24 2.9221476970600607-122 + rne18ap_to_mg22 8.7339936863597313-119 + rne18ap_to_mg24 8.7339936863597313-119 + rne19pg_to_mg22 1.9681988619045489D-49 + rne19pg_to_mg24 1.9681988619045489D-49 rne20ap_aux 0.0000000000000000D+00 rne20ap_to_mg24 0.0000000000000000D+00 rne20gp_aux 0.0000000000000000D+00 @@ -888,16 +888,16 @@ rni56gprot_to_fe52 -1.0000000000000000D+00 rni56gprot_to_fe54 -1.0000000000000000D+00 rni56ng_to_fe54 -1.0000000000000000D+00 - rni57na_aux 0.0000000000000000D+00 + rni57na_aux 2.1998591723647793D+06 ro16_to_n14 3.2893588509057685D-40 ro16_to_o17 3.2893588509057685D-40 ro16ap_aux 0.0000000000000000D+00 ro16ap_to_ne20 0.0000000000000000D+00 ro16gp_aux 0.0000000000000000D+00 ro16gp_to_c12 0.0000000000000000D+00 - ro17_to_o18 0.0000000000000000D+00 - rp31pa_aux 0.0000000000000000D+00 - rp31pg_aux 0.0000000000000000D+00 + ro17_to_o18 2.0334669089529022D-40 + rp31pa_aux 4.5316890996502531D-72 + rp31pg_aux 3.1039228145612459D-65 rpep_to_he3 1.2376709082622805D-27 rpp_to_he3 1.9001962005810273D-23 rprot_neutg_aux 4.5224899785473506D+04 @@ -906,8 +906,8 @@ rs32ap_to_ar36 0.0000000000000000D+00 rs32gp_aux 0.0000000000000000D+00 rs32gp_to_si28 0.0000000000000000D+00 - rsc43pa_aux 0.0000000000000000D+00 - rsc43pg_aux 0.0000000000000000D+00 + rsc43pa_aux 4.0726684295068816D-79 + rsc43pg_aux 2.2722098637423061D-78 rsi28ap_aux 0.0000000000000000D+00 rsi28ap_to_s32 0.0000000000000000D+00 rsi28gp_aux 0.0000000000000000D+00 @@ -917,54 +917,54 @@ rti44gp_aux 0.0000000000000000D+00 rti44gp_to_ca40 0.0000000000000000D+00 rv47pa_aux 0.0000000000000000D+00 - rv47pg_aux 0.0000000000000000D+00 + rv47pg_aux 1.2282139359429897D-70 r_h1_h1_wk_h2 1.9001962005810273D-23 r_h1_h1_ec_h2 1.2376709082622805D-27 - rn14ag_to_ne22 0.0000000000000000D+00 + rn14ag_to_ne22 5.7595359837054342D-97 rf19pa_aux 3.0412796782211088D-41 r_b8_wk_he4_he4 9.0019114358434449D-01 rmn51pa_aux 0.0000000000000000D+00 rfe54gn_aux 0.0000000000000000D+00 rco55pa_aux 0.0000000000000000D+00 rco55prota_aux 0.0000000000000000D+00 - rn14ag_to_o18 0.0000000000000000D+00 + rn14ag_to_o18 5.7595359837054342D-97 r_h1_he3_wk_he4 2.7547694540397324D-29 r_be7_wk_li7 6.9000396298341936D-09 - r_al27_pg_si28 0.0000000000000000D+00 + r_al27_pg_si28 1.3424791010808044D-58 r_si28_gp_al27 0.0000000000000000D+00 r_si28_ap_p31 0.0000000000000000D+00 - r_p31_pa_si28 0.0000000000000000D+00 - r_p31_pg_s32 0.0000000000000000D+00 + r_p31_pa_si28 4.5316890996502531D-72 + r_p31_pg_s32 3.1039228145612459D-65 r_s32_gp_p31 0.0000000000000000D+00 r_s32_ap_cl35 0.0000000000000000D+00 - r_cl35_pa_s32 0.0000000000000000D+00 - r_cl35_pg_ar36 0.0000000000000000D+00 + r_cl35_pa_s32 1.2491308643655965D-76 + r_cl35_pg_ar36 4.3247045064840110D-67 r_ar36_gp_cl35 0.0000000000000000D+00 r_ar36_ap_k39 0.0000000000000000D+00 - r_k39_pa_ar36 0.0000000000000000D+00 - r_k39_pg_ca40 0.0000000000000000D+00 + r_k39_pa_ar36 1.1033282777003295D-78 + r_k39_pg_ca40 7.2279010928964974D-72 r_ca40_gp_k39 0.0000000000000000D+00 r_ca40_ap_sc43 0.0000000000000000D+00 - r_sc43_pa_ca40 0.0000000000000000D+00 - r_sc43_pg_ti44 0.0000000000000000D+00 + r_sc43_pa_ca40 4.0726684295068816D-79 + r_sc43_pg_ti44 2.2722098637423061D-78 r_ti44_gp_sc43 0.0000000000000000D+00 r_ti44_ap_v47 0.0000000000000000D+00 r_v47_pa_ti44 0.0000000000000000D+00 - r_v47_pg_cr48 0.0000000000000000D+00 + r_v47_pg_cr48 1.2282139359429897D-70 r_cr48_gp_v47 0.0000000000000000D+00 - r_cr48_ap_mn51 0.0000000000000000D+00 + r_cr48_ap_mn51 6.1613814457736093-230 r_mn51_pa_cr48 0.0000000000000000D+00 - r_mn51_pg_fe52 0.0000000000000000D+00 + r_mn51_pg_fe52 2.5610235199081346D-90 r_fe52_gp_mn51 0.0000000000000000D+00 - r_fe52_ap_co55 0.0000000000000000D+00 + r_fe52_ap_co55 2.4948190097338029-243 r_co55_pa_fe52 0.0000000000000000D+00 - r_co55_pg_ni56 0.0000000000000000D+00 + r_co55_pg_ni56 1.3620833659952599D-95 r_ni56_gp_co55 0.0000000000000000D+00 - r_fe52_ng_fe53 0.0000000000000000D+00 + r_fe52_ng_fe53 7.5610631459147874D+06 r_fe53_gn_fe52 0.0000000000000000D+00 - r_fe53_ng_fe54 0.0000000000000000D+00 + r_fe53_ng_fe54 1.1905625028580947D+07 r_fe54_gn_fe53 0.0000000000000000D+00 - r_fe54_pg_co55 0.0000000000000000D+00 + r_fe54_pg_co55 3.3614445565734901D-93 r_co55_gp_fe54 0.0000000000000000D+00 r_he3_ng_he4 2.5565521689834746D+01 r_he4_gn_he3 0.0000000000000000D+00 @@ -972,17 +972,17 @@ r_h2_gn_h1 0.0000000000000000D+00 r_h2_pg_he3 8.5789830644283769D-07 r_he3_gp_h2 0.0000000000000000D+00 - r_fe54_ng_fe55 0.0000000000000000D+00 + r_fe54_ng_fe55 1.8827867900769722D+08 r_fe55_gn_fe54 0.0000000000000000D+00 - r_fe55_ng_fe56 0.0000000000000000D+00 + r_fe55_ng_fe56 1.1450537918056973D+03 r_fe56_gn_fe55 0.0000000000000000D+00 r_fe54_ap_co57 0.0000000000000000D+00 - r_co57_pa_fe54 0.0000000000000000D+00 - r_fe56_pg_co57 0.0000000000000000D+00 + r_co57_pa_fe54 6.4177508709520918-111 + r_fe56_pg_co57 2.7509929131992614D-93 r_co57_gp_fe56 0.0000000000000000D+00 - r_c12_c12_to_h1_na23 0.0000000000000000D+00 + r_c12_c12_to_h1_na23 1.0572981387735790-221 r_he4_ne20_to_c12_c12 0.0000000000000000D+00 - r_c12_c12_to_he4_ne20 0.0000000000000000D+00 + r_c12_c12_to_he4_ne20 1.4271003292515259-221 r_he4_mg24_to_c12_o16 0.0000000000000000D+00 r_al26-1_to_al26-2 0.0000000000000000D+00 r_al26-2_to_al26-1 0.0000000000000000D+00