anchor.diff:  
      ! Calculate root translational and angular accelerations
      u%RootMotion%TranslationAcc(:,:) = 0.0D0
  -!   u%RootMotion%RotationAcc(:,:) = 0.0D0
  +   u%RootMotion%RotationAcc(:,:) = 0.0D0
   !   u%RootMotion%TranslationAcc(:,1) = MATMUL(BD_Tilde(u%RootMotion%RotationVel(:,1)), &
   !               MATMUL(BD_Tilde(u%RootMotion%RotationVel(:,1)),temp_rr))
      ! END Calculate root translational and angular accelerations
  
positive.code: SUBROUTINE BD_InputSolve( t, u,  p, ErrStat, ErrMsg)
   
     USE BeamDyn
     USE BeamDyn_Types
  
     REAL(DbKi),                     INTENT(IN   ):: t
  
negative.code:    u%DistrLoad%Moment(:,:) = 0.0D0
  
  !   DO i=1,p%ngp*p%elem_total+2
  !       u%DistrLoad%Force(1,i) = 1.0D+04
  !   ENDDO
  
  END SUBROUTINE BD_InputSolve
  