import numpy as np
[docs]def initSKIDAttribs(obj):
obj.locs.SKID_DIR = '%s/%s/SKID' % (obj.locs.snapdir,
obj.rto.GROUPDIRPREFIX)
obj.locs.SKID_FILE = '%s/gals_%04d_z%0.3f' % (obj.locs.SKID_DIR,
obj.snapnum,
obj.redshift)
obj.locs.SKID_STATFILE = '%s.stat' % obj.locs.SKID_FILE
obj.locs.SKID_GRPFILE = '%s.grp' % obj.locs.SKID_FILE
"""
##!!! SKID UNIT CONVERSIONS !!!##
CM_PER_MPC = obj.constants.cm_per_mpc
UnitLength_in_cm = obj.constants.UnitLength_in_cm
UnitMass_in_g = obj.constants.UnitMass_in_g
UnitVelocity_in_cm_per_s = obj.constants.UnitVelocity_in_cm_per_s
UnitLength_in_cm /= obj.h
UnitMass_in_g /= obj.h
L = obj.boxsize_h * 1.0e-3
UnitTime_in_s = UnitLength_in_cm / UnitVelocity_in_cm_per_s
UnitDensity_in_cgs = UnitMass_in_g / UnitLength_in_cm**3
UnitEnergy_in_cgs = UnitMass_in_g * UnitLength_in_cm**2 / UnitTime_in_s**2
unit_Time = np.sqrt(8.*np.pi/3.) * CM_PER_MPC / (100.*obj.h*UnitVelocity_in_cm_per_s)
unit_Density = 1.8791e-29*obj.h**2
unit_Length = L*CM_PER_MPC
unit_Mass = unit_Density * unit_Length**3
unit_Velocity = unit_Length / unit_Time
#print("COSMO PARAMS: L=%g Mpc, h=%g, Omega=%g" % (L,obj.h,obj.O0))
#print("UNITS: T=%g rho=%g L=%g M=%g v=%g" %
# (unit_Time,unit_Density,unit_Length,unit_Mass,unit_Velocity))
## convert back to GADGET units,
## convert gadget_mass -> Msun,
## then divide out that damned little h
obj.conversionFactors.TIPSY_M_convert = unit_Mass / UnitMass_in_g
obj.conversionFactors.TIPSY_M_convert *= 1.0e10 / obj.h
obj.conversionFactors.TIPSY_L_convert = unit_Length / UnitLength_in_cm
obj.conversionFactors.TIPSY_L_convert /= obj.h
obj.conversionFactors.TIPSY_V_convert = unit_Velocity / UnitVelocity_in_cm_per_s
obj.conversionFactors.TIPSY_RHO_convert = (UnitLength_in_cm**3 * unit_Density) / UnitMass_in_g
obj.conversionFactors.TIPSY_RHO_convert*= 1.0e10 * obj.h**2
##!!!##
"""