anchor.diff:            write(*,*)
         endif
         write(*,*)
  +     1 '-az                     Displacement output -> (az hdsp Z)'
  +      if (str.eq.'long') then
  +          write(*,*)
  +          write(*,*)
  +     1 '    Format: stlo stla stdp AZIM U_HOR(m) u_Z(m)'
  +          write(*,*)
  +     1 '    Vertical displacement is positive up.'
  +          write(*,*)
  +     1 '    Azimuth is clockwise from north.'
  +          write(*,*)
  +          write(*,9999)
  +          write(*,*)
  +      endif
  +      write(*,*)
        1 '-strain STRNFILE        Strain matrix (EE NN ZZ EN EZ NZ)'
         if (str.eq.'long') then
             write(*,*)
  
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: C Compute displacements, strains, and stresses from shear dislocations.
  C----
        IMPLICIT NONE
        REAL*8 pi,r2d
        PARAMETER (pi=4.0d0*atan(1.0d0),r2d=1.8d2/pi)
        CHARACTER*40 staf,trgf,dspf,stnf,stsf,norf,shrf,coulf
        INTEGER i,nsta,ntrg,progout
        REAL*8 stlo,stla,stdp,trgstr,trgdip,trgrak,frict,vp,vs,dens
  