CTITLEMERXX  -- COMPUTES SYNTHETIC AND ERROR                            00000010
C***********************************************************************
C                 COPYRIGHT ATLANTIC RICHFIELD COMPANY 1991            *
C***********************************************************************
CA       AUTHOR                 D. D. THOMPSON                          00000020
CA       DESIGNER               D. D. THOMPSON                          00000030
CA       LANGUAGE               FORTRAN 77                              00000040
CA       SYSTEM                 IBM & CRAY                              00000041
CA       WRITTEN                1972                                    00000050
C        REVISED                MO-DA-YR                                00000060
C        REVISED                05-20-85 TWH. ADAPTED TO IBM & CRAY.    00000061
C                                                                       00000070
CA                                                                      00000080
CA                                                                      00000090
CA       CALL MERXX (E, XX, X, Z, KR, KW, KX, NZERO)                    00000100
CA             OUTPUT E      = RETURNED SUM OF SQUARES OF TRACE    R4   00000110
CA                             ERRORS                                   00000120
CA             OUTPUT XX     = SYNTHETIC TRACE                     R4   00000130
CA             INPUT  X      = TRACE INPUT                         R4   00000140
CA             INPUT  Z      = FIRST KR SAMPLES ARE INPUT          R4   00000150
CA                             REFLECTION COEFFICIENTS.                 00000160
CA                             Z(KR+1) . . . Z(KR+KW) IS INPUT          00000170
CA                             WAVELET                                  00000180
CA             INPUT  KR     = NUMBER OF REFLECTION COEFFICIENTS   I4   00000190
CA             INPUT  KW     = NUMBER OF WAVELET SAMPLES           I4   00000200
CA             INPUT  KX     = NUMBER OF TRACE SAMPLES             I4   00000210
CA             INPUT  NZERO  = INPUT FLAG                          I4   00000220
CA                             0 - ASSUME SYNTHETIC IS ZERO AND         00000230
CA                                 NEGLECT Z                            00000240
CA                             1 - COMPUTE SYNTHETIC FROM Z             00000250
CA                                                                      00000260
CA                                                                      00000270
CA              THIS ROUTINE COMPUTES THE SYNTHETIC TRACE FOR AN        00000280
CA          INPUT ESTIMATE OF WAVELET AND RC'S AND COMPUTES THE         00000290
CA          SUM OF SQUARES OF THE ERROR BETWEEN THIS SYNTHETIC AND      00000300
CA          THE INPUT TRACE.                                            00000310
CA                                                                      00000320
C                                                                       00000330
C    SUBROUTINES CALLED: NONE                                           00000340
C                                                                       00000350
C     EJECT                                                             00000360
C                                                                       00000370
      SUBROUTINE MERXX(E,XX,X,Z,KR,KW,KX,NZERO)                         00000380
CAEND                                                                   00000390
C                                                                       00000400
      DIMENSION Z(1),XX(1),X(1)                                         00000410
      IA=KW+1                                                           00000420
      IF((KW/2)*2.NE.KW)IA=KW                                           00000430
      E=0.                                                              00000440
C                                                                       00000450
          DO                                                        30  00000460
     *    I=1,KX                                                        00000470
          A=0.                                                          00000480
          IF(NZERO.EQ.0) GO TO                                      20  00000490
          J1=1                                                          00000500
          IF((I/2)*2.NE.I) J1=2                                         00000510
          IB=IA+I                                                       00000520
C                                                                       00000530
              DO                                                    10  00000540
     *        J=J1,KW,2                                                 00000550
C                                                                       00000560
   10         A=A+Z(KR+J)*Z((IB-J)/2)                                   00000570
C                                                                       00000580
   20     F=A-X(I)                                                      00000590
          XX(I)=A                                                       00000600
C                                                                       00000610
   30     E=E+F*F                                                       00000620
C                                                                       00000630
      RETURN                                                            00000640
      END                                                               00000650