def temperature(r): MP_OVER_KB = 121.148 HYDROGEN_MASSFRAC = 0.76 meanweight_n = 4.0 / (1 + 3*HYDROGEN_MASSFRAC) meanweight_i = 4.0 / (3 + 5*HYDROGEN_MASSFRAC) integral = integrate.quad(opt.T_integrand, r, np.inf, args=(M_gas, a_gas, M_dm, a_dm, gamma_gas, gamma_dm), full_output=-1) result = integral[0] / opt.dehnen_density(r, M_gas, a_gas, gamma_gas) temp_i = MP_OVER_KB * meanweight_i * result temp_n = MP_OVER_KB * meanweight_n * result if(temp_i > 1.0e4): return temp_to_internal_energy(temp_i) else: return temp_to_internal_energy(temp_n)
def temperature(r): MP_OVER_KB = 121.148 HYDROGEN_MASSFRAC = 0.76 meanweight_n = 4.0 / (1 + 3 * HYDROGEN_MASSFRAC) meanweight_i = 4.0 / (3 + 5 * HYDROGEN_MASSFRAC) integral = integrate.quad(opt.T_integrand, r, np.inf, args=(M_gas, a_gas, M_dm, a_dm, gamma_gas, gamma_dm), full_output=-1) result = integral[0] / opt.dehnen_density(r, M_gas, a_gas, gamma_gas) temp_i = MP_OVER_KB * meanweight_i * result temp_n = MP_OVER_KB * meanweight_n * result if (temp_i > 1.0e4): return temp_to_internal_energy(temp_i) else: return temp_to_internal_energy(temp_n)
def set_temperatures(coords_gas): U_grid = np.zeros((N_rho, Nz)) U = np.zeros(N_gas) # Constantless temperature, will be used in the circular # velocity determination for the gas. T_cl_grid = np.zeros((N_rho, Nz)) MP_OVER_KB = 121.148 HYDROGEN_MASSFRAC = 0.76 meanweight_n = 4.0 / (1 + 3 * HYDROGEN_MASSFRAC) meanweight_i = 4.0 / (3 + 5 * HYDROGEN_MASSFRAC) disk_temp = 10000 U.fill(temp_to_internal_energy(disk_temp)) if(disk_temp >= 1.0e4): T_cl_grid.fill(disk_temp / MP_OVER_KB / meanweight_i) else: T_cl_grid.fill(disk_temp / MP_OVER_KB / meanweight_n) return U, T_cl_grid