Fix cell temperature bug for multipole. Results change a bit.

This commit is contained in:
Paul Romano 2016-08-30 21:49:22 -05:00
parent a95266ab22
commit fbebf7bf70
2 changed files with 5 additions and 5 deletions

View file

@ -334,11 +334,11 @@ contains
else
! Determine temperature
!!$ if (temperature_method == TEMPERATURE_MULTIPOLE) then
!!$ kT = p % sqrtkT**2
!!$ else
if (temperature_method == TEMPERATURE_MULTIPOLE) then
kT = p % sqrtkT**2
else
kT = nuc % kTs(micro_xs(i_nuclide) % index_temp)
!!$ end if
end if
! Perform collision physics for elastic scattering
call elastic_scatter(i_nuclide, nuc % reactions(1), kT, &

View file

@ -1,5 +1,5 @@
k-combined:
1.457760E+00 1.119656E-02
1.425673E+00 1.779969E-02
Cell
ID = 11
Name =