anchor.diff:        ! --- 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
  
positive.code:    CASE (3) ! Member-based model: coefficients set using member-specific coefficient tables
         do i = 1, member%NElements + 1
           ! Pull member  end-node data from the tables and then linearly interpolate it onto the interior member nodes    
           s = (i-1) / member%NElements
           if ( member%tMG(i) > 0.0_ReKi ) then
              member%Cd    (i) = CoefMembers(MmbrCoefIDIndx)%MemberCdMG1*(1-s) + CoefMembers(MmbrCoefIDIndx)%MemberCdMG2*s
              member%Ca    (i) = CoefMembers(MmbrCoefIDIndx)%MemberCaMG1*(1-s) + CoefMembers(MmbrCoefIDIndx)%MemberCaMG2*s
  
negative.code:             exit
           end if
        end do
     else
        member%i_floor = 0 ! lower end is at or above the seabed
     end if
  
  
  