Replica Fields

If option R is chosen, OASN writes the replica field toa binary, sequential file with the name input.rpo . The file is opened using the following statements

C        *****  OPEN RPO FILE...Note that the logical unit 14 must 
C        be assigned a filename external to the program
C        In Unix: setenv FOR014 input.rpo.
C
         LUN=14
         CALL GETENV('FOR014',RPOFILE)
         OPEN  (UNIT       =  LUN        ,
     -          FILE       =  RPOFILE    ,
     -          STATUS     = 'UNKNOWN'   ,
     -          FORM       = 'UNFORMATTED')

The rpo file will first have a header for identification in terms of title, frequency sampling, array geometry and replica scanning space. The header is written by the following code, and should clearly be read accordingly. The parameters are defined in Tables 9 and 10.

        :
      CHARACTER*80  TITLE
        :
        :
c *** WRITE HEADER
       WRITE (LUN) TITLE
       WRITE (LUN) NRCV, NFREQ
c >>> DELFRQ is the frequency increment (FREQ2 - FREQ1)/(NFREQ-1)
       WRITE (LUN) FREQ1, FREQ2, DELFRQ
c >>> Replica sampling
       WRITE (LUN) ZMINR, ZMAXR, NZR
       WRITE (LUN) XMINR, XMAXR, NXR
       WRITE (LUN) YMINR, YMAXR, NYR
c >>> Array element data
       DO 10 IRCV=1,NRCV
        WRITE(LUN) X(IRCV),Y(IRCV),Z(IRCV),ITYP(IRCV),GAIN(IRCV)
 10    CONTINUE

The complex replicas follow frequency by frequency, written with the following loop structure

        :
      COMPLEX REPLIC(NRCV,NYR,NXR,NZR,NFREQ)
        :
        :
       DO 20 IFREQ=1,NFREQ
        DO 20 IZR=1,NZR
         DO 20 IXR=1,NXR
          DO 20 IYR=1,NYR
           DO 20 IRCV=1,NRCV
            WRITE (LUN) REPLIC(IRCV,IYR,IXR,IZR,IFREQ)
20     CONTINUE