C       Nc - number of centers, Np - number of dm particles
C       Nrad - number of shells for halos
C
      PARAMETER (Nc  =15000)            ! Max number of initial centers
      PARAMETER (Np  = Nmaxpart) ! Max Nparticles, make extra 
      PARAMETER (Nrad=75)               ! Number of shells around a halo
      PARAMETER (Nrad_lim =75)
      COMMON /CAT/ Xp(Np),Yp(Np),Zp(Np),VXp(Np),VYp(Np),VZp(Np)
      COMMON /MXM/ Xm(Nc),Ym(Nc),Zm(Nc),Amc(Nc),Rmc(Nc),Vrm(Nc),
     +              Vcmax(Nc),Vmx(Nc),Vmy(Nc),Vmz(Nc),
     +               m1(Nc),m2(Nc),m3(Nc),m4(Nc),
     +               aJ(Nc),aLambda1(Nc),aLambda2(Nc)
      COMMON /MXC/ Wxc(Nc,Nrad),Wyc(Nc,Nrad),Wzc(Nc,Nrad),
     +             Vrmc(Nc,Nrad),Nbc(Nc,Nrad),Rrc(Nc,Nrad),
     +              NbcM(Nc,Nrad),NbcG(Nc,Nrad),NbcS(Nc,Nrad)
     +              ,Agx(Nc,Nrad),Agy(Nc,Nrad),Agz(Nc,Nrad)