<jupyter_start><jupyter_text>pyIsoDep Copyright (c) <NAME>  A scenario with multiple cross section sets   Description Example on how to execute a full depletion calculations.  Load modules<jupyter_code>import numpy as np
from pyIsoDep.functions.maindepletionsolver import MainDepletion
from pyIsoDep.functions.generatedata import TransmutationData
from pyIsoDep.functions.postprocessresults import Results
from pyIsoDep.functions.read_csv import ReadCsv<jupyter_output><empty_output><jupyter_text>Data Generation Stage<jupyter_code>FY_WGT = 0.0  # determines the fission yield wieghting
VOL = 332097.8  # cm^3

xsTable, fields = ReadCsv("./bootstrap.csv")
bootstrap = TransmutationData(libraryFlag=True, wgtFY=1.0)
bootstrap.ReadData(ID=xsTable[:, 0], sig_f=xsTable[:, 3], sig_c=xsTable[:, 2],
                   sig_alpha=xsTable[:, 4], flagBarns=True)

xsTable, fields = ReadCsv("./tempramp.csv")
tempramp = TransmutationData(libraryFlag=True, wgtFY=1.0)
tempramp.ReadData(ID=xsTable[:, 0], sig_f=xsTable[:, 3], sig_c=xsTable[:, 2],
                  sig_alpha=xsTable[:, 4], flagBarns=True)

xsTable, fields = ReadCsv("./fullthrust.csv")
fullthrust = TransmutationData(libraryFlag=True, wgtFY=1.0)
fullthrust.ReadData(ID=xsTable[:, 0], sig_f=xsTable[:, 3], sig_c=xsTable[:, 2],
                    sig_alpha=xsTable[:, 4], flagBarns=True)<jupyter_output><empty_output><jupyter_text>Depletion Stage  Execute Depletion<jupyter_code>dep = MainDepletion([0.0, 5.5, 30.0], bootstrap, tempramp, fullthrust)
# define metadata (steps, flux, and so on)
power = 1E+6*np.array([16.545, 118.49, 272.52, 330.22, 272.52, 214.82, 118.49])
dt = np.array([5.5, 24.5, 7., 1800., 7., 180., 40.])
dep.SetDepScenario(power=power, timeUnits="seconds", timesteps=dt)
# set initial composition
dep.SetInitialComposition(xsTable[:, 0], xsTable[:, 1], vol=VOL)
# solve the Bateman equations
dep.SolveDepletion(method="cram", xsinterp=False)<jupyter_output><empty_output><jupyter_text>Post depletion analysis<jupyter_code>dep.DecayHeat()
dep.Radiotoxicity()
dep.Activity()
dep.Mass()<jupyter_output><empty_output><jupyter_text>Post-process results<jupyter_code>res = Results(dep)
res.getvalues("totalQt")
res.plot("Nt", timeUnits="seconds", markers=['--^', '--o'], isotopes=[531350, 922360],
         ylabel="Atomic density, #/b/cm")
res.plot("totalQt", timeUnits="seconds", norm=1E+6, ylabel="Decay heat, MW")
res.plot("flux", timeUnits="seconds", markers="--ro", ylabel="Flux, n/cm2/s")<jupyter_output><empty_output>
{'fraction_non_alphanumeric': 0.15316760224538895, 'fraction_numerical': 0.04210104250200481, 'mean_word_length': 6.447761194029851, 'pattern_counts': {'":': 0, '<': 22, '<?xml version=': 0, '>': 22, 'https://': 0, 'lorem ipsum': 0, 'www.': 0, 'xml': 0}, 'pii_count': 0, 'substrings_counts': 0, 'word_list_counts': {'cursed_substrings.json': 0, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
{'chain_length': '5', 'dir': 'jupyter-structured-clean-dedup', 'id': '86575', 'license': 'MIT', 'n_tokens_mistral': 988, 'n_tokens_neox': 904, 'n_words': 175, 'path': 'examples/multiple_xs_states.ipynb', 'repo_name': 'MattKrecicki/PYTHON-ISOTOPIC-DEPLETION-PACKAGE', 'size': '2499'}