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

fix to halo catalog handling

Browse source
This commit is contained in:
P.M. Sutter 2012-12-04 21:21:30 -06:00
parent 3cb2a990da
commit f155564d6d
1 changed files with 2 additions and 1 deletions
Download patch file Download diff file Expand all files Collapse all files

3
python_tools/pipeline_source/prepareCatalogs.in.py
View file

@ -356,7 +356,7 @@ if (args.halos or args.all) and haloFileBase != "":
numPart += 1 numPart += 1
inFile.close() inFile.close()
sampleName = prefix+"halos_min"+str(minHaloMass)+"_z"+redshift sampleName = prefix+"halos_min"+str(minHaloMass)+"_z"+fileNums[iRedshift]
outFile = open(catalogDir+"/"+sampleName+".dat", 'w') outFile = open(catalogDir+"/"+sampleName+".dat", 'w')
outFile.write("%f\n" %(lbox)) outFile.write("%f\n" %(lbox))
@ -367,6 +367,7 @@ if (args.halos or args.all) and haloFileBase != "":
inFile = open(dataFile, 'r') inFile = open(dataFile, 'r')
for (iHalo,line) in enumerate(inFile): for (iHalo,line) in enumerate(inFile):
if iHalo < haloFileNumComLines: continue
line = line.split(haloFileColSep) line = line.split(haloFileColSep)
if minHaloMass == "none" or float(line[haloFileMCol]) > minHaloMass: if minHaloMass == "none" or float(line[haloFileMCol]) > minHaloMass:
x = float(line[haloFileXCol]) x = float(line[haloFileXCol])