anchor.diff: -      PROGRAM flt2displacement
  +      PROGRAM fltxy2displacement
   C----
   C Calculate the north, east and vertical displacements at stations
   C due to shear dislocations in an isotropic halfspace.
   C
   C To run, requires the following files:
   C     faults.txt: list of ruptures
  -C     stations.txt: list of receiver lon, lat, and depth (km)
  +C     stations.txt: list of receiver position (km), and depth (km)
   C     structure.txt: vp (m/s) vs (m/s) density (kg/m^3)
   C
   C Produces the output file:
  -C     disp.out: stlo, stla, uN, uE, uZ
  +C     disp.out: x, y, uN, uE, uZ
  +C
  +C MODIFICATIONS
  +C   2014-02-24: Added help/usage option
  +C
   C----
         IMPLICIT NONE
         REAL*8 pi,d2r
         PARAMETER (pi=4.0d0*datan(1.0d0),d2r=pi/1.8d2)
         
  +      LOGICAL verbose
         CHARACTER*20 srcfile,stafile,haffile
         INTEGER flt,nflt,maxflt
         PARAMETER (maxflt=50)
  
positive.code:           write (*,8888)
            read *,haffile
            if (haffile.eq.'quit') stop
            goto 13
        else
             write(*,9999)
  
negative.code:           write (*,8887)
            read *,stafile
            if (stafile.eq.'quit') stop
            if (stafile.eq.'help') call usage()
            goto 12
        else
             write(*,9999)
  