npart_tile(bi,bj) = npart_read
npart_read = MIN( npart_read, max_npart_tile )
DO ip=1,npart_read
INTEGER npart_read
npart_read = 0
npart_read = NINT(tmp(1))
irecord = npart_read+ip+1
_RL npart_read,npart_times
npart_read = tmp(1)
npart_read = 0.
tmp(1) = DBLE(npart_tile(bi,bj))+npart_read
_RL npart_read, npart_times