examples: updated OUTP examples

Author: BitterB0NG0

examples/example_read_F1F8out.py

@@ -0,0 +1,16 @@
+"""
+Examples for reading INP cells using ``from_mcnp``.
+"""
+
+import pathlib
+
+import pymcnp
+
+# Reading INP.
+path = pathlib.Path(__file__).parent / 'files' / 'inp' / 'F1F8.i'
+inp = pymcnp.Inp.from_file(path)
+
+# Reading cell.
+cell = inp.cells[0]
+
+print(cell)

examples/example_read_energy_and_time.py

@@ -0,0 +1,16 @@
+"""
+Examples for reading INP materials using ``from_mcnp``.
+"""
+
+import pathlib
+
+import pymcnp
+
+# Reading INP.
+path = pathlib.Path(__file__).parent / 'files' / 'inp' / 'F1F8.i'
+inp = pymcnp.Inp.from_file(path)
+
+# Reading material.
+material = inp.data[0]
+
+print(material)

examples/example_read_outp.py

@@ -0,0 +1,16 @@
+"""
+Examples for reading INP surfaces using ``from_mcnp``.
+"""
+
+import pathlib
+
+import pymcnp
+
+# Reading INP.
+path = pathlib.Path(__file__).parent / 'files' / 'inp' / 'F1F8.i'
+inp = pymcnp.Inp.from_file(path)
+
+# Reading surface.
+surface = inp.surfaces[0]
+
+print(surface)

examples/example_read_time.py

@@ -0,0 +1,29 @@
+"""
+Example reading OUTP type-1 tallies.
+"""
+
+import pathlib
+
+import pymcnp
+import matplotlib.pyplot as plt
+
+TALLY = '1'
+SURFACE = '2.1'
+BIN_WIDTH = 5
+
+# Reading OUTP.
+path = pathlib.Path(__file__).parent / 'files' / 'outp' / 'F1.o'
+outp = pymcnp.Outp.from_file(path)
+
+# Reading tallies.
+tallies = outp.to_dataframe()
+tally = tallies[TALLY]
+surface = tally.loc[tally['surface'] == SURFACE]
+
+# Plotting tally.
+plt.figure()
+plt.bar(surface['bins'], surface['counts'], width=BIN_WIDTH)
+plt.xlabel('Bins')
+plt.ylabel('Counts')
+plt.title(f'Type-1 Tally: Bins vs Counts')
+plt.show()

examples/get_cell.py

@@ -0,0 +1,27 @@
+"""
+Example reading OUTP type-2 tallies.
+"""
+
+import pathlib
+
+import pymcnp
+import matplotlib.pyplot as plt
+
+TALLY = '2'
+SURFACE = '8'
+BIN_WIDTH = 5
+
+# Reading tallies.
+path = pathlib.Path(__file__).parent / 'files' / 'outp' / 'png_tmp.o'
+outp = pymcnp.Outp.from_file(path)
+tallies = outp.to_dataframe()
+tally = tallies[TALLY]
+surface = tally.loc[tally['surface'] == SURFACE]
+
+# Plotting tally.
+plt.figure()
+plt.bar(surface['bins'], surface['counts'], width=BIN_WIDTH)
+plt.xlabel('Bins')
+plt.ylabel('Counts')
+plt.title(f'Type-2 Tally: Bins vs Counts')
+plt.show()

examples/get_material.py

@@ -0,0 +1,27 @@
+"""
+Example reading OUTP type-4 tallies.
+"""
+
+import pathlib
+
+import pymcnp
+import matplotlib.pyplot as plt
+
+TALLY = '14'
+CELL = '12'
+BIN_WIDTH = 5
+
+# Reading tallies.
+path = pathlib.Path(__file__).parent / 'files' / 'outp' / 'F4.o'
+outp = pymcnp.Outp.from_file(path)
+tallies = outp.to_dataframe()
+tally = tallies[TALLY]
+cell = tally.loc[tally['cell'] == CELL]
+
+# Plotting tally.
+plt.figure()
+plt.bar(cell['bins'], cell['counts'], width=BIN_WIDTH)
+plt.xlabel('Bins')
+plt.ylabel('Counts')
+plt.title(f'Type-4 Tally: Bins vs Counts')
+plt.show()