anchor.diff:          end do
       end do
       
  +    ! Check: The eigenvector corresponding to the zero eigenvalue should
  +    ! be identical to eqRatios
  +    !do j = 1, nIsom+nProd
  +    !    if (abs(X(j,nIsom+nProd) - eqRatios(j)) / eqRatios(j) > 0.001) then
  +    !        return
  +    !    end if
  +    !end do
  +
       ! Invert the chemical eigenvector matrix
       Xinv = X
       allocate(iPiv(1:nIsom+nProd))
  
positive.code:         return
      elseif (info < 0) then
          msg = 'Illegal argument passed to DSYEV.'
          return
      end if
  
  
negative.code:     
      ! The indices of the master equation matrix
      integer, dimension(1:nGrains,1:nIsom) :: indices
  
      ! Intermediate matrices
      real(8), dimension(:,:), allocatable :: Mcoll, Mrxn, M, X, Xinv
      real(8), dimension(:), allocatable :: S, Sinv, V0, V
  