Skip to content

Commit 52e48cc

Browse files
rjfarmerrjfarmer
committed
Add missing r_he4_ap_li7
When we switched r_h1_li7_to_he4_he4 to r_li7_pa_he4 we lost the reverse rate. Reported-by Hannah Brinkman
1 parent c263191 commit 52e48cc

File tree

5 files changed

+16
-1
lines changed

5 files changed

+16
-1
lines changed

‎data/rates_data/reactions.list‎

Lines changed: 1 addition & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -355,6 +355,7 @@ r_h2_pg_he3 1 h1 + 1 h2 => 1 he3
355355
r_he3_ag_be7 1 he3 + 1 he4 => 1 be7 2 he3 he4 pp he4(he3,g)be7
356356
r_he3_be7_to_h1_h1_he4_he4 1 he3 + 1 be7 => 2 h1 + 2 he4 2 he3 be7 pp he3 + be7 => 2 p + 2 he4
357357
r_li7_pa_he4 1 li7 + 1 h1 => 2 he4 2 h1 li7 pp li7(p,a)he4
358+
r_he4_ap_li7 1 he4 + 1 he4 => 1 li7 + 1 h1 2 he4 he4 other he4(a,p)li7
358359
r_n13_ap_o16 1 n13 + 1 he4 => 1 h1 + 1 o16 2 he4 n13 cno n13(a,p)o16
359360
r_n13_gp_c12 1 n13 => 1 h1 + 1 c12 1 photo n13(g,p)c12
360361
r_n13_pg_o14 1 n13 + 1 h1 => 1 o14 2 h1 n13 cno n13(p,g)o14

‎rates/private/rates_names.f90‎

Lines changed: 1 addition & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -87,6 +87,7 @@ subroutine set_reaction_names
8787
reaction_Name(ir_h1_h1_he4_to_he3_he3) = 'r_h1_h1_he4_to_he3_he3'
8888
reaction_Name(ir_he4_he4_he4_to_c12) = 'r_he4_he4_he4_to_c12'
8989
reaction_Name(ir_li7_pa_he4) = 'r_li7_pa_he4'
90+
reaction_Name(ir_he4_ap_li7) = 'r_he4_ap_li7'
9091
reaction_Name(ir_mg24_ag_si28) = 'r_mg24_ag_si28'
9192
reaction_Name(ir_mg24_ap_al27) = 'r_mg24_ap_al27'
9293
reaction_Name(ir_mg24_ga_ne20) = 'r_mg24_ga_ne20'

‎rates/private/raw_rates.f90‎

Lines changed: 3 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -224,6 +224,9 @@ subroutine set_raw_rate(ir, temp, tf, raw_rate, ierr)
224224
case(ir_li7_pa_he4) ! li7(p, a)he4
225225
call do1(rate_li7pa_nacre)
226226

227+
case(ir_he4_ap_li7) ! li7(p, a)he4
228+
call do1_reverse(rate_li7pa_nacre)
229+
227230
case(ir_n13_gp_c12) ! n13(g, p)c12
228231
call do1_reverse(rate_c12pg_nacre)
229232

‎rates/public/rates_def.f90‎

Lines changed: 2 additions & 1 deletion
Original file line numberDiff line numberDiff line change
@@ -362,7 +362,8 @@ end subroutine interpolate_weak_table
362362
integer, parameter :: ir_h1_h1_he4_to_he3_he3 = ir_he3_he3_to_h1_h1_he4+1
363363
integer, parameter :: ir_he4_he4_he4_to_c12 = ir_h1_h1_he4_to_he3_he3+1
364364
integer, parameter :: ir_li7_pa_he4 = ir_he4_he4_he4_to_c12+1
365-
integer, parameter :: ir_mg24_ag_si28 = ir_li7_pa_he4+1
365+
integer, parameter :: ir_he4_ap_li7 = ir_li7_pa_he4+1
366+
integer, parameter :: ir_mg24_ag_si28 = ir_he4_ap_li7+1
366367
integer, parameter :: ir_mg24_ap_al27 = ir_mg24_ag_si28+1
367368
integer, parameter :: ir_mg24_ga_ne20 = ir_mg24_ap_al27+1
368369
integer, parameter :: ir_n13_ap_o16 = ir_mg24_ga_ne20+1

‎rates/test/test_output‎

