anchor.diff:                                      temp_Force)
      ENDIF
   
  -   CALL BD_MotionTensor(p%GlbRot,p%GlbPos,MoTens,1)
      temp6(:) = 0.0D0
      IF(p%analysis_type .EQ. 2) THEN
  -!       temp6(:) = temp_ReactionForce(1:6)
          temp6(1:6) = -OS_tmp%Acc(1:6)
  -       temp6(:) = MATMUL(TRANSPOSE(MoTens),temp6)
  -       y%ReactionForce%Force(1:3,1) = temp6(1:3)
  -       y%ReactionForce%Moment(1:3,1) = temp6(4:6)
  +       y%ReactionForce%Force(1:3,1) = MATMUL(p%GlbRot,temp6(1:3))
  +       y%ReactionForce%Moment(1:3,1) = MATMUL(p%GlbRot,temp6(4:6))
      ENDIF
      DO i=1,p%node_total
          temp_id = (i-1)*p%dof_node
          temp6(:) = 0.0D0
          temp6(:) = temp_Force(temp_id+1:temp_id+6)
  -       temp6(:) = MATMUL(TRANSPOSE(MoTens),temp6)
  -       y%BldForce%Force(1:3,i) = temp6(1:3)
  -       y%BldForce%Moment(1:3,i) = temp6(4:6)
  -   ENDDO
  -
  -!   CALL BD_CopyOtherState(OS_tmp, OtherState, MESH_NEWCOPY, ErrStat, ErrMsg)
  -!   OtherState%Acc(1:3) = u_tmp%RootMotion%TranslationAcc(1:3,1)
  -!   OtherState%Acc(4:6) = u_tmp%RootMotion%RotationAcc(1:3,1)
  -!   OtherState%Xcc(:) = OtherState%Acc(:)
  +       y%BldForce%Force(1:3,i) = MATMUL(p%GlbRot,temp6(1:3))
  +       y%BldForce%Moment(1:3,i) = MATMUL(p%GlbRot,temp6(4:6))
  +   ENDDO
  +
   
      CALL BD_DestroyInput(u_tmp, ErrStat, ErrMsg)
      CALL BD_DestroyContState(x_tmp, ErrStat, ErrMsg )
  
positive.code:    CALL BD_CopyOtherState(OtherState, OS_tmp, MESH_NEWCOPY, ErrStat, ErrMsg)
     CALL BD_CopyInput(u, u_tmp, MESH_NEWCOPY, ErrStat, ErrMsg)
  
     CALL BD_MotionTensor(p%GlbRot,p%GlbPos,MoTens,0)
     DO i=1,p%elem_total
         DO j=1,p%node_elem
             temp_id = ((i-1)*(p%node_elem-1)+j-1)*p%dof_node
  
negative.code:    CALL BD_CrvExtractCrv(InitInp%GlbRot,temp_glbrot,ErrStat2,ErrMsg2)
        CALL SetErrStat( ErrStat2, ErrMsg2, ErrStat, ErrMsg, RoutineName )
     CALL AllocAry(p%GlbPos,3,'Global position vector',ErrStat2,ErrMsg2)
        CALL SetErrStat( ErrStat2, ErrMsg2, ErrStat, ErrMsg, RoutineName )
     CALL AllocAry(p%GlbRot,3,3,'Global rotation tensor',ErrStat2,ErrMsg2)
        CALL SetErrStat( ErrStat2, ErrMsg2, ErrStat, ErrMsg, RoutineName )
     p%GlbPos(:) = 0.0D0
     p%GlbRot(:,:) = 0.0D0
     p%GlbPos(1:3) = InitInp%GlbPos(1:3)
     p%GlbRot(1:3,1:3) = InitInp%GlbRot(1:3,1:3)
  
     CALL AllocAry(SP_Coef,InputFileData%kp_total-1,4,4,'Spline coefficient matrix',ErrStat2,ErrMsg2)
        CALL SetErrStat( ErrStat2, ErrMsg2, ErrStat, ErrMsg, RoutineName )
     SP_Coef(:,:,:) = 0.0D0
  
     temp_id = 0
  