anchor.diff:        prog = 0
         progtag = 0
         call progbar(prog,prog100,progtag)
  -      if (.not.trgmatchsta) read(24,*) trgstr,trgdip,trgrak,frict
  +      if (trgmatchsta.eq.0) read(24,*) trgstr,trgdip,trgrak,frict
     101 read (22,*,end=107) stlo,stla,stdp
  -          if (trgmatchsta) read(24,*) trgstr,trgdip,trgrak,frict
  +          if (trgmatchsta.eq.1) read(24,*) trgstr,trgdip,trgrak,frict
             stdp = 1.0d3*stdp
             do 103 i = 1,3
                 do 102 j = 1,3
  
positive.code:   104             continue
    105         continue
    106     continue
  C
            call strain2stress(stress,enet,vp,vs,dens)
            call coulomb(coul,norml,shear,stress,trgstr,trgdip,trgrak,
       1                 frict)
  
negative.code: 
        RETURN
        END
  
  C----------------------------------------------------------------------C
  
        SUBROUTINE mattr(matout,matin)
        IMPLICIT none
  