anchor.diff:          end do
       end do
   
  -    ! Generate new M containing only chemically-significant eigenmodes via
  -    ! multiplication of the above matrices
  -    do i = 1, nRows
  -        do j = 1, nRows
  -            M(i,j) = sum(X(i,:) * (V * Xinv(:,j)))
  -        end do
  -    end do
  -
       ! Generate set of initial condition vectors (one per isomer and reactant)
       allocate( eqDist(1:nRows, 1:nIsom+nReac) )
       do j = 1, nIsom + nReac
  
positive.code: !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
  ! 
  !   MEMURN - Master Equation Model for Unimolecular Reaction Networks
  ! 
  !   Copyright (c) 2009 by Josh Allen (chemejosh@gmail.com)
  ! 
  !   Permission is hereby granted, free of charge, to any person obtaining a
  !   copy of this software and associated documentation files (the 'Software'),
  !   to deal in the Software without restriction, including without limitation
  !   the rights to use, copy, modify, merge, publish, distribute, sublicense,
  !   and/or sell copies of the Software, and to permit persons to whom the
  !   Software is furnished to do so, subject to the following conditions:
  ! 
  !   The above copyright notice and this permission notice shall be included in
  !   all copies or substantial portions of the Software.
  ! 
  !   THE SOFTWARE IS PROVIDED 'AS IS', WITHOUT WARRANTY OF ANY KIND, EXPRESS OR
  !   IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY,
  !   FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO EVENT SHALL THE
  
negative.code:     integer, dimension(1:nGrains,1:nIsom) :: indices
  
      ! Intermediate matrices
      real(8), dimension(:,:), allocatable :: Mcoll, Mrxn, M, X, Xinv, eqDist, dXij
      real(8), dimension(:), allocatable :: S, Sinv, V0, V, C
  
      ! BLAS and LAPACK temporary variables
      real(8), dimension(:), allocatable :: work
      integer info
  
      ! Indices
      integer i, j, n, r, h, index, count
  
      ! Construct accounting matrix
  	! Row is grain number, column is well number, value is index into matrix
  