anchor.diff:  !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
  -! 
  +!
   !   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
  
positive.code: !   LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING
  !   FROM, OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER
  !   DEALINGS IN THE SOFTWARE.
  ! 
  !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
  
  subroutine estimateRateCoefficients_CSE(T, P, E, Mcoll0, E0, densStates, &
  
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
  