def setup_class(cls):
libfile = os.getenv("BRIGHT_DATA") + '/LWR.h5'
BriPy.load_isos2track_hdf5(libfile)
cls.r1g = Reactor1G(default_rp, 'r1g')
cls.r1g.loadLib(libfile)
cls.r1g.IsosIn = MassStream({922350: 0.5, 922380: 0.5})
cls.r1g.doCalc()
def test_setParams(self):
lf = os.getenv("BRIGHT_DATA") + "/LWR.h5"
BriPy.load_isos2track_hdf5(lf)
rp = BriPy.LWRDefaults()
rp.BUt = 50.0
lwr = LightWaterReactor1G(lf, rp, 'lwr')
lwr.doCalc(MassStream({922350: 0.05, 922380:0.95}))
lwr.setParams()
assert_equal(lwr.ParamsIn["BUd"], 0.0)
assert_equal(lwr.ParamsOut["BUd"], lwr.BUd)
assert_equal(lwr.ParamsIn["U"], lwr.InU.mass)
assert_equal(lwr.ParamsOut["U"], lwr.OutU.mass)
assert_equal(lwr.ParamsIn["TRU"], lwr.InTRU.mass)
assert_equal(lwr.ParamsOut["TRU"], lwr.OutTRU.mass)
assert_equal(lwr.ParamsIn["ACT"], lwr.InACT.mass)
assert_equal(lwr.ParamsOut["ACT"], lwr.OutACT.mass)
assert_equal(lwr.ParamsIn["LAN"], lwr.InLAN.mass)
assert_equal(lwr.ParamsOut["LAN"], lwr.OutLAN.mass)
assert_equal(lwr.ParamsIn["FP"], 1.0 - lwr.InACT.mass - lwr.InLAN.mass)
assert_equal(lwr.ParamsOut["FP"], 1.0 - lwr.OutACT.mass - lwr.OutLAN.mass)
def test_setParams(self):
lf = os.getenv("BRIGHT_DATA") + "/FR.h5"
BriPy.load_isos2track_hdf5(lf)
rp = BriPy.FRDefaults()
rp.BUt = 140.0
fr = FastReactor1G(lf, rp, 'fr')
fr.doCalc(MassStream({922350: 0.30, 922380: 0.70}))
fr.setParams()
assert_equal(fr.ParamsIn["BUd"], 0.0)
assert_equal(fr.ParamsOut["BUd"], fr.BUd)
assert_equal(fr.ParamsIn["TRUCR"], 0.0)
assert_equal(fr.ParamsOut["TRUCR"], fr.TruCR)
assert_equal(fr.ParamsIn["P_NL"], 0.0)
assert_equal(fr.ParamsOut["P_NL"], fr.P_NL)
assert_equal(fr.ParamsIn["U"], fr.InU.mass)
assert_equal(fr.ParamsOut["U"], fr.OutU.mass)
assert_equal(fr.ParamsIn["TRU"], fr.InTRU.mass)
assert_equal(fr.ParamsOut["TRU"], fr.OutTRU.mass)
assert_equal(fr.ParamsIn["ACT"], fr.InACT.mass)
assert_equal(fr.ParamsOut["ACT"], fr.OutACT.mass)
assert_equal(fr.ParamsIn["LAN"], fr.InLAN.mass)
assert_equal(fr.ParamsOut["LAN"], fr.OutLAN.mass)
assert_equal(fr.ParamsIn["FP"], 1.0 - fr.InACT.mass - fr.InLAN.mass)
assert_equal(fr.ParamsOut["FP"], 1.0 - fr.OutACT.mass - fr.OutLAN.mass)
def test_load_isos2track_hdf5_4(self):
old_isos = BriPy.isos2track()
BriPy.isos2track([80160])
BriPy.load_isos2track_hdf5('isos.h5', '/NotIsos', True)
assert_equal(BriPy.isos2track(), [10010, 922350, 922380])
BriPy.isos2track(old_isos)
def test_load_isos2track_hdf5_1(self):
old_isos = BriPy.isos2track()
BriPy.isos2track([80160])
BriPy.load_isos2track_hdf5('isos.h5')
assert_equal(BriPy.isos2track(), [10010, 80160, 922350, 922380])
BriPy.isos2track(old_isos)
import os
import sys
import time
import tables as tb
import BriPy as bp
from facility_info import get_reactor_id_map
# Load isos2track
data_dir = os.getenv("BRIGHT_DATA")
lwr_data = data_dir + "/LWR.h5"
bp.load_isos2track_hdf5(lwr_data)
# We are going to capture the output, so don't write it out automatically.
bp.write_text(False)
bp.write_hdf5(False)
# Generate a used fuel table description
UsedFuel = {
'assembly_id': tb.StringCol(itemsize=15, pos=0),
'discharge_k': tb.Float64Col(pos=1),
'mass': tb.Float64Col(pos=2),
}
nbuf = len(UsedFuel)
i2t = bp.isos2track()
for i in xrange(len(i2t)):
UsedFuel[bp.isoname.zzaaam_2_LLAAAM(i2t[i])] = tb.Float64Col(pos=nbuf + i)
def setup_class(cls):
libfile = os.getenv("BRIGHT_DATA") + '/LWR.h5'
BriPy.load_isos2track_hdf5(libfile)
cls.r1g = Reactor1G()
cls.r1g.loadLib(libfile)
