anchor.diff:                    ktrg = 2
                 endif
             elseif (tag(1:5).eq.'-disp') then
  -              kdsp = 1
  +              if (kdsp.eq.0) kdsp = 1
                 i = i + 1
                 call getarg(i,dspf)
  +          elseif (tag(1:3).eq.'-az') then
  +              kdsp = 2
             elseif (tag(1:7).eq.'-strain') then
                 kstn = 1
                 i = i + 1
  
positive.code:           write(*,*)
        endif
        write(*,*)
       1 '-strain STRNFILE        Strain matrix (EE NN ZZ EN EZ NZ)'
        if (str.eq.'long') then
            write(*,*)
  
negative.code:               read(13,*) trgstr,trgdip,trgrak,frict
            endif
  C         Compute displacements
            if (kdsp.ne.0) 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.and.kdsp.eq.1) then
                    write(21,9990) stlo,stla,stdp*1d-3,uE,uN,uZ
                elseif (long.eq.1.and.kdsp.eq.1) then
                    write(21,9993) stlo,stla,stdp*1d-3,uE,uN,uZ
                elseif (kdsp.eq.2) then
                    az = atan2(uE,uN)*r2d
                    hdsp = dsqrt(uE*uE+uN*uN)
                    write(21,9990) stlo,stla,stdp*1d-3,az,hdsp,uz
                endif
            endif
  C         Compute (3x3) strain matrix
  