anchor.diff:           CALL BDyn_CopyConstrState (BDyn_ConstraintState, BDyn_ConstraintState_pred, 0, Errstat, ErrMsg)
   
            CALL BDyn_CopyDiscState   (BDyn_DiscreteState,   BDyn_DiscreteState_pred,   0, Errstat, ErrMsg)
  -
  +         
            CALL BeamDyn_UpdateStates( t_global, n_t_global, BDyn_Input, BDyn_InputTimes, BDyn_Parameter, &
                                      BDyn_ContinuousState_pred, &
                                      BDyn_DiscreteState_pred, BDyn_ConstraintState_pred, &
  
positive.code:    tmp_vector = 0.     
  
     ! Point mesh: Force 
     u%PointMesh%Force(:,1)  = tmp_vector
  
  
  END SUBROUTINE BDyn_InputSolve
  !----------------------------------------------------------------------------------------------------------------------------------
  
negative.code:    t_initial = 0.0D0
     t_final   = 2.0D0
  
     pc_max = 1  ! Number of predictor-corrector iterations; 1 corresponds to an explicit calculation where UpdateStates 
                 ! is called only once  per time step for each module; inputs and outputs are extrapolated in time and 
                 ! are available to modules that have an implicit dependence on other-module data
  
  