anchor.diff:  
               ! near wake region (linearized)
               else if (deltax0 > -x0) then
  -
  +    
                   ! horizontal spread
                   call sigmay_func(ky, deltax0, rotorDiameter(turb), yaw(turb), sigmay)
                   
  
positive.code:     
  end subroutine sigmay_func
      
  ! calculates the vertical spread of the wake at a given distance from the onset of far 
  ! wake condition
  subroutine sigmaz_func(kz, deltax0, rotor_diameter, sigmaz)
  
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
  
  