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

Added more output on verbosity level 10 (for traces).

Browse source
This commit is contained in:
Paul Romano 2011-12-05 12:39:59 -05:00
parent f1f443dadc
commit f95cb888d7
2 changed files with 16 additions and 7 deletions
Download patch file Download diff file Expand all files Collapse all files

4
src/geometry.f90
View file

@ -441,7 +441,9 @@ contains
lat => lattices(p % coord % lattice)
if (verbosity >= 10 .or. trace) then
message = " Crossing lattice " // int_to_str(lat % id)
message = " Crossing lattice " // trim(int_to_str(lat % id)) // &
". Current position (" // trim(int_to_str(p % coord % lattice_x)) &
// "," // trim(int_to_str(p % coord % lattice_y)) // ")"
call write_message()
end if

19
src/physics.f90
View file

@ -14,7 +14,7 @@ module physics
use particle_header, only: Particle, LocalCoord
use random_lcg, only: prn
use search, only: binary_search
use string, only: int_to_str
use string, only: int_to_str, real_to_str
use tally, only: score_tally, score_surface_current
implicit none
@ -412,6 +412,13 @@ contains
! Sample nuclide/reaction for the material the particle is in
call sample_reaction(p, MT)
! Display information about collision
if (verbosity >= 10 .or. trace) then
message = " " // trim(reaction_name(MT)) // ". Energy = " // &
trim(real_to_str(p % E * 1e6_8)) // " eV."
call write_message()
end if
! check for very low energy
if (p % E < 1.0e-100_8) then
p % alive = .false.
@ -419,11 +426,6 @@ contains
call warning()
end if
! Score collision estimator tallies for any macro tallies -- this is done
! after a collision has occurred rather than before because we need
! information on the outgoing energy for any tallies with an outgoing energy
! filter
! Check if particle scattered or fissioned
if (survival_biasing) then
fissioned = .false.
@ -433,6 +435,11 @@ contains
scattered = is_scatter(MT)
end if
! Score collision estimator tallies for any macro tallies -- this is done
! after a collision has occurred rather than before because we need
! information on the outgoing energy for any tallies with an outgoing energy
! filter
if (tallies_on) then
call score_tally(p, scattered, fissioned)
call score_surface_current(p)