anchor.diff:     subroutine PackConvectingPoints()
         ! Counting total number of control points that convects
         nCPs = CountCPs(p, m%nNW, nFWEff)
  -      m%Uind=0.0_ReKi ! very important due to side effects of ui_seg
         m%CPs=-999.9_ReKi
         ! Packing
         iHeadP=1
  
positive.code:    if (DEV_VERSION) then
        print'(A,I0,A,I0,A,I0)','Convection - nSeg:',nSeg,' - nSegP:',nSegP, ' - nCPs:',nCPs
     endif
     call ui_seg( 1, nCPs, m%CPs, 1, nSeg, nSeg, nSegP, m%SegPoints, m%SegConnct, m%SegGamma, p%RegFunction, m%SegEpsilon, m%Uind)
     call UnPackInducedVelocity()
  
  contains
  
negative.code:    integer(IntKi) :: iW, nSeg, nSegP, nCPs, iHeadP
     integer(IntKi) :: nFWEff  ! Number of farwake panels that are free at current tmie step
     logical        :: bMirror ! True if we mirror the vorticity wrt ground
     ! TODO new options
     integer(IntKi) :: nPart
     integer(IntKi) :: nPartPerSeg
     integer(IntKi) :: UIMethod
     integer(IntKi) :: RegFunctionPart
     real(ReKi) :: BranchFactor, BranchSmall
     type(T_Tree)   :: Tree
     real(ReKi), dimension(:,:), allocatable :: PartPoints !< Particle points
     real(ReKi), dimension(:,:), allocatable :: PartAlpha  !< Particle circulation
     real(ReKi), dimension(:)  , allocatable :: PartEpsilon !< Regularization parameter
     ErrStat= ErrID_None
     ErrMsg =''
  
  