anchor.diff:        ! Check if we are fully initialized
      if ( OtherState%firstPass ) then
         ! TODO: This entire block needs to be reviewed
  +                 
  +      ! Define a temparary version that is the horizontal component of u%xhat_disk (used to initialize all downwind wake planes).
  +      xxdisk = (/u%xhat_disk(1), 0.0, 0.0/)
  +      yydisk = (/0.0, u%xhat_disk(2), 0.0/)
  +      xydisknorm = TwoNorm(xxdisk + yydisk)
  +      tmp_xhat_disk = ( xxdisk + yydisk ) / xydisknorm
  +       
         do i = 0, p%NumPlanes - 1
            x_plane = u%Vx_rel_disk*real(i,ReKi)*real(p%DT,ReKi)
          
  
positive.code:    integer(intKi)                               :: ErrStat2
     character(ErrMsgLen)                         :: ErrMsg2
     character(*), parameter                      :: RoutineName = 'WD_CalcOutput'
     real(ReKi)                                   :: correction(3)
     real(ReKi)                                   :: x_plane
     errStat = ErrID_None
     errMsg  = ""
  
negative.code:          end if  
        do j = 1,p%NumRadii-1
              ! NOTE: xd%Vr_wake(0,:) was initialized to 0 and remains 0.
           xd%Vr_wake(j,i) = real(  j-1,ReKi)*(  xd%Vr_wake(j-1,i)  )/real(j,ReKi) &
              !  Vx_wake is for the                         [n+1]       ,      [n+1]        ,      [n]          , and    [n]        increments             
                           - real(2*j-1,ReKi)*p%dr * (  xd%Vx_wake(j,i) + xd%Vx_wake(j-1,i) - xd%Vx_wake(j,i-1) - xd%Vx_wake(j-1,i-1)  ) / ( real(4*j,ReKi) * dx )
        end do  
     end do ! i = 1,p%NumPlanes-1 
   
  