anchor.diff:  
               end if
               if (wtVelocity(turbI) .lt. 0.0_dp) then 
  -				wtVelocity(turbI) = 0.0_dp
  -			end if
  +                wtVelocity(turbI) = 0.0_dp
  +            end if
           end do
       end do
   
  
positive.code:             deltax0 = x - x0
  
              ! far wake region
              if (deltax0 > 0.0_dp) then
              
                  ! horizontal spread
                  call sigmay_func(ky, deltax0, rotorDiameter(turb), yaw(turb), sigmay)
  
negative.code:                 wtVelocity(turbI) = wtVelocity(turbI) - deltav
  
              ! near wake region (linearized)
              else if (x > 0.0_dp) then
      
                  ! horizontal spread
                  call sigmay_func(ky, deltax0, rotorDiameter(turb), yaw(turb), sigmay)
  