anchor.diff:          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)
  
positive.code:         real(8), dimension(:), allocatable, intent(inout) :: Elist
  
          integer nE, done, maxIndex
          real(8) Emin, Emax0, Emax, dE, mult, maxValue, tol
          type(Isomer) isom
  
          integer i, r
  
negative.code:             end if
          end do
  
          ! Average energy transferred in a deactivating collision
          dEdown = net%bathGas%dEdown
  
  