anchor.diff:  ! implementation of the Bastankhah and Porte Agel (BPA) wake model for analysis
   subroutine porteagel_analyze(nTurbines, turbineXw, turbineYw, turbineZ, &
                                rotorDiameter, Ct, wind_speed, &
  -                             yawDeg, ky, kz, alpha, beta, I, wtVelocity)
  +                             yawDeg, ky, kz, alpha, beta, I, wake_combination_method, &
  +                             TI_calculation_method, wtVelocity)
   
       ! independent variables: turbineXw turbineYw turbineZ rotorDiameter
       !                        Ct yawDeg
  
positive.code:                 wsArray(loc) = wsArray(loc) - deltav
  
              end if
          end do
      end do
  
  
negative.code: 
      ! in
      integer, intent(in) :: nTurbines
      integer, intent(in) :: wake_combination_method, TI_calculation_method
      real(dp), dimension(nTurbines), intent(in) :: turbineXw, turbineYw, turbineZ
      real(dp), dimension(nTurbines), intent(in) :: rotorDiameter, yawDeg
      real(dp), dimension(nTurbines), intent(in) :: Ct
  