C       Nc - number of centers, Np - number of dm particles
C       Nrad - number of shells for halos
C
      PARAMETER (Nm  =-1, Nb=260)  ! define size of linker mesh
      PARAMETER (Nc  =900000)            ! Max number of initial centers
      PARAMETER (Ncc  =Nc/15)    
C      PARAMETER (Np  =INT(1.05*NROW**3)) ! Max Nparticles
      PARAMETER (Np  = Nmaxpart) ! Max Nparticles, make extra 
      PARAMETER (Nrad=45)              ! Number of shells around a halo
      PARAMETER (Nrad_lim =45)
      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),
     +             iRadius(Nc) 
      COMMON /GRD/ Label(Nm:Nb,Nm:Nb,Nm:Nb),Lst(Np)
      COMMON /MXC/ Wxc(Ncc,Nrad),Wyc(Ncc,Nrad),Wzc(Ncc,Nrad),
     +             Vrmc(Ncc,Nrad),Nbc(Ncc,Nrad),Rrc(Ncc,Nrad),
     +              NbcM(Ncc,Nrad),NbcG(Ncc,Nrad),NbcS(Ncc,Nrad)
     +              ,Nrad0(Ncc,Nrad),Wrad0(Ncc,Nrad),Wrad2(Ncc,Nrad)
c     +              ,Agx(Ncc,Nrad),Agy(Ncc,Nrad),Agz(Ncc,Nrad)
c      COMMON /PNT/ Xc,Yc,Zc,a0,Radius,RadSr1,RadSr2,
      COMMON /PNT/ RadSr1,RadSr2,
     +             Rad(Nrad),Rad2(Nrad),Cell
c      COMMON /STS/ Vx1(Nrad),Vy1(Nrad),Vz1(Nrad),Vxx(Nrad),Nob1(Nrad),
c     .             Rob1(Nrad),NobM(Nrad),NobG(Nrad),NobS(Nrad)
c   	COMMON / ESCP/ Vescape2(Nrad)