anchor.diff:        end if
   
   ! 6/26/08 gmagoon: call EdgeFlux if AUTOFLAG is 1 to
  -! determine EDGEFLAG (-1 if flux threshhold has not been met,
  +! determine EDGEFLAG and EDGEFLAGCORE (-1 if flux threshhold has not yet been met,
   ! positive integer otherwise)
   	EDGEFLAG = -1
  +	EDGEFLAGCORE = -1
  +	!initialize flag CORESTORED to -1 indicating that the output values (when core inclusion threshold is exceeded) have not been stored yet; this will be set to 1 once the values are stored
  +	CORESTORED = -1
   	!initialize MAXRATIO to a vector of all zeroes;
   	!this variable will track the maxiumum
   	!of rate(i)/Rchar for each edge species i, with respect to time
  
positive.code:     &     YPRIME_OUTPT(NSTATE), TIME_OUTPT,  &
      &     RWORK_OUTPT(51 + 9*NEQ + NSTATE**2)
  
        INTEGER PRUNEVEC(ESPECIES),PRUNEVEC_OUTPT(ESPECIES)
  
        DOUBLE PRECISION RPAR(REACTIONSIZE+THIRDBODYREACTIONSIZE+ &
       &     TROEREACTIONSIZE+LINDEREACTIONSIZE+NSTATE-1), del(neq), CJ, &
  
negative.code: 	! the threshhold and update the MAXRATIO vector
  	FLOOP: DO I=1, ESPECIES
  		RATIO = RATE(I)/FLUXRC
  		!if we reach the error tolerance (aka core tolerance),
  		!set EDGEFLAGCORE to the species ID
  		!if we also reach the termination tolerance, set EDGEFLAG to the species ID
  		IF (RATIO .GE. CORETHRESH) THEN
  		    EDGEFLAGCORE = I
  		    IF(RATIO .GE. THRESH) THEN
  			EDGEFLAG = I
  		    END IF
  		END IF
  		!update the MAXRATIO vector
  		IF(RATIO .GT. MAXRATIO(I)) THEN
  