Skip to content

Commit

Permalink
fix Urca compilation (#62)
Browse files Browse the repository at this point in the history
we renamed the rate keys
  • Loading branch information
zingale authored Nov 28, 2023
1 parent 260d2f5 commit e9abc84
Showing 1 changed file with 4 additions and 4 deletions.
8 changes: 4 additions & 4 deletions urca/spherical/init_1d.H
Original file line number Diff line number Diff line change
Expand Up @@ -37,8 +37,8 @@ AMREX_INLINE void fopt_urca_23(eos_t& eos_state, Real& fopt, Real& r_ecap, Real&
constexpr int do_T_derivatives{0};
evaluate_rates<do_T_derivatives, rate_t>(burn_state, rates);

r_ecap = rates.screened_rates(k_na23_to_ne23);
r_beta = rates.screened_rates(k_ne23_to_na23);
r_ecap = rates.screened_rates(k_Na23_to_Ne23);
r_beta = rates.screened_rates(k_Ne23_to_Na23);

// calculate e-capture and beta decay rates are in equilibrium
Real xr_ecap = burn_state.xn[ina23] * r_ecap;
Expand Down Expand Up @@ -103,8 +103,8 @@ AMREX_INLINE void composition_equilibrium(eos_t& eos_state){
constexpr int do_T_derivatives{0};
evaluate_rates<do_T_derivatives, rate_t>(burn_state, rates);

Real r_ecap = rates.screened_rates(k_na23_to_ne23);
Real r_beta = rates.screened_rates(k_ne23_to_na23);
Real r_ecap = rates.screened_rates(k_Na23_to_Ne23);
Real r_beta = rates.screened_rates(k_Ne23_to_Na23);

eos_state.xn[ine23] = problem_rp::na_ne_23/(1.0_rt + r_beta/r_ecap);
eos_state.xn[ina23] = problem_rp::na_ne_23 - eos_state.xn[ine23];
Expand Down

0 comments on commit e9abc84

Please sign in to comment.