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:             end do
          end do
          if (invalidRate == 1) then
              write (1,fmt='(A)') 'ERROR: One ore more rate coefficients not properly estimated.'
              write (1,*) 'Temperature =', T, 'K, Pressure =', P, 'Pa, Rates ='
              do i = 1, nIsom+nReac+nProd
  
negative.code:         integer, intent(in) :: method
          real(8), dimension(1:nIsom+nReac+nProd,1:nIsom+nReac+nProd), intent(out) :: K
  
          real(8), dimension(1:nIsom,1:nGrains) :: densStates
          real(8), dimension(1:nIsom+nReac+nProd) :: Eres
          real(8), dimension(1:nIsom) :: E0
  