Merge pull request #1844 from paulromano/surface-tally-bug

Disallow non-current scores with a surface filter
openmc-dev · Jun 8, 2021 · f69149b · f69149b
2 parents f629ea0 + c6d2376
commit f69149b
Showing 1 changed file with 4 additions and 1 deletion.
diff --git a/src/tallies/tally.cpp b/src/tallies/tally.cpp
@@ -507,8 +507,11 @@ void Tally::set_scores(const vector<std::string>& scores)
   if (type_ == TallyType::SURFACE || type_ == TallyType::MESH_SURFACE) {
     if (scores_.size() != 1)
       fatal_error("Cannot tally other scores in the same tally as surface "
-        "currents");
+        "currents.");
   }
+  if ((surface_present || meshsurface_present) && scores_[0] != SCORE_CURRENT)
+    fatal_error("Cannot tally score other than 'current' when using a surface "
+      "or mesh-surface filter.");
 }
 
 void