anchor.diff:          real(8), dimension(:), allocatable, intent(inout) :: Elist
   
           integer nE, done, maxIndex
  -        real(8) Emin, Emax0, Emax, dE, mult, maxValue, tol
  +        real(8) Emin, Emax0, Emax, dE, mult, maxValue, maxIsomerE0, maxReactionE0, tol
           type(Isomer) isom
   
           integer i, r
  
positive.code:         do i = 1, nIsom+nReac+nProd
              Eres(i) = isomer_getActiveSpaceEnergy(i, net%reactions)
          end do
  
          ! Isomerization, dissociation, and association microcanonical rate
          ! coefficients, respectively
          do r = 1, size(net%reactions)
  
negative.code:             
              call estimateratecoefficients_rs(T, P, Elist, Mcoll, densStates, E0, Eres, &
                  Kij, Fim, Gnj, dEdown, nIsom, nReac, nProd, nGrains, K, msg)
  
          end if
  
          ! Check for validity of phenomenological rate coefficients
  