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:    !DO n_t_global = 0, 1000
  
  
        CALL BDyn_InputSolve( BDyn_Input(1), BDyn_Output(1), BDyn_Parameter, ErrStat, ErrMsg)
  
  
  !     CALL BeamDyn_CalcOutput( t_global, BDyn_Input(1), BDyn_Parameter, BDyn_ContinuousState, BDyn_DiscreteState, &
  
negative.code: !      WRITE(*,*) t_global, BDyn_ContinuousState%q(15)
        WRITE(QiDisUnit,6000) t_global,BDyn_ContinuousState%q(BDyn_Parameter%dof_total-5),BDyn_ContinuousState%q(BDyn_Parameter%dof_total-4),&
                             &BDyn_ContinuousState%q(BDyn_Parameter%dof_total-3),BDyn_ContinuousState%q(BDyn_Parameter%dof_total-2),&
                             &BDyn_ContinuousState%q(BDyn_Parameter%dof_total-1),BDyn_ContinuousState%q(BDyn_Parameter%dof_total)
  
        ! Save all final variables 
  
  