anchor.diff:        DO J=1,p%NMOutputs     
            DO I=1,p%MOutLst(J)%NOutLoc   
            
  -            ! These indices are in the DistribMesh index system, not the overall nodes index system, so distribToNodeIndx() mapping needs to be performed if you want to index into the nodes array or wave kinematics arrays
  -            ! But, all of the D_* arrays are already using the DistribMesh index system, so we are OK for those
  -            
  +             
               m1 = p%MOutLst(J)%Marker1(I)
               m2 = p%MOutLst(J)%Marker2(I)
               s  = p%MOutLst(J)%s      (I)
  
positive.code:    ErrMsg  = ''
     
     ! This only works because the member nodes are equally spaced
     
     i1 = floor(member%NElements*d)  ! 0<= d <= 1.0
     
     if ( i1 == member%NElements ) then
        ! special case when d = 1.0_ReKi
        i2 = i1 + 1            
        s  = 1.0_ReKi
     else 
        s  = member%NElements*d - i1  ! here i1 is zero-based
        i1 = i1 + 1                   ! Shift to one-based indices
        i2 = i1 + 1
     end if
     
  
negative.code:             ! Need to search mesh for the two markers which surround the requested output location and then store those marker indices and compute the
              ! scale factor based on how far they are from the requested output location.
              ! Since this is being done on markers and not nodes, the subroutine must be called after the Morison_Init() subroutine is called
  
              CALL GetNeighboringNodes(p%Members(memberIndx), p%MOutLst(I)%NodeLocs(J),  m1, m2, s, ErrStat, ErrMsg)
                    
              p%MOutLst(I)%Marker1(J) = m1
              p%MOutLst(I)%Marker2(J) = m2 ! The 2nd marker indx which is used to
              p%MOutLst(I)%s(J)       = s ! linear interpolation factor     
  