[wzrzainh]

In text output, I don't find particles velocity or position. Did I not turn on particles tracking? In addition, I used this subroutine to Counting Particle Concentration, it dosent work.




!------ Count particles in grid cells for mesh block nbl
!
where sclr(ijk,1) = number of particles in cell ijk
!---- initialize scalar array #1 to zero
if(nsc.gt.0) then
do ijk=ijklim(nbl,1),ijklim(nbl,2)
sclr(ijk,1)=zero
enddo
!---- loop over particles in mesh block nbl to make the count
do k1=1,itnps
if(ipmbl(k1).ne.nbl) cycle ! skip particles that do not belong
!
ijk_str=ijkprt(k1)
!
!
!
to current mesh block
! extract structured index of the cell where
particle k1 is located
ijk=ijk_str2unstr(ijk_str) ! convert structured cell index to fluid
sub-domain index
sclr(ijk,1)=sclr(ijk,1)+one ! increment the particle count
in cell ijk
enddo
endif