anchor.diff:              ! Get BeamDyn's contribution: (note if p_FAST%CompElast /= Module_BD, SizeJac_ED_SD_HD_BD(5) = 0)
               !............................... 
               DO nb=1,p_FAST%nBeams
  +               CALL BD_CopyOutput(y_BD(nb),y_BD_perturb(nb),MESH_UPDATECOPY,ErrStat2,ErrMsg2)
  +                  CALL SetErrStat( ErrStat2, ErrMsg2, ErrStat, ErrMsg, RoutineName  )
  +            END DO
  +                        
  +            DO nb=1,p_FAST%nBeams
                  
                  ! make sure we perturb the outputs from only the current blade (overwrite the previous perturbation with y_BD values)
                  if (nb > 1) then               
  
positive.code:    REAL(ReKi)                        , INTENT(IN   ) :: u_in(:)
     REAL(ReKi)                        , INTENT(  OUT) :: U_Resid(:)
  
     INTEGER(IntKi)                                    :: i                      ! counter for ice leg loop
     
     !..................
     ! Set mooring line and ice inputs (which don't have acceleration fields and aren't used elsewhere in this routine, thus we're using the actual inputs (not a copy) 
  
negative.code:       
        MeshMapData%u_ED_HubPtLoad%Force  = 0.0_ReKi
        MeshMapData%u_ED_HubPtLoad%Moment = 0.0_ReKi
        IF ( p_FAST%CompElast == Module_BD ) THEN
           
           ! Transfer ED motions to BD inputs:
           call Transfer_ED_to_BD_tmp( y_ED2, MeshMapData, ErrStat2, ErrMsg2 )   ! sets MeshMapData%u_BD_RootMotion(:)
              CALL SetErrStat(ErrStat2,ErrMsg2, ErrStat, ErrMsg, RoutineName)
              
           ! Transfer BD loads to ED hub input:
              ! we're mapping loads, so we also need the sibling meshes' displacements:
           do i=1,p_FAST%nBeams            
              CALL Transfer_Point_to_Point( y_BD2(i)%ReactionForce, MeshMapData%u_ED_HubPtLoad_2, MeshMapData%BD_P_2_ED_P(i), ErrStat2, ErrMsg2, MeshMapData%u_BD_RootMotion(i), y_ED2%HubPtMotion) !u_BD_RootMotion and y_ED2%HubPtMotion contain the displaced positions for load calculations
                 CALL SetErrStat(ErrStat2,ErrMsg2, ErrStat, ErrMsg, RoutineName)
                 
              MeshMapData%u_ED_HubPtLoad%Force  = MeshMapData%u_ED_HubPtLoad%Force  + MeshMapData%u_ED_HubPtLoad_2%Force  
  