Update grid.cc

artis-mcrt · Mar 7, 2025 · f7ea42f · f7ea42f
1 parent 5dd43a4
commit f7ea42f
Showing 1 changed file with 7 additions and 6 deletions.
diff --git a/grid.cc b/grid.cc
@@ -2543,13 +2543,14 @@ auto get_totmassradionuclide(const int z, const int a) -> double {
         d_coordminboundary[d] = -1;
       }
 
-      const double t_coordmaxboundary = ((pktposgridcoord[d] - (pktvelgridcoord[d] * tstart)) /
-                                         (cellcoordmax[d] - (pktvelgridcoord[d] * globals::tmin)) * globals::tmin) -
-                                        tstart;
-
-      d_coordmaxboundary[d] = CLIGHT_PROP * t_coordmaxboundary;
-      if ((pktvelgridcoord[d] - (pktposgridcoord[d] / tstart)) < 0) {
+      if (pktvelgridcoord[d] < (pktposgridcoord[d] / tstart)) {
         d_coordmaxboundary[d] = -1;
+      } else {
+        const double t_coordmaxboundary = ((pktposgridcoord[d] - (pktvelgridcoord[d] * tstart)) /
+                                           (cellcoordmax[d] - (pktvelgridcoord[d] * globals::tmin)) * globals::tmin) -
+                                          tstart;
+
+        d_coordmaxboundary[d] = CLIGHT_PROP * t_coordmaxboundary;
       }
     }
   } else {