anchor.diff:  C Check for input receiver, half-space, and target parameter files
   C If resolving stresses, verify NTRG=1 or NTRG=NSTA
   C----
  +      if (verbos.ge.1) then
  +          write(0,*) 'CHECKING STATION, HALFSPACE, AND TARGET INPUTS'
  +      endif
         call chkhafstatrg(haff,staf,trgf,nsta,ntrg)
  +      if (verbos.ge.1) then
  +          write(0,*) 'OTHER INPUTS DEFINED'
  +          write(0,*)
  +      endif
         call readhaff(haff,vp,vs,dens)
   
   C----
  -C Read in station coordinates
  +C Read in station coordinates and see if they make sense
   C----
         call chksta(sta,staf,nsta,xy)
   
   C----
   C Compute displacements, strains, and stresses
   C----
  +      if (verbos.ge.1) then
  +          write(0,*) 'CALCULATING DEFORMATION'
  +      endif
  +      if (verbos.ge.2) then
  +          if (kdsp.ge.1) write(0,*) '    CALCULATING DISPLACEMENTS'
  +          if (kstn.ge.1) write(0,*) '    CALCULATING STRAINS'
  +          if (ksts.ge.1) write(0,*) '    CALCULATING STRESSES'
  +          if (knor.ge.1) write(0,*) '    CALCULATING NORMAL STRESSES'
  +          if (kshr.ge.1) write(0,*) '    CALCULATING SHEAR STRESSES'
  +          if (kcou.ge.1) write(0,*) '    CALCULATING COULOMB STRESSES'
  +      endif
         call calcdefm(staf,trgf,nsta,ntrg,nflt,evlo,evla,evdp,str,dip,rak,
        1              dx,dy,slip,flttyp,vp,vs,dens,dspf,stnf,stsf,norf,
        2              shrf,coulf,xy,prog,long)
  +      if (verbos.ge.1) then
  +          write(0,*) 'FINISHED WITH DEFORMATION CALCULATION'
  +          write(0,*)
  +      endif
   
   C----
   C Compute deformation for volume sources
  
positive.code:       COMMON /ICHECK/ kffm,kfsp,kflt,kmag,ksta,khaf,ktrg,kvol
        COMMON /OCHECK/ kdsp,kstn,ksts,knor,kshr,kcou,kgmt
  
  C----
  C Parse command line
  C----
        call gcmdln(ffmf,fspf,fltf,magf,haff,staf,trgf,dspf,stnf,stsf,
       1            norf,shrf,coulf,gmtf,flttyp,auto,incr,adist,xy,prog,
       2            long,emprel,volf,sthr,verbos)
        if (verbos.eq.1) then
            write(0,*) 'FFMF:   ',trim(ffmf)
            write(0,*) 'FSPF:   ',trim(fspf)
            write(0,*) 'FLTF:   ',trim(fltf)
  
negative.code: C----------------------------------------------------------------------C
  
        SUBROUTINE readsrc(ffmf,fspf,fltf,magf,nflt,evlo,evla,evdp,str,
       1                   dip,rak,dx,dy,slip,hylo,hyla,emprel,haff,sthr,
       2                   verbos)
  C----
  C Read shear dislocation parameters from fault source files. Parameter
  C FMAX defines maximum number of shear dislocations that can be stored.
  