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:                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) +
  
negative.code:      $                 (TIME - PREVTIME)/2
                    PREVREACTIONFLUX(I) = CURRENTREACTIONFLUX(I)
                 END DO
  c               write(*,*) time
                 PREVTIME = TIME
                 if (time .le. 1e6) then
                    tout = 200*time
                 end if
                 go to 2
              END IF
           END IF
  