CC---PROGRAM TO GENERATE AMPLITUDE FILES
      SUBROUTINE ABQMAIN
C
      IMPLICIT REAL*8(A-H,O-Z)
      DIMENSION P2(84),TIME(84),P3(84)
      DATA P2/30.,32.5,35.,37.5,40.,40.5,40.,40.5,41.,41.5,42.,
     1 42.5,43.,43.5,44.,44.5,45.,45.5,46.,46.188,46.,45.75,45.5,
     2 45.25,45.,44.75,44.5,44.25,44.,43.75,43.5,43.25,43.,42.75,
     3 42.5,42.25,42.,41.75,41.5,41.25,41.,40.75,40.5,40.25,
     4 40.,39.,38.,37.,36.,34.,32.,30.,27.5,25.,20.,15.,10.,5.,0.,
     6 -5.,-10.,-14.,-17.,-20.,-21.,-22.,-23.,-24.,-25.,-26.,
     7 -27.,-28.,-29.,-30.,-31.,-32.,-33.,-34.,-35.,-36.,-37.,-38.,
     8 -39.,-40./
      OPEN(UNIT=16,FILE='PLATE3.AMP')
      OPEN(UNIT=17,FILE='PLATE2.AMP')
      DO 100 I=1,84
      XI=I
      TIME(I)=XI/100.
      IF(I-7) 1,2,2
    1 P3(I)=0.
      GO TO 40
    2 ARG=P2(I)*P2(I)-4.*(P2(I)*P2(I)-40.*40.)
      ARG=SQRT(ARG)
      IF(I.GT.20) GO TO 20
      P3(I)=(P2(I)-ARG)/2.
      GO TO 40
   20 CONTINUE
      P3(I)=(P2(I)+ARG)/2.
   40 CONTINUE
      II=MOD(I,4)
      IF(II.NE.0) GO TO 100
      WRITE(17,99) TIME(I-3),P2(I-3),TIME(I-2),P2(I-2),TIME(I-1),
     1P2(I-1),TIME(I),P2(I)
      WRITE(16,99) TIME(I-3),P3(I-3),TIME(I-2),P3(I-2),TIME(I-1),
     1P3(I-1),TIME(I),P3(I)
C     WRITE(6,90) TIME(I-3),P2(I-3),TIME(I-2),P2(I-2),TIME(I-1),
C    1P2(I-1),TIME(I),P2(I),TIME(I-3),P3(I-3),TIME(I-2),P3(I-2),
C    2TIME(I-1),P3(I-1),TIME(I),P3(I)
   99 FORMAT(8(F10.6,','))
C  90 FORMAT(3X,8F10.6)
  100 CONTINUE
      REWIND 17
      REWIND 16
      STOP
      END

