pro ReadParts




;base='/mybase/'
;ext='myext.halos'

fname=strcompress(base+ext,/remove_all)



nparts=0L


print ,'Reading : ',fname
close, 1
openr, 1, fname, /f77_unformatted





readu,1 , nparts
print,'nparts : ',nparts

PartPos=fltarr(3,nparts)
PartVel=intarr(3,nparts)


readu, 1, PartPos

; e.g. SGZ(10th part)=PartPos(2,10)
; NOTE : the particles are NOT ordered !


readu, 1, PartVel


; e.g. vx(2nd part)=PartVel(0,2)
; NOTE : the particles are NOT ordered !


close, 1

print, 'Read ...'

print,'min,max PartPos : ',min(PartPos(*,*)),max(PartPos(*,*))
print,'min,max PartPos : ',min(PartVel(*,*)),max(PartVel(*,*))


end