anchor.diff:  			dE = dE0
           end if
   
  +        if (allocated(Elist)) deallocate(Elist)
           allocate(Elist(1:nGrains))
           do r = 1, nGrains
               Elist(r) = Emin + (r - 1) * dE
  
positive.code:         integer, intent(in) :: method
          real(8), dimension(1:nIsom+nReac+nProd,1:nIsom+nReac+nProd), intent(out) :: K
  
          real(8), dimension(1:nGrains) :: densStates0
          real(8), dimension(1:nIsom,1:nGrains) :: densStates
          real(8), dimension(1:nIsom+nReac+nProd) :: Eres
          real(8), dimension(1:nIsom) :: E0
  
negative.code:         integer index, i, r
  
          ! Initialize density of states arrays
          if (.not. allocated(isom%densStates)) then
              allocate(isom%densStates(1:nGrains))
          end if
          do r = 1, nGrains
              isom%densStates(r) = 0.0
              densStates(r) = 0.0
  