CTITLEARMPPW -- MULTIPLY COMPLEX ARRAY BY SQRT OF REAL ARRAY 00010001 C*********************************************************************** C COPYRIGHT ATLANTIC RICHFIELD COMPANY 1991 * C*********************************************************************** CA AUTHOR STU NELAN 00020000 CA DESIGNER STU NELAN 00030000 CA LANGUAGE FORTRAN 00040000 CA SYSTEM IBM AND CRAY 00050001 CA WRITTEN 05-30-86 00060000 C REVISED 07-11-86 ESN. FOR USE ON THE IBM. 00070001 CA 00080000 CA 00090000 CA CALL ARMPPW (Z, W, P, N) 00100000 CA INPUT Z = COMPLEX ARRAY C4 00110000 CA OUTPUT W = COMPLEX ARRAY C4 00120000 CA INPUT P = REAL ARRAY R4 00130000 CA INPUT N = NUMBER OF ELEMENTS IN EACH ARRAY Z, W, I4 00140000 CA AND P 00150000 CA 00160000 CA 00170000 CA ARMPPW MULTIPLIES THE COMPLEX ARRAY Z BY THE REAL ARRAY P 00180000 CA AND STORES THE RESULT IN THE COMPLEX ARRAY W. 00190000 CA IE. IF Z(J) = (X+IY) AND PS = SQRT ( P(J) ) 00200000 CA THEN W(J) = (PS*X+PS*IY) 00210000 CA Z AND W MAY BE THE SAME ARRAY. 00220000 CA 00230000 CA THE ARRAY P IS NORMALLY A POWER SPECTRUM 00240000 CAEND 00250000 SUBROUTINE ARMPPW (Z, W, P, N) 00260000 C 00270000 REAL P(N) 00280000 COMPLEX W(N) 00290000 COMPLEX Z(N) 00300000 C 00310000 DO 10 I = 1, N 00320000 10 W(I) = Z(I) * SQRT(P(I)) 00330000 C 00340000 RETURN 00350000 C 00360000 END 00370000