From 57505522aa33b05fbb7e093b919b04248693fc95 Mon Sep 17 00:00:00 2001 From: Luke Shingles Date: Sun, 16 Feb 2025 10:54:43 +0000 Subject: [PATCH] Fixup call to calculate_nt_ionization_ratecoeff() --- nonthermal.cc | 3 +-- 1 file changed, 1 insertion(+), 2 deletions(-) diff --git a/nonthermal.cc b/nonthermal.cc index 120fc0503..5deef8245 100644 --- a/nonthermal.cc +++ b/nonthermal.cc @@ -1190,10 +1190,9 @@ auto nt_ionization_ratecoeff_wfapprox(const int modelgridindex, const int elemen // coefficient i.e. multiply this by ion population to get a rate of ionizations per second Do not call during packet // propagation, as the y vector may not be in memory! IMPORTANT: we are dividing by the shell potential, not the // valence potential here! To change this set assumeshellpotentialisvalence to true -auto calculate_nt_ionization_ratecoeff(const int modelgridindex, const int element, const int ion, +auto calculate_nt_ionization_ratecoeff(const int nonemptymgi, const int element, const int ion, const bool assumeshellpotentialisvalence, const std::array &yfunc) -> double { - const auto nonemptymgi = grid::get_nonemptymgi_of_mgi(modelgridindex); std::array cross_section_vec{}; auto gsl_cross_section_vec = gsl_vector_view_array(cross_section_vec.data(), SFPTS).vector; std::array cross_section_vec_allshells{};