2. cooling uv heating

test

./pycool_init

1. plot number of ions as a function of Temprature
2. for a given z and n_H

./pycool_NvsT

1. plot cooling functions as a function of temperature
2. for a given z and n_H

./pycool_LHvsT

A FAIRE

• ploter les valeurs de Gamma pour le flux UV, comparer avec les tables de Hardt & Madau

idee:

--> faire des tables hdf5

  dep : Redshift, n_H

combiner :
  - wiersma pour > 4
  - heating (Katz) T<4
  - cooling (Maio) T<4

... utiliser touts les ions pour ameliorer mu

2. cooling with metals (pNbody default)

check init

./pycool_init_cooling_with_metals # ok, init with metals

compute cooling for a given entropy, density, metalicity

./pycool_compute_cooling_with_metals

plot lambda

./pycool_plot_cooling_with_metals --coolingfile cooling.dat

plot single lambda

./pycool_cooling_with_metals_plot_cooling_vs_Temperature

2d map

./pycool_cooling_with_metals_plot_cooling_vs_TemperatureDensity

2d map : cooling time

./pycool_cooling_with_metals_plot_CoolingTime_vs_TemperatureDensity

2. cooling from hdf5 tables ()

check init

./pycool_init_cooling_from_hdf5

get lamda

./pycool_get_cooling_from_hdf5

plot lambda

./pycool_plot_cooling_from_hdf5 # Lambda, as a fct of temperature

./pycool_plot_cooling_from_hdf5_bis # Lambda, as a fct of redshift

2. cooling Wiersma

./pycool_cooling_wiersma_plot_cooling_vs_Temperature ./pycool_cooling_wiersma_plot_cooling_vs_Redshift

2d map

./pycool_cooling_wiersma_plot_cooling_vs_TemperatureDensity --Redshift 10.0

2d map : cooling time

./pycool_cooling_wiersma_plot_CoolingTime_vs_TemperatureDensity --Redshift 0