Commit ff3429b6 authored by April Novak's avatar April Novak
Browse files

Use consistent naming convention for rescaled pebble coordinates.

parent 8ba1cc54
...@@ -91,7 +91,7 @@ ...@@ -91,7 +91,7 @@
type = NearestPointAverage type = NearestPointAverage
variable = temp variable = temp
execute_on = 'timestep_end' # Evaluate on timestep end so OpenMC can get updated values next timestep execute_on = 'timestep_end' # Evaluate on timestep end so OpenMC can get updated values next timestep
points_file = 'pebble_centers.txt' points_file = 'pebble_centers_rescaled.txt'
[] []
[] []
......
...@@ -31,7 +31,7 @@ ...@@ -31,7 +31,7 @@
tally_type = 'cell' tally_type = 'cell'
pebble_cell_level = 1 pebble_cell_level = 1
mesh_template = 'sphere.e' mesh_template = 'sphere.e'
centers_file = 'pebble_centers.txt' centers_file = 'pebble_centers_rescaled.txt'
volumes = '14.137166941154067' volumes = '14.137166941154067'
[] []
......
...@@ -76,7 +76,7 @@ ...@@ -76,7 +76,7 @@
type = NearestPointAverage type = NearestPointAverage
variable = temp variable = temp
execute_on = 'initial timestep_end' # Evaluate on timestep end so OpenMC can get updated values next timestep execute_on = 'initial timestep_end' # Evaluate on timestep end so OpenMC can get updated values next timestep
points_file = 'pebble_centers.txt' points_file = 'pebble_centers_rescaled.txt'
[] []
[] []
......
...@@ -50,9 +50,9 @@ tolerance = 0.00 # Tolerance for pebbles univer ...@@ -50,9 +50,9 @@ tolerance = 0.00 # Tolerance for pebbles univer
radius_pebble_flibe = radius_pebble_outer+tolerance radius_pebble_flibe = radius_pebble_outer+tolerance
# Pebble centers coordinates (x,y,z) # Pebble centers coordinates (x,y,z)
print("Reading pseudo-random pebble centers from file pebble_centers.txt") print("Reading pseudo-random pebble centers from file pebble_centers_rescaled.txt")
pebble_centers = np.loadtxt('pebble_centers.txt', delimiter=' ') pebble_centers = np.loadtxt('pebble_centers_rescaled.txt', delimiter=' ')
print("File pebble_centers.txt reading completed") print("File pebble_centers_rescaled.txt reading completed")
# Vessel Dimensions # Vessel Dimensions
vessel_x, vessel_y = (0.0, 0.0) vessel_x, vessel_y = (0.0, 0.0)
......
...@@ -91,7 +91,7 @@ ...@@ -91,7 +91,7 @@
type = NearestPointAverage type = NearestPointAverage
variable = temp variable = temp
execute_on = 'timestep_end' # Evaluate on timestep end so OpenMC can get updated values next timestep execute_on = 'timestep_end' # Evaluate on timestep end so OpenMC can get updated values next timestep
points_file = 'pebble_centers.txt' points_file = 'pebble_centers_rescaled.txt'
[] []
[] []
......
...@@ -61,12 +61,10 @@ radius_pebble_central = radius_pebble_outer - 0.1 # UCBTH-14-002, Table 2-1; Cis ...@@ -61,12 +61,10 @@ radius_pebble_central = radius_pebble_outer - 0.1 # UCBTH-14-002, Table 2-1; Cis
tolerance = 0.00 # Tolerance for pebbles universe filling tolerance = 0.00 # Tolerance for pebbles universe filling
radius_pebble_flibe = radius_pebble_outer+tolerance radius_pebble_flibe = radius_pebble_outer+tolerance
# Pebble centers coordinates (x,y,z) # Pebble centers coordinates (x,y,z)
print("Reading pseudo-random pebble centers from file pebble_centers.txt") print("Reading pseudo-random pebble centers from file pebble_centers_rescaled.txt")
pebble_centers = np.loadtxt('list_pebbles_1568.csv', delimiter=',', skiprows=1)[:, 1:] pebble_centers = np.loadtxt('pebble_centers_rescaled.txt', delimiter=',')
print("File pebble_centers.txt reading completed") print("File pebble_centers_rescaled.txt reading completed")
rescaled_file = 'pebble_centers.txt' rescaled_file = 'pebble_centers_rescaled.txt'
np.savetxt(rescaled_file, pebble_centers.reshape((-1,3)))
print("Saved rescaled centers to " + rescaled_file)
# Vessel # Vessel
extra_thickness = 12.125 # change this value to set keff=1 extra_thickness = 12.125 # change this value to set keff=1
x_vessel = 0.0 # x position vessel x_vessel = 0.0 # x position vessel
......
...@@ -30,7 +30,7 @@ ...@@ -30,7 +30,7 @@
power = 156800 # 1000 per pebble power = 156800 # 1000 per pebble
tally_type = 'mesh' tally_type = 'mesh'
mesh_template = 'sphere.e' mesh_template = 'sphere.e'
centers_file = 'pebble_centers.txt' centers_file = 'pebble_centers_rescaled.txt'
volumes = '14.137166941154067' volumes = '14.137166941154067'
[] []
......
...@@ -76,7 +76,7 @@ ...@@ -76,7 +76,7 @@
type = NearestPointAverage type = NearestPointAverage
variable = temp variable = temp
execute_on = 'initial timestep_end' # Evaluate on timestep end so OpenMC can get updated values next timestep execute_on = 'initial timestep_end' # Evaluate on timestep end so OpenMC can get updated values next timestep
points_file = 'pebble_centers.txt' points_file = 'pebble_centers_rescaled.txt'
[] []
[] []
......
...@@ -50,9 +50,9 @@ tolerance = 0.00 # Tolerance for pebbles univer ...@@ -50,9 +50,9 @@ tolerance = 0.00 # Tolerance for pebbles univer
radius_pebble_flibe = radius_pebble_outer+tolerance radius_pebble_flibe = radius_pebble_outer+tolerance
# Pebble centers coordinates (x,y,z) # Pebble centers coordinates (x,y,z)
print("Reading pseudo-random pebble centers from file pebble_centers.txt") print("Reading pseudo-random pebble centers from file pebble_centers_rescaled.txt")
pebble_centers = np.loadtxt('list_pebbles_1568.csv', delimiter=',', skiprows=1)[:, 1:] pebble_centers = np.loadtxt('pebble_centers_rescaled.txt', delimiter=',')
print("File pebble_centers.txt reading completed") print("File pebble_centers_rescaled.txt reading completed")
# Vessel Dimensions # Vessel Dimensions
vessel_x, vessel_y = (0.0, 0.0) vessel_x, vessel_y = (0.0, 0.0)
......
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment