From 727ec909112d4bd5d6c02a492385496667ce41f5 Mon Sep 17 00:00:00 2001 From: Luke Shingles Date: Sat, 8 Mar 2025 16:54:00 +0000 Subject: [PATCH] Fix pos_component_vel_relative_to_flow --- grid.cc | 4 ++-- 1 file changed, 2 insertions(+), 2 deletions(-) diff --git a/grid.cc b/grid.cc index 44ae9c452..c84b7e606 100644 --- a/grid.cc +++ b/grid.cc @@ -2451,12 +2451,12 @@ auto get_totmassradionuclide(const int z, const int a) -> double { // packet pos below min const double boundaryposmin = grid::get_cellcoordmin(cellindex, d) / globals::tmin * tstart; delta = pktposgridcoord[d] - boundaryposmin; - isoutside_thisside = pos_component_vel_relative_to_flow && (pktposgridcoord[d] < (boundaryposmin - 10.)); + isoutside_thisside = !pos_component_vel_relative_to_flow && (pktposgridcoord[d] < (boundaryposmin - 10.)); } else { // packet pos above max const double boundaryposmax = cellcoordmax[d] / globals::tmin * tstart; delta = pktposgridcoord[d] - boundaryposmax; - isoutside_thisside = !pos_component_vel_relative_to_flow && (pktposgridcoord[d] > (boundaryposmax + 10.)); + isoutside_thisside = pos_component_vel_relative_to_flow && (pktposgridcoord[d] > (boundaryposmax + 10.)); } if (isoutside_thisside) {