anchor.diff:       CALL BD_CalcOutput( t_global, BD_Input(2), BD_Parameter, BD_ContinuousState, BD_DiscreteState, &
                                BD_ConstraintState, &
                                BD_OtherState,  BD_Output(2), ErrStat, ErrMsg)
  -!     IF(BD_Parameter%analysis_type .EQ. 2 .AND. n_t_global .EQ. 0) THEN
  -!         CALL BD_InitAcc( t_global, BD_Input(1), BD_Parameter, &
  -!               BD_ContinuousState,BD_OtherState,ErrStat,ErrMsg)
  -!!      WRITE(*,*) 'Initial Acc'
  -!!      WRITE(*,*) BD_OtherState%acc(:)
  -!!      WRITE(*,*) 'Initial Xcc'
  -!!      WRITE(*,*) BD_OtherState%xcc(:)
  -!     ENDIF
  -WRITE(*,*) 'TEST'
  +
        IF(BD_Parameter%analysis_type .EQ. 1 .AND. n_t_global .EQ. 1) EXIT 
   
        CALL BD_UpdateStates( t_global, n_t_global, BD_Input, BD_InputTimes, BD_Parameter, &
  
positive.code:    ! Initialize some variables:
     ErrStat = ErrID_None
     ErrMsg  = ""
     
     CALL GetNewUnit(UnIn,ErrStat2,ErrMsg2)
        CALL SetErrStat( ErrStat2, ErrMsg2, ErrStat, ErrMsg, RoutineName )
  
negative.code:    ! Initialization of glue-code time-step variables
     ! -------------------------------------------------------------------------
     CALL GET_COMMAND_ARGUMENT(1,DvrInputFile)
     CALL GetRoot(DvrInputFile,RootName)
     CALL BD_ReadDvrFile(DvrInputFile,t_initial,t_final,dt_global,BD_InitInput,&
                            ErrStat,ErrMsg)
     BD_InitInput%RootName  = TRIM(BD_Initinput%InputFile)
  