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

Got cmfd working, not sure about with feedback just yet

Browse source
This commit is contained in:
Adam Nelson 2015-11-07 14:23:20 -05:00
parent 39445f059f
commit c571a4b126
3 changed files with 17 additions and 8 deletions
Download patch file Download diff file Expand all files Collapse all files

6
examples/xml/multigroup/c5g7/2d/cmfd.xml
View file

@ -1,9 +1,9 @@
<?xml version="1.0"?>
<cmfd>
<begin>1</begin>
<feedback>true</feedback>
<solver>jfnk</solver>
<begin>2</begin>
<feedback>false</feedback>
<solver>power</solver>
<mesh>
<lower_left> 0.0 0.0 -100.0</lower_left>
<upper_right> 64.26 64.26 100.0</upper_right>

4
examples/xml/multigroup/c5g7/2d/settings.xml
View file

@ -33,7 +33,7 @@
<!-- Establish statepoints to aid in examining convergence -->
<state_point>
<interval> 2000 </interval>
<interval> 2500 </interval>
<source_write>true</source_write>
</state_point>
@ -43,7 +43,7 @@
<tallies>true</tallies>
</output>
<run_cmfd>false</run_cmfd>
<run_cmfd>true</run_cmfd>
<cross_sections>../data.xml</cross_sections>

15
src/tally.F90
View file

@ -966,8 +966,11 @@ contains
! For scattering production, we need to use the pre-collision
! weight times the multiplicity as the estimate for the number of
! neutrons exiting a reaction with neutrons in the exit channel
score = p % wgt * nuc % get_xs(p % g, 'f_mu', p % last_g, &
score = p % wgt
if (i_nuclide > 0) then
score = score * nuc % get_xs(p % g, 'f_mu', p % last_g, &
p % last_uvw, p % mu)
end if
case (SCORE_NU_SCATTER_PN)
@ -980,8 +983,11 @@ contains
! For scattering production, we need to use the pre-collision
! weight times the multiplicity as the estimate for the number of
! neutrons exiting a reaction with neutrons in the exit channel
score = p % wgt * nuc % get_xs(p % g, 'f_mu', p % last_g, &
score = p % wgt
if (i_nuclide > 0) then
score = score * nuc % get_xs(p % g, 'f_mu', p % last_g, &
p % last_uvw, p % mu)
end if
case (SCORE_NU_SCATTER_YN)
@ -994,8 +1000,11 @@ contains
! For scattering production, we need to use the pre-collision
! weight times the multiplicity as the estimate for the number of
! neutrons exiting a reaction with neutrons in the exit channel
score = p % wgt * nuc % get_xs(p % g, 'f_mu', p % last_g, &
score = p % wgt
if (i_nuclide > 0) then
score = score * nuc % get_xs(p % g, 'f_mu', p % last_g, &
p % last_uvw, p % mu)
end if
case (SCORE_TRANSPORT)