Commit 89551d92 authored by Gurecky, William's avatar Gurecky, William
Browse files

update gold nodal_p7_power3d nodal inp file

parent 2559f6ef

SubKit/build/SquareLatticeLWR_Nodal.py

+4 −4
Original line number Diff line number Diff line
@@ -1075,8 +1075,8 @@ class SquareLatticeLWR_Nodal(CoreBuilder.Core): 


               # interpolate to the z_bounds

               asm_axial_power_ctf = asm_axial_power_interpolant(me.z_power_grid)

               avg_axpwr = np.trapz(asm_axial_power_ctf, x=me.z_power_grid) / \

                     (np.max(me.z_power_grid) - np.min(me.z_power_grid))

               avg_axpwr = np.max((np.trapz(asm_axial_power_ctf, x=me.z_power_grid) / \

                     (np.max(me.z_power_grid) - np.min(me.z_power_grid)), 1e-12))



               me.model.addAxialPowerShape(powID, me.z_power_grid, asm_axial_power_ctf / avg_axpwr)

               rods = me.rodsInAssem[assemID]
@@ -1174,8 +1174,8 @@ class SquareLatticeLWR_Nodal(CoreBuilder.Core): 


                  # interpolate to the z_bounds

                  node_axial_power_ctf = node_axial_power_interpolant(me.z_power_grid)

                  avg_axpwr = np.trapz(node_axial_power_ctf, x=me.z_power_grid) / \

                        (np.max(me.z_power_grid) - np.min(me.z_power_grid))

                  avg_axpwr = np.max((np.trapz(node_axial_power_ctf, x=me.z_power_grid) / \

                        (np.max(me.z_power_grid) - np.min(me.z_power_grid)), 1e-12))



                  # set axial power profile in the ctf model

                  me.model.addAxialPowerShape(powID, me.z_power_grid, node_axial_power_ctf / avg_axpwr)

tests/nodal_p7_pow3d/deck.inp.gold

+79419 −41591

File changed.

Preview size limit exceeded, changes collapsed.