anchor.diff:  
           real(8) dE, dEdown
           character(len=128) :: msg
  -
  -        integer invalidRate
   
           integer i, j, r
   
  
positive.code:         real(8), dimension(1:nT,1:nP,1:nIsom+nReac+nProd,1:nIsom+nReac+nProd), intent(out) :: K
  
          real(8) T, P
          integer i, u, v, w
  
          do u = 1, nT
  
  
negative.code:                 Kij, Fim, Gnj, dEdown, nIsom, nReac, nProd, nGrains, K, msg)
  
          end if
  
          if (msg /= '') then
              write (*,fmt='(A)') 'ERROR: One or more rate coefficients not properly estimated. See fame.log for details.'
              write (*,fmt='(A)') 'The message returned was:', msg
              write (1,fmt='(A)') 'ERROR: One or more rate coefficients not properly estimated.'
  