Commit 54ad9f5119b6ddae48b70ccdb31eaaae477089f4
1 parent
0917f1d1
Exists in
master
feat correction
Showing
1 changed file
with
4 additions
and
4 deletions
Show diff stats
four_levels_model.py
... | ... | @@ -236,13 +236,13 @@ class HeatingGas: |
236 | 236 | self.pah_cross_dc = self.pah_cross_dc[self.energy_range] |
237 | 237 | |
238 | 238 | ''' yield from anion to neutral ''' |
239 | - part = np.where(self.energy_negative_charged >= 0)[0] | |
239 | + part = np.where(self.energy_negative_charged <=13.6)[0] | |
240 | 240 | |
241 | 241 | self.detachment_yield = np.full(len(part), 1) |
242 | 242 | #the anion being unstable, any energy is enough to detach the electron |
243 | 243 | |
244 | 244 | ''' yield from neutral to the first photoionization ''' |
245 | - first_part = np.where(self.energy_neutral<self.ip_neutral)[0] | |
245 | + first_part = np.where(self.energy_neutral<=self.ip_neutral)[0] | |
246 | 246 | |
247 | 247 | second_part = np.where( (self.energy_neutral>=self.ip_neutral) & (self.energy_neutral<(self.ip_neutral+9.2) ) )[0] |
248 | 248 | third_part = np.where((self.energy_neutral>self.ip_neutral+9.2))[0] |
... | ... | @@ -367,8 +367,8 @@ class HeatingGas: |
367 | 367 | #erg s-1 /(molecule of size n_c) |
368 | 368 | |
369 | 369 | ''' heating rate of the molecule itself ''' |
370 | - heating_pah_a = np.trapz(photo_absorption_a,\ | |
371 | - self.energy) #erg s-1 | |
370 | + heating_pah_a = np.trapz(photo_absorption_a[energy_range_power_absorbed],\ | |
371 | + self.energy[energy_range_power_absorbed]) #erg s-1 | |
372 | 372 | heating_pah_n = np.trapz(photo_absorption_n[energy_range_power_absorbed],\ |
373 | 373 | self.energy[energy_range_power_absorbed]) #erg s-1 |
374 | 374 | heating_pah_c = np.trapz(photo_absorption_c[energy_range_power_absorbed],\ | ... | ... |