def calculate_Vr2_cm_per_s_2(input_data) :
d = input_data
Vr_cm_per_s_2 = {}
for hid in d['halo_ids'] :
Vr_cm_per_s_2[hid] = make_profile(x=d['Rscaled'][hid],
y=(d['V_r'][hid] * km2cm)**2)
return Vr_cm_per_s_2
def calculate_Vr2_cm_per_s_2(input_data):
d = input_data
Vr_cm_per_s_2 = {}
for hid in d['halo_ids']:
Vr_cm_per_s_2[hid] = make_profile(x=d['Rscaled'][hid],
y=(d['V_r'][hid] * km2cm)**2)
return Vr_cm_per_s_2
def calculate_Tmw_cm_per_s_2(input_data) :
d = input_data
Tmw_cm_per_s_2 = {}
for hid in d['halo_ids'] :
Tmw_erg = d['T_mw'][hid] * kb
Tmw_cm_per_s_2[hid] = make_profile(x=d['Rscaled'][hid],
y=Tmw_erg/mu/mp)
return Tmw_cm_per_s_2
def calculate_Tmw_cm_per_s_2(input_data):
d = input_data
Tmw_cm_per_s_2 = {}
for hid in d['halo_ids']:
Tmw_erg = d['T_mw'][hid] * kb
Tmw_cm_per_s_2[hid] = make_profile(x=d['Rscaled'][hid],
y=Tmw_erg / mu / mp)
return Tmw_cm_per_s_2
def calculate_Tnt_cm_per_s_2(input_data) :
d = input_data
Tnt_cm_per_s_2 = {}
for hid in d['halo_ids'] :
Tnt_erg = keV2erg*calculate_T_nonthermal(sigr=d['sigr'][hid],vr=d['vr'][hid],
sigt=d['sigt'][hid],vt=d['vt'][hid])
Tnt_cm_per_s_2[hid] = make_profile(x=d['Rscaled'][hid],
y=Tnt_erg/mu/mp)
return Tnt_cm_per_s_2
def calculate_Tnt_cm_per_s_2(input_data):
d = input_data
Tnt_cm_per_s_2 = {}
for hid in d['halo_ids']:
Tnt_erg = keV2erg * calculate_T_nonthermal(sigr=d['sigr'][hid],
vr=d['vr'][hid],
sigt=d['sigt'][hid],
vt=d['vt'][hid])
Tnt_cm_per_s_2[hid] = make_profile(x=d['Rscaled'][hid],
y=Tnt_erg / mu / mp)
return Tnt_cm_per_s_2
def calculate_Ttot_cm_per_s_2(input_data) :
d = input_data
Ttot_cm_per_s_2 = {}
for hid in d['halo_ids'] :
Ttot_erg = keV2erg*calculate_Ttot(sigr=d['sigr'][hid],vr=d['vr'][hid],
sigt=d['sigt'][hid],vt=d['vt'][hid],
Tmw=d['T_mw'][hid])
Ttot_cm_per_s_2[hid] = make_profile(x=d['Rscaled'][hid],
y=Ttot_erg/mu/mp)
return Ttot_cm_per_s_2
def calculate_Ttot_cm_per_s_2(input_data):
d = input_data
Ttot_cm_per_s_2 = {}
for hid in d['halo_ids']:
Ttot_erg = keV2erg * calculate_Ttot(sigr=d['sigr'][hid],
vr=d['vr'][hid],
sigt=d['sigt'][hid],
vt=d['vt'][hid],
Tmw=d['T_mw'][hid])
Ttot_cm_per_s_2[hid] = make_profile(x=d['Rscaled'][hid],
y=Ttot_erg / mu / mp)
return Ttot_cm_per_s_2