anchor.diff:        TYPE(MD_ContinuousStateType)                   :: dxdt    ! time derivatives of continuous states (initialized in CalcContStateDeriv)
         INTEGER(IntKi)                                 :: I       ! counter
         INTEGER(IntKi)                                 :: J       ! counter
  -      REAL(ReKi)                                     :: t2      ! real version of t (double)
   
         INTEGER(IntKi)                                 :: ErrStat2   ! Error status of the operation
         CHARACTER(ErrMsgLen)                           :: ErrMsg2    ! Error message if ErrStat2 /= ErrID_None
   
   
         ! below updated to make sure outputs are current (based on provided x and u)  - similar to what's in UpdateStates
  -
  -      t2 = real(t, ReKi)
   
         ! go through fairleads and apply motions from driver
         DO I = 1, p%NFairs
  
positive.code: 
     !========================================================================================================
     SUBROUTINE TimeStep ( t, dtStep, u, utimes, p, x, xd, z, other, m, ErrStat, ErrMsg )
        REAL(ReKi)                     , INTENT(INOUT)      :: t
        REAL(ReKi)                     , INTENT(IN   )      :: dtStep     ! how long to advance the time for
        TYPE( MD_InputType )           , INTENT(INOUT)      :: u(:)       ! INTENT(IN   )
        REAL(DbKi)                     , INTENT(IN   )      :: utimes(:)  ! times corresponding to elements of u(:)?
  
negative.code: 
           ! step 1
  
           CALL MD_Input_ExtrapInterp(u, utimes, u_interp, t          , ErrStat, ErrMsg)   ! interpolate input mesh to correct time (t)
        
           CALL MD_CalcContStateDeriv( t, u_interp, p, x, xd, z, other, m, dxdt, ErrStat, ErrMsg )
           DO J = 1, Nx
  