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

added zero flux trap

Browse source
This commit is contained in:
Bryan Herman 2012-01-14 16:54:19 -05:00
parent 9d3689c650
commit 5c5d120220
1 changed files with 15 additions and 0 deletions
Download patch file Download diff file Expand all files Collapse all files

15
src/cmfd_execute.F90
View file

@ -129,6 +129,21 @@ contains
flux = t % scores(score_index,1) % val
cmfd % flux(h,i,j,k) = flux
! detect zero flux
if ((flux - 0.0D0) < 1.0D-10) then
if (.not. allocated(cmfd%coremap)) then
write(*,*) 'Fatal: detected zero flux without coremap'
stop
else
write(*,*) 'Warning: detected zero flux at:',i,j,k
flux = 99999.0D0
if (.not. cmfd%coremap(i,j,k) == 99999) then
write(*,*) 'Fatal: need to check core map with zero flux'
stop
end if
end if
end if
! get total rr and convert to total xs
cmfd % totalxs(h,i,j,k) = t % scores(score_index,2) % val / flux