anchor.diff:  !  This way, when RK4 is called using ExtrapInterp, it will grab the EXACT answers that you defined at the time
   !  step endpionts and midpoint.
   
  -      CALL BD_InputSolve( t_global,                   BD_Input(1), BD_InputTimes(1), BD_Parameter, BD_ContinuousState, ErrStat, ErrMsg)
  -      CALL BD_InputSolve( t_global + dt_global,       BD_Input(2), BD_InputTimes(2), BD_Parameter, BD_ContinuousState, ErrStat, ErrMsg)
  -      CALL BD_InputSolve( t_global + 2.0D0*dt_global, BD_Input(3), BD_InputTimes(3), BD_Parameter, BD_ContinuousState, ErrStat, ErrMsg)
  +      CALL BD_InputSolve( t_global,                   BD_Input(1), BD_InputTimes(1), BD_Parameter, ErrStat, ErrMsg)
  +      CALL BD_InputSolve( t_global + dt_global,       BD_Input(2), BD_InputTimes(2), BD_Parameter, ErrStat, ErrMsg)
  +      CALL BD_InputSolve( t_global + 2.0D0*dt_global, BD_Input(3), BD_InputTimes(3), BD_Parameter, ErrStat, ErrMsg)
   !--------------------------------------------
   ! Compute initial condition given root motion
   !--------------------------------------------
  
positive.code: 
  
  
     SUBROUTINE BD_InputSolve( t, u, ut, p, x, ErrStat, ErrMsg)
   
     USE BeamDyn_SP
     USE BeamDyn_Types
  
negative.code: 
  
     END SUBROUTINE BD_InputSolve
  
  
     SUBROUTINE BD_InitialCondition(u,p,x,ErrStat,ErrMsg)
  
  