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 extra output about lattices upon errors.

Browse source
This commit is contained in:
Paul Romano 2011-12-01 13:05:31 -05:00
parent 392fecb7db
commit 86009101de
1 changed files with 9 additions and 7 deletions
Download patch file Download diff file Expand all files Collapse all files

16
src/geometry.f90
View file

@ -483,13 +483,13 @@ contains
type(Lattice), pointer :: lat
type(Universe), pointer :: univ
lat => lattices(p % lattice)
if (verbosity >= 10) then
message = " Crossing lattice"
message = " Crossing lattice " // int_to_str(lat % id)
call write_message()
end if
lat => lattices(p % lattice)
u = p % uvw(1)
v = p % uvw(2)
@ -550,12 +550,14 @@ contains
i_x = p % index_x
i_y = p % index_y
if (i_x < 1 .or. i_x > lat % n_x) then
message = "Reached edge of lattice (" // trim(int_to_str(i_x)) // &
"," // trim(int_to_str(i_y)) // ")."
message = "Reached edge of lattice " // trim(int_to_str(lat % id)) // &
" at position (" // trim(int_to_str(i_x)) // "," // &
trim(int_to_str(i_y)) // ")."
call fatal_error()
elseif (i_y < 1 .or. i_y > lat % n_y) then
message = "Reached edge of lattice (" // trim(int_to_str(i_x)) // &
"," // trim(int_to_str(i_y)) // ")."
message = "Reached edge of lattice " // trim(int_to_str(lat % id)) // &
" at position (" // trim(int_to_str(i_x)) // "," // &
trim(int_to_str(i_y)) // ")."
call fatal_error()
end if