anchor.diff:          ek(ix) = one
           el(ix) = zero
         else if(abs(el(ix)+two) <= pieni) then
  -        dki(ix,2) = ed(ix)
  +        dki(ix,2) = -ed(ix)
           dki(ix,3) = ek(ix)
  +        print *, "heerea"
           ed(ix) = one
           ek(ix) = one
           el(ix) = zero
  
positive.code:     strack(i) = smiv(i)*c1m24
    case(11)
      strack(i) = dki(ix,1)/dki(ix,3)
    case(12)
      strack(i) = dki(ix,1)
    case(13)
      strack(i) = dki(ix,2)/dki(ix,3)
    case(14)
      strack(i) = dki(ix,2)
    case default
  
negative.code: 
  
  ! ================================================================================================ !
  !  Routines related to the initialisation of elements
  !  Need to be in a separate file due to dependency issues
  