Commit aca0789aef8d7d56eff1a89edfdc5f687ef4352a
1 parent
c8415783
Exists in
master
add commentaries
Showing
1 changed file
with
4 additions
and
2 deletions
Show diff stats
four_levels_model.py
| @@ -146,7 +146,8 @@ class HeatingGas: | @@ -146,7 +146,8 @@ class HeatingGas: | ||
| 146 | ''' others parameters to be returned ''' | 146 | ''' others parameters to be returned ''' |
| 147 | self.g_0, self.distance, self.wavelength, self.wavelength_intensity, self.energy_intensity, self.energy, self.ISRF, RF_list = rf.radiation_field(self.filename, self.star_radius, self.parsec, self.ISRF, RF_list=False) | 147 | self.g_0, self.distance, self.wavelength, self.wavelength_intensity, self.energy_intensity, self.energy, self.ISRF, RF_list = rf.radiation_field(self.filename, self.star_radius, self.parsec, self.ISRF, RF_list=False) |
| 148 | self.gamma = ( self.g_0 * np.sqrt(self.t_gas) ) / self.n_e #ionization parameter | 148 | self.gamma = ( self.g_0 * np.sqrt(self.t_gas) ) / self.n_e #ionization parameter |
| 149 | - | 149 | + #in cm3 K^1/2 |
| 150 | + | ||
| 150 | ''' Ionization Potential (IP) estimation ''' | 151 | ''' Ionization Potential (IP) estimation ''' |
| 151 | a = (self.n_c/468)**(1/3) #molecule diameter, in nm | 152 | a = (self.n_c/468)**(1/3) #molecule diameter, in nm |
| 152 | self.ip_neutral = 3.9 + one_in_4_pi_eps_0 * ( ( z_0 + (1/2) ) * (ev**2/a) + ( z_0 + 2 ) * (ev**2/a) *(0.03/a) ) * (1/ev) | 153 | self.ip_neutral = 3.9 + one_in_4_pi_eps_0 * ( ( z_0 + (1/2) ) * (ev**2/a) + ( z_0 + 2 ) * (ev**2/a) *(0.03/a) ) * (1/ev) |
| @@ -390,5 +391,6 @@ class HeatingGas: | @@ -390,5 +391,6 @@ class HeatingGas: | ||
| 390 | ''' gas heating ''' | 391 | ''' gas heating ''' |
| 391 | self.total_gas_heating = total_injected_power * (self.fc_pah/self.n_c) * 2.7e-4 | 392 | self.total_gas_heating = total_injected_power * (self.fc_pah/self.n_c) * 2.7e-4 |
| 392 | #2.7e-4 : elemental abundance of C relative to H (Tielens 2021) | 393 | #2.7e-4 : elemental abundance of C relative to H (Tielens 2021) |
| 393 | - | 394 | + #total gas heating in erg s-1 H-1 |
| 395 | + | ||
| 394 | return self.total_gas_heating, self.heating_efficiency, self.g_0 , self.gamma, self.t_gas, self.n_e, self.n_c | 396 | return self.total_gas_heating, self.heating_efficiency, self.g_0 , self.gamma, self.t_gas, self.n_e, self.n_c |
| 395 | \ No newline at end of file | 397 | \ No newline at end of file |