anchor.diff:  
                   ! Add difference between isomer ground-state energy and highest
                   ! transition state or isomer energy
  -                maxValue = 0.0
  +                maxIsomerE0 = 0.0
  +                maxReactionE0 = 0.0
                   do r = 1, size(net%isomers)
  -                    if (net%isomers(r)%E0 > maxValue) maxValue = net%isomers(r)%E0
  +                    if (net%isomers(r)%E0 > maxIsomerE0) maxIsomerE0 = net%isomers(r)%E0
                   end do
                   do r = 1, size(net%reactions)
  -                    if (net%reactions(r)%E0 > maxValue) maxValue = net%reactions(r)%E0
  +                    if (net%reactions(r)%E0 > maxReactionE0) maxReactionE0 = net%reactions(r)%E0
                   end do
  -                Emax = Elist(r) + maxValue - Emax0
  +                Emax = Elist(r) + maxIsomerE0 + maxReactionE0 - Emin
   
                   ! Round Emax up to nearest integer
                   Emax = ceiling(Emax)
  
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:         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)
  