C   
C
      PARAMETER (Np  = Nmaxpart) ! Max Nparticles, make extra 
      PARAMETER (Nrad=80)              ! Number of shells
      PARAMETER (Nrad2=2*Nrad+1) ! length of an array for rad
      PARAMETER (Nzp =20)              ! Number of bins in Z
      PARAMETER (Nbp =18)              ! Number of bins in phi
      PARAMETER (Nbt =Nrad2*Nbp)              !  length of an array for phi
      COMMON /CAT/ Xp(Np),Yp(Np),Zp(Np),VXp(Np),VYp(Np),VZp(Np)
      COMMON /CATB/ GXp(Np),GYp(Np),GZp(Np),Fip(Np),!Jp(Np),
     +                        XL(Np),YL(Np),ZL(Np),
     +                         Idisk,Ihalo
      COMMON /DISK/ Nradd(-Nrad:Nrad),Vfi(-Nrad:Nrad),Vrad(-Nrad:Nrad),
     +             Vzz(-Nrad:Nrad),Vfi2(-Nrad:Nrad),
     +             Vrad2(-Nrad:Nrad),Vzz2(-Nrad:Nrad),Radd(-Nrad:Nrad),
     +             Nzz(Nzp,Nrad),Wradd(-Nrad:Nrad)
      COMMON /DISKp/ Nradp(-Nrad:Nrad),Vfip(-Nrad:Nrad),
     +             Vradp(-Nrad:Nrad),Vzzp(-Nrad:Nrad),Vfi2p(-Nrad:Nrad),
     +             Vrad2p(-Nrad:Nrad),Vzz2p(-Nrad:Nrad),Radp(-Nrad:Nrad)
      COMMON /BAR/ Nradb(-Nrad:Nrad,Nbp),Vfib(-Nrad:Nrad,Nbp),
     +             Vradb(-Nrad:Nrad,Nbp),Vzzb(-Nrad:Nrad,Nbp),
     +             Vfi2b(-Nrad:Nrad,Nbp),Vrad2b(-Nrad:Nrad,Nbp),
     +             Vzz2b(-Nrad:Nrad,Nbp),Radb(-Nrad:Nrad,Nbp)
      COMMON /HALO/ Nradh(-Nrad:Nrad),Vfih(-Nrad:Nrad),
     +             Vradh(-Nrad:Nrad),Vfih2(-Nrad:Nrad),
     +             Vradh2(-Nrad:Nrad),Radh(-Nrad:Nrad)