anchor.diff:              aeg(i,ih,2) = hs/hi
               aeg(i,ih,3) = hs*hi
               aeg(i,ih,4) = hc
  -            aeg(i,ih,5) = ((wf/afok)*(one-hc))/dpsq ! hr06
  -            aeg(i,ih,6) = (((-one*wf)/hi)*hs)/dpsq  ! hr06
  +            aeg(i,ih,5) = ((wf/afok)*(one-hc))/dpsq
  +            aeg(i,ih,6) = (((-one*wf)/hi)*hs)/dpsq
               ih = ih+1
               if(ih.gt.2) ih = 1
               hi = sqrt(abs(ek(i)/dpd))
  
positive.code:             co = cos_mb(fi)
              aeg(i,ih,1) = co
              aeg(i,ih,2) = si/hi
              aeg(i,ih,3) = (-one*si)*hi                      ! hr06
              aeg(i,ih,4) = co
              aeg(i,ih,5) = (((-one*wf)/afok)*(one-co))/dpsq  ! hr06
              aeg(i,ih,6) = (((-one*wf)/hi)*si)/dpsq          ! hr06
              ih = ih+1
              if(ih.gt.2) ih=1
              hi = sqrt(abs(ek(i)/dpd))
  
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
        rho=(one/ed(i))*dpsq
  