2. used chemical parameter files ############
   • for dSph tests (dSph-1324-8new !!! Si->Ba : Fe,Mg,O,Si,Metals) ../scripts/pychem_generate_parameters -p chemistryparameters/chimieparam.yr.py -o chimie.yr.dat
     1. hdf5 version ../scripts/pychem_generate_hdf5_parameters -p chemistryparameters/chimieparam.yr.py -o chimie.yr.h5

• for dSph with new code (compatible with paper simulations, only 4 elts, Fe,Mg,O,Metals) ../scripts/pychem_generate_parameters -p chemistryparameters/chimieparam.yr.4.py -o chimie.yr.4.dat

• for dSph2 ("Fe","Mg","O","Ba","Metals") ../scripts/pychem_generate_parameters -p chemistryparameters/chimieparam.yr.2.py -o chimie.yr.2.dat

• salpeter IMF ../scripts/pychem_generate_parameters -p chemistryparameters/chimieparam.salpeter.py -o chimie.salpeter.dat
• 2 different IMF ../scripts/pychem_generate_parameters -p chemistryparameters/chimieparam.yr.py -o chimie.yr.dat.0 ../scripts/pychem_generate_parameters -p chemistryparameters/chimieparam.yr.py.1 -o chimie.yr.dat.1

• yields from Francois (Mon Apr 18 13:11:44 CEST 2011) ../scripts/pychem_generate_parameters -p chemistryparameters/chimieparam.francois.py -o chimie.francois.dat

./plot_SNII_yields.py --log xy --legend --xmin 0.05 --xmax 100 --ymin 1e-4 --ymax 1 --elts=Ej,Ejnp,Fe,Mg,O chimie.yr.dat chimie.francois.dat ./plot_SNII_interated_yields.py --log xy --legend --xmin 0.05 --xmax 100 --ymin 1e-4 --ymax 1 --elts=Ej,Ejnp,Fe,Mg,O chimie.yr.dat chimie.francois.dat

• chimie.francois-met.dat = chimie.francois.dat with Metals from chimie.yr.dat
• yields from Francois + SNIa from Francois ../scripts/pychem_generate_parameters -p chemistryparameters/chimieparam.francois-002.py -o chimie.francois-002.dat

• yields from Woosley (Tue May 3 15:18:43 CEST 2011) ../scripts/pychem_generate_parameters -p chemistryparameters/chimieparam.woosley.py -o chimie.woosley.dat

( * Element Solar Abundances ( becomes ( Element Solar Mass Abundances (=Mx/MH) ( ( ../scripts/pychem_generate_parameters -p chemistryparameters/chimieparam.yr.py -o chimie.yr-new.dat

• truncate SNII tables at 12.9 instead of 10 (Thu Apr 25 15:29:19 CEST 2013) ../scripts/pychem_generate_parameters -p chemistryparameters/chimieparam.yr-25042013.py -o chimie.yr-25042013.dat

• AGB ../scripts/pychem_generate_parameters-31052013 -p chemistryparameters/chimieparam.yr-31052013.py -o chimie.yr-31052013.dat ../scripts/pychem_generate_hdf5_parameters-31052013 -p chemistryparameters/chimieparam.yr-31052013.py -o chimie.yr-31052013.h5
• with Cobalt ../scripts/pychem_generate_hdf5_parameters -p chemistryparameters/chimieparam-AGB+OMgSFeZnSrYBaEuCo-20062014.h5.py -o chemistry-AGB+OMgSFeZnSrYBaEuCo-20062014.h5
• with Ca, Ni, Mn, Cr ../scripts/pychem_generate_hdf5_parameters -p chemistryparameters/chimieparam-AGB+OMgSFeZnSrYBaEuCoCaNiMnCr-12022015.h5.py -o chemistry-AGB+OMgSFeZnSrYBaEuCoCaNiMnCr-12022015.h5

• for agora
  
  cp chemistryparameters/chimieparam-AGB+OMgSFeZnSrYBaEu-16072013.h5.py chemistryparameters/chimieparam-agora-24092015.h5.py cp tables/SolarAbundances/Grevesse98.txt tables/SolarAbundances/Grevesse98+agora.txt ../scripts/pychem_generate_hdf5_parameters -p chemistryparameters/chimieparam-agora-24092015.h5.py -o chemistry-agora-24092015.h5

• (Thu Dec 22 17:11:21 CET 2016) would like to re-create chemistry-AGB+OMgSFeZnSrYBaEu-16072013.h5
  
  

../scripts/pychem_generate_hdf5_parameters -p chemistryparameters/chimieparam-AGB+OMgSFeZnSrYBaEu_TOPH-16072013.h.py -o chemistry-AGB+OMgSFeZnSrYBaEu_TOPH-16072013.h5

# change the stars mass (ensure the same num of SNII compared to a Kroupa IMF)
../scripts/pychem_generate_hdf5_parameters -p chemistryparameters/chimieparam-AGB+OMgSFeZnSrYBaEu_TOPHv2-16072013.h.py -o chemistry-AGB+OMgSFeZnSrYBaEu_TOPHv2-16072013.h5

	# checks

../scripts/pychem_sample_IMF_and_properties chemistry-AGB+OMgSFeZnSrYBaEu_TOPHv2-16072013.h5 --mstar 152200 ../scripts/pychem_plot_SNII_yields --log xy --xmin 0.05 --xmax 1000 --ymin 1e-10 --ymax 1 chemistry-AGB+OMgSFeZnSrYBaEu_TOPHv2-16072013.h5

        ../scripts/pychem_plot_lifetime --log xy --xmin 0.05 --xmax 1000 --ymin 1e-1 --ymax 1e10 chemistry-AGB+OMgSFeZnSrYBaEu_TOPHv2-16072013.h5
        ../scripts/pychem_SSP_discret_evolution3 chemistry-AGB+OMgSFeZnSrYBaEu_TOPHv2-16072013.h5 --mstar 1522 --xmin 0      --xmax 10    --timeunit Myr -o NSNII,CumNSNII,StellarMass --Z=0  --dt .001

2. open and read the chemistry parameter file ############

../scripts/pychem_read chimie.yr.dat

2. plot IMF ############

../scripts/pychem_plot_IMF --log xy chimie.yr.dat chimie.salpeter.dat --legend

2. sample IMF ############

../scripts/pychem_sample_IMF chimie.yr.dat --mstar 1e4

../scripts/pychem_sample_IMF_and_properties chimie.yr.dat --mstar 1e4

2. lifetime ############

Plot the stars lifetime as a function of mass and metallicity

../scripts/pychem_plot_lifetime --legend --log xy --xmin 0.05 --xmax 100 --ymin 1e-1 --ymax 1e10 chimie.yr.dat

../scripts/pychem_lifetime chimie.yr.dat.1

2. Helium core ############

./pychem_plot_HeliumCore.py tables/HeliumCore.dat tables/AGBs/OriginalData/karakas-HeliumCore-0.0200.txt

2. yields with Z ############

only for DYIN

../scripts/pychem_plot_DYIN_yields_withZ --log xy --xmin 0.05 --xmax 100 --ymin 1e-10 --ymax 1 -Z 1e-1 chemistry.simple.h5

total (single) ejection

../scripts/pychem_plot_TOTAL_yields_withZ --log xy --legend --xmin 0.05 --xmax 100 --ymin 1e-4 --ymax 1 --NSNII=1 --NSNIa=0 --NDYIN=1 chemistry.simple.h5 ../scripts/pychem_plot_TOTAL_yields_withZ --log xy --legend --xmin 0.05 --xmax 100 --ymin 1e-4 --ymax 1 --NSNII=1 --NSNIa=0 --NDYIN=0 chemistry.simple.h5 ../scripts/pychem_plot_TOTAL_yields_withZ --log xy --legend --xmin 0.05 --xmax 100 --ymin 1e-4 --ymax 1 --NSNII=0 --NSNIa=1 --NDYIN=0 chemistry.simple.h5 ../scripts/pychem_plot_TOTAL_yields_withZ --log xy --legend --xmin 0.05 --xmax 100 --ymin 1e-4 --ymax 1 --NSNII=1 --NSNIa=1 --NDYIN=1 chemistry.simple.h5

total (integrated) ejection

../scripts/pychem_plot_TOTAL_integrated_yields_withZ --log xy --xmin 0.05 --xmax 100 chemistry.simple.h5

2. yields ############

Plot the Mass fraction of ejected elements due to the explotion of one SNII of mass m

../scripts/pychem_plot_SNII_yields --log xy --legend --xmin 0.05 --xmax 100 --ymin 1e-4 --ymax 1 chimie.yr.dat ../scripts/pychem_plot_DYIN_yields --log xy --legend --xmin 0.05 --xmax 100 --ymin 1e-4 --ymax 1 chimie.yr.dat ../scripts/pychem_plot_SNII_yields --log xy --legend --xmin 0.05 --xmax 100 --ymin 1e-4 --ymax 1 --elts=Ej,Ejnp,Fe56,Mg24,O16,Metals chimie.yr.dat

compute in addition, Mg/Fe ratio

../scripts/pychem_plot_SNII_yields --log y --legend --xmin 8 --xmax 30 --ymin 1e-4 --ymax 1 --MgFe --elts=Mg,Fe chimie.yr.dat

total (single) ejection

../scripts/pychem_plot_TOTAL_yields --log xy --legend --xmin 0.05 --xmax 100 --ymin 1e-4 --ymax 1 --NSNII=1 --NSNIa=0 --NDYIN=1 chimie.yr.dat ../scripts/pychem_plot_TOTAL_yields --log xy --legend --xmin 0.05 --xmax 100 --ymin 1e-4 --ymax 1 --NSNII=1 --NSNIa=0 --NDYIN=0 chimie.yr.dat ../scripts/pychem_plot_TOTAL_yields --log xy --legend --xmin 0.05 --xmax 100 --ymin 1e-4 --ymax 1 --NSNII=0 --NSNIa=1 --NDYIN=0 chimie.yr.dat ../scripts/pychem_plot_TOTAL_yields --log xy --legend --xmin 0.05 --xmax 100 --ymin 1e-4 --ymax 1 --NSNII=1 --NSNIa=1 --NDYIN=1 chimie.yr.dat

total (integrated) ejection

../scripts/pychem_plot_TOTAL_integrated_yields --log xy --xmin 0.05 --xmax 100 chimie.yr.dat

Plot the integrated yields (!! the zero depends on the zero of the function...)

../scripts/pychem_plot_SNII_interated_yields --log xy --legend --xmin 0.05 --xmax 100 --ymin 1e-4 --ymax 1 chimie.yr.dat

Estimation of the energy lost by stars

../scripts/pychem_mass_ejection chimie.yr.dat.1

2. chemical evolution ############

chemical enrichment due to a SSP

../scripts/pychem_SSP_evolution chimie.yr.dat --xmin 0 --xmax 30 --y Mg --timeunit Myr -o Y,ESN ../scripts/pychem_SSP_evolution chimie.yr.dat --xmin 0.001 --xmax 13 --y Mg --timeunit Gyr -o Y,ESN

../scripts/pychem_SSP_evolution chimie.yr.dat --xmin 0.001 --xmax 20 --y Mg --timeunit Gyr -o Y,ESN --log x ../scripts/pychem_SSP_evolution chimie.yr.dat --xmin 1 --xmax 20000 --y Mg --timeunit Myr -o Y,ESN --log x

../scripts/pychem_SSP_evolution chimie.yr.dat --xmin 1 --xmax 20000 --y Mg --timeunit Myr -o D

evolution of the stellar mass

../scripts/pychem_SSP_evolution chimie.yr.dat --xmin 1 --xmax 20000 --y Mg --timeunit Myr -o MSTAR --log x --ymin 0 --ymax 1e-5 --NumberOfTables 2 --DefaultTable 0 ../scripts/pychem_SSP_evolution chimie.yr.dat --xmin 1 --xmax 20000 --y Mg --timeunit Myr -o MSTAR --log x --ymin 0 --ymax 1e-5 --NumberOfTables 2 --DefaultTable 1

(il faudrait un pgm style SSP_evolution.py mais ou l'on plot ce qui nous interesse. ou alors ecrire les output de SSP_evolution.py et le plotter avec un autre pgm...)

using different tables

../scripts/pychem_SSP_evolution chimie.yr.dat --xmin 1 --xmax 20000 --y Mg --timeunit Myr -o Y,ESN --log x --NumberOfTables 2 ../scripts/pychem_SSP_evolution chimie.yr.dat --xmin 1 --xmax 20000 --y Mg --timeunit Myr -o Y,ESN --log x --NumberOfTables 2 --DefaultTable 1

using the SN trick

../scripts/pychem_SSP_evolution chimie.yr.dat --mstar 1e7 --xmin 0 --xmax 50 --y Mg --timeunit Myr -o NSNII,CumNSNII --mstar 3e4 --dt .1 --discsn ../scripts/pychem_SSP_evolution chimie.yr.dat --mstar 1e7 --xmin 0 --xmax 3000 --y Mg --timeunit Myr -o NSNIa,CumNSNIa --mstar 3e4 --dt .1 --discsn

../scripts/pychem_SSP_discret_evolution2 chimie.yr.dat --mstar 1e7 --xmin 0 --discsn --x Fe --y Mg --timeunit Myr --dt .01 --xmax 3000 -o CumNSNII,CumNSNIa,CumNDYIN,CumEjMass,StellarMass,Y,Yg,D

check SN monte carlo

../scripts/pychem_SSP_discret_evolution2 chimie.yr.dat --mstar 105 --xmin 0 --discsn --x Fe --y Mg --timeunit Myr --dt .01 --xmax 3000 -o CumNSNII,CumNSNIa,CumNDYIN,CumEjMass