Update update_grid.cc

artis-mcrt · May 1, 2024 · 67086f0 · 67086f0
1 parent 895a2e8
commit 67086f0
Showing 1 changed file with 8 additions and 4 deletions.
diff --git a/update_grid.cc b/update_grid.cc
@@ -1253,11 +1253,15 @@ void update_grid(FILE *estimators_file, const int nts, const int nts_prev, const
     // generate a loop over all mgi
     if constexpr (DETAILED_BF_ESTIMATORS_ON) {
       for (int mgi = 0; mgi < grid::get_npts_model(); mgi++) {
-        double estimator_normfactor =
-            1 / (grid::get_modelcell_assocvolume_tmin(mgi) * pow(globals::timesteps[nts_prev].mid / globals::tmin, 3)) /
-            globals::timesteps[nts_prev].width / globals::nprocs;
+        const int assoc_cells = grid::get_numassociatedcells(mgi);
+        if (assoc_cells > 0) {
+          double estimator_normfactor =
+              1 /
+              (grid::get_modelcell_assocvolume_tmin(mgi) * pow(globals::timesteps[nts_prev].mid / globals::tmin, 3)) /
+              globals::timesteps[nts_prev].width / globals::nprocs;
 
-        radfield::normalise_bf_estimators(mgi, estimator_normfactor / H);
+          radfield::normalise_bf_estimators(mgi, estimator_normfactor / H);
+        }
       }
     }