added fix neutron balance option

This commit is contained in:
Bryan Herman 2014-03-02 19:42:15 -05:00
parent b5346c5045
commit b3236e434f

View file

@ -172,7 +172,7 @@ contains
cmfd_downscatter = .true.
end if
! Run an adjoint calc
! Reset dhat parameters
if (check_for_node(doc, "dhat_reset")) then
call get_node_value(doc, "dhat_reset", temp_str)
call lower_case(temp_str)
@ -180,6 +180,14 @@ contains
dhat_reset = .true.
end if
! Enforce perfect neutorn balance
if (check_for_node(doc, "balance")) then
call get_node_value(doc, "balance", temp_str)
call lower_case(temp_str)
if (trim(temp_str) == 'true' .or. trim(temp_str) == '1') &
cmfd_fix_balance = .true.
end if
! Set the solver type
if (check_for_node(doc, "solver")) &
call get_node_value(doc, "solver", cmfd_solver_type)