anchor.diff:       + 2X,'*             December 1991                               *'/
        + 2X,'***********************************************************')
   C
  -      WRITE(*,200)
  -  200 FORMAT(4X,'Name of Input File =')
  -      READ(*,10) FIN
  -C
  -      WRITE(*,300)
  -  300 FORMAT(4X,'Name of Output File #1 (input, peak values .. etc) =')
  -      READ(*,10) FOUT
  -C
  -      WRITE(*,400)
  -  400 FORMAT(4X,'Name of Output File #2 (time histories .. etc) =')
  -      READ(*,10) PUNCH
  -C
  -   10 FORMAT(A32)
  -C
  -      OPEN(5,FILE=FIN,STATUS='OLD')
  -      OPEN(6,FILE=FOUT,STATUS='NEW')
  -      OPEN(7,FILE=PUNCH,STATUS='NEW')
  -C
  -      WRITE(6,100)
  +      write(*,200)
  +  200 format(4X,'Name of Input File =')
  +      read(*,10) FIN
  +C
  +      write(*,300)
  +  300 format(4X,'Name of Output File #1 (input, peak values .. etc) =')
  +      read(*,10) FOUT
  +C
  +      write(*,400)
  +  400 format(4X,'Name of Output File #2 (time histories .. etc) =')
  +      read(*,10) PUNCH
  +C
  +   10 format(A32)
  +C
  +      open(5,FILE=FIN,STATUS='OLD')
  +      open(6,FILE=FOUT,STATUS='NEW')
  +      open(7,FILE=PUNCH,STATUS='NEW')
  +C
  +      write(6,100)
         WW = .0624
         GT = 32.2
         MAMAX=4096
  
positive.code:       NS = NAA + 2*MAMAX
        NINV = NS + NAX/8 + 1
        NTOT = NINV + NAX/8 + 1
        IF (SKO .LT. .000001)   SKO = .45
        WRITE(6,2000) MAMAX, NTOT
  C
        CALL SHAKIT(X(1), X(NAX), X(NAA), X(NS), X(NINV))
  C
        STOP
  C ****************************************************
   1000 FORMAT(I5, F10.0)
   2000 FORMAT( 45H  MAX. NUMBER OF TERMS IN FOURIER TRANSFORM =   I10/
       1          45H  NECESSARY LENGTH OF BLANK COMMON X        =   I10)
        END
  C********************************************************************
        SUBROUTINE EARTHQ(X,AX,S,INV)
  C***********************************************************************
  C
  C   THIS ROUTINE  READS THE MOTION IN THE TIME DOMAIN, ADDS TRAILING
  
negative.code: C
  C***********************************************************************
  C
        character*6 TITLE
        character*30 FINPEQ
        character*80 HEAD
        character*12 FMAT
        complex X, AX
        dimension XR(8),X(300),AX(3,270),S(70),INV(70)
        common /EQ/ MFOLD,MA2,TITLE(5),DT, MA , MMA, DF, MX
        common/FRCUT/ NCUT,NZERO
        common /JISCK/ JIS,FINPEQ
  C
        PI2 = 6.283185307
        read(5,1001) NV, MA, DT, FINPEQ, FMAT
        read(5,1004) XF,XMAX,FMAX,NHEAD, NPL
        if (FMAX .LT. .001) FMAX = 100000.
        if (FMAT .EQ. '            ') FMAT = '(8F9.6,I7)'
        MA2=2
      2 if (MA2.GE.MA) go to 3
        MA2=MA2*2
        go to 2
      3 MA=MA2
  C  ....................................................................
        write(6,2012) FINPEQ, NV, MA, NHEAD, NPL, DT, FMAT
        write (*,2026) FINPEQ,FMAT
        open(8,FILE=FINPEQ,STATUS='OLD')
        write (6,2021)
        do 4 I=1,NHEAD
        read(8,2022) HEAD
      4 write(6,2022) HEAD
  C      write (6,2023)
        MMA = NV + NV / 10
        if (MMA .GT. MA) MMA=MA
        MA2 = MA + 2
        MFOLD = MA2/2
        MFOLD = MFOLD + 1
  