OxyYield = (np.sum(Stellar_oxygen*S['m']) + np.sum(Gas_oxygen*G['m'])) / np.sum(S['m'])
EmpiricalYield = Total_metal_mass / np.sum(S['m'])
a = G['header'][2]
z = 1.0/a - 1.0
redshift = G['header'][3]
boxsize = G['header'][9]
omega_matter = G['header'][10]
omega_L = G['header'][11]
h = G['header'][12]
redshiftstring = "{0:.3f}".format(redshift)
Hubble = hubble_param(a, omega_matter, h)
PhysTemp, PhysRho = SF.convertTemp(G['u'], G['ne'], G['rho'], h)
halostats = SF.find_halo_now(halo_to_do, a, therod=use_fixed_halos)
print 'halo stats', halostats
haloN = halostats[1]
#read halo catalog, assign halo properties
Rvir = halostats[11]
Vsig = halostats[10]
haloX = halostats[2]
haloY = halostats[3]
haloZ = halostats[4]
haloVX =halostats[5]
haloVY = halostats[6]
haloVZ = halostats[7]
M = halostats[8]
Mstar = halostats[13]
Syoung = stellar_age_cut Rstars = where_do_we_go * h / a RstarsPhys = where_do_we_go workcount = 0 Hubble = hubble_param(a, omega_matter, h) ########### #Work function ########### C = SF.read_halo_catalog(Nsnapstring, redshiftstring) while (workcount < len(halo_to_do)): halostats = SF.find_halo_now(halo_to_do[workcount], a, therod=use_fixed_halos) haloN = halostats[1] Rvir = halostats[11] haloX = halostats[2] haloY = halostats[3] haloZ = halostats[4] Mstar = halostats[13] Vsig = halostats[10] Mgas = halostats[12] M = halostats[4] Rvir, Vsig, M = SF.check_Rvir_growth(halo_to_do[workcount], a, Rvir, Vsig, M, therod=use_fixed_halos) if (adaptive_Rstar): Rstars = Inner_SF_thresh*Rvir RstarsPhys = Rstars * a / h
