anchor.diff:      &     YPRIME_OUTPT(NSTATE), TIME_OUTPT,  &
       &     RWORK_OUTPT(51 + 9*NEQ + NSTATE**2)
   
  -      INTEGER PRUNEVEC(ESPECIES),PRUNEVEC_OUTPT(ESPECIES)
  +      INTEGER PRUNEVEC(ESPECIES)
   
         DOUBLE PRECISION RPAR(REACTIONSIZE+THIRDBODYREACTIONSIZE+ &
        &     TROEREACTIONSIZE+LINDEREACTIONSIZE+NSTATE-1), del(neq), CJ, &
  
positive.code: 	! the threshhold and update the MAXRATIO vector
  	FLOOP: DO I=1, ESPECIES
  		RATIO = RATE(I)/FLUXRC
  		!if we reach the termination tolerance,
  		!set EDGEFLAG to a positive value
  		!(in particular, the species ID)
  		IF(RATIO .GE. THRESH) THEN
  		   EDGEFLAG = I
  		END IF
  		!update the MAXRATIO vector
  		IF(RATIO .GT. MAXRATIO(I)) THEN
  
negative.code: 			    !IWORK_OUTPT = IWORK
  			    !RWORK_OUTPT = RWORK
  			    TOTALREACTIONFLUX_OUTPT = TOTALREACTIONFLUX
  			    CORESTORED=1
  			END IF
  		END IF
                 
  