crp/OpenMC
1
0
Fork 0
mirror of https://github.com/openmc-dev/openmc.git synced 2026-07-28 22:26:08 -04:00
Code Issues Projects Releases Packages Wiki Activity Actions 5

added residual output file

Browse source
This commit is contained in:
Bryan Herman 2012-04-10 09:52:31 -07:00
parent 5570281cee
commit 96557cf729
1 changed files with 9 additions and 2 deletions
Download patch file Download diff file Expand all files Collapse all files

11
src/cmfd_output.F90
View file

@ -74,8 +74,8 @@ contains
interactions = cmfd % totalxs(g,i,j,k) * cmfd % flux(g,i,j,k)
! get scattering and fission
scattering = 0.0
fission = 0.0
scattering = 0.0_8
fission = 0.0_8
GROUPH: do h = 1,ng
scattering = scattering + cmfd % scattxs(h,g,i,j,k) * &
@ -88,6 +88,13 @@ contains
! compute residual
res = leakage + interactions - scattering - (ONE/keff)*fission
write(670,*) 'Group:',g
write(670,*) 'Leakage:',leakage
write(670,*) 'Interactions:',interactions
write(670,*) 'Scattering:',scattering
write(670,*) 'Fission:',fission
write(670,*) 'k-eff:',keff
write(670,*) 'Balance:',res
! normalize by flux
res = res/cmfd%flux(g,i,j,k)