C---------------------------- parameters for PMhalosM --------
C      iDebug      = 1 gives long output for debuging
C                  = 0 makes output shorter
C      iHelp       = 1 gives short description of parameters
C                            and creates a file, which stores
C                            all your answers. Use it for
C                            rerunning the code.
C
C
C       Nc - number of centers, Np - number of dm particles
C       Nrad - number of shells for halos
C
      PARAMETER (iDebug  =     0)  ! debugging parameteter
      PARAMETER (iHelp   =     0)  ! help
      PARAMETER (Nm      =    -20)  ! define size of linker mesh
      PARAMETER (Nb      =    260)  ! define size of linker mesh
      PARAMETER (Nc      =   2000000)  ! Max number of initial centers
      PARAMETER (Nrad    =     45)  ! Number of shells for halo profile
      PARAMETER (Cell    =    1.0)  ! Cell-size for Linker List in grid units
      PARAMETER (Ncc     =  Nc/5)    
      PARAMETER (Np      = Nmaxpart) ! Max Nparticles
      PARAMETER (Nrad_lim =Nrad)
      COMMON /CAT/ Wp(6,Np)
      COMMON /MXM/ Xm(Nc),Ym(Nc),Zm(Nc),Amc(Nc),Rmc(Nc),Vrm(Nc),
     +             iRadius(Nc),RmaxV(Nc) 
      COMMON /GRD/ Label(Nm:Nb,Nm:Nb,Nm:Nb),Lst(Np)
      COMMON /MXC/ Wxc(Ncc,0:Nrad),Wyc(Ncc,0:Nrad),Wzc(Ncc,0:Nrad),
     +       Vrmc(Ncc,0:Nrad),Nbc(Ncc,0:Nrad),Rrc(Ncc,0:Nrad),
     +       NbcM(Ncc,0:Nrad),NbcG(Ncc,0:Nrad),NbcS(Ncc,0:Nrad)
c     +       Nrad0(Ncc,0:Nrad),Wrad0(Ncc,0:Nrad),Wrad2(Ncc,0:Nrad)
c     +      ,Agx(Ncc,0:Nrad),Agy(Ncc,0:Nrad),Agz(Ncc,0:Nrad)
      COMMON /PNT/ RadSr1,RadSr2,
     +             Rad(Nrad),Rad2(Nrad)