anchor.diff:     u%RootMotion%RotationVel(:,:) = 0.0D0
      u%RootMotion%RotationVel(1,1) = IniVelo(5)
      u%RootMotion%RotationVel(2,1) = IniVelo(6)
  -   u%RootMotion%RotationVel(3,1) = IniVelo(1)
  +   u%RootMotion%RotationVel(3,1) = IniVelo(4)
      u%RootMotion%TranslationVel(:,:) = 0.0D0
      u%RootMotion%TranslationVel(:,1) = MATMUL(BD_Tilde(real(u%RootMotion%RotationVel(:,1),BDKi)),temp_rr)
      ! END Calculate root translational and angular velocities
  
positive.code:    CALL BD_CrvMatrixR(temp_vec,temp33,ErrStat,ErrMsg)
     temp_rr(:) = MATMUL(u%RootMotion%Orientation(:,:,1),temp_r0)
     u%RootMotion%Orientation(:,:,1) = TRANSPOSE(u%RootMotion%Orientation(:,:,1))
  !   IF(t .GT. 0.2) THEN
  !       temp3(:) = 0.0D0
  !       temp3(3) = 0.0175
  !       u%RootMotion%Orientation(:,:,1) = EulerConstruct(temp3)
  !   ENDIF
     u%RootMotion%TranslationDisp(:,:)  = 0.0D0
     u%RootMotion%TranslationDisp(:,1) = temp_rr(:) - temp_r0(:)
     ! END Calculate root displacements and rotations
  
negative.code: !     reading of inputs.
  !------------------------------------------------------------------------------------
     USE BeamDyn
     USE BeamDyn_Subs
     USE BeamDyn_Types
     USE NWTC_Library
  
  