From b7aaf333114a5cbf8241f05448043e739e1a4e98 Mon Sep 17 00:00:00 2001
From: church89 <l.chierici89@gmail.com>
Date: Mon, 27 Mar 2023 11:30:00 +0200
Subject: [PATCH] add plots

---
 msre_power_history.py | 23 +++++++++++++++++++++++
 1 file changed, 23 insertions(+)

diff --git a/msre_power_history.py b/msre_power_history.py
index 6ae0574..0f91754 100644
--- a/msre_power_history.py
+++ b/msre_power_history.py
@@ -256,9 +256,32 @@ else:
 tally.scores.append(heating_score)
 tallies = openmc.Tallies([tally])
 
+#Plots
+colors = {salt:'yellow', graphite:'black', inor: 'grey', helium: 'cyan', inconel: 'grey',
+              bush: 'blue', ss316: 'grey', concrete: 'brown', shield: 'red', insulation: 'green',
+              sandwater: 'lightgreen', steel: 'grey'}
+plot = openmc.Plot.from_geometry(geometry)
+plot.basis = 'xy'
+plot.width = (150,150)
+plot.pixels = (200,200)
+plot.origin = (0,0,150)
+plot.color_by = 'material'
+plot.colors = colors
+openmc.plot_inline(plot)
+
+plot.basis = 'xz'
+plot.width = (150,200)
+plot.pixels = (150,200)
+plot.origin = (0,-5,150)
+plot.color_by = 'material'
+plot.colors = colors
+openmc.plot_inline(plot)
+
 # Depletion settings
 model = openmc.model.Model(geometry,mats,settings,tallies)
 
+results=model.run()
+
 op = openmc.deplete.CoupledOperator(model,
     normalization_mode = "energy-deposition",
     chain_file='/home/lorenzo/ca_depletion_chains/ENDF-B-VIII.0_chain_msr.xml')