c     Particle properties:
c     hei=n_HeI/n_H, heiii=n_HeIII/n_H, H=hi+hii
c     if hi<0, then     n_HII=-hi * n_H (in old TT version)
c               else    n_Hi =+hi * n_H
      real*4 rm,r,v,hi,hei,heiii,dn,h,e,zmetal
      integer*4 itype
      common /part1/ rm(Nmax+4)
      common /part2/ r(3,Nmax+4)
      common /part3/ v(3,Nmax+4)
		common /abun/  hi(nmax),hei(nmax),heiii(nmax)
      common /dens/  dn(Nmax+4),h(Nmax+4),e(Nmax+4),zmetal(Nmax+4)
      common /types/ itype(Nmax)