Add support for TRC macrobody (#19)

README.md

@@ -29,7 +29,7 @@ Some geometry features are not currently supported:
 
 - Periodic boundary conditions
 - `X`, `Y`, and `Z` surfaces with 1 or 3 coordinate pairs
-- `RHP`, `REC`, `TRC`, `ELL`, `WED`, and `ARB` macrobodies
+- `RHP`, `REC`, `ELL`, `WED`, and `ARB` macrobodies
 - Hexagonal lattices
 - One-dimensional lattices
 - Two-dimensional lattices with basis other than x-y

src/openmc_mcnp_adapter/openmc_conversion.py

@@ -15,6 +15,7 @@
     RightCircularCylinder as RCC,
     RectangularParallelepiped as RPP,
     OrthogonalBox as BOX,
+    ConicalFrustum as TRC,
 )
 from openmc.model import surface_composite
 
@@ -45,6 +46,11 @@
         4: ('ymin', False),
         5: ('zmax', True),
         6: ('zmin', False)
+    },
+    TRC: {
+        1: ('cone', False),
+        2: ('plane_top', False),
+        3: ('plane_bottom', True),
     }
 }
 
@@ -356,6 +362,12 @@ def flip_sense(surf):
                 surf = BOX(v, a1, a2, a3)
             else:
                 surf = BOX(v, a1, a2)
+        elif s['mnemonic'] == 'trc':
+            v = coeffs[:3]
+            h = coeffs[3:6]
+            r1 = coeffs[6]
+            r2 = coeffs[7]
+            surf = TRC(v, h, r1, r2)
         else:
             raise NotImplementedError('Surface type "{}" not supported'
                                       .format(s['mnemonic']))