Upped particles per batch and added entropy trackingHEADmaster

1 files changed, 11 insertions, 3 deletions

diff --git a/krusty.py b/krusty.py
index 9e0ea06..a29dcb7 100755
--- a/krusty.py
+++ b/krusty.py
@@ -5,6 +5,7 @@ import argparse
 import hashlib
 import sys
 import time
+import numpy as np
 
 # KRUSTY Model in OpenMC
 # C.Moore, <connor.moore@psi.ch>
@@ -264,7 +265,7 @@ geo_obj.export_to_xml()
 
 ## Settings ##
 set_obj = openmc.Settings()
-set_obj.particles = 20000
+set_obj.particles = 100000
 set_obj.batches = 4500
 set_obj.inactive = 500
 set_obj.temperature = {"method": "interpolation"}
@@ -272,12 +273,19 @@ set_obj.temperature = {"method": "interpolation"}
 set_src = openmc.IndependentSource()
 set_src.space = openmc.stats.Point((0.0,0.0,30.0))
 set_src.space = openmc.stats.CylindricalIndependent(
-        r = openmc.stats.Uniform(geo_air_inner_cyl.r,geo_fuel_cyl.r),
-        z = openmc.stats.Uniform(geo_z_core_bot.z0,geo_z_core_top.z0),
+        r = openmc.stats.Uniform(0.0, geo_rad_cyl[-1].r),
+        z = openmc.stats.Uniform(geo_z_abs_bot.z0,geo_z_abs_top.z0),
         phi = openmc.stats.Uniform(0.0,2*math.pi),
         origin = (0.0,0.0,0.0)
         )
 
+set_entropy_mesh = openmc.RegularMesh()
+set_entropy_mesh.lower_left = (-geo_rad_cyl[-1].r, -geo_rad_cyl[-1].r, geo_z_abs_bot.z0)
+set_entropy_mesh.upper_right = (geo_rad_cyl[-1].r, geo_rad_cyl[-1].r, geo_z_abs_top.z0)
+set_entropy_mesh.dimension = (30,30,30)
+
+set_obj.entropy_mesh=set_entropy_mesh
+
 #set_obj.verbosity=10
 set_obj.export_to_xml()