Merge pull request #1 from wbinventor/multipole-option

Added user option to use multipole cross sections

.gitignore

@@ -98,4 +98,9 @@ ENV/
 *.png
 
 # Temp files
-*.xml~
+*.xml~
+
+# Cobalt
+*.output
+*.error
+*.cobaltlog

2x2-periodic/build-xml.py

@@ -170,6 +170,10 @@ openmc_geometry.export_to_xml()
 
 #### Create OpenMC "settings.xml" file
 
+# Query the user on whether to use multipole cross sections
+multipole = input('Use multipole cross sections? (y/n): ').lower()
+multipole = True if multipole == 'y' else False
+
 # Construct uniform initial source distribution over fissionable zones
 lower_left = two_by_two.bounds[:3]
 upper_right = two_by_two.bounds[3:]
@@ -184,6 +188,10 @@ settings_file.ptables = True
 settings_file.output = {'tallies': False}
 settings_file.source = source
 settings_file.sourcepoint_write = False
+
+if multipole:
+    settings_file.temperature = {'multipole': True, 'tolerance': 1000}
+
 settings_file.export_to_xml()
 
 
@@ -276,4 +284,4 @@ capture_rates.scores = ['absorption', 'fission']
 tallies_file = openmc.Tallies([fission_rates, capture_rates])
 cell_mgxs_lib.add_to_tallies_file(tallies_file, merge=True)
 mat_mgxs_lib.add_to_tallies_file(tallies_file, merge=True)
-tallies_file.export_to_xml()
+tallies_file.export_to_xml()

2x2-reflector/build-xml.py

@@ -178,6 +178,10 @@ openmc_geometry.export_to_xml()
 
 #### Create OpenMC "settings.xml" file
 
+# Query the user on whether to use multipole cross sections
+multipole = input('Use multipole cross sections? (y/n): ').lower()
+multipole = True if multipole == 'y' else False
+
 # Construct uniform initial source distribution over fissionable zones
 lower_left = reflector.bounds[:3]
 upper_right = reflector.bounds[3:]
@@ -196,6 +200,10 @@ settings_file.ptables = True
 settings_file.output = {'tallies': False}
 settings_file.source = source
 settings_file.sourcepoint_write = False
+
+if multipole:
+    settings_file.temperature = {'multipole': True, 'tolerance': 1000}
+
 settings_file.export_to_xml()
 
 
@@ -288,4 +296,4 @@ capture_rates.scores = ['absorption', 'fission']
 tallies_file = openmc.Tallies([fission_rates, capture_rates])
 cell_mgxs_lib.add_to_tallies_file(tallies_file, merge=True)
 mat_mgxs_lib.add_to_tallies_file(tallies_file, merge=True)
-tallies_file.export_to_xml()
+tallies_file.export_to_xml()

assembly/build-xml.py

@@ -90,6 +90,10 @@ openmc_geometry.export_to_xml()
 
 #### Create OpenMC "settings.xml" file
 
+# Query the user on whether to use multipole cross sections
+multipole = input('Use multipole cross sections? (y/n): ').lower()
+multipole = True if multipole == 'y' else False
+
 # Construct uniform initial source distribution over fissionable zones
 lower_left = fuel_assembly.bounds[:3]
 upper_right = fuel_assembly.bounds[3:]
@@ -104,6 +108,10 @@ settings_file.ptables = True
 settings_file.output = {'tallies': False}
 settings_file.source = source
 settings_file.sourcepoint_write = False
+
+if multipole:
+    settings_file.temperature = {'multipole': True, 'tolerance': 1000}
+
 settings_file.export_to_xml()
 
 
@@ -196,4 +204,4 @@ capture_rates.scores = ['absorption', 'fission']
 tallies_file = openmc.Tallies([fission_rates, capture_rates])
 cell_mgxs_lib.add_to_tallies_file(tallies_file, merge=True)
 mat_mgxs_lib.add_to_tallies_file(tallies_file, merge=True)
-tallies_file.export_to_xml()
+tallies_file.export_to_xml()

pin-cell/build-xml.py

@@ -90,6 +90,10 @@ openmc_geometry.export_to_xml()
 
 #### Create OpenMC "settings.xml" file
 
+# Query the user on whether to use multipole cross sections
+multipole = input('Use multipole cross sections? (y/n): ').lower()
+multipole = True if multipole == 'y' else False
+
 # Construct uniform initial source distribution over fissionable zones
 lower_left = pin_geometry.bounds[:2] + [-10.]
 upper_right = pin_geometry.bounds[3:5] + [10.]
@@ -104,6 +108,10 @@ settings_file.ptables = True
 settings_file.output = {'tallies': False}
 settings_file.source = source
 settings_file.sourcepoint_write = False
+
+if multipole:
+    settings_file.temperature = {'multipole': True, 'tolerance': 1000}
+
 settings_file.export_to_xml()
 
 
@@ -153,4 +161,4 @@ mgxs_lib.build_library()
 # Create a "tallies.xml" file for the MGXS Library
 tallies_file = openmc.Tallies()
 mgxs_lib.add_to_tallies_file(tallies_file, merge=True)
-tallies_file.export_to_xml()
+tallies_file.export_to_xml()