Lines changed: 9 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -417,6 +417,7 @@
417417
r_h1_h1_he4_to_he3_he3 0.0000000000000000D+00
418418
r_he4_he4_he4_to_c12 0.0000000000000000D+00
419419
r_li7_pa_he4 1.1481972998876799D-26
420+
r_he4_ap_li7 0.0000000000000000D+00
420421
r_mg24_ag_si28 0.0000000000000000D+00
421422
r_mg24_ap_al27 0.0000000000000000D+00
422423
r_mg24_ga_ne20 0.0000000000000000D+00
@@ -727,6 +728,7 @@
727728
r_h1_h1_he4_to_he3_he3 0.0000000000000000D+00
728729
r_he4_he4_he4_to_c12 0.0000000000000000D+00
729730
r_li7_pa_he4 2.8536849590579373D-15
731+
r_he4_ap_li7 0.0000000000000000D+00
730732
r_mg24_ag_si28 0.0000000000000000D+00
731733
r_mg24_ap_al27 0.0000000000000000D+00
732734
r_mg24_ga_ne20 0.0000000000000000D+00
@@ -1037,6 +1039,7 @@
10371039
r_h1_h1_he4_to_he3_he3 0.0000000000000000D+00
10381040
r_he4_he4_he4_to_c12 6.0425276466603135D-71
10391041
r_li7_pa_he4 1.3036488512272843D-07
1042+
r_he4_ap_li7 0.0000000000000000D+00
10401043
r_mg24_ag_si28 0.0000000000000000D+00
10411044
r_mg24_ap_al27 0.0000000000000000D+00
10421045
r_mg24_ga_ne20 0.0000000000000000D+00
@@ -1347,6 +1350,7 @@
13471350
r_h1_h1_he4_to_he3_he3 0.0000000000000000D+00
13481351
r_he4_he4_he4_to_c12 2.1189548985111926D-46
13491352
r_li7_pa_he4 1.7169556787367157D-02
1353+
r_he4_ap_li7 0.0000000000000000D+00
13501354
r_mg24_ag_si28 3.4094273819191269-162
13511355
r_mg24_ap_al27 9.1263133445433935-274
13521356
r_mg24_ga_ne20 0.0000000000000000D+00
@@ -1657,6 +1661,7 @@
16571661
r_h1_h1_he4_to_he3_he3 0.0000000000000000D+00
16581662
r_he4_he4_he4_to_c12 2.0403192412842943D-24
16591663
r_li7_pa_he4 4.3342837593883942D+01
1664+
r_he4_ap_li7 0.0000000000000000D+00
16601665
r_mg24_ag_si28 1.7640238996806857D-49
16611666
r_mg24_ap_al27 2.1146812954972443D-91
16621667
r_mg24_ga_ne20 0.0000000000000000D+00
@@ -1967,6 +1972,7 @@
19671972
r_h1_h1_he4_to_he3_he3 1.1772013589467652-212
19681973
r_he4_he4_he4_to_c12 7.9590873420121651D-13
19691974
r_li7_pa_he4 7.8177297732306224D+03
1975+
r_he4_ap_li7 2.2679834766347775-273
19701976
r_mg24_ag_si28 2.6469697259346481D-15
19711977
r_mg24_ap_al27 8.0770504592160291D-29
19721978
r_mg24_ga_ne20 3.1486672973522443-150
@@ -2277,6 +2283,7 @@
22772283
r_h1_h1_he4_to_he3_he3 1.2168141674471324D-69
22782284
r_he4_he4_he4_to_c12 3.4041066123320119D-10
22792285
r_li7_pa_he4 2.2984785945795634D+05
2286+
r_he4_ap_li7 4.0535520962769382D-82
22802287
r_mg24_ag_si28 8.0870253378865792D-04
22812288
r_mg24_ap_al27 3.1845928303516034D-06
22822289
r_mg24_ga_ne20 1.0831093286103609D-39
@@ -2587,6 +2594,7 @@
25872594
r_h1_h1_he4_to_he3_he3 2.9036145678610530D-23
25882595
r_he4_he4_he4_to_c12 2.2461964456422845D-10
25892596
r_li7_pa_he4 2.1991028257065220D+06
2597+
r_he4_ap_li7 1.3073398877915164D-20
25902598
r_mg24_ag_si28 1.4662488667811619D+01
25912599
r_mg24_ap_al27 9.0039206583785949D+02
25922600
r_mg24_ga_ne20 1.6781713325994766D-02
@@ -2897,6 +2905,7 @@
28972905
r_h1_h1_he4_to_he3_he3 1.2872091694058283D-07
28982906
r_he4_he4_he4_to_c12 2.1320831532048911D-11
28992907
r_li7_pa_he4 1.6571622519109823D+07
2908+
r_he4_ap_li7 4.4349794705459109D+00
29002909
r_mg24_ag_si28 6.9791363732282798D+02
29012910
r_mg24_ap_al27 3.0549748551998041D+06
29022911
r_mg24_ga_ne20 3.2362267060968735D+10

0 commit comments

Comments
 (0)