Bug fix for Polygon 'yz' basis (#3259)

openmc-dev · Jan 10, 2025 · c226c78 · c226c78
1 parent 51f0e6f
commit c226c78
Showing 1 changed file with 3 additions and 3 deletions.
diff --git a/openmc/model/surface_composite.py b/openmc/model/surface_composite.py
@@ -1252,11 +1252,11 @@ def _get_convex_hull_surfs(self, qhull):
             else:
                 op = operator.neg
                 if basis == 'xy':
-                    surf = openmc.Plane(a=dx, b=dy, d=-c)
+                    surf = openmc.Plane(a=dx, b=dy, c=0.0, d=-c)
                 elif basis == 'yz':
-                    surf = openmc.Plane(b=dx, c=dy, d=-c)
+                    surf = openmc.Plane(a=0.0, b=dx, c=dy, d=-c)
                 elif basis == 'xz':
-                    surf = openmc.Plane(a=dx, c=dy, d=-c)
+                    surf = openmc.Plane(a=dx, b=0.0, c=dy, d=-c)
                 else:
                     y0 = -c/dy
                     r2 = dy**2 / dx**2