anchor.diff:      !print *, "initialized TIturbs: ", TIturbs
       ky_local = ky
       kz_local = kz
  +    
  +    !print *, "ky_local: ", ky_local
  +    !print *, "kz_local: ", kz_local
  +    !print *, "TIturbs init: ", TIturbs
   
       do, d=1, nTurbines
       
  
positive.code:     logical, intent(in) :: calc_k_star
  
      ! local (General)
      real(dp), dimension(nTurbines) :: yaw, TIturbs, k_star, wtVelocity
      real(dp) :: x0, deltax0, deltay, theta_c_0, sigmay, sigmaz, wake_offset, discontinuity_point
      real(dp) :: x, deltav, deltaz, sigmay_dp, sigmaz_dp, deltax0_dp, deficit_sum, tol
      real(dp) :: LocalRotorPointY, LocalRotorPointZ, point_velocity, point_z, point_velocity_with_shear
      real(dp), dimension(nTurbines) :: ky_local, kz_local
      Integer :: u, d, turb, turbI, loc, p
  
negative.code:     real(dp), dimension(nRotorPoints), intent(in) :: RotorPointsY, RotorPointsZ
  
      ! local (General)
      real(dp), dimension(nTurbines) :: yaw, TIturbs
      real(dp) :: x0, deltax0, deltay, theta_c_0, sigmay, sigmaz, wake_offset, k_star
      real(dp) :: x, deltav, deltaz, sigmay_dp, sigmaz_dp, deltax0_dp, deficit_sum
      real(dp) :: ky_local, kz_local, tol, discontinuity_point
      real(dp) :: LocalRotorPointY, LocalRotorPointZ, point_velocity, point_z, point_velocity_with_shear
  