Commit 70f1d98f authored by Wysocki, Aaron's avatar Wysocki, Aaron
Browse files

Setting gt radpw to 0 and ax shape to 0s; CTF reg test giving timeout

parent d635907a

SubKit/build/SquareLatticeLWR_Nodal.py

+10 −3
Original line number Diff line number Diff line
@@ -1162,12 +1162,13 @@ class SquareLatticeLWR_Nodal(CoreBuilder.Core): 
        # renormalize

        p = list(np.array(p) / np.average(p))



        me.model.addAxialPowerShape('allrods', z, p)

        me.model.addAxialPowerShape('fuelrods', z, p)

#        me.model.addAxialPowerShape('guidetubes', z, p)



        radPowers = []

        weights = []

        for obj in me.model.solidObjects.values():

            obj.powID = 'allrods'

#            obj.powID = 'allrods'

            # Need to weight the rod power when normalizing radial power shape

            assert(obj.symtype == 1)

            if radialPower:
@@ -1185,6 +1186,9 @@ class SquareLatticeLWR_Nodal(CoreBuilder.Core): 
                                    objType = me.model.solidTypeIDs[solidID]

                                    if isinstance(me.model.solidTypes[objType], Solid.FuelRod):

                                        obj.power = radialPower[i][j]

                                        obj.powID = 'fuelrods'

#                                    else:

#                                        obj.powID = 'guidetubes'

                else:

                    raise RuntimeError(

                        "Radial function not supported in setCorePowerShape")
@@ -1210,7 +1214,10 @@ class SquareLatticeLWR_Nodal(CoreBuilder.Core): 
#            mults.append(obj.mult)

        radP = np.array(radP) / np.average(np.array(radP), weights=mults)

        for i, obj in enumerate(me.model.solidObjects.values()):

            if mults[i] > 0:

                obj.power = radP[i]

            else:

                obj.power = 1.0



    def setCorePowerShape3D(me, powerDist, wgts=None, coreStart=0.0, interp_type="nearest"):

        """