anchor.diff:           !.......................................................................................
   
               ! Orientation of rotor centerline, normal to disk
  -         farm%m%AllOuts(RtAxsXT(i_turb)) = farm%FWrap(i_turb)%y%xHat_Disk(1)
  -         farm%m%AllOuts(RtAxsYT(i_turb)) = farm%FWrap(i_turb)%y%xHat_Disk(2)
  -         farm%m%AllOuts(RtAxsZT(i_turb)) = farm%FWrap(i_turb)%y%xHat_Disk(3)
  +         farm%m%AllOuts(RtAxsXT(nt)) = farm%FWrap(nt)%y%xHat_Disk(1)
  +         farm%m%AllOuts(RtAxsYT(nt)) = farm%FWrap(nt)%y%xHat_Disk(2)
  +         farm%m%AllOuts(RtAxsZT(nt)) = farm%FWrap(nt)%y%xHat_Disk(3)
            
               ! Center position of hub, m
  -         farm%m%AllOuts(RtPosXT(i_turb)) = farm%FWrap(i_turb)%y%p_hub(1)
  -         farm%m%AllOuts(RtPosYT(i_turb)) = farm%FWrap(i_turb)%y%p_hub(2)
  -         farm%m%AllOuts(RtPosZT(i_turb)) = farm%FWrap(i_turb)%y%p_hub(3)
  +         farm%m%AllOuts(RtPosXT(nt)) = farm%FWrap(nt)%y%p_hub(1)
  +         farm%m%AllOuts(RtPosYT(nt)) = farm%FWrap(nt)%y%p_hub(2)
  +         farm%m%AllOuts(RtPosZT(nt)) = farm%FWrap(nt)%y%p_hub(3)
            
               ! Rotor diameter, m
  -         farm%m%AllOuts(RtDiamT(i_turb)) = farm%FWrap(i_turb)%y%D_rotor
  +         farm%m%AllOuts(RtDiamT(nt)) = farm%FWrap(nt)%y%D_rotor
            
               ! Nacelle-yaw error at the wake planes, deg 
  -         farm%m%AllOuts(YawErrT(i_turb)) = farm%FWrap(i_turb)%y%YawErr*R2D
  +         farm%m%AllOuts(YawErrT(nt)) = farm%FWrap(nt)%y%YawErr*R2D
            
               ! Ambient turbulence intensity of the wind at the rotor disk, percent
  -            ! TODO: Is this really in percent form? GJH 3/21/2017
  -         farm%m%AllOuts(TIAmbT(i_turb))  = farm%AWAE%y%TI_amb(i_turb)*100.0_ReKi
  +         farm%m%AllOuts(TIAmbT(nt))  = farm%AWAE%y%TI_amb(nt)*100.0_ReKi 
            
               ! Rotor-disk-averaged ambient wind speed (normal to disk, not including structural motion, local induction or wakes from upstream turbines), m/s
  -         farm%m%AllOuts(RtVAmbT(i_turb)) = farm%AWAE%y%Vx_wind_disk(i_turb)
  +         farm%m%AllOuts(RtVAmbT(nt)) = farm%AWAE%y%Vx_wind_disk(nt)
            
               ! Rotor-disk-averaged relative wind speed (normal to disk, including structural motion and wakes from upstream turbines, but not including local induction), m/s
  -         farm%m%AllOuts(RtVRelT(i_turb)) = farm%FWrap(i_turb)%y%DiskAvg_Vx_Rel
  +         farm%m%AllOuts(RtVRelT(nt)) = farm%FWrap(nt)%y%DiskAvg_Vx_Rel
            
               ! Azimuthally averaged thrust force coefficient (normal to disk), distributed radially, -
            do ir = 1, farm%p%NOutRadii
  -            farm%m%AllOuts(CtTN(ir, i_turb)) = farm%FWrap(i_turb)%y%AzimAvg_Ct(farm%p%OutRadii(ir)+1)  ! y%AzimAvg_Ct is a 1-based array but the user specifies 0-based node indices, so we need to add 1
  +            farm%m%AllOuts(CtTN(ir, nt)) = farm%FWrap(nt)%y%AzimAvg_Ct(farm%p%OutRadii(ir)+1)  ! y%AzimAvg_Ct is a 1-based array but the user specifies 0-based node indices, so we need to add 1
            end do
            
            !.......................................................................................
  
positive.code:    ! close output file
     !.......................................................................................
     call Farm_EndOutput( farm, ErrStat2, ErrMsg2 )
        call SetErrStat(ErrStat2, ErrMsg2, ErrStat, ErrMsg, 'T'//trim(num2lstr(i_turb))//':'//RoutineName)
  end subroutine FARM_End
  !----------------------------------------------------------------------------------------------------------------------------------
  SUBROUTINE Transfer_FAST_to_WD(farm)
     type(All_FastFarm_Data),  INTENT(INOUT) :: farm                            !< FAST.Farm data  
  
     integer(intKi)  :: i_turb
     
     DO i_turb = 1,farm%p%NumTurbines   
        farm%WD(i_turb)%u%xhat_disk      = farm%FWrap(i_turb)%y%xHat_Disk       ! Orientation of rotor centerline, normal to disk
        farm%WD(i_turb)%u%p_hub          = farm%FWrap(i_turb)%y%p_hub           ! Center position of hub, m
        farm%WD(i_turb)%u%D_rotor        = farm%FWrap(i_turb)%y%D_rotor         ! Rotor diameter, m
        farm%WD(i_turb)%u%Vx_rel_disk    = farm%FWrap(i_turb)%y%DiskAvg_Vx_Rel  ! Rotor-disk-averaged relative wind speed (ambient + deficits + motion), normal to disk, m/s
        farm%WD(i_turb)%u%Ct_azavg       = farm%FWrap(i_turb)%y%AzimAvg_Ct      ! Azimuthally averaged thrust force coefficient (normal to disk), distributed radially, -
        farm%WD(i_turb)%u%YawErr         = farm%FWrap(i_turb)%y%YawErr          ! Nacelle-yaw error at the wake planes, rad   
     END DO
     
  END SUBROUTINE Transfer_FAST_to_WD
  
negative.code:       !--------------------
        ! 1. call WD_CO and transfer y_WD to u_AWAE        
     
     DO nt = 1,farm%p%NumTurbines
        
        call WD_CalcOutput( t, farm%WD(nt)%u, farm%WD(nt)%p, farm%WD(nt)%x, farm%WD(nt)%xd, farm%WD(nt)%z, &
                       farm%WD(nt)%OtherSt, farm%WD(nt)%y, farm%WD(nt)%m, ErrStat2, ErrMsg2 )         
           call SetErrStat(ErrStat2, ErrMsg2, ErrStat, ErrMsg, 'T'//trim(num2lstr(nt))//':'//RoutineName)
                 
     END DO
     if (ErrStat >= AbortErrLev) return
  