C---------------------------- parameters for PMselect.f --------
C
C       Np    - max number of halos
C       Nrmax - mx number of shells for a halo
C
      PARAMETER (Np    = 300 000)  
      PARAMETER (Nrmax =   150) !Number of shells for  profile
      PARAMETER (Riso  =    2.) ! radius for isolated halo in Rvirs
      PARAMETER (aMiso =    2.) ! mass ratio for isolation
      PARAMETER (errLimit = 0.6) ! error of NFW fit
      COMMON /CAT/ Xp(Np),Yp(Np),Zp(Np),VXp(Np),VYp(Np),VZp(Np)
      COMMON /ADD/ aM(Np),Radh(Np),Vrms(Np),Vcirc(Np),
     +             Npart(Np),Rmax(Np),Conc(Np),Nlines(Np),errFit(Np)
      COMMON /PRF/ Rpr(Nrmax,Np),Npr(Nrmax,Np),
     +             Vcpr(Nrmax,Np),RHOpr(Nrmax,Np),Vrpr(Nrmax,Np),
     +             bM(Nrmax,Np),Dpr(Nrmax,Np),Vdpr(Nrmax,Np)
      COMMON/SEL/ aMass, dM, X0, Y0, Z0, dR, Cc, dCc