c      IMPLICIT double precision (A-H,O-Z)
      PARAMETER (Nmax  =1000000)
      PARAMETER (Nmaxp  =30000000)
      PARAMETER (Nshells  = 300)
      PARAMETER (Nshell2  = 2*Nshells+2)
      PARAMETER (hubble = 100.d0)
      PARAMETER (Box      = 120.d0)
      PARAMETER ( aMasslimit = 1.d10)
      PARAMETER ( Radlimit = 1000.d0) ! radial shell in projection
      PARAMETER ( Radlimit2 =0.d0)
      PARAMETER ( RadLmax = 1.5d0)   ! isolation creterion
      PARAMETER ( Vlimit = 1000.d0)
      PARAMETER ( Dlimit1 = 0.0d0)    ! take satellites inside 3D shell
      PARAMETER ( Dlimit2 = 100.5d0)
      PARAMETER ( Vclimit = 30.d0)   ! do not take halos below this
      PARAMETER ( Vcl = 700.)       ! Vc limit on clusters
      PARAMETER ( Dcl = 5.)       ! Distance limit to clusters
      PARAMETER ( Vratio = 2.)
      PARAMETER ( pi   = 3.14159265)
      COMMON /DATA/ xc(Nmax),yc(Nmax),zc(Nmax),
     &                     vxc(Nmax),vyc(Nmax),vzc(Nmax),      
     &                     vc(Nmax),vrms(Nmax),aMass(Nmax),
     &                     Rvir(Nmax)
       COMMON /DATAP/ xp(Nmaxp),yp(Nmaxp),zp(Nmaxp),
     &                     vxp(Nmaxp),vyp(Nmaxp),vzp(Nmaxp)      
      COMMON /SET/Nh,Next,Np,NpBox
      Character Line*79,LineZ*45
      COMMON /iSol/Lox(Nmax),Loy(Nmax),Loz(Nmax)
      COMMON /ClusterList/ ListCl(Nmax),Nclusters