CTITLEGATHER -- GATHERS A VECTOR FROM A SOURCE VECTOR                   00010000
C***********************************************************************
C                 COPYRIGHT ATLANTIC RICHFIELD COMPANY 1991            *
C***********************************************************************
CA    AUTHOR        MIKE MAROLDA                                        00020000
CA    DESIGNER      MIKE MAROLDA                                        00030000
CA    LANGUAGE      FORTRAN                                             00040000
CA    SYSTEM        IBM (SEE CRAY SYSTEM LIBRARY)                       00050000
CA    WRITTEN       01-24-86                                            00060000
C     REVISED       XX-XX-XX III.  ...                                  00070000
C     REVISED       07-28-86 JMP.  RECOMPILE FOR VECTORIZATION.         00070100
CA                                                                      00080000
CA                                                                      00090000
CA    CALL GATHER (N, A, B, INDEX)                                      00100000
CA         INPUT      N = NUMBER OF ELEMENTS IN A AND INDEX         I4  00110000
CA         OUTPUT     A = REAL ARRAY OF LENGTH N                    R4  00120000
CA         INPUT      B = REAL ARRAY                                R4  00130000
CA         INPUT  INDEX = INDEX VECTOR INTO B OF LENGTH N           I4  00140000
CA                                                                      00150000
CA                                                                      00160000
CA       GATHER OUTPUTS VECTOR A WHICH IS GATHER FROM A(INDEX(I))       00170000
CAEND                                                                   00180000
      SUBROUTINE GATHER ( N, A, B, INDEX )                              00190000
C                                                                       00200000
      DIMENSION  A   (1)                                                00210000
      DIMENSION  B   (1)                                                00220000
      DIMENSION  INDEX (1)                                              00230000
C                                                                       00240000
      DO 100 J=1,N                                                      00250000
      A(J) = B(INDEX(J))                                                00260000
  100 CONTINUE                                                          00270000
C                                                                       00280000
      RETURN                                                            00290000
C                                                                       00300000
      END                                                               00310000