anchor.diff:  
   
      END SUBROUTINE BD_InputSolve
  +
   
      SUBROUTINE BD_InitialCondition(u,p,x,ErrStat,ErrMsg)
   
  
positive.code: !  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)
  !--------------------------------------------
  ! Compute initial condition given root motion
  !--------------------------------------------
  
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
  