1 SUBROUTINE VRADF2 (MP,IDO,L1,CC,CH,MDIMC,WA1) 2C 3C VRFFTPK, VERSION 1, AUGUST 1985 4C 5 DIMENSION CH(MDIMC,IDO,2,L1) ,CC(MDIMC,IDO,L1,2) , 6 1 WA1(IDO) 7 DO 101 K=1,L1 8 DO 1001 M=1,MP 9 CH(M,1,1,K) = CC(M,1,K,1)+CC(M,1,K,2) 10 CH(M,IDO,2,K) = CC(M,1,K,1)-CC(M,1,K,2) 11 1001 CONTINUE 12 101 CONTINUE 13 IF (IDO-2) 107,105,102 14 102 IDP2 = IDO+2 15 DO 104 K=1,L1 16 DO 103 I=3,IDO,2 17 IC = IDP2-I 18 DO 1003 M=1,MP 19 CH(M,I,1,K) = CC(M,I,K,1)+(WA1(I-2)*CC(M,I,K,2)- 20 1 WA1(I-1)*CC(M,I-1,K,2)) 21 CH(M,IC,2,K) = (WA1(I-2)*CC(M,I,K,2)-WA1(I-1)* 22 1 CC(M,I-1,K,2))-CC(M,I,K,1) 23 CH(M,I-1,1,K) = CC(M,I-1,K,1)+(WA1(I-2)*CC(M,I-1,K,2)+ 24 1 WA1(I-1)*CC(M,I,K,2)) 25 CH(M,IC-1,2,K) = CC(M,I-1,K,1)-(WA1(I-2)*CC(M,I-1,K,2)+ 26 1 WA1(I-1)*CC(M,I,K,2)) 27 1003 CONTINUE 28 103 CONTINUE 29 104 CONTINUE 30 IF (MOD(IDO,2) .EQ. 1) RETURN 31 105 DO 106 K=1,L1 32 DO 1006 M=1,MP 33 CH(M,1,2,K) = -CC(M,IDO,K,2) 34 CH(M,IDO,1,K) = CC(M,IDO,K,1) 35 1006 CONTINUE 36 106 CONTINUE 37 107 RETURN 38 END 39