anchor.diff:       1 '               -sta STAFILE -trg TRGFILE|S/D/R/F ',
        2                               '[-haf HAFFILE] [-xy]'
         write(*,*)
  -     1 '               -disp DSPFILE -strain STNFILE -stress ',
  +     1 '               -disp DSPFILE [-az] -strain STNFILE -stress ',
        2                             'STSFILE'
         write(*,*)
        1 '               -normal NORFILE -shear SHRFILE -coul COULFILE'
  
positive.code:               read(13,*) trgstr,trgdip,trgrak,frict
            endif
  C         Compute displacements
            if (kdsp.eq.1) then
                call calcdisp(uN,uE,uZ,stlo,stla,stdp,nflt,evlo,evla,evdp,
       1                      str,dip,rak,dx,dy,slip,vp,vs,dens,flttyp,xy)
                if (long.eq.0) then
                    write(21,9990) stlo,stla,stdp*1d-3,uE,uN,uZ
                else
                    write(21,9993) stlo,stla,stdp*1d-3,uE,uN,uZ
                endif
            endif
  C         Compute (3x3) strain matrix
  
negative.code:       REAL*8 evlo(FMAX),evla(FMAX),evdp(FMAX),str(FMAX),dip(FMAX),
       1       rak(FMAX),dx(FMAX),dy(FMAX),slip(FMAX)
        REAL*8 uN,uE,uZ,strain(3,3),stress(3,3),norml,shear,coul
        REAL*8 az,hdsp
        INTEGER kffm,kflt,kmag,ksta,khaf,ktrg
        INTEGER kdsp,kstn,ksts,knor,kshr,kcou,kgmt
        COMMON /ICHECK/ kffm,kflt,kmag,ksta,khaf,ktrg
        COMMON /OCHECK/ kdsp,kstn,ksts,knor,kshr,kcou,kgmt
  C     Open files specified for output
        if (kdsp.ne.0) open(unit=21,file=dspf,status='unknown')
        if (kstn.eq.1) open(unit=22,file=stnf,status='unknown')
        if (ksts.eq.1) open(unit=23,file=stsf,status='unknown')
        if (knor.eq.1) open(unit=24,file=norf,status='unknown')
  