How to use the pynbody.units.ratio function in pynbody
To help you get started, we’ve selected a few pynbody examples, based on popular ways it is used in public projects.
Secure your code as it's written. Use Snyk Code to scan source code in minutes - no build needed - and fix issues immediately.
pynbody / pynbody / pynbody / halo / ahf.py View on Github 
if not os.path.exists(groupfinder):
raise RuntimeError("Path to AHF (%s) is invalid" % groupfinder)
if (os.path.basename(groupfinder) == 'AHFstep'):
isAHFstep = True
else:
isAHFstep = False
# build units file
if isAHFstep:
f = open('tipsy.info', 'w')
f.write(str(sim.properties['omegaM0']) + "\n")
f.write(str(sim.properties['omegaL0']) + "\n")
f.write(str(sim['pos'].units.ratio(
units.kpc, a=1) / 1000.0 * sim.properties['h']) + "\n")
f.write(
str(sim['vel'].units.ratio(units.km / units.s, a=1)) + "\n")
f.write(str(sim['mass'].units.ratio(units.Msol)) + "\n")
f.close()
# make input file
f = open('AHF.in', 'w')
f.write(sim._filename + " " + str(typecode) + " 1\n")
f.write(sim._filename + "\n256\n5\n5\n0\n0\n0\n0\n")
f.close()
else:
# make input file
f = open('AHF.in', 'w')
lgmax = np.min([int(2 ** np.floor(np.log2(
1.0 / np.min(sim['eps'])))), 131072])
# hardcoded maximum 131072 might not be necessary
print >>f, config_parser.get('AHFCatalogue', 'Config', vars={raise RuntimeError("Path to AHF (%s) is invalid" % groupfinder)
if (os.path.basename(groupfinder) == 'AHFstep'):
isAHFstep = True
else:
isAHFstep = False
# build units file
if isAHFstep:
f = open('tipsy.info', 'w')
f.write(str(sim.properties['omegaM0']) + "\n")
f.write(str(sim.properties['omegaL0']) + "\n")
f.write(str(sim['pos'].units.ratio(
units.kpc, a=1) / 1000.0 * sim.properties['h']) + "\n")
f.write(
str(sim['vel'].units.ratio(units.km / units.s, a=1)) + "\n")
f.write(str(sim['mass'].units.ratio(units.Msol)) + "\n")
f.close()
# make input file
f = open('AHF.in', 'w')
f.write(sim._filename + " " + str(typecode) + " 1\n")
f.write(sim._filename + "\n256\n5\n5\n0\n0\n0\n0\n")
f.close()
else:
# make input file
f = open('AHF.in', 'w')
lgmax = np.min([int(2 ** np.floor(np.log2(
1.0 / np.min(sim['eps'])))), 131072])
# hardcoded maximum 131072 might not be necessary
print >>f, config_parser.get('AHFCatalogue', 'Config', vars={
'filename': sim._filename,lgmax = np.min([int(2 ** np.floor(np.log2(
1.0 / np.min(sim['eps'])))), 131072])
# hardcoded maximum 131072 might not be necessary
print >>f, config_parser.get('AHFCatalogue', 'Config', vars={
'filename': sim._filename,
'typecode': typecode,
'gridmax': lgmax
})
print >>f, config_parser.get('AHFCatalogue', 'ConfigTipsy', vars={
'omega0': sim.properties['omegaM0'],
'lambda0': sim.properties['omegaL0'],
'boxsize': sim['pos'].units.ratio('Mpc a h^-1', **sim.conversion_context()),
'vunit': sim['vel'].units.ratio('km s^-1 a', **sim.conversion_context()),
'munit': sim['mass'].units.ratio('Msol h^-1', **sim.conversion_context()),
'eunit': 0.03 # surely this can't be right?
})
f.close()
if (not os.path.exists(sim._filename)):
os.system("gunzip " + sim._filename + ".gz")
# determine parallel possibilities
if os.path.exists(groupfinder):
# run it
os.system(groupfinder + " AHF.in")
returngroupfinder = ahf
break
if not os.path.exists(groupfinder):
raise RuntimeError("Path to AHF (%s) is invalid" % groupfinder)
if (os.path.basename(groupfinder) == 'AHFstep'):
isAHFstep = True
else:
isAHFstep = False
# build units file
if isAHFstep:
f = open('tipsy.info', 'w')
f.write(str(sim.properties['omegaM0']) + "\n")
f.write(str(sim.properties['omegaL0']) + "\n")
f.write(str(sim['pos'].units.ratio(
units.kpc, a=1) / 1000.0 * sim.properties['h']) + "\n")
f.write(
str(sim['vel'].units.ratio(units.km / units.s, a=1)) + "\n")
f.write(str(sim['mass'].units.ratio(units.Msol)) + "\n")
f.close()
# make input file
f = open('AHF.in', 'w')
f.write(sim._filename + " " + str(typecode) + " 1\n")
f.write(sim._filename + "\n256\n5\n5\n0\n0\n0\n0\n")
f.close()
else:
# make input file
f = open('AHF.in', 'w')
lgmax = np.min([int(2 ** np.floor(np.log2(
1.0 / np.min(sim['eps'])))), 131072])groupfinder = ahf
break
if not os.path.exists(groupfinder):
raise RuntimeError("Path to AHF (%s) is invalid" % groupfinder)
if (os.path.basename(groupfinder) == 'AHFstep'):
isAHFstep = True
else:
isAHFstep = False
# build units file
if isAHFstep:
f = open('tipsy.info', 'w')
f.write(str(sim.properties['omegaM0']) + "\n")
f.write(str(sim.properties['omegaL0']) + "\n")
f.write(str(sim['pos'].units.ratio(
units.kpc, a=1) / 1000.0 * sim.properties['h']) + "\n")
f.write(
str(sim['vel'].units.ratio(units.km / units.s, a=1)) + "\n")
f.write(str(sim['mass'].units.ratio(units.Msol)) + "\n")
f.close()
# make input file
f = open('AHF.in', 'w')
f.write(sim._filename + " " + str(typecode) + " 1\n")
f.write(sim._filename + "\n256\n5\n5\n0\n0\n0\n0\n")
f.close()
else:
# make input file
f = open('AHF.in', 'w')
lgmax = np.min([int(2 ** np.floor(np.log2(
1.0 / np.min(sim['eps'])))), 131072])raise RuntimeError("Path to AHF (%s) is invalid" % groupfinder)
if (os.path.basename(groupfinder) == 'AHFstep'):
isAHFstep = True
else:
isAHFstep = False
# build units file
if isAHFstep:
f = open('tipsy.info', 'w')
f.write(str(sim.properties['omegaM0']) + "\n")
f.write(str(sim.properties['omegaL0']) + "\n")
f.write(str(sim['pos'].units.ratio(
units.kpc, a=1) / 1000.0 * sim.properties['h']) + "\n")
f.write(
str(sim['vel'].units.ratio(units.km / units.s, a=1)) + "\n")
f.write(str(sim['mass'].units.ratio(units.Msol)) + "\n")
f.close()
# make input file
f = open('AHF.in', 'w')
f.write(sim._filename + " " + str(typecode) + " 1\n")
f.write(sim._filename + "\n256\n5\n5\n0\n0\n0\n0\n")
f.close()
else:
# make input file
f = open('AHF.in', 'w')
lgmax = np.min([int(2 ** np.floor(np.log2(
1.0 / np.min(sim['eps'])))), 131072])
# hardcoded maximum 131072 might not be necessary
print >>f, config_parser.get('AHFCatalogue', 'Config', vars={
'filename': sim._filename,lgmax = np.min([int(2 ** np.floor(np.log2(
1.0 / np.min(sim['eps'])))), 131072])
# hardcoded maximum 131072 might not be necessary
print >>f, config_parser.get('AHFCatalogue', 'Config', vars={
'filename': sim._filename,
'typecode': typecode,
'gridmax': lgmax
})
print >>f, config_parser.get('AHFCatalogue', 'ConfigTipsy', vars={
'omega0': sim.properties['omegaM0'],
'lambda0': sim.properties['omegaL0'],
'boxsize': sim['pos'].units.ratio('Mpc a h^-1', **sim.conversion_context()),
'vunit': sim['vel'].units.ratio('km s^-1 a', **sim.conversion_context()),
'munit': sim['mass'].units.ratio('Msol h^-1', **sim.conversion_context()),
'eunit': 0.03 # surely this can't be right?
})
f.close()
if (not os.path.exists(sim._filename)):
os.system("gunzip " + sim._filename + ".gz")
# determine parallel possibilities
if os.path.exists(groupfinder):
# run it
os.system(groupfinder + " AHF.in")
returnpynbody
Light-weight astronomical N-body/SPH analysis for python
Package Health Score
78 / 100Popular pynbody functions
- pynbody.analysis
- pynbody.analysis.angmom
- pynbody.analysis.angmom.faceon
- pynbody.analysis.profile.Profile
- pynbody.analysis.profile.Profile.profile_property
- pynbody.array
- pynbody.array.SimArray
- pynbody.config_parser
- pynbody.config_parser.get
- pynbody.family
- pynbody.load
- pynbody.snapshot.gadgethdf.GadgetHDFSnap.derived_quantity
- pynbody.snapshot.gadgethdf.SubFindHDFSnap.derived_quantity
- pynbody.units
- pynbody.units.ratio
- pynbody.units.Unit
- pynbody.units.UnitsException
- pynbody.util
- pynbody.util.open_
- pynbody.util.skip_fortran