anchor.diff:  END FUNCTION InterpWrappedStpLogical ! ( XVal, XAry, YAry, Ind, AryLen )
   
   SUBROUTINE DistrBuoyancy( densWater, R, tMG, dRdz, Z, C, g, F_B  ) 
  +   ! This calculates the distributed buoyancy forces and moments on a given node
   
      REAL(ReKi),         INTENT ( IN    )  :: densWater
  -   REAL(ReKi),         INTENT ( IN    )  :: R
  -   REAL(ReKi),         INTENT ( IN    )  :: tMG
  -   REAL(ReKi),         INTENT ( IN    )  :: dRdz
  -   REAL(ReKi),         INTENT ( IN    )  :: Z
  +   REAL(ReKi),         INTENT ( IN    )  :: R                      ! Radius of node [m]
  +   REAL(ReKi),         INTENT ( IN    )  :: tMG                    ! Thickness of marine growth (adds to radius) [m]
  +   REAL(ReKi),         INTENT ( IN    )  :: dRdz                   ! Rate of change in radius with length at node [-]
  +   REAL(ReKi),         INTENT ( IN    )  :: Z                      ! z elevation of node [m] (not currently used)
      REAL(ReKi),         INTENT ( IN    )  :: C(3,3)
      REAL(ReKi),         INTENT ( IN    )  :: g
  -   REAL(ReKi),         INTENT (   OUT )  :: F_B(6)
  +   REAL(ReKi),         INTENT (   OUT )  :: F_B(6)                 ! Distributed force and moment vector [N/m and N-m/m]
      
      REAL(DbKi)                           :: Reff,ReffSq,ReffCub,f1,f2,f3
      
      REAL(DbKi) :: CC(3,3)
      
      CC    = REAL(C,DbKi)
  -   Reff  = REAL(R + tMG,DbKi)
  +   Reff  = REAL(R + tMG,DbKi)                                      ! Effective radius after adding marine growth
     
      
      
  
positive.code:    END IF
     D_dragConst = 0.0
     
     ! Loop over nodes to create all loads on the resulting markers except for the buoyancy loads
     ! 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.
     
     count = 1 
     
  
negative.code:             
            END IF      ! IF ( .NOT. nodes(I)%PropPot )
              
            ! This is always computed, but may be zereod out for any given timestep during the CalcOutput work   <---@mhall: what is this?
           CALL DistrDragConst( densWater, nodes(I)%Cd, nodes(I)%R, nodes(I)%tMG, D_dragConst(count) ) 
           
           IF ( nodes(I)%FillFlag ) THEN
  