added option to reset dhats

This commit is contained in:
Bryan Herman 2014-03-02 19:39:16 -05:00
parent f1436d942f
commit b5346c5045
3 changed files with 9 additions and 7 deletions

View file

@ -309,8 +309,6 @@ contains
if (associated(m)) nullify(m)
#ifdef CMFD_DEBUG
open(file='cmfd_src_' // trim(to_str(current_batch)) // '.dat', unit=100)
open(file='source_bank_' // trim(to_str(current_batch)) // '.dat', unit=101)
open(file='openmc_src_' // trim(to_str(current_batch)) // '.dat', unit=102)
open(file='totalxs1_' // trim(to_str(current_batch)) // '.dat', unit=103)
open(file='totalxs2_' // trim(to_str(current_batch)) // '.dat', unit=104)
@ -347,8 +345,6 @@ contains
((cmfd % current(4*l-2,2,i,j,1) - &
cmfd % current(4*l-3,2,i,j,1)))
end do
write(100,*) cmfd % cmfd_src(1,i,j,1)
write(101,*) cmfd % sourcecounts(1,i,j,1)
write(102,*) cmfd % openmc_src(1,i,j,1)
write(103,*) cmfd % totalxs(1,i,j,1)
write(104,*) cmfd % totalxs(2,i,j,1)
@ -370,8 +366,6 @@ contains
write(120,*) leak2
end do
end do
close(100)
close(101)
close(102)
close(103)
close(104)

View file

@ -243,7 +243,7 @@ contains
open(file='cmfd_src_' // trim(to_str(current_batch)) // '.dat', unit=100)
do i = 1, nx
do j = 1,ny
write(100,*) cmfd % sourcecounts(1,i,j,1)
write(100,*) cmfd % cmfd_src(1,i,j,1)
end do
end do
close(100)

View file

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