anchor.diff:        END
         SUBROUTINE VRADF5 (MP,IDO,L1,CC,CH,MDIMC,WA1,WA2,WA3,WA4)
         implicit double precision (a-h, o-z)
  -C
  +C     
   C     VRFFTPK, VERSION 1, AUGUST 1985
  -C
  +C     
         DIMENSION  CC(MDIMC,IDO,L1,5)    ,CH(MDIMC,IDO,5,L1)     ,
  -     1           WA1(IDO)     ,WA2(IDO)     ,WA3(IDO)     ,WA4(IDO)
  +     1     WA1(IDO)     ,WA2(IDO)     ,WA3(IDO)     ,WA4(IDO)
         ARG=2.d0*PIMACH(1.0d0)/5.d0
         TR11=COS(ARG)
         TI11=SIN(ARG)
  
positive.code:          DO 103 I=3,IDO,2
              IC = IDP2-I
              DO 1003 M=1,MP
              CH(M,I-1,1,K) = ((WA1(I-2)*CC(M,I-1,K,2)+WA1(I-1)*
       1       CC(M,I,K,2))+(WA3(I-2)*CC(M,I-1,K,4)+WA3(I-1)*
       1       CC(M,I,K,4)))+(CC(M,I-1,K,1)+(WA2(I-2)*CC(M,I-1,K,3)+
       1       WA2(I-1)*CC(M,I,K,3)))
              CH(M,IC-1,4,K) = (CC(M,I-1,K,1)+(WA2(I-2)*CC(M,I-1,K,3)+
       1       WA2(I-1)*CC(M,I,K,3)))-((WA1(I-2)*CC(M,I-1,K,2)+
       1       WA1(I-1)*CC(M,I,K,2))+(WA3(I-2)*CC(M,I-1,K,4)+
       1       WA3(I-1)*CC(M,I,K,4)))
              CH(M,I,1,K) = ((WA1(I-2)*CC(M,I,K,2)-WA1(I-1)*
       1       CC(M,I-1,K,2))+(WA3(I-2)*CC(M,I,K,4)-WA3(I-1)*
       1       CC(M,I-1,K,4)))+(CC(M,I,K,1)+(WA2(I-2)*CC(M,I,K,3)-
       1       WA2(I-1)*CC(M,I-1,K,3)))
              CH(M,IC,4,K) = ((WA1(I-2)*CC(M,I,K,2)-WA1(I-1)*
       1       CC(M,I-1,K,2))+(WA3(I-2)*CC(M,I,K,4)-WA3(I-1)*
       1       CC(M,I-1,K,4)))-(CC(M,I,K,1)+(WA2(I-2)*CC(M,I,K,3)-
       1       WA2(I-1)*CC(M,I-1,K,3)))
              CH(M,I-1,3,K) = ((WA1(I-2)*CC(M,I,K,2)-WA1(I-1)*
       1       CC(M,I-1,K,2))-(WA3(I-2)*CC(M,I,K,4)-WA3(I-1)*
       1       CC(M,I-1,K,4)))+(CC(M,I-1,K,1)-(WA2(I-2)*CC(M,I-1,K,3)+
       1       WA2(I-1)*CC(M,I,K,3)))
              CH(M,IC-1,2,K) = (CC(M,I-1,K,1)-(WA2(I-2)*CC(M,I-1,K,3)+
       1       WA2(I-1)*CC(M,I,K,3)))-((WA1(I-2)*CC(M,I,K,2)-WA1(I-1)*
       1       CC(M,I-1,K,2))-(WA3(I-2)*CC(M,I,K,4)-WA3(I-1)*
       1       CC(M,I-1,K,4)))
              CH(M,I,3,K) = ((WA3(I-2)*CC(M,I-1,K,4)+WA3(I-1)*
       1       CC(M,I,K,4))-(WA1(I-2)*CC(M,I-1,K,2)+WA1(I-1)*
       1       CC(M,I,K,2)))+(CC(M,I,K,1)-(WA2(I-2)*CC(M,I,K,3)-
       1       WA2(I-1)*CC(M,I-1,K,3)))
              CH(M,IC,2,K) = ((WA3(I-2)*CC(M,I-1,K,4)+WA3(I-1)*
       1       CC(M,I,K,4))-(WA1(I-2)*CC(M,I-1,K,2)+WA1(I-1)*
       1       CC(M,I,K,2)))-(CC(M,I,K,1)-(WA2(I-2)*CC(M,I,K,3)-WA2(I-1)*
       1       CC(M,I-1,K,3)))
   1003       CONTINUE
    103    CONTINUE
    104 CONTINUE
  
negative.code:       END
        SUBROUTINE VRADF3 (MP,IDO,L1,CC,CH,MDIMC,WA1,WA2)
        implicit double precision (a-h, o-z)
  C     
  C     VRFFTPK, VERSION 1, AUGUST 1985
  C     
        DIMENSION   CH(MDIMC,IDO,3,L1)  ,CC(MDIMC,IDO,L1,3)     ,
       1     WA1(IDO)     ,WA2(IDO)
        ARG=2.d0*PIMACH(1.0d0)/3.d0
        TAUR=COS(ARG)
        TAUI=SIN(ARG)
        DO 101 K=1,L1
           DO 1001 M=1,MP
              cc23 = (CC(M,1,K,2)+CC(M,1,K,3))
              CH(M,1,1,K) = CC(M,1,K,1)+cc23
              CH(M,1,3,K) = TAUI*(CC(M,1,K,3)-CC(M,1,K,2))
              CH(M,IDO,2,K) = CC(M,1,K,1)+TAUR*cc23
   1001    CONTINUE
   101  CONTINUE
        IF (IDO .EQ. 1) RETURN
        IDP2 = IDO+2
        DO K=1,L1
           DO I=3,IDO,2
              IC = IDP2-I
              DO M=1,MP
                 p1 = WA1(I-1)*CC(M,I,K,2) + WA1(I-2)*CC(M,I-1,K,2)
                 p2 = WA2(I-1)*CC(M,I,K,3) + WA2(I-2)*CC(M,I-1,K,3)
                 p3 = WA1(I-2)*CC(M,I,K,2) - WA1(I-1)*CC(M,I-1,K,2)
                 p4 = WA2(I-2)*CC(M,I,K,3) - WA2(I-1)*CC(M,I-1,K,3)
                 p1pp2 = p1 + p2
                 p3pp4 = p3 + p4
                 tp1pp2 = taur*p1pp2
                 tp3pp4 = TAUR*p3pp4
                 tp3mp4 = TAUI*(p3-p4)
                 tp2mp1 = TAUI*(p2-p1)
                 cc12 = CC(M,I-1,K,1)+tp1pp2
                 cc34 = CC(M,I,K,1)+tp3pp4
                 CH(M,I-1,1,K)  = CC(M,I-1,K,1)+p1pp2
                 CH(M,I,1,K)    = CC(M,I,K,1)+p3pp4
                 CH(M,I-1,3,K)  = cc12+tp3mp4
                 CH(M,IC-1,2,K) = cc12-tp3mp4
                 CH(M,I,3,K)    = cc34+tp2mp1
                 CH(M,IC,2,K)   = tp2mp1-cc34
              enddo
           enddo
        enddo
        RETURN
        END
        SUBROUTINE VRADF4 (MP,IDO,L1,CC,CH,MDIMC,WA1,WA2,WA3)
  