anchor.diff:  	!IWORK_OUTPT = IWORK
   	!RWORK_OUTPT = RWORK
   	TOTALREACTIONFLUX_OUTPT = TOTALREACTIONFLUX
  -	PRUNEVEC_OUTPT = PRUNEVEC
         END IF
         IWORK_OUTPT = IWORK
         RWORK_OUTPT = RWORK
  
positive.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
        IF ((AUTOFLAG .NE. 1) .OR. (HIGHESTRATIO .LT. CORETHRESH)) THEN
  	ITER_OUTPT=ITER
  	TIME_OUTPT=TIME
  	Y_OUTPT = Y
  
negative.code:             
  ! 6/26/08 gmagoon: adding subroutines EDGEFLUX and RCHAR
  
  ! EDGEFLUX determines whether core or termination threshhold flux has been
  ! reached by any species; if so, EDGEFLAGCORE and EDGEFLAG, respectively, are set to a
  ! positive value corresponding to the offending species ID;
  ! otherwise, these values remain at whatever they were when passed to the function
  ! EDGEFLUX also calls RCHAR
  ! when edgeflux exceeds core threshold, edgeflagcore is set to the offending species ID
  ! if edgeflux also exceeds the termination threshold flux, edgeflag is set to the offending species ID
  
  SUBROUTINE EDGEFLUX(EDGEFLAG,EDGEFLAGCORE,Y,YPRIME,THRESH,CORETHRESH,ESPECIES, &
      		 &     EREACTIONSIZE,NEREAC,NEPROD,IDEREAC,IDEPROD, &
  		&      KVEC, MAXRATIO,PRUNEVEC,NSTATE)
  	IMPLICIT NONE
  	INTEGER EDGEFLAG, EDGEFLAGCORE, ESPECIES, EREACTIONSIZE, NSTATE, I, J, K
  	DOUBLE PRECISION THRESH,CORETHRESH,FLUXRC, RFLUX, Y(NSTATE), YPRIME(NSTATE)
  	DOUBLE PRECISION, DIMENSION(:), ALLOCATABLE :: RATE
  	INTEGER NEREAC(EREACTIONSIZE),NEPROD(EREACTIONSIZE)
          INTEGER IDEREAC(EREACTIONSIZE,3), IDEPROD(EREACTIONSIZE,3)
  