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:    TYPE(BD_InputType),             INTENT(INOUT):: u
     REAL(DbKi),                     INTENT(INOUT):: ut
     TYPE(BD_ParameterType),         INTENT(IN   ):: p
     INTEGER(IntKi),                 INTENT(  OUT):: ErrStat     ! Error status of the operation
     CHARACTER(*),                   INTENT(  OUT):: ErrMsg      ! Error message if ErrStat /= ErrID_None
     ! local variables
  