tboudreaux/GridFire
1
0
Fork 0
Code Issues Pull Requests Actions Packages Projects Releases Wiki Activity

refactor(misc): removed debugging output

Browse Source
This commit is contained in:
Emily Boudreaux
2025-11-18 08:25:17 -05:00
parent 05c49d8dcf
commit d5d767a21c
2 changed files with 1 additions and 18 deletions
Download Patch File Download Diff File Expand all files Collapse all files

17
utils/reaclib/format.py
View File

@@ -361,28 +361,11 @@ def determine_reaction_type(reactants: List[str],
ejectiles: List[str] = [] ejectiles: List[str] = []
debug = False
if reactants == ['b8'] and products == ['be8']:
debug = True
BETA_PLUS_THRESHOLD_ENERGY = 1.022 # MeV BETA_PLUS_THRESHOLD_ENERGY = 1.022 # MeV
# ----------------------------------------------------------------------- # -----------------------------------------------------------------------
# 1. Charged-lepton bookkeeping (|ΔZ| = 1) ------------------------------ # 1. Charged-lepton bookkeeping (|ΔZ| = 1) ------------------------------
# ----------------------------------------------------------------------- # -----------------------------------------------------------------------
if debug:
print("============")
print("Reactant Species: ", reactantSpecies)
print("Product Species: ", productSpecies)
print("Target Species: ", targetSpecies)
print("Residual Species: ", residualSpecies)
print("Nuclear Projectiles: ", nuclearProjectiles)
print("Nuclear Ejectiles: ", nuclearEjectiles)
print("aReact, aProd: ", aReact, aProd)
print("zReact, zProd: ", zReact, zProd)
print("nReact, nProd: ", nReact, nProd)
print("dA, dZ, dN: ", dA, dZ, dN)
print("qValue: ", qValue)
if dZ == -1 and chapter == 1: if dZ == -1 and chapter == 1:
ejectiles.append("e-") # β- decay ejectiles.append("e-") # β- decay

2
validation/pynucastro/GridFireEquiv/evolve.py
View File

@@ -50,8 +50,8 @@ from fourdst.atomic import species, Species
from datetime import datetime from datetime import datetime
from typing import List, Dict, Set, Tuple from typing import List, Dict, Set, Tuple
symbols : list[str] = ["H-1", "He-3", "He-4", "C-12", "N-14", "O-16", "Ne-20", "Mg-24"]
X : list[float] = [0.708, 2.94e-5, 0.276, 0.003, 0.0011, 9.62e-3, 1.62e-3, 5.16e-4] X : list[float] = [0.708, 2.94e-5, 0.276, 0.003, 0.0011, 9.62e-3, 1.62e-3, 5.16e-4]
symbols : list[str] = ["H-1", "He-3", "He-4", "C-12", "N-14", "O-16", "Ne-20", "Mg-24"]
comp = Composition() comp = Composition()