def TBOL(nu,S,Snu):
#calculate integrals for I0 and I1
SUM_Snu = -trapizum.trap(nu,Snu)
SUM_S = -trapizum.trap(nu,S)
V_m = SUM_Snu/SUM_S
TBOL = V_m*bol
return TBOL
def LBOL(d,nu,nu_smm,S,S_smm):
L = 3.839E26
#calculate integrals for I0 and I1
SUM_S = -trapizum.trap(nu,S)
#SubMM integral
SMM = -trapizum.trap(nu_smm,S_smm)
#luminosity calc.
L_bol = (4*np.pi*(d**(2.))*SUM_S)/L
L_SMM = (4*np.pi*(d**(2.))*SMM)/L
L_R = (L_SMM/L_bol)*100
Lumin = [L_bol,L_SMM,L_R]
return Lumin