def yso_class(nu_star, fnu_star, nu, fnu, nu_fnu, t_star, r_star):
# Protostellar bolometric luminosity
L_star = fs.Lbol_calc(nu_star, fnu_star)
# Model bolometric luminosity
L_bol = fs.Lbol_calc(nu, fnu)
# Ratio of model sub-mm to bolometric luminosity
L_ratio = fs.L_ratio_calc(nu, fnu)
# Compute bolometric temperature for model SED
T_bol = fs.Tbol_calc(nu, fnu, nu_fnu)
# Print values to terminal
print("\nProtostellar luminosity (SED) is: {0} L_sol\n".format(L_star))
print("Bolometric luminosity is: {0} L_sol\n".format(L_bol))
print("Ratio of Bolometric and sub-mm "+ \
"luminosity is {0} %\n".format(L_ratio) )
print("Bolometric temperature is: {0} K\n".format(T_bol))
return
nu_q, fnu_q, nu_fnu_q = fs.SED_read(f_sed_q)
#
flam_int_o = interp1d(wavs_o,flam_o, kind="linear")
flam_int_q = interp1d(wavs_q,flam_q, kind="linear")
for w in range(len(wavs_compare)):
flam_ratio[cut][s][c][0][i][w] = \
flam_int_o(wavs_compare[w]) / flam_int_q(wavs_compare[w])
# Model bolometric luminosity/temperature - quiescent
if (s == 0):
#
wavs_q *= cs.cm_per_micron
#
L_tmp = fs.Lbol_calc(nu_q, fnu_q)
L_bol[cut][0][c][0][i].append( L_tmp )
#
T_tmp = fs.Tbol_calc(nu_q, fnu_q, nu_fnu_q)
T_bol[cut][0][c][0][i].append( T_tmp )
# Model bolometric luminosity/temperature - outbursting
wavs_o *= cs.cm_per_micron
#
L_tmp = fs.Lbol_calc(nu_o, fnu_o)
L_bol[cut][1 + s][c][0][i].append( L_tmp )
#
T_tmp = fs.Tbol_calc(nu_o, fnu_o, nu_fnu_o)
T_bol[cut][1 + s][c][0][i].append( T_tmp )
#
# Read SED data
fsed = dat_dir + theta[i] + "/spectrum6L_ISRF_" + str(
best_inclin) + "i.out"
wav[i], flam[i], lam_flam[i], nu, fnu, nu_fnu = fs.SED_read(fsed)
for j in range(len(lam_flam[i])):
lam_flam[i][j] /= dist**(2.0)
if bolometrics:
# Model bolometric luminosity
L_bol = fs.Lbol_calc(nu, fnu)
# Model bolometric temperature
T_bol = fs.Tbol_calc(nu, fnu, nu_fnu)
# Model sub-mm to bolometric luminosity ratio
L_ratio = fs.L_ratio_calc(nu, fnu)
# Print results to terminal
print("\nFor theta = {0} deg., {1} model: \n".format(best_cavity, \
leg_names[i]))
print("Bolometric luminosity is: {0} L_sol\n".format(L_bol))