anchor.diff:     !DO n_t_global = 0, 1000
   
   
  -      CALL BDyn_InputSolve( BDyn_Input(1), BDyn_Output(1), BDyn_Parameter, ErrStat, ErrMsg)
   
   
   !     CALL BeamDyn_CalcOutput( t_global, BDyn_Input(1), BDyn_Parameter, BDyn_ContinuousState, BDyn_DiscreteState, &
  
positive.code: 
        ! extrapolate inputs and outputs to t + dt; will only be used by modules with an implicit dependence on input data.
  
        CALL BDyn_Input_ExtrapInterp(BDyn_Input, BDyn_InputTimes, u1, t_global + dt_global, ErrStat, ErrMsg)
  
  !     CALL BDyn_Output_ExtrapInterp(BDyn_Output, BDyn_OutputTimes, y1, t_global + dt_global, ErrStat, ErrMsg)
  
        ! Shift "window" of the Mod1_Input and Mod1_Output
  
        DO i = BDyn_interp_order, 1, -1
           CALL BDyn_CopyInput (BDyn_Input(i),  BDyn_Input(i+1), MESH_UPDATECOPY, Errstat, ErrMsg)
           CALL BDyn_CopyOutput (BDyn_Output(i),  BDyn_Output(i+1),  MESH_UPDATECOPY, Errstat, ErrMsg)
           BDyn_InputTimes(i+1) = BDyn_InputTimes(i)
           BDyn_OutputTimes(i+1) = BDyn_OutputTimes(i)
        ENDDO
  
        CALL BDyn_CopyInput (u1,  BDyn_Input(1),  MESH_UPDATECOPY, Errstat, ErrMsg)
        CALL BDyn_CopyOutput (y1,  BDyn_Output(1),  MESH_UPDATECOPY, Errstat, ErrMsg)
        BDyn_InputTimes(1) = t_global + dt_global
        BDyn_OutputTimes(1) = t_global + dt_global
  
        ! Shift "window" of the Mod1_Input and Mod1_Output
  
  
negative.code:    ErrMsg  = ''
  
     ! gather point forces and line forces
  
     ! Point mesh: Force 
     u%PointMesh%TranslationDisp(:,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
  !----------------------------------------------------------------------------------------------------------------------------------
  