search

dfielding14 / MultiphaseGalacticWind

Simple code demonstrating how to solve for the steady state structure of a multiphase galactic wind
4 stars 1 forks source link

Cooling time calculation prefactor correction #2

Closed dutta-alankar closed 2 years ago

dutta-alankar commented 2 years ago

Line 111 in MultiphaseGalacticWind.py return 1.5 * (muH/mu)**2 * kb * T / ( nH_actual * Lambda_z0((np.log10(nH),np.log10(T), metallicity))) needs to be changed to return 1.5 * (muH/mu) * kb * T / ( nH_actual * Lambda_z0((np.log10(nH),np.log10(T), metallicity))) Gas adiabatic index gamma is hard coded!

dfielding14 commented 2 years ago

Good catch!