Commit 11edf982 authored by Guang's avatar Guang Committed by Guang
Browse files

cut AGN module to lam>5nm

parent 4e345650
......@@ -124,6 +124,15 @@ class Fritz2006(SedModule):
self.l_agn_agn = np.trapz(self.fritz2006.lumin_agn,
x=self.fritz2006.wave)
# Apply wavelength cut
lam_cut = 5.
lam_use_idxs = self.fritz2006.wave>=lam_cut
self.fritz2006.wave = self.fritz2006.wave[lam_use_idxs]
self.fritz2006.lumin_agn = self.fritz2006.lumin_agn[lam_use_idxs]
self.fritz2006.lumin_scatt = self.fritz2006.lumin_scatt[lam_use_idxs]
self.fritz2006.lumin_therm = self.fritz2006.lumin_therm[lam_use_idxs]
self.fritz2006.lumin_intrin_agn = self.fritz2006.lumin_intrin_agn[lam_use_idxs]
# Apply polar-dust obscuration
# We define various constants necessary to compute the model
self.c = cst.c * 1e9
......
......@@ -54,7 +54,7 @@ class SKIRTOR2016(SedModule):
('oa', (
'cigale_list(options=10 & 20 & 30 & 40 & 50 & 60 & 70 & 80)',
"Angle measured between the equatorial plan and edge of the torus. "
"Half-opening angle of the dust-free cone is 90-oa"
"Half-opening angle of the dust-free cone is 90-oa. "
"Possible values are: 10, 20, 30, 40, 50, 60, 70, and 80",
40
)),
......@@ -120,10 +120,18 @@ class SKIRTOR2016(SedModule):
self.oa, self.R, self.Mcl,
self.i)
# Apply wavelength cut
lam_cut = 5.
lam_use_idxs = self.SKIRTOR2016.wave>=lam_cut
self.SKIRTOR2016.wave = self.SKIRTOR2016.wave[lam_use_idxs]
self.SKIRTOR2016.disk = self.SKIRTOR2016.disk[lam_use_idxs]
self.SKIRTOR2016.dust = self.SKIRTOR2016.dust[lam_use_idxs]
self.SKIRTOR2016.intrin_iso_disk = self.SKIRTOR2016.intrin_iso_disk[lam_use_idxs]
# Modify the model so that the central source is isotropic
# Calculate the coverting factors disk emission
# (only for type-1 case, since disk emission in type-2 is zero)
# (since dust emission is normalized to unity; it remains same)
# (dust emission is normalized to unity, so it remains same)
if self.i<=(90-self.oa):
self.SKIRTOR2016.disk*= (4*np.sin(self.oa*np.pi/180)**2 + 3*np.sin(self.oa*np.pi/180)) / \
(12*np.cos(self.i*np.pi/180)**2 + 6*np.cos(self.i *np.pi/180))
......
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