anchor.diff:  
           real(8) dE, dEdown
           character(len=128) :: msg
  -
  -        integer invalidRate
   
           integer i, j, r
   
  
positive.code:                 ! Determine phenomenological rate coefficients
                  call network_applyApproximateMethod(net, nIsom, nReac, nProd, Elist, nGrains, T, P, method, K(u,v,:,:))
  
  
              end do
  
          end do
  
      end subroutine
  
  
negative.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, j, u, v, w
  
          integer invalidRate, invalidRates
  
          do u = 1, nT
  
  