crp/OpenMC
1
0
Fork 0
mirror of https://github.com/openmc-dev/openmc.git synced 2026-07-27 13:45:36 -04:00
Code Issues Projects Releases Packages Wiki Activity Actions

disallow negative resonance scattering energy bounds

Browse source
This commit is contained in:
Jon Walsh 2017-04-04 07:46:39 -07:00
parent 7f588e878b
commit 008f2ec220
1 changed files with 3 additions and 3 deletions
Download patch file Download diff file Expand all files Collapse all files

6
src/physics.F90
View file

@ -875,9 +875,9 @@ contains
wgt = wcf * wgt
case (RES_SCAT_DBRC, RES_SCAT_ARES)
E_red = sqrt((awr * E) / kT)
E_low = (((E_red - FOUR)**2) * kT) / awr
E_up = (((E_red + FOUR)**2) * kT) / awr
E_red = sqrt(awr * E / kT)
E_low = max(ZERO, E_red - FOUR)**2 * kT / awr
E_up = (E_red + FOUR)**2 * kT / awr
! find lower and upper energy bound indices
! lower index