Commit fb1cde23 authored by Médéric Boquien's avatar Médéric Boquien

Eliminate the unphysical wavy structure of the nebular continuum in the...

Eliminate the unphysical wavy structure of the nebular continuum in the mid-infrared by pre-interpolating it on a refined wavelength grid in log-log space.
parent e5a80f58
Pipeline #2026 skipped with stage
......@@ -39,6 +39,7 @@
- Make sure that upper limits on physical properties are correctly taken into account. (Médéric Boquien)
- Improve the sanitation of input data so that upper limits of extensive properties are not eliminated from the quantities to be fitted when upper limits are activated. (Médéric Boquien)
- Ensure that the radio module still works with recent versions of numpy. (Médéric Boquien & Laure Ciesla, reported by Wenjia Zhou)
- The unphysical wavy structure of the nebular continuum in the mid-infrared has been eliminated. (Médéric Boquien)
### Optimised
- Slight speedup of the computation of the likelihood from the χ² using a multiplication rather than a division. (Médéric Boquien)
- Speedup of the computation of the χ² by ~10% taking the opposite of a scalar rather than of an array. (Médéric Boquien)
......
......@@ -760,7 +760,7 @@ def build_nebular(base):
# Convert wavelength from Å to nm
wave_lines *= 0.1
wave_cont = cont[:1600, 0] * 0.1
wave_cont = cont[:3729, 0] * 0.1
# Get the list of metallicities
metallicities = np.unique(lines[:, 1])
......@@ -790,7 +790,7 @@ def build_nebular(base):
# Import continuum
for idx, metallicity in enumerate(metallicities):
spectra = cont[1600 * idx: 1600 * (idx+1), :]
spectra = cont[3729 * idx: 3729 * (idx+1), :]
for logU, spectrum in zip(np.around(np.arange(-4., -.9, .1), 1),
spectra.T):
models_cont.append(NebularContinuum(metallicity, logU, wave_cont,
......
This diff is collapsed.
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