anchor.diff:                  wtVelocity(turbI) = wtVelocity(turbI) - deltav
   
               ! near wake region (linearized)
  -            else if (deltax0 > -x0) then
  +            else if (x > 0.0_dp) then
       
                   ! horizontal spread
                   call sigmay_func(ky, deltax0, rotorDiameter(turb), yaw(turb), sigmay)
  
positive.code:     deltav = (((deltav0m - deltavs)/x0) * x + deltavs) *              &
          exp(-0.5_dp * (deltay / sigmay) ** 2) *                       &
          exp(-0.5_dp * (deltaz / sigmaz) ** 2)
  	! deltav = deltav0m *              &
  !         exp(-0.5_dp * (deltay / sigmay) ** 2) *                       &
  !         exp(-0.5_dp * (deltaz / sigmaz) ** 2)
      
  end subroutine deltav_near_wake_lin_func
  
negative.code:             deltax0 = x - x0
  
              ! far wake region
              if (x > x0) then
              
                  ! horizontal spread
                  call sigmay_func(ky, deltax0, rotorDiameter(turb), yaw(turb), sigmay)
  