anchor.diff:     END IF
      D_dragConst = 0.0
      
  -   ! Loop over nodes to create all loads on the resulting markers except for the buoyancy loads
  +   ! Loop over nodes to create all loads on the resulting markers except for the buoyancy loads  <---@mhall: what does "create all loads" mean?
      ! For the buoyancy loads, loop over the elements and then apply one half of the resulting value
      ! to each of the interior element nodes but the full value to an end node.  This means that an internal member node will receive 1/2 of its
      ! load from element A and 1/2 from element B.  If it is the end of a member it will simply receive
  -   ! the element A load.
  +   ! the element A load.                                                                         <---@mhall: shouldn't end nodes only receive half too?
      
      count = 1 
      
  
positive.code:                                   InitOut%Morison_Rad,  ErrStat, ErrMsg )
                                      
                                   
  
        
                                   
       IF ( ErrStat > ErrID_None ) RETURN
  
negative.code: !END SUBROUTINE SplitElementsForWtr
  !====================================================================================================
  SUBROUTINE SubdivideMembers( numNodes, nodes, numElements, elements, ErrStat, ErrMsg )
     ! This subdivides all of the (already-split) Morison elements according to each element's maximum desired 
     ! element length (MDivSize), adding resuling new nodes and elements to the respective master arrays.
     
     INTEGER,                  INTENT ( INOUT )  :: numNodes
     TYPE(Morison_NodeType),   INTENT ( INOUT )  :: nodes(:)   
     INTEGER,                  INTENT ( INOUT )  :: numElements
  