anchor.diff:                TEMP=ALPHA*Z(I,K)+BETA*Z(I,KP)
                 Z(I,KP)=ALPHA*Z(I,KP)-BETA*Z(I,K)
     120         Z(I,K)=TEMP
  +              print *, '   k, alpha, beta, tot', k, alpha, beta, tot
             END IF
             K=K-1
             GOTO 100
  
positive.code:       KK=IACT(ICON)
        DO 90 I=1,N
     90 DXNEW(I)=A(I,KK)
        TOT=0.0d0
        K=N
    100 IF (K .GT. NACT) THEN
            SP=0.0d0
            SPABS=0.0d0
  
negative.code:   110     SPABS=SPABS+DABS(TEMP)
            ACCA=SPABS+0.1d0*DABS(SP)
            ACCB=SPABS+0.2d0*DABS(SP)
            print *, '   sp, spabs, acca, accb', sp, spabs, acca, accb
            IF ((SPABS .GE. ACCA) .OR. (ACCA .GE. ACCB)) SP=0.0D0
            IF (TOT .EQ. 0.0D0) THEN
                TOT=SP
                print *, '   simple tot = ', tot
            ELSE
                KP=K+1
                TEMP=DSQRT(SP*SP+TOT*TOT)
  