fixed determination of number of energy groups from cmfd input file

This commit is contained in:
Bryan Herman 2011-12-13 13:48:05 -05:00
parent 0e886bbf76
commit 4120d974ea

View file

@ -44,7 +44,7 @@ contains
! get number of energy groups
if (len_trim(mesh_ % energy) > 0) then
call split_string(mesh_ % energy, words, n_words)
ng = n_words
ng = n_words - 1
end if
cmfd % indices(4) = ng ! sets energy group dimension
@ -52,9 +52,11 @@ contains
cmfd % albedo = mesh_ % albedo
! get acceleration map
allocate(cmfd % coremap(cmfd % indices(1), cmfd % indices(2), &
& cmfd % indices(3)))
cmfd % coremap = reshape(mesh_ % map,(cmfd % indices(1:3)))
if (associated(mesh_ % map)) then
allocate(cmfd % coremap(cmfd % indices(1), cmfd % indices(2), &
& cmfd % indices(3)))
cmfd % coremap = reshape(mesh_ % map,(cmfd % indices(1:3)))
end if
! create tally objects
call create_cmfd_tally()