Commit ec286843 authored by Guang's avatar Guang Committed by Guang
Browse files

bug fixed for AGN disk SED shape

parent 5c4765e7
......@@ -747,11 +747,11 @@ def build_skirtor2016(base):
disk_new_unnormed = wl_uncut.copy()
# 2nd part: 50-125 nm
wave_idxs = (wl_uncut>=50) & (wl_uncut<125)
# Here, the factor 50**0.2 is to re-nomalize for continuity
norm_fac = 50**0.2
disk_new_unnormed[wave_idxs] = disk_new_unnormed[wave_idxs]**0.8 * norm_fac
# Here, norm_fac is to re-nomalize for continuity
norm_fac = 50**1.2
disk_new_unnormed[wave_idxs] = disk_new_unnormed[wave_idxs]**-0.2 * norm_fac
# 3rd part: 125-1e4 nm
norm_fac *= 125**2.3
norm_fac *= 125**1.3
wave_idxs = (wl_uncut>=125) & (wl_uncut<1e4)
disk_new_unnormed[wave_idxs] = disk_new_unnormed[wave_idxs]**-1.5 * norm_fac
# 4th part: >1e4 nm
......
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