Plot a reactor facemap

Load a test reactor from the test suite and plot a dummy power distribution from it. You can plot any block parameter.

Out:

=========== Case Information ===========
=========== Input File Information ===========
=========== Reactor Cycle Information ===========
=========== Adding Assemblies to <Core: core id:140034066137536> ===========
=========== Verifying Assembly Configurations ===========
=========== Summarizing Source of Material Data for <Core: core id:140034066137536> ===========
=========== Applying Geometry Modifications ===========
=========== Adding Assemblies to <Core: sfp id:140034078157408> ===========
=========== Summarizing Source of Material Data for <Core: sfp id:140034078157408> ===========
=========== Applying Geometry Modifications ===========
=========== Interface Stack Summary  ===========
===========  Triggering Init Event ===========
=========== 01 - main                           Init            ===========
=========== 02 - fissionProducts                Init            ===========
=========== 03 - xsGroups                       Init            ===========
=========== 04 - fuelHandler                    Init            ===========
=========== 05 - history                        Init            ===========
=========== 06 - database                       Init            ===========
=========== 07 - memoryProfiler                 Init            ===========
=========== 08 - snapshot                       Init            ===========
===========  Completed Init Event ===========

from armi.reactor.tests import test_reactors
from armi.utils import plotting
from armi import configure, runLog

# configure ARMI
configure(permissive=True)

operator, reactor = test_reactors.loadTestReactor()
reactor.core.growToFullCore(None)
# set dummy power
for b in reactor.core.getBlocks():
    x, y, z = b.spatialLocator.getGlobalCoordinates()
    b.p.pdens = x ** 2 + y ** 2 + z ** 2

plotting.plotFaceMap(reactor.core, param="pdens", labelFmt="{0:.1e}")