#############################################
#             Nyx THERMAL GRID              #
# Contact: J. Onorbe jose.onnorbe@gmail.com #
#############################################

HM12 = Haardt & Madau 2012
FG = Faucher-Giguere et al. 2009 (last update available in the web: Deceember 2011)
HM01 = Haardt & Madau 2001
HM96 = Haardt & Madau 1996

A and B parameters stand for factors used to modify the HM12 photoheating
rates (a la Bolton et al., Becker et al., etc):
    heating_rates= A * Delta**(B) * heatin_rates_HM12
where Delta is the overdensity.
 
###### Thermal Models 

T000: HM12 (so A=1 B=0)

T000: HM12 (so A=1 B=0)
T001: A=2.0 B=0.0
T002: A=3.0 B=0.0
T003: A=0.5 B=0.0
#T004: A=0.3 B=0.0
T005: A=1.0 B=-0.5
T006: A=2.0 B=-0.5
T007: A=3.0 B=-0.5
T008: A=0.5 B=-0.5
#T009: A=0.3 B=-0.5
T010: A=1.0 B=-1.0
T011: A=2.0 B=-1.0
T012: A=3.0 B=-1.0
T013: A=0.5 B=-1.0
#T014: A=0.3 B=-1.0 
T015: A=4.0 B=-1.0 
#T016: A=4.0 B=0.0
T020: A=0.5 B=0.25
T021: A=1.5 B=0.25
T022: A=2.5 B=0.25
T023: A=1.5 B=-0.25
T024: A=2.5 B=-0.25
T025: A=1.5 B=-0.75
T026: A=2.5 B=-0.75
T027: A=0.5 B=0.5
T028: A=1.5 B=0.5
T029: A=2.5 B=0.5
#T030: A=0.5 B=-0.25
#T031: A=0.5 B=-0.75
#T032: A=0.25 B=0.0
#T033: A=0.25 B=-0.5
#T034: A=0.25 B=-0.25
#T034: A=0.25 B=0.25 
#T036: A=0.25 B=0.5


T100: A=1.0 B=0.0 starts at z=7
T103: A=0.5 B=0.0 starts at z=7
T170: A=1.0 B=0.0 starts at z=9
T190: A=1.0 B=0.0 starts at z=11

T999: FG
T998: HM01
T997: HM96
T996: HM12 NO helium 

T632: Late HI reionization model  (standard heating) Use 632 instead of T692
T693: Late HI reionization model (He warm)
T694: Late HI reionization model (He cold) 
T695: Late HI reionization model (He hot)
T696: Late HI reionization model (H very cold, He cold)
T697: Late HI reionization model (H very cold, He standard)
T698: Late HI reionization model (H very cold, He hot)
T699: Late HI reionization model (H warm, He standard)
T690: Late HI reionization model (H hot, He standard)


T682: Early HI reionization model (standard heating)
T687: Early HI reionization model (H hot) 
T688: Early HI reionization model (H hot He hot)

T671: Middle HI reionization (standard heating) No HeII model
T672: Middle HI reionization model (standard heating)
T673: Middle HI reionization model ( He warm) 
T674: Middle HI reionization model (He cold)
T675: Middle HI reionization model (He hot)

T665: Middle HI reionization model (H very cold, He cold)
T666: Middle HI reionization model (H hot He hot)

T652: Very late reion (standard heating)
T657: Very late reion (Ultra cold H and He)

T641: Fast z=6.6 HI reionization (standard heating)

T649: Reionization 1 sigma above Planck 2016
T643: Late HI reionization. Very Early HeII reionization, standard heating (z=5) (e.g.Madau+2015)

T633: Late HI reionization, He z=2.8 (standard heating)
T635: Late HI reionization, He z=2.8 (He hot)