anchor.diff:     ErrMsg  = ''
   
      ! gather point forces and line forces
  -   tmp_vector = 0.     
   
      ! Point mesh: Force 
      u%PointMesh%TranslationDisp(:,1)  = tmp_vector
  -   u%PointMesh%TranslationVel(:,1)  = tmp_vector
  -   u%PointMesh%TranslationAcc(:,1)  = tmp_vector
  +   u%PointMesh%TranslationVel(:,1)   = tmp_vector
  +   u%PointMesh%TranslationAcc(:,1)   = tmp_vector
  +
  +   u%PointMesh%TranslationDisp(3,1)  = +0.1*sin(t+dt)
  +   u%PointMesh%TranslationVel(3,1)   = -0.1*cos(t+dt)
  +   u%PointMesh%TranslationAcc(3,1)   = -0.1*sin(t+dt)
  +
  +   ut(1) = t+dt
  +
  +   u%PointMesh%TranslationDisp(:,2)  = tmp_vector
  +   u%PointMesh%TranslationVel(:,2)   = tmp_vector
  +   u%PointMesh%TranslationAcc(:,2)   = tmp_vector
  +
  +   u%PointMesh%TranslationDisp(3,2)  = +0.1*sin(t+0.5*dt)
  +   u%PointMesh%TranslationVel(3,2)   = -0.1*cos(t+0.5*dt)
  +   u%PointMesh%TranslationAcc(3,2)   = -0.1*sin(t+0.5*dt)
  +   ut(2) = t+0.5*dt
  +
  +   u%PointMesh%TranslationDisp(:,3)  = tmp_vector
  +   u%PointMesh%TranslationVel(:,3)  = tmp_vector
  +   u%PointMesh%TranslationAcc(:,3)  = tmp_vector
  +
  +   u%PointMesh%TranslationDisp(3,3)  = +0.1*sin(t)
  +   u%PointMesh%TranslationVel(3,3)   = -0.1*cos(t)
  +   u%PointMesh%TranslationAcc(3,3)   = -0.1*sin(t)
  +   ut(3) = t
   
   END SUBROUTINE BDyn_InputSolve
   !----------------------------------------------------------------------------------------------------------------------------------
  
positive.code: 
  
  
  SUBROUTINE BDyn_InputSolve( u, y, p, ErrStat, ErrMsg)
   
     USE BeamDyn
     USE BeamDyn_Types
  
     ! Module1 Derived-types variables; see Registry_Module1.txt for details
  
     TYPE(BDyn_InputType),           INTENT(INOUT) :: u
     TYPE(BDyn_OutputType),          INTENT(IN   ) :: y
     TYPE(BDyn_ParameterType),       INTENT(IN   ) :: p
  
  
     INTEGER(IntKi),                 INTENT(  OUT)  :: ErrStat     ! Error status of the operation
  
negative.code:    !DO n_t_global = 0, 1000
  
  
  
  
  !     CALL BeamDyn_CalcOutput( t_global, BDyn_Input(1), BDyn_Parameter, BDyn_ContinuousState, BDyn_DiscreteState, &
  