guilhem_lavaux/vide_public
0
0
Fork 0
mirror of https://bitbucket.org/cosmicvoids/vide_public.git synced 2025-07-04 15:21:11 +00:00
Code Issues Projects Releases Packages Wiki Activity Actions

added utility to read off particle velocities

Browse source
This commit is contained in:
P.M. Sutter 2014-01-27 20:26:27 -06:00
parent d3554b0732
commit b466730c42
1 changed files with 24 additions and 0 deletions
Download patch file Download diff file Expand all files Collapse all files

24
python_tools/void_python_tools/partUtil/partUtil.py
View file

@ -105,6 +105,30 @@ def loadPart(workDir, sampleDir, sample):
return partData, boxLen, volNorm, isObservationData
# -----------------------------------------------------------------------------
def loadPartVel(workDir, sampleDir, sample):
#print " Loading particle velocities..."
sys.stdout.flush()
infoFile = workDir+"/"+sampleDir+"/zobov_slice_"+sample.fullName+".par"
File = NetCDFFile(infoFile, 'r')
isObservation = getattr(File, 'is_observation')
if isObservation:
print "No velocities for observations!"
return -1
vx = File.variables['vel_x'][0:]
vy = File.variables['vel_y'][0:]
vx = File.variables['vel_z'][0:]
File.close()
partVel = np.column_stack((vx,vy,vz))
return partVel
# -----------------------------------------------------------------------------
def shiftPart(inPart, newCenter, periodicLine, boxLen):