anchor.diff:        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
  +      if (SKO .LT. .000001)   SKO = .45
  +      write(6,2000) MAMAX, NTOT
  +C
  +      callSHAKIT(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/
  + 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
  +      end
   C********************************************************************
  -      SUBROUTINE EARTHQ(X,AX,S,INV)
  +      subroutine EARTHQ(X,AX,S,INV)
   C***********************************************************************
   C
   C   THIS ROUTINE  READS THE MOTION IN THE TIME DOMAIN, ADDS TRAILING
  
positive.code:      + 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)
        WW = .0624
        GT = 32.2
        MAMAX=4096
  
negative.code: C $NOFLOATCALLS
  C $NODEBUG
  C ...................................................................
        program SHAKE91
        character*32 FIN,FOUT,PUNCH
        common X(25620)
        common /TIME/ T(9)
        common /WGK/  WW, GT, SKO
  C ...................................................................
        write(*,100)
    100 format(2X,'*****************************************************'/
       + 2X,'* SHAKE  --   A COMPUTER PROGRAM FOR EARTHQUAKE RESPONSE  *'/
       + 2X,'*             ANALYSIS OF HORIZONTALLY LAYERED SITES      *'/
       + 2X,'*             by: Per B. Schnabel & John Lysmer -- 1970   *'/
  