anchor.diff:  REAL, DIMENSION(nnod, ng) :: ft, ftx1, fty1, ftz1, ftx2, fty2, ftz2  ! Flux at previous time step
   REAL, DIMENSION(nnod, ng) :: af                                      ! adjoint flux
   REAL, DIMENSION(nnod, ng) :: sigrp                                   ! Temporary sigr
  -    
  +
   REAL :: rho
   REAL :: t1, t2
   REAL :: tpow1, tpow2
  
positive.code: 
      WRITE(ounit,'(I4, F10.3, F10.4, ES15.4, 12F9.2)') step, t2, rho/tbeta, &
      tpow2/tpow1, (bpos(n), n = 1, nb)
      
      IF (maxi) THEN
          WRITE(*,'(I4, F10.3, F10.4, ES15.4, A35)') step, t2, rho/tbeta, &
          tpow2/tpow1, 'OUTER ITERATION DID NOT CONVERGE'
      ELSE
          WRITE(*,'(I4, F10.3, F10.4, ES15.4)') step, t2, rho/tbeta, &
          tpow2/tpow1
      END IF  
  
      IF (maxi) tranw = .TRUE.    
      
  END DO
  
  ! Second Time Step
  
negative.code:         WRITE(*,'(I4, F10.3, F10.4, ES15.4)') step, t2, rho/tbeta, &
          xppow
      END IF
  
      IF (maxi) tranw = .TRUE.
  
      IF (step>1000) THEN
  