anchor.diff:  
      ! calculate h_floor if seabed-piercing
      member%h_floor = 0.0_ReKi
  -   member%i_floor = 0
  +   member%i_floor = member%NElements+1  ! Default to entire member is below the seabed
      if (Za < -WtrDepth) then
         do i= 2, member%NElements+1
            Za = InitInp%Nodes(member%NodeIndx(i))%Position(3)
  
positive.code:             exit
           end if
        end do
     end if
  
  
  
negative.code:       ! --- external buoyancy loads: ends ---
  
        if ( .not. mem%PropPot ) then
           pos1    = u%Mesh%TranslationDisp(:, mem%NodeIndx(1))   + u%Mesh%Position(:, mem%NodeIndx(1))
           pos2    = u%Mesh%TranslationDisp(:, mem%NodeIndx(N+1)) + u%Mesh%Position(:, mem%NodeIndx(N+1))
           z1 = pos1(3)
           z2 = pos2(3)
           if (mem%i_floor == 0) then  ! both ends above or at seabed
              if (z2<= 0.0_ReKi) then
                 ! Compute loads on both ends
  