anchor.diff:  
       if(ix <= 0) then
         ix = -1*ix
  -      jj = mel(ix)+1         !initial idx
  -      dj = -1                !step
  -    endif
  +      jj = mel(ix)+1 ! initial idx
  +      dj = -1        ! step
  +    end if
       jm = mel(ix)
   
   !-- Loop over elements inside the block
  -    do 150 j=1,jm
  +    do j=1,jm
         jj = jj+dj       ! Subelement index of current sub=element
         jk = mtyp(ix,jj) ! Single-element index of the current sub-element
         if(ithick == 1 .and. kz(jk) /= 0) goto 120
  
positive.code:       select case(kzz)
        case(1)
  !--VERTICAL DIPOLE
          ekk=ekk*c1e3
  #include "include/kickl01v.f90"
  #include "include/kickq01v.f90"
  !--SKEW QUADRUPOLE
  
negative.code: #include "include/kicklxxv.f90"
  !--SKEW 14-POLE
        case(7)
          ekk = ekk*c1m15
  #include "include/kicksho.f90"
  #include "include/kicksho.f90"
  #include "include/kicksho.f90"
  