anchor.diff:  	 &     LINDEREACTIONRATEARRAY, LINDEREACTIONARRAY
   
         INTEGER  INFO(30), LIW, LRW, IWORK(41 + NSTATE), IPAR(1),IDID, iter, &
  -     &     IMPSPECIES, conc, ITER_OUTPT, IWORK_OUTPT(41 + NSTATE), COPYQ
  +     &     IMPSPECIES, conc, ITER_OUTPT, IWORK_OUTPT(41 + NSTATE)
         DOUBLE PRECISION Y(NSTATE), YPRIME(NSTATE), Time, TOUT, RTOL, ATOL &
        &     , RWORK(51 + 9*NSTATE + NSTATE**2), TARGETCONC, Y_OUTPT(NSTATE) &
        &     , YPRIME_OUTPT(NSTATE), TIME_OUTPT  &
        &     , RWORK_OUTPT(51 + 9*NSTATE + NSTATE**2)
  -     INTEGER PRUNEVEC(ESPECIES),PRUNEVEC_OUTPT(ESPECIES)
  +     INTEGER PRUNEVEC(ESPECIES)
   
         DOUBLE PRECISION RPAR(REACTIONSIZE+THIRDBODYREACTIONSIZE+ &
        &     TROEREACTIONSIZE+LINDEREACTIONSIZE),&
  
positive.code: 	! 4/24/08 gmagoon: call EdgeFlux if AUTOFLAG is 1 to
  	! determine EDGEFLAG 
  		IF (AUTOFLAG.eq.1) THEN
  			CALL EDGEFLUX(EDGEFLAG, Y, YPRIME,THRESH,ESPECIES, &
      		 &     EREACTIONSIZE,NEREAC,NEPROD,IDEREAC,IDEPROD,KVEC, &
  		 &     MAXRATIO,PRUNEVEC,NSTATE)
  		 	!update HIGHESTRATIO and store relevant variables that will be printed
  			!later; note that we cannot exit the do loop once the IF statement is
  			!caught because, even though output variables will be stored properly,
  			!the HIGHESTRATIO will not necessarily have the highest value
  			COPYQ=0
  			DO I=1, ESPECIES
  			    IF (MAXRATIO(I) .GT. HIGHESTRATIO) THEN
  				HIGHESTRATIO = MAXRATIO(I)
  				COPYQ=1 !mark this ODE time step for copying into output variables
  			    END IF
  			END DO
  			IF (COPYQ .EQ. 1) THEN
  			    ITER_OUTPT=ITER
  			    TIME_OUTPT=TIME
  			    Y_OUTPT = Y
  
negative.code:       !core inclusion threshhold, and we therefore want to report the
        !final state back to the Java code so it will recognize that the
        !target time/conversion has been reached
        !in either case, CORESTORED should be -1
        IF (CORESTORED .EQ. -1) THEN
  	ITER_OUTPT=ITER
  	TIME_OUTPT=TIME
  	Y_OUTPT = Y
  