anchor.diff:  
       ! 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
  
positive.code:                 wsArray(loc) = wsArray(loc) - deltav
  
              end if
          end do
      end do
  
  
negative.code: ! 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, wake_combination_method, &
                               TI_calculation_method, wtVelocity)
  
      ! independent variables: turbineXw turbineYw turbineZ rotorDiameter
      !                        Ct yawDeg
  