anchor.diff:       &     CURRENTReactionFlux(REACTIONSIZE+THIRDBODYREACTIONSIZE+ &
        &     TROEREACTIONSIZE+LINDEREACTIONSIZE),&
        &     PREVReactionFlux(REACTIONSIZE+THIRDBODYREACTIONSIZE+&
  -     &     TROEREACTIONSIZE+LINDEREACTIONSIZE), PREVTIME
  +     &     TROEREACTIONSIZE+LINDEREACTIONSIZE), PREVTIME, &
  +     &     TOTALREACTIONFLUX_OUTPT(REACTIONSIZE+THIRDBODYREACTIONSIZE+ &
  +     &     TROEREACTIONSIZE+LINDEREACTIONSIZE)
   
         EXTERNAL RES, JAC, GETFLUX
       ! 5/13/08 gmagoon: vars for timing
  
positive.code: 	DOUBLE PRECISION, DIMENSION(:), ALLOCATABLE :: RATE
  	INTEGER NEREAC(EREACTIONSIZE),NEPROD(EREACTIONSIZE)
          INTEGER IDEREAC(EREACTIONSIZE,3), IDEPROD(EREACTIONSIZE,3)
          DOUBLE PRECISION KVEC(EREACTIONSIZE), MAXRATIO(ESPECIES)
  	
  	ALLOCATE(RATE(ESPECIES))
  	! initialize array to have all zeroes
  
negative.code:      &     TROEREACTIONARRAY, ConstantConcentration, &
  	 &     LINDEREACTIONRATEARRAY, LINDEREACTIONARRAY
  
        INTEGER  INFO(30), LIW, LRW, IWORK(1541), IPAR(1),IDID, iter, &
       &     IMPSPECIES, conc, ITER_OUTPT, IWORK_OUTPT(1541)
        DOUBLE PRECISION Y(NSTATE), YPRIME(NSTATE), Time, TOUT, RTOL, ATOL &
       &     , RWORK(51+9*SPMAX+SPMAX**2), TARGETCONC, Y_OUTPT(NSTATE) &
       &     , YPRIME_OUTPT(NSTATE), TIME_OUTPT  &
       &     , RWORK_OUTPT(51+9*SPMAX+SPMAX**2)
  
        DOUBLE PRECISION RPAR(REACTIONSIZE+THIRDBODYREACTIONSIZE+ &
       &     TROEREACTIONSIZE+LINDEREACTIONSIZE),&
  