anchor.diff:           DO I=1,6
                           
               ! We are now combining the dynamic pressure term into the inertia term
  -            m%L_F_I(I,J) = InterpWrappedStpReal ( REAL(Time, SiKi), p%WaveTime(:), p%L_F_I(:,I,J), &
  -                                    m%LastIndWave, p%NStepWave + 1       ) 
  +            m%L_F_I(I,J) = InterpolateWithSlope(InterpolationSlope, m%LastIndWave, p%L_F_I(:,I,J))
               
               IF (I < 4 ) THEN
         
  
positive.code:       CHARACTER(*),                      INTENT(  OUT)  :: ErrMsg      !< Error message if ErrStat /= ErrID_None
  
        REAL(ReKi)                                        :: F_D(6), F_DP(6), D_F_I(3), kvec(3), v(3),  vf(3), vrel(3), vmag
        INTEGER                                           :: I, J, K, nodeIndx
        REAL(ReKi)                                        :: elementWaterState
        REAL(ReKi)                                        :: AllOuts(MaxMrsnOutputs)  ! TODO: think about adding to OtherState
        REAL(ReKi)                                        :: qdotdot(6) ,qdotdot2(3)     ! The structural acceleration of a mesh node
  
negative.code:       REAL(ReKi)                                        :: D_AM_M(6,6)
        REAL(ReKi)                                        :: nodeInWater
        REAL(ReKi)                                        :: D_dragConst     ! The distributed drag factor
        REAL(ReKi)                                        :: InterpolationSlope 
        
           ! Initialize ErrStat
           
        ErrStat = ErrID_None         
        ErrMsg  = ""               
        
        InterpolationSlope = GetInterpolationSlope(Time, p, m, IntWrapIndx)
        
           ! Compute outputs here:
        
  