anchor.diff:                 GO TO 1
               else if (idid .eq. 1) then
                  iter = iter+1
  -               
  +
                  CALL REACTION_FLUX(Y, CURRENTREACTIONFLUX, RPAR)
                 DO I=1,REACTIONSIZE+THIRDBODYREACTIONSIZE+TROEREACTIONSIZE
                     TOTALREACTIONFLUX(I) = TOTALREACTIONFLUX(I) +
  
positive.code:      $                 (TIME - PREVTIME)/2
                    PREVREACTIONFLUX(I) = CURRENTREACTIONFLUX(I)
                 END DO
                 PREVTIME = TIME
                 tout = 200*time
                 go to 2
              END IF
           END IF
  
negative.code:       ELSE
   2       IF (Y(IMPSPECIES) .GE. TARGETCONC*Y(NSTATE))
       $        THEN
  
              CALL DDASPK(RES, NEQ, Time, Y, YPRIME, TOUT, INFO, RTOL, 
       $           Atol,IDID,RWORK, LRW, IWORK, LIW, RPAR, IPAR, JAC,PSOL,
       $           SENPAR, G_RES)
  