Compare commits

..

2 Commits

2 changed files with 9 additions and 3 deletions

View File

@ -62,6 +62,12 @@ class Particles:
col[key]-=1 col[key]-=1
return cls.from_array(pp,col,time,period,select_col=select_col) return cls.from_array(pp,col,time,period,select_col=select_col)
@classmethod
def from_ucf(cls,file,period,select_col=None):
from .ucf import read_particles
(pp,col,time) = read_particles(file,step=1,verbosity=False,debug=False)
return cls.from_array(pp,col,time,period,select_col=select_col)
@classmethod @classmethod
def from_position(cls,x,y,z,time,period): def from_position(cls,x,y,z,time,period):
assert x.ndim==1 and y.ndim==1 and z.ndim==1,\ assert x.ndim==1 and y.ndim==1 and z.ndim==1,\

6
ucf.py
View File

@ -864,9 +864,9 @@ def colmap_from_flags(irank,ihybrid,idem,iscal):
col['tyc'] = ioffset; ioffset+=1 col['tyc'] = ioffset; ioffset+=1
col['tzc'] = ioffset; ioffset+=1 col['tzc'] = ioffset; ioffset+=1
if iscal>0: if iscal>0:
for ii in range(0,iscal): for ii in range(iscal):
col['s'+str(ii)] = ioffset; ioffset+=1 col['s'+str(ii+1)] = ioffset; ioffset+=1
col['q'+str(ii)] = ioffset; ioffset+=1 col['q'+str(ii+1)] = ioffset; ioffset+=1
return col return col
def grid_chunk(chunk,gridg): def grid_chunk(chunk,gridg):