anchor.diff:  	& EREACTIONSIZE, LINDEMAX
         INTEGER NEREAC(EREACTIONSIZE),NEPROD(EREACTIONSIZE)
         INTEGER IDEREAC(EREACTIONSIZE,3), IDEPROD(EREACTIONSIZE,3)
  -      DOUBLE PRECISION KVEC(EREACTIONSIZE),MAXRATIO(ESPECIES)
  +      DOUBLE PRECISION KVEC(EREACTIONSIZE),MAXRATIO(ESPECIES), &
  +	& HIGHESTRATIO
   
         PARAMETER (SPMAX = 1500, RMAX = 100000, TBRMAX=100, TROEMAX=100, &
   	  &    LINDEMAX=100)
  
positive.code:      &     CURRENTReactionFlux(REACTIONSIZE+THIRDBODYREACTIONSIZE+ &
       &     TROEREACTIONSIZE+LINDEREACTIONSIZE),&
       &     PREVReactionFlux(REACTIONSIZE+THIRDBODYREACTIONSIZE+&
       &     TROEREACTIONSIZE+LINDEREACTIONSIZE), PREVTIME
  
        EXTERNAL RES, JAC, GETFLUX
      ! 5/13/08 gmagoon: vars for timing
  
negative.code: 	DOUBLE PRECISION, DIMENSION(:), ALLOCATABLE :: RATE
  	INTEGER NEREAC(EREACTIONSIZE),NEPROD(EREACTIONSIZE)
          INTEGER IDEREAC(EREACTIONSIZE,3), IDEPROD(EREACTIONSIZE,3)
          DOUBLE PRECISION KVEC(EREACTIONSIZE), RATIO, MAXRATIO(ESPECIES)
  	
  	ALLOCATE(RATE(ESPECIES))
  	! initialize array to have all zeroes
  