crp/OpenMC
1
0
Fork 0
mirror of https://github.com/openmc-dev/openmc.git synced 2026-07-29 06:35:48 -04:00
Code Issues Projects Releases Packages Wiki Activity Actions 5

Fix differential fission eout tallies

Browse source
This commit is contained in:
Sterling Harper 2016-01-23 22:13:41 -05:00
parent c7350a3da4
commit 5634993bf9
2 changed files with 10 additions and 6 deletions
Download patch file Download diff file Expand all files Collapse all files

8
src/tally.F90
View file

@ -1046,8 +1046,12 @@ contains
! determine score based on bank site weight and keff
score = keff * fission_bank(n_bank - p % n_bank + k) % wgt
! Add derivative information for differenetial tallies.
if (allocated(t % deriv)) score = score * t % deriv % flux_deriv
! Add derivative information for differenetial tallies. Note that the
! i_nuclide and atom_density arguments do not matter since this is an
! analog estimator.
if (allocated(t % deriv)) then
call apply_derivative_to_score(p, t, 0, ZERO, SCORE_NU_FISSION, score)
end if
! determine outgoing energy from fission bank
E_out = fission_bank(n_bank - p % n_bank + k) % E

8
tests/test_diff_tally/results_true.dat
View file

@ -83,15 +83,15 @@ d_material,d_nuclide,d_variable,score,mean,std. dev.
3,,density,nu-fission,2.40e-02,6.37e-02
3,,density,nu-fission,0.00e+00,0.00e+00
3,,density,nu-fission,0.00e+00,0.00e+00
1,,density,nu-fission,-2.85e-02,2.57e-03
1,,density,nu-fission,-6.16e-02,5.51e-03
1,,density,nu-fission,3.45e-03,8.06e-04
1,,density,nu-fission,7.16e-03,3.66e-03
1,,density,nu-fission,0.00e+00,0.00e+00
1,,density,nu-fission,0.00e+00,0.00e+00
1,O-16.71c,nuclide_density,nu-fission,1.70e-01,2.18e-01
1,O-16.71c,nuclide_density,nu-fission,6.30e-01,5.11e-01
1,O-16.71c,nuclide_density,nu-fission,0.00e+00,0.00e+00
1,O-16.71c,nuclide_density,nu-fission,0.00e+00,0.00e+00
1,U-235.71c,nuclide_density,nu-fission,-2.00e+02,1.68e+01
1,U-235.71c,nuclide_density,nu-fission,-4.76e+02,1.41e+01
1,U-235.71c,nuclide_density,nu-fission,1.35e+02,7.94e+00
1,U-235.71c,nuclide_density,nu-fission,2.47e+02,1.29e+01
1,U-235.71c,nuclide_density,nu-fission,0.00e+00,0.00e+00
1,U-235.71c,nuclide_density,nu-fission,0.00e+00,0.00e+00