diff --git a/src/cmfd_data.F90 b/src/cmfd_data.F90
index a2d27c3a0f..e84801aa43 100644
--- a/src/cmfd_data.F90
+++ b/src/cmfd_data.F90
@@ -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) 
diff --git a/src/cmfd_execute.F90 b/src/cmfd_execute.F90
index a1b8853042..f89f8c3c46 100644
--- a/src/cmfd_execute.F90
+++ b/src/cmfd_execute.F90
@@ -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) 
diff --git a/src/cmfd_input.F90 b/src/cmfd_input.F90
index 3808f2bdf2..3309d2d130 100644
--- a/src/cmfd_input.F90
+++ b/src/cmfd_input.F90
@@ -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)