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: 		! Convert density of states from (cm^-1)^-1 to mol/J
  		densStates = densStates / conv
  
          deallocate(vib, rot, hind)
  
      end subroutine
  
      !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
  