anchor.diff:                      if ((x > 0.0_dp) .and. (TI_calculation_method > 0)) then
                           !print *, "turbI, turb: ", turbI, turb
                           ! calculate TI value at each turbine
  -                        print *, "turb, turbI: ", turb, turbI
  +!                         print *, "turb, turbI: ", turb, turbI
                           call added_ti_func(TI, Ct(turb), x, ky_local, rotorDiameter(turb), & 
                                              & rotorDiameter(turbI), deltay, turbineZ(turb), &
                                              & turbineZ(turbI), TIturbs(turb), &
  
positive.code:     ! initialize TI of all turbines to free-stream value
      !print *, "start TIturbs: ", TIturbs
      TIturbs = TI
      !print *, "initialized TIturbs: ", TIturbs
      ky_local = ky
      kz_local = kz
      
  
negative.code:         ! get index of downstream turbine
          turbI = sorted_x_idx(d) + 1
          
          do, p=1, nRotorPoints
          
              ! initialize the TI_area_ratio to 0.0 for each turbine
              ti_area_ratio = 0.0_dp
      
              ! initialize deficit summation term to zero
              deficit_sum = 0.0_dp
  