anchor.diff:  C     starting a new iteration (or integration) step.
   C     'resid_tol' is the steady-state convergence tolerance:
   C     steady state when max(|dT/dt|,|dC/dt|)<resid_tol
  -C     This should give you an idea of where to insert the common block 
  -C     in your CFD solver 
  +C     This should give you an idea of where to insert the common block
  +C     in your CFD solver
   C     resid_tol should ONLY be changed by adapchem routine
   C *********************************************************************
   C
  
positive.code:                CTB(IFAL(N)) = 1.0
              ELSE
                 PR = RKLOW * C(KFAL(N)) / RKFT(IFAL(N))
  	       ispecies(kfal(n)) = 1
              ENDIF
  C
              PCOR = PR / (1.0 + PR)
  
negative.code: C
        DO 75 N = 1, NRLT
           I = IRLT(N)
         if(nur(i).eq.0) goto 75
           TFAC = RPLT(1,N)/T**(1.0/3.0) + RPLT(2,N)/T**(2.0/3.0)
           RKRT(I) = RKRT(I) * EXP(TFAC)
           EQK(I) = RKFT(I)/RKRT(I)
  