anchor.diff:              co = cos_mb(fi)
               a(i,ih,1) = co
               a(i,ih,2) = si/hi
  -            a(i,ih,3) = (-one*si)*hi ! hr06
  +            a(i,ih,3) = (-one*si)*hi
               a(i,ih,4) = co
   
           case (9) ! EDGE FOCUSSING
               rhoi = ed(i)/dpsq
  -            fok  = rhoi*tan_mb((el(i)*rhoi)*half) ! hr06
  +            fok  = rhoi*tan_mb((el(i)*rhoi)*half)
               a(i,1,1) = one
               a(i,1,2) = zero
               a(i,1,3) = fok
               a(i,1,4) = one
               a(i,2,1) = one
               a(i,2,2) = zero
  -            a(i,2,3) = -one*fok ! hr06
  +            a(i,2,3) = -one*fok
               a(i,2,4) = one
   
           end select
  
positive.code:           al(2,l,j,i) = el(i)
            al(3,l,j,i) = zero
            al(4,l,j,i) = one
            as(6,l,j,i) = ((-one*rv)*al(2,l,j,i))/c2e3 ! hr05
          end do
          as(1,1,j,i) = (el(i)*(one-rv))*c1e3          ! hr05
          cycle
        end if
        ih=0
  
negative.code:           al(2,l,j,i) = el(i)
            al(3,l,j,i) = zero
            al(4,l,j,i) = one
            as(6,l,j,i) = ((-one*rv)*al(2,l,j,i))/c2e3
          end do
          as(1,1,j,i) = (el(i)*(one-rv))*c1e3
          cycle
        end if
        afok=abs(fok)
  