Commit 1f2a11c4 authored by Guang's avatar Guang Committed by Guang
Browse files

bug fixed for polar-dust extinction

parent 27eb17bc
......@@ -142,10 +142,9 @@ class Fritz2006(SedModule):
self.c = cst.c * 1e9
lambda_0 = 200e3
# Calculate the extinction (SMC)
k_lam = 1.39*(self.fritz2006.wave/1e3)**-1.2-0.38
# The analytical formula is from Bongiorno+2012
k_lam = 1.39*(self.fritz2006.wave/1e3)**-1.2
A_lam = k_lam * self.EBV
# Set negative value to zero
A_lam[A_lam<0] = 0
# The extinction factor, flux_1/flux_0
ext_fac = 10**(A_lam/-2.5)
# Calculate the new AGN SED shape after extinction
......
......@@ -125,10 +125,9 @@ class SKIRTOR2016(SedModule):
self.c = cst.c * 1e9
lambda_0 = 200e3
# Calculate the extinction (SMC)
k_lam = 1.39*(self.SKIRTOR2016.wave/1e3)**-1.2-0.38
# The analytical formula is from Bongiorno+2012
k_lam = 1.39*(self.SKIRTOR2016.wave/1e3)**-1.2
A_lam = k_lam * self.EBV
# Set negative value to zero
A_lam[A_lam<0] = 0
# The extinction factor, flux_1/flux_0
ext_fac = 10**(A_lam/-2.5)
# Calculate the new AGN SED shape after extinction
......
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment