anchor.diff:                  a(i,1,2) = rho*si
                   a(i,1,3) = zero
                   a(i,1,4) = one
  -                a(i,1,5) = ((-one*rho)*(one-co))/dpsq         ! hr06
  -                a(i,1,6) = ((-one*two)*tan_mb(fok*half))/dpsq ! hr06
  +                a(i,1,5) = ((-one*rho)*(one-co))/dpsq
  +                a(i,1,6) = ((-one*two)*tan_mb(fok*half))/dpsq
                   ! VERTICAL
                   a(i,2,1) = one-gl
                   a(i,2,2) = el(i)
  -                a(i,2,3) = (-one*g)*(two-gl) ! hr06
  +                a(i,2,3) = (-one*g)*(two-gl)
                   a(i,2,4) = a(i,2,1)
               case (4)
                   ! HORIZONTAL
                   a(i,1,1) = co
                   a(i,1,2) = rho*si
  -                a(i,1,3) = (-one*si)/rho              ! hr06
  +                a(i,1,3) = (-one*si)/rho
                   a(i,1,4) = co
  -                a(i,1,5) = ((-one*rho)*(one-co))/dpsq ! hr06
  -                a(i,1,6) = (-one*si)/dpsq             ! hr06
  +                a(i,1,5) = ((-one*rho)*(one-co))/dpsq
  +                a(i,1,6) = (-one*si)/dpsq
                   ! VERTICAL
                   a(i,2,1) = one
                   a(i,2,2) = el(i)
  
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
        rho=(one/ed(i))*dpsq
  
negative.code:                 fokq = ek(i)
              else
                  ih   = 1
                  fokq = -one*ek(i)
              end if
              wf  = ed(i)/dpsq
              fok = fokq/dpd-wf**2
              if(abs(fok).le.pieni) then
                  do l=1,2
                      aeg(i,l,1) = one
  