def read_patch(self, file, verbose=0):
''' Read particles in a patch '''
name = file + '.peb'
f = FortranFile(name, 'rb')
fmt = f.readInts()
if verbose > 2:
print('ioformat:', fmt[0])
dim = f.readInts()
dim = dim[0:2]
a = f.readReals()
a = a.reshape(dim[1], dim[0]).transpose()
idv = f.readInts()
idx = idv[0:dim[0]]
f.close()
dict = {}
for i in range(size(idx)):
id = idx[i]
dict[id] = {
'p': a[i, 0:3],
'v': a[i, 0:3],
't': a[i, 6],
's': a[i, 7],
'w': a[i, 8]
}
return idx, dict
