anchor.diff:   9991 format(3F12.4,X,6E16.6)
    9992 format(3F12.4,X,1E16.6)
    9993 format(3F12.4,X,3F16.8)
  - 9994 format(2F12.4,X,3F12.4)
         RETURN
         END
   
  
positive.code:           write(*,*)
        endif
        write(*,*)
       1 '-noz                    Output format: X Y VAL (instead of X ',
       2                               'Y Z VAL)'
        if (str.eq.'long') then
            write(*,*)
            write(*,9999)
            write(*,*)
        endif
        write(*,*)
       1 '-long                   Increase displacement output precision'
        if (str.eq.'long') then
            write(*,*)
  
negative.code:               call calcstn(strain,stlo,stla,stdp,nflt,
       1                     evlo,evla,evdp,str,dip,rak,dx,dy,slip,
       2                     vp,vs,dens,flttyp,xy)
                if (kstn.eq.1) then
                    write(22,9991) stlo,stla,stdp*1d-3,
       1                           strain(1,1),strain(2,2),strain(3,3),
       2                           strain(1,2),strain(1,3),strain(2,3)
                endif
            endif
  C         Compute (3x3) stress matrix
            if (ksts.ge.1) then
                call stn2sts(stress,strain,vp,vs,dens)
                if (ksts.eq.1) then
                    write(23,9991) stlo,stla,stdp*1d-3,
       1                           stress(1,1),stress(2,2),stress(3,3),
       2                           stress(1,2),stress(1,3),stress(2,3)
                elseif (ksts.eq.3) then
                    write(23,9991) stlo,stla,stdp*1d-3,estrs(stress)
                endif
            endif
  C         Resolve stresses onto target fault planes
            if (kshr.ge.1.or.knor.eq.1.or.kcou.eq.1) then
                call coulomb(coul,norml,shear,stress,trgstr,trgdip,trgrak,
       1                     frict)
                if (kshr.ge.1) then
                    write(25,9992) stlo,stla,stdp*1d-3,shear
                endif
                if (knor.eq.1) then
                    write(24,9992) stlo,stla,stdp*1d-3,norml
                endif
                if (kcou.eq.1) then
                    write(26,9992) stlo,stla,stdp*1d-3,coul
                endif
            endif
  C         Update progress indicator
  