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

Call np.trapz() once instead of twice.

parent f6449ef8
...@@ -81,8 +81,8 @@ class BC03(SedModule): ...@@ -81,8 +81,8 @@ class BC03(SedModule):
# compute the energy absorbed by the dust before ionising gas. # compute the energy absorbed by the dust before ionising gas.
wave = self.ssp.wavelength_grid wave = self.ssp.wavelength_grid
w = np.where(wave <= 91.1) w = np.where(wave <= 91.1)
lum_lyc_young = np.trapz(spec_young[w], wave[w]) lum_lyc_young, lum_lyc_old = np.trapz([spec_young[w], spec_old[w]],
lum_lyc_old = np.trapz(spec_old[w], wave[w]) wave[w])
sed.add_module(self.name, self.parameters) sed.add_module(self.name, self.parameters)
......
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