anchor.diff:            endif
   
                    !hr08
  -            strack(i)=crad*ptnfac(ix)
  -
  +strack(i)=crad*ptnfac(ix)
   if(ibbc.eq.0) then
     crkveb_d=parbe(ix,5)
     cikveb_d=parbe(ix,6)
  
positive.code:       xbb_d=sigmanq(2,imbb(i))*xrb_d
        zbb_d=sigmanq(1,imbb(i))*zrb_d
        call errf(zbb_d,xbb_d,cbzb_d,cbxb_d)
      else if(ibtyp.eq.1) then
        r2b_d=two*(sigman2(2,imbb(i))-sigman2(1,imbb(i)))
        rb_d=sqrt(r2b_d)
  
negative.code:       r2b_d=two*(sigman2(2,imbb(i))-sigman2(1,imbb(i)))
        rb_d=sqrt(r2b_d)
        rkb_d=(strack(i)*pisqrt)/rb_d                            !hr03
        xrb_d=abs(crkveb_d)/rb_d
        zrb_d=abs(cikveb_d)/rb_d
        call errf(zrb_d,xrb_d,crzb_d,crxb_d)
  