CTITLEARCONR -- CORRELATION 00010001 C*********************************************************************** C COPYRIGHT ATLANTIC RICHFIELD COMPANY 1991 * C*********************************************************************** CA AUTHOR L. B. LIN 00020000 CA DESIGNER L. B. LIN 00030000 CA LANGUAGE FORTRAN 00040000 CA SYSTEM IBM AND CRAY 00050001 CA WRITTEN 10/21/85 00060000 C REVISED 06-05-86 ESN. DERIVED FROM ARCON. 00070000 CA 00080000 CA 00090000 CA CALL ARCONR (F, NF, A, B, N, NRM, IR) 00100000 CA INPUT F = REAL FILTER ARRAY R4 00110000 CA INPUT NF = NUMBER OF ELEMENTS IN F I4 00120000 CA INPUT A = REAL DATA ARRAY R4 00130000 CA OUTPUT B = REAL DATA ARRAY R4 00140000 CA INPUT N = NUMBER OF ELEMENTS IN B (OUTPUT POINTS) I4 00150000 CA INPUT NRM = NORMALIZATION KEY, IF NON-ZERO THE OUTPUT I4 00160000 CA WILL BE DIVIDED BY THE NUMBER OF FILTER 00170000 CA POINTS 00180000 CA INPUT IR = PARAMETER TO INDICATE THAT F WILL I4 00190000 CA BE INDEXED IN REVERSE 00200000 CA 00210000 CA 00220000 CA ARCONR PERFORMS CORRELATION ON THE INPUT REAL DATA ARRAY A WITH00230000 CA THE REAL FILTER ARRAY F. THE N OUTPUT POINTS ARE STORED IN THE00240000 CA REAL ARRAY B. NORMALIZATION BY THE NUMBER OF FILTER POINTS IS 00250000 CA PERFORMED ON THE OUTPUT IF NRM IS NON-ZERO. 00260000 CA 00270000 CAEND 00280000 SUBROUTINE ARCONR (F, NF, A, B, N, NRM, IR) 00290000 C 00300000 REAL A (1) 00310000 REAL B (N) 00320000 REAL F (NF) 00330000 C 00340000 DO 20 I = 1, N 00350000 B(I) = 0.0 00360000 C 00370000 DO 10 J = 1, NF 00380000 10 B(I) = B(I) + F(NF+1-J) * A(I+J-1) 00390000 C 00400000 20 CONTINUE 00410000 C 00420000 IF (NRM .NE. 0) THEN 00430000 CONS = NF 00440000 CALL ARDVFC (B, B, CONS, N) 00450000 ENDIF 00460000 C 00470000 RETURN 00480000 C 00490000 END 00500000