anchor.diff:        JSTEP=NT
         JDIF=NTV2
         S(JDIF)=SIN(THETA)
  -      DO 660 L=2,MT
  +      do 660 L=2,MT
         THETA=THETA/2.
         JSTEP2=JSTEP
         JSTEP=JDIF
  
positive.code:       SC= COS(DEL)
        SI=0.
        CO =1.0
        DO 50 I=1,NT
        K6=NTOT2-2*I+5
        CIRE= A(2*I-1) + A(K6)
        CIIM=A(2*I)-A(K6+1)
        CNIRE=(-SI*(A(2*I)+A(K6+1))+CO*(A(2*I-1)-A(K6)))
        IF(SI)62,61,62
     62 CNIIM=(A(2*I-1)-A(K6)-CO*CNIRE)/SI
        GO TO 63
     61 CNIIM=0.
     63 A(2*I-1)=CIRE
        A(2*I)=CIIM
  
negative.code:       N1=2**M1
        N2=2**M2
        N3=2**M3
        if (IFSET) 50,50,70
  50    NX=N1*N2*N3
        FN=NX
        do 60 I=1,NX
        A(2*I-1)=A(2*I-1)/FN
  60    A(2*I)=-A(2*I)/FN
  70    NP(1)=N1*2
        NP(2)=NP(1)*N2
        NP(3)=NP(2)*N3
        do 330 ID=1,3
        IL=NP(3)-NP(ID)
        IL1=IL+1
        MI=M(ID)
        if (MI) 330,330,80
  80    IDIF=NP(ID)
        KBIT=NP(ID)
        MEV=2*(MI/2)
        if (MI-MEV) 120,120,90
  90    KBIT=KBIT/2
        KL=KBIT-2
        do 100 I=1,IL1,IDIF
        KLAST=KL+I
        do 100 K=I,KLAST,2
        KD=K+KBIT
        T=A(KD)
        A(KD)=A(K)-T
  