anchor.diff:                 return
               end if
         
  -         y%p_plane     (:,i) = u%p_hub(:) + x_plane*u%xhat_disk(:) + correction
  -         y%xhat_plane  (:,i) = u%xhat_disk(:)
  +         y%p_plane     (:,i) = u%p_hub(:) + x_plane*tmp_xhat_disk(:) + correction
  +
  +         if ( i == 0 ) then
  +            y%xhat_plane   (:,i) = u%xhat_disk(:)
  +         else
  +            y%xhat_plane   (:,i) = tmp_xhat_disk(:)
  +         end if         
            
               ! NOTE: Since we are in firstPass=T, then xd%Vx_wake is already set to zero, so just pass that into WakeDiam
            y%D_wake(i)  =  WakeDiam( p%Mod_WakeDiam, p%NumRadii, p%dr, p%r, xd%Vx_wake(:,i), u%Vx_wind_disk, u%D_rotor, p%C_WakeDiam)
  
positive.code:       ! Check if we are fully initialized
     if ( OtherState%firstPass ) then
        ! TODO: This entire block needs to be reviewed
        do i = 0, p%NumPlanes - 1
           x_plane = u%Vx_rel_disk*real(i,ReKi)*real(p%DT,ReKi)
         
  
negative.code:    integer(intKi)                               :: ErrStat2
     character(ErrMsgLen)                         :: ErrMsg2
     character(*), parameter                      :: RoutineName = 'WD_CalcOutput'
     real(ReKi)                                   :: correction(3), xxdisk(3), yydisk(3), xydisknorm, tmp_xhat_disk(3)
     real(ReKi)                                   :: x_plane
     errStat = ErrID_None
     errMsg  = ""
  