anchor.diff:                                                 wake_combination_method, deficit_sum)
                   
                       if ((x > 0.0_dp) .and. (TI_calculation_method > 0)) then
  +                    
  +                        ! save ti_area_ratio value to new memory location to avoid aliasing
  +                        TI_area_ratio_temp = TI_area_ratio
  +                        TI_dst_tmp = TIturbs(turbI)
                   
                           ! calculate TI value at each turbine
                           call added_ti_func(TI, Ct(turb), x, ky_local(turb), rotorDiameter(turb), & 
                                              & rotorDiameter(turbI), deltay, turbineZ(turb), &
                                              & turbineZ(turbI), TIturbs(turb), &
  -                                           & TI_calculation_method, ti_area_ratio, &
  +                                           & TI_calculation_method, TI_area_ratio_temp, 
  +                                           & TI_dst_tmp, TI_area_ratio, &
                                              & TIturbs(turbI))
                       end if
                   
  
positive.code: ! combines wakes using various methods
  subroutine added_ti_func(TI, Ct_ust, x, k_star_ust, rotor_diameter_ust, rotor_diameter_dst, & 
                          & deltay, wake_height, turbine_height, TI_ust, &
                          & TI_calculation_method, ti_area_ratio, TI_dst)
                                   
      implicit none
          
  
negative.code:     ! in
      real(dp), intent(in) :: Ct_ust, x, k_star_ust, rotor_diameter_ust, rotor_diameter_dst
      real(dp), intent(in) :: deltay, wake_height, turbine_height, TI_ust, TI
      real(dp), intent(in) :: TI_area_ratio_in, TI_dst_in
      integer, intent(in) :: TI_calculation_method
      
      ! local
      real(dp) :: axial_induction_ust, beta, epsilon, sigma, wake_diameter, wake_overlap
      real(dp) :: TI_added, TI_tmp, rotor_area_dst, TI_area_ratio_tmp
      real(dp), parameter :: pi = 3.141592653589793_dp
      
      ! out  
      real(dp), intent(out) :: TI_dst, TI_area_ratio
      
      ! intrinsic functions
      intrinsic sqrt
  