anchor.diff:  
   
   
  -SUBROUTINE BDyn_InputSolve( u, y, p, ErrStat, ErrMsg)
  +!SUBROUTINE BDyn_InputSolve( u, y, p, ErrStat, ErrMsg)
  +SUBROUTINE BDyn_InputSolve( t, dt, u, ut, ErrStat, ErrMsg)
    
      USE BeamDyn
      USE BeamDyn_Types
   
      ! Module1 Derived-types variables; see Registry_Module1.txt for details
   
  -   TYPE(BDyn_InputType),           INTENT(INOUT) :: u
  -   TYPE(BDyn_OutputType),          INTENT(IN   ) :: y
  -   TYPE(BDyn_ParameterType),       INTENT(IN   ) :: p
  +   REAL(ReKi),                     INTENT(IN   ) :: t
  +   REAL(ReKi),                     INTENT(IN   ) :: dt
  +   TYPE(BDyn_InputType),           INTENT(INOUT) :: u(:)
  +   REAL(DbKi),                     INTENT(INOUT) :: ut(:)
   
   
      INTEGER(IntKi),                 INTENT(  OUT)  :: ErrStat     ! Error status of the operation
  
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:    !DO n_t_global = 0, 1000
  
  
  
  
  !     CALL BeamDyn_CalcOutput( t_global, BDyn_Input(1), BDyn_Parameter, BDyn_ContinuousState, BDyn_DiscreteState, &
  