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
      PARAMETER (bar_max = 6.5e-3)      ! max bar radius in Mpc
      PARAMETER (bar_min = 0.5e-30)      ! min bar radius in Mpc
      PARAMETER (dZ_max  = 3.5e-3)     ! max dZ for bar search
      PARAMETER (dZ_slice= 1.5e-30)     ! max dZ for slice
      PARAMETER (Nmlink   = 0)         ! size of linker list
      PARAMETER (Nblink   = 100)
      PARAMETER (Nsearch  = 70)
      PARAMETER (Nbarbins = 70)      ! number of bins for bar profile
      PARAMETER (pi =3.14159265) 
      COMMON /DISK/ Vfi(-Nrad:Nrad),Nradd(-Nrad:Nrad),Vrad(-Nrad:Nrad),
     +             Vzz(-Nrad:Nrad),Vfi2(-Nrad:Nrad),
     +             Vrad2(-Nrad:Nrad),Vzz2(-Nrad:Nrad),Radd(-Nrad:Nrad)
      COMMON /DISKp/ Vfip(-Nrad:Nrad),Nradp(-Nrad:Nrad),
     +             Vradp(-Nrad:Nrad),Vzzp(-Nrad:Nrad),Vfi2p(-Nrad:Nrad),
     +             Vrad2p(-Nrad:Nrad),Vzz2p(-Nrad:Nrad),Radp(-Nrad:Nrad)
      COMMON /BAR/ Vfib(-Nrad:Nrad,Nbp),Nradb(-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),
     +      RadBar(0:Nbarbins),eccentr(0:Nbarbins),nBarin(0:Nbarbins),
     +      amSin(0:Nbarbins),amCos(0:Nbarbins),barangle(0:Nbarbins)
      COMMON /BULGE/ Nradg(-Nrad:Nrad),Vfig(-Nrad:Nrad),
     +        Vradg(-Nrad:Nrad),Vfig2(-Nrad:Nrad),Nradu(-Nrad:Nrad),
     +             Vradg2(-Nrad:Nrad),Radg(-Nrad:Nrad)
      COMMON /HALO/ Nradh(-Nrad:Nrad),Vfih(-Nrad:Nrad),
     +             Vradh(-Nrad:Nrad),Vfih2(-Nrad:Nrad),
     +             Vradh2(-Nrad:Nrad),Radh(-Nrad:Nrad)
      COMMON/BarLength/bar_angle,bar_length,bar_amp ! params of Bar
      COMMON/DensA/Label(Nmlink:Nblink,Nmlink:Nblink,Nmlink:Nblink)
      COMMON/DensB/Dens(Np),
     +           Lst(Np),Cell,dStep,RadDens(Nsearch)
      COMMON /CATx/ Xp(Np)
      COMMON /CATy/ Yp(Np)
      COMMON /CATz/ Zp(Np)
      COMMON /CAT1/ VXp(Np)
      COMMON /CAT2/ VYp(Np)
      COMMON /CAT3/ VZp(Np)
      COMMON /GENERAL/time
      double precision time