anchor.diff:          tpow2/tpow1
       END IF
   
  -    IF (maxi) tranw = .TRUE.    
  +    IF (maxi) tranw = .TRUE.
   
   END DO
   
  
positive.code: 
      ! Calculate reactivity
      CALL react(af, sigrp, rho)
      
      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'
  
negative.code: REAL :: tpow1, tpow2
  INTEGER :: n, i, j, g, imax, step
  
  LOGICAL :: maxi   ! Maximum Outer Iteration Reached?
  
  ! Allocate precusor density
  ALLOCATE (c0(nf,nnod),cx1(nf,nnod),cy1(nf,nnod),cz1(nf,nnod))
  