anchor.diff:                   , ErrStat  = ErrStat2         &
                    , ErrMess  = ErrMsg2          )
         CALL SetErrStat( ErrStat2, ErrMsg2, ErrStat, ErrMsg, RoutineName )
  -      
  +      if (ErrStat>=AbortErrLev) RETURN
         ! initialization (not necessary)
      !y%BldForce%Force(:,:)  = 0.0_BDKi
      !y%BldForce%Moment(:,:) = 0.0_BDKi
         
      !.................................
  -   ! y%ReactionForce
  +   ! y%ReactionForce (for coupling with ElastoDyn)
      !.................................
         
      CALL MeshCopy( SrcMesh   = u%RootMotion     &
  
positive.code:    character(*),                 intent(  out)  :: ErrMsg            ! Error message if ErrStat /= ErrID_None
  
        ! local variables
     real(ReKi)                                   :: DCM(3,3)
     real(BDKi)                                   :: TmpDCM(3,3)
     real(BDKi)                                   :: TmpPos(3)
     real(BDKi)                                   :: temp_POS(3)
     real(BDKi)                                   :: temp_CRV(3)
     
  
negative.code: end subroutine SetInitOut
  !-----------------------------------------------------------------------------------------------------------------------------------
  subroutine Init_y( p, u, y, ErrStat, ErrMsg)
  ! this routine initializes the outputs, y, that are used in the BeamDyn interface for coupling in the FAST framework.
  
     type(BD_ParameterType),       intent(inout)  :: p                 ! Parameters  ! intent(out) only because it changes p%NdIndx
     type(BD_InputType),           intent(inout)  :: u                 ! Inputs
  