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:         
          ! Check if this is the max and use it if it is
          if (ti_area_ratio_tmp > ti_area_ratio) then
  !            print *, "ti_area_ratio_tmp > ti_area_ratio"
             TI_dst = TI_tmp
             ti_area_ratio = ti_area_ratio_tmp
          end if
  