anchor.diff:                 PREVREACTIONFLUX(I) = CURRENTREACTIONFLUX(I)
               END DO
               PREVTIME = TIME
  +            if (conc .eq. 1) then
  +               tout = 200*time
  +            end if
   c            write(*,*) time, y(impspecies), targetConc*y(nstate)
               GO TO 1
            END IF
  
positive.code: 
  
        INTEGER  INFO(30), LIW, LRW, IWORK(*), IPAR(1),IDID, iter,
       $     IMPSPECIES
        DOUBLE PRECISION Y(NSTATE), YPRIME(NSTATE), Time, TOUT, RTOL, ATOL
       $     , RWORK(*), TARGETCONC
  
  
negative.code: 
        EXTERNAL RES, JAC, GETFLUX
  
        conc = 1
        IPAR(1)= NSTATE
  
        LIW = 41 + NSTATE
  