anchor.diff:      step = step + 1
       t1 = t2
       t2 = REAL(i)*tstep1
  +
  +    IF (i > 1) THEN
  +       omeg = LOG(f0 / ft) / tstep1
  +    ELSE
  +       omeg = 0.
  +    END IF
   
       ! Rod bank changes
       DO n = 1, nb
  
positive.code:   tbeta = tbeta + iBeta(j)
  END DO
  
  ! ReCalculate forward flux
  ! CALL outer4(0)
  
  ! Calculate reactivity
  CALL react(af, sigr, rho)
  
  
negative.code:     step = step + 1
      t1 = t2
      t2 = tdiv + REAL(i)*tstep2
  
      omeg = LOG(f0 / ft) / tstep2
  
      ! Rod bank changes
      DO n = 1, nb
  