anchor.diff:                elseif (noz.eq.1.and.kdsp.eq.2) then
                     az = atan2(uE,uN)*r2d
                     hdsp = dsqrt(uE*uE+uN*uN)
  -                  write(21,*) stlo,stla,az,hdsp,uz
  +                  write(21,9994) stlo,stla,az,hdsp,uz
                 endif
             endif
   C         Compute (3x3) strain matrix
  
positive.code:               kstn = 1
                i = i + 1
                call getarg(i,stnf)
            elseif (tag(1:7).eq.'-stress') then
                if (kstn.eq.0) kstn = 2
                ksts = 1
  
negative.code:       stress(2,1) = stress(1,2)
        stress(3,1) = stress(1,3)
        stress(3,2) = stress(2,3)
        RETURN
        END
  
  C----------------------------------------------------------------------C
  
        FUNCTION estrs(sts)
        IMPLICIT none
        REAL*8 estrs,sts(3,3),s12,s13,s23
        s12 = sts(1,1)-sts(2,2)
        s13 = sts(1,1)-sts(3,3)
        s23 = sts(2,2)-sts(3,3)
        estrs = (1.0d0/6.0d0)*(s12*s12+s13*s13+s23*s23)
       1          + sts(1,2)*sts(1,2)+sts(1,3)*sts(1,3)+sts(2,3)*sts(2,3)
        RETURN
        END
  
  