anchor.diff:  af = f0       ! Save adjoint flux to af
   
   ALLOCATE(beta(nf), C(nf))
  +ALLOCATE(xvdum(nnod,ng))
   
   END SUBROUTINE adj_calc
   
  
positive.code: END DO
  
  xA = Integrate(vdum2) / F2    ! Calculate neutron generation time
  
  ! Calculate reactivity (rho)
  vdum2 = 0.; vdum3 = 0.; vdum4 = 0.
  
negative.code: 
  REAL, DIMENSION(nnod,ng) ::  osigr, onuf
  REAL, DIMENSION(nnod,ng,ng) :: osigs
  
  REAL :: A, rho
  
  