From 6fd5188890859cd069d22c27e8255cd7ab6147ce Mon Sep 17 00:00:00 2001
From: Paul Romano <paul.k.romano@gmail.com>
Date: Wed, 24 Jul 2019 22:16:54 -0500
Subject: [PATCH] Remove neutron shields, which don't exist for NuScale model

---
 smr/smr/core.py | 67 +------------------------------------------------
 1 file changed, 1 insertion(+), 66 deletions(-)

diff --git a/smr/smr/core.py b/smr/smr/core.py
index 4d80fbf..aef9444 100644
--- a/smr/smr/core.py
+++ b/smr/smr/core.py
@@ -133,75 +133,10 @@ def core_geometry(num_rings, num_axial, depleted):
                    +surfs['lower bound'] & -surfs['upper bound'])
     root_univ.add_cell(cell)
 
-    # Neutron shield panels
-    cell = openmc.Cell(name='neutron shield panel NW')
-    cell.fill = mats['SS']
-    cell.region = (+surfs['core barrel OR'] & -surfs['neutron shield OR'] &
-                   +surfs['neutron shield NWbot SEtop'] &
-                   -surfs['neutron shield NWtop SEbot'] &
-                   +surfs['lower bound'] & -surfs['upper bound'])
-    root_univ.add_cell(cell)
-
-    cell = openmc.Cell(name='neutron shield panel N')
-    cell.fill = mats['H2O']
-    cell.region = (+surfs['core barrel OR'] & -surfs['neutron shield OR'] &
-                   +surfs['neutron shield NWtop SEbot'] &
-                   -surfs['neutron shield NEtop SWbot'] &
-                   +surfs['lower bound'] & -surfs['upper bound'])
-    root_univ.add_cell(cell)
-
-    cell = openmc.Cell(name='neutron shield panel SE')
-    cell.fill = mats['SS']
-    cell.region = (+surfs['core barrel OR'] & -surfs['neutron shield OR'] &
-                   -surfs['neutron shield NWbot SEtop'] &
-                   +surfs['neutron shield NWtop SEbot'] &
-                   +surfs['lower bound'] & -surfs['upper bound'])
-    root_univ.add_cell(cell)
-
-    cell = openmc.Cell(name='neutron shield panel E')
-    cell.fill = mats['H2O']
-    cell.region = (+surfs['core barrel OR'] & -surfs['neutron shield OR'] &
-                   +surfs['neutron shield NWbot SEtop'] &
-                   +surfs['neutron shield NEbot SWtop'] &
-                   +surfs['lower bound'] & -surfs['upper bound'])
-    root_univ.add_cell(cell)
-
-    cell = openmc.Cell(name='neutron shield panel NE')
-    cell.fill = mats['SS']
-    cell.region = (+surfs['core barrel OR'] & -surfs['neutron shield OR'] &
-                   +surfs['neutron shield NEbot SWtop'] &
-                   -surfs['neutron shield NEtop SWbot'] &
-                   +surfs['lower bound'] & -surfs['upper bound'])
-    root_univ.add_cell(cell)
-
-    cell = openmc.Cell(name='neutron shield panel S')
-    cell.fill = mats['H2O']
-    cell.region = (+surfs['core barrel OR'] & -surfs['neutron shield OR'] &
-                   -surfs['neutron shield NWtop SEbot'] &
-                   +surfs['neutron shield NEtop SWbot'] &
-                   +surfs['lower bound'] & -surfs['upper bound'])
-    root_univ.add_cell(cell)
-
-    cell = openmc.Cell(name='neutron shield panel SW')
-    cell.fill = mats['SS']
-    cell.region = (+surfs['core barrel OR'] & -surfs['neutron shield OR'] &
-                   -surfs['neutron shield NEbot SWtop'] &
-                   +surfs['neutron shield NEtop SWbot'] &
-                   +surfs['lower bound'] & -surfs['upper bound'])
-    root_univ.add_cell(cell)
-
-    cell = openmc.Cell(name='neutron shield panel W')
-    cell.fill = mats['H2O']
-    cell.region = (+surfs['core barrel OR'] & -surfs['neutron shield OR'] &
-                   -surfs['neutron shield NWbot SEtop'] &
-                   -surfs['neutron shield NEbot SWtop'] &
-                   +surfs['lower bound'] & -surfs['upper bound'])
-    root_univ.add_cell(cell)
-
     # Downcomer
     cell = openmc.Cell(name='downcomer')
     cell.fill = mats['H2O']
-    cell.region = (+surfs['neutron shield OR'] & -surfs['RPV IR'] &
+    cell.region = (+surfs['core barrel OR'] & -surfs['RPV IR'] &
                    +surfs['lower bound'] & -surfs['upper bound'])
     root_univ.add_cell(cell)