npart_read = tmp(1)
npart_tile(bi,bj) = INT(npart_read)
do ip=1,INT(npart_read)
if (.not. globalFile) _GLOBAL_SUM_R8(npart_read,myThid)
write(errormessageunit,*) 'sum npart_read: ',npart_read
_RL npart_read, npart_dist
irecord=npart_read+ip+1
_RL npart_read,npart_times
npart_read = 0.
tmp(1) = DBLE(npart_tile(bi,bj))+npart_read