Skip to content

Commit

Permalink
Update rpkt.cc
Browse files Browse the repository at this point in the history
  • Loading branch information
lukeshingles committed Oct 26, 2024
1 parent 088070f commit 23b1b6f
Showing 1 changed file with 5 additions and 5 deletions.
10 changes: 5 additions & 5 deletions rpkt.cc
Original file line number Diff line number Diff line change
Expand Up @@ -1135,11 +1135,11 @@ void calculate_chi_rpkt_cont(const double nu_cmf, Rpkt_continuum_absorptioncoeff
chi_rpkt_cont.total = chi_rpkt_cont.ffescat + chi_rpkt_cont.bf + chi_rpkt_cont.ffheat;

if (!std::isfinite(chi_rpkt_cont.total)) {
printout("[fatal] calculate_chi_rpkt_cont: resulted in non-finite chi_rpkt_cont.total ... abort\n");
printout("[fatal] es %g, ff %g, bf %g\n", chi_rpkt_cont.ffescat, chi_rpkt_cont.ffheat, chi_rpkt_cont.bf);
printout("[fatal] nbfcontinua %d\n", globals::nbfcontinua);
printout("[fatal] in cell %d with density %g\n", modelgridindex, grid::get_rho(modelgridindex));
printout("[fatal] pkt.nu_cmf %g\n", nu_cmf);
printoutf("[fatal] calculate_chi_rpkt_cont: resulted in non-finite chi_rpkt_cont.total ... abort\n");
printoutf("[fatal] es %g, ff %g, bf %g\n", chi_rpkt_cont.ffescat, chi_rpkt_cont.ffheat, chi_rpkt_cont.bf);
printoutf("[fatal] nbfcontinua %d\n", globals::nbfcontinua);
printoutf("[fatal] in cell %d with density %g\n", modelgridindex, grid::get_rho(modelgridindex));
printoutf("[fatal] pkt.nu_cmf %g\n", nu_cmf);
if (std::isfinite(chi_rpkt_cont.ffescat)) {
chi_rpkt_cont.ffheat = 0.;
chi_rpkt_cont.bf = 0.;
Expand Down

0 comments on commit 23b1b6f

Please sign in to comment.