anchor.diff:     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
  
positive.code:    CALL ReadCom(UnIn,DvrInputFile,'Section Header: Initial Velocity Parameter',ErrStat2,ErrMsg2,UnEc)
        CALL SetErrStat( ErrStat2, ErrMsg2, ErrStat, ErrMsg, RoutineName )
     InitInputData%RootVel(:)   = 0.0D0
     CALL ReadVar(UnIn,DvrInputFile,InitInputData%RootVel(1),"InitInputData%IniRootVel(1)", "velocity vector X",ErrStat2,ErrMsg2,UnEc)
        CALL SetErrStat( ErrStat2, ErrMsg2, ErrStat, ErrMsg, RoutineName )  
     CALL ReadVar(UnIn,DvrInputFile,InitInputData%RootVel(2),"InitInputData%IniRootVel(2)", "velocity vector Y",ErrStat2,ErrMsg2,UnEc)
        CALL SetErrStat( ErrStat2, ErrMsg2, ErrStat, ErrMsg, RoutineName )
     CALL ReadVar(UnIn,DvrInputFile,InitInputData%RootVel(3),"InitInputData%IniRootVel(3)", "velocity vector Z",ErrStat2,ErrMsg2,UnEc)
        CALL SetErrStat( ErrStat2, ErrMsg2, ErrStat, ErrMsg, RoutineName )
     CALL ReadVar(UnIn,DvrInputFile,InitInputData%RootVel(4),"InitInputData%IniRootVel(4)", "angular velocity vector X",ErrStat2,ErrMsg2,UnEc)
        CALL SetErrStat( ErrStat2, ErrMsg2, ErrStat, ErrMsg, RoutineName ) 
     CALL ReadVar(UnIn,DvrInputFile,InitInputData%RootVel(5),"InitInputData%IniRootVel(5)", "angular velocity vector Y",ErrStat2,ErrMsg2,UnEc)
        CALL SetErrStat( ErrStat2, ErrMsg2, ErrStat, ErrMsg, RoutineName ) 
     CALL ReadVar(UnIn,DvrInputFile,InitInputData%RootVel(6),"InitInputData%IniRootVel(6)", "angular velocity vector Z",ErrStat2,ErrMsg2,UnEc)
        CALL SetErrStat( ErrStat2, ErrMsg2, ErrStat, ErrMsg, RoutineName ) 
     if (ErrStat >= AbortErrLev) then
         call cleanup()
         return
  
negative.code: !     reading of inputs.
  !------------------------------------------------------------------------------------
     USE BeamDyn
     USE BeamDyn_Subs
     USE BeamDyn_Types
     USE NWTC_Library
  
  