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

Fixed entropy mesh width calculation

Browse source
This commit is contained in:
Nick Horelik 2014-10-06 10:33:22 -04:00
parent bae255680a
commit 2c433d3061
2 changed files with 10 additions and 4 deletions
Download patch file Download diff file Expand all files Collapse all files

8
src/eigenvalue.F90
View file

@ -536,12 +536,12 @@ contains
m % n_dimension = 3
allocate(m % dimension(3))
m % dimension = n
! determine width
m % width = (m % upper_right - m % lower_left) / m % dimension
end if
! allocate and determine width
allocate(m % width(3))
m % width = (m % upper_right - m % lower_left) / m % dimension
! allocate p
allocate(entropy_p(1, m % dimension(1), m % dimension(2), &
m % dimension(3)))

6
src/input_xml.F90
View file

@ -3135,6 +3135,12 @@ contains
call fatal_error()
end if
if (.not. allocated(entropy_mesh % dimension)) then
message = "No dimension specified on entropy mesh for " //
"meshlines on plot " // trim(to_str(pl % id))
call fatal_error()
end if
pl % meshlines_mesh => entropy_mesh
case ('tally')