anchor.diff:     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 =''
   
  
positive.code:    ! Implementation 
     integer(IntKi), parameter :: iNWStart=2 !< Index in r%NW where the near wake start (if >1 then the Wing panels are included in r_NW)
     integer(IntKi), parameter :: FWnSpan=1  !< Number of spanwise far wake panels ! TODO make it an input later
     logical       , parameter :: DEV_VERSION=.FALSE.
  contains
  
  !==========================================================================
  
negative.code:    subroutine PackConvectingPoints()
        ! Counting total number of control points that convects
        nCPs = CountCPs(p, m%nNW, nFWEff)
        m%CPs=-999.9_ReKi
        ! Packing
        iHeadP=1
  