anchor.diff:      sigmaz = rotor_diameter * (kz * deltax0 / rotor_diameter + 1.0_dp / sqrt(8.0_dp))
       
   end subroutine sigmaz_func
  +
   
   ! calculates the horizontal distance from the wake center to the hub of the turbine making
   ! the wake
  
positive.code:                             
  
      ! determine the onset location of far wake
  	x0 = rotor_diameter * (cos(yaw) * (1.0_dp + sqrt(1.0_dp - Ct)) / &
  								(sqrt(2.0_dp) * (alpha * I + beta * &
  												& (1.0_dp - sqrt(1.0_dp - Ct)))))
  
  end subroutine x0_func
  
  ! calculates the wake angle at the onset of far wake conditions
  subroutine theta_c_0_func(yaw, Ct, theta_c_0)
  
negative.code:     ! local (General)
      real(dp), dimension(nTurbines) :: yaw
      real(dp) :: x0, deltax0, deltay, theta_c_0, sigmay, sigmaz, wake_offset
      real(dp) :: x, deltav, deltav0m, deltaz, sigmay0, sigmaz0
      Integer :: turb, turbI
      real(dp), parameter :: pi = 3.141592653589793_dp
  
  