anchor.diff:        kshr = 0
         kcou = 0
         kgmt = 0
  +C Parse command line
         narg = iargc()
         if (narg.eq.0) call usage('!! Error: No command line arguments '//
        1                                     'specified')
  
positive.code: 
  C----------------------------------------------------------------------C
  
        SUBROUTINE calcdefm(stafile,trgfile,nsta,ntrg,nflt,evlo,evla,evdp,
       1                    str,dip,rak,dx,dy,slip,flttyp,vp,vs,dens,
       2                    dspfile,stnfile,stsfile,norfile,shrfile,
       3                    coufile,xy)
  C----
  C Compute displacements, strains, and stresses from shear dislocations.
  C----
        IMPLICIT NONE
        CHARACTER*40 stafile,trgfile,dspfile,stnfile,stsfile,norfile,
       4             shrfile,coufile
        INTEGER i,nsta,ntrg,progout
        REAL*8 stlo,stla,stdp,trgstr,trgdip,trgrak,frict,vp,vs,dens
        INTEGER nflt,FMAX,flttyp,xy
        PARAMETER (FMAX=1500)
        REAL*8 evlo(FMAX),evla(FMAX),evdp(FMAX),str(FMAX),dip(FMAX),
       1       rak(FMAX),dx(FMAX),dy(FMAX),slip(FMAX)
  
negative.code:      1 '    Format: stlo stla s_EE(Pa) s_NN(Pa) s_ZZ(Pa) s_EN(Pa) ',
       2                         's_EZ(Pa) s_NZ(Pa)'
            write(*,*)
       1 '    The stress matrix is symmetric; the output contains ',
       2       'the six independent'
            write(*,*)
       1 '    components. In an isotropic, elastic material (summation ',
  