from visualize import visualize #get some bestfit-value from 'Eris' from bestfit_param_omega.find_bestfit_param_v0.set_param_from_eris import * #get values and filename from document from read_parameter_space import read_param #Set cmdline arguments try: from sys import argv n = int(argv[1]) except IndexError: n = 10 #Run calculations of GCE with 'Omega' timesteps = n loa_inflow_vals, save_name = read_param( "inflow_rate") #get values from 'parameter_space.txt' print "Using inflow values: ", loa_inflow_vals bestfit_mgal = 1.0e+10 loa_omega_inst = [ omega(special_timesteps=timesteps, in_out_control=True, inflow_rate=inflow, imf_type=bestfit_imf_type, sfh_array=bestfit_sfh_array, ns_merger_on=bestfit_ns_merger_on, nsmerger_table=bestfit_nsmerger_table, mgal=bestfit_mgal) for inflow in loa_inflow_vals ] #omega-instances with new inflow_rate loa_omega_names = [ "$\dot{M}_{in}$=%1.2e" % inflow for inflow in loa_inflow_vals ]
from visualize import visualize #get some bestfit-value from 'Eris' from bestfit_param_omega.find_bestfit_param_v0.set_param_from_eris import * #get values and filename from document from read_parameter_space import read_param #Set cmdline arguments try: from sys import argv n = int(argv[1]) except IndexError: n = 10 #Run calculations of GCE with 'Omega' timesteps = n loa_mergefrac_vals, save_name = read_param( "f_merger") #get values from 'parameter_space.txt' bestfit_mgal = 4.0e+10 loa_omega_inst = [ omega(special_timesteps=timesteps, f_merger=mergefrac, imf_type=bestfit_imf_type, sfh_array=bestfit_sfh_array, ns_merger_on=bestfit_ns_merger_on, nsmerger_table=bestfit_nsmerger_table, mgal=bestfit_mgal) for mergefrac in loa_mergefrac_vals ] #omega-instances with new outflow loa_omega_names = [ "$F_{merge}\in$%1.2e" % mergefrac for mergefrac in loa_mergefrac_vals ] #visualize masses and sfr with 'visualize'
from visualize import visualize #get some bestfit-value from 'Eris' from bestfit_param_omega.find_bestfit_param_v0.set_param_from_eris import * #get values and filename from document from read_parameter_space import read_param #Set cmdline arguments try: from sys import argv n = int(argv[1]) except IndexError: n = 10 #Run calculations of GCE with 'Omega' timesteps = n loa_binaryfrac_vals, save_name = read_param( "f_binary") #get values from 'parameter_space.txt' bestfit_mgal = 4.0e+10 loa_omega_inst = [ omega(special_timesteps=timesteps, f_binary=binaryfrac, imf_type=bestfit_imf_type, sfh_array=bestfit_sfh_array, ns_merger_on=bestfit_ns_merger_on, nsmerger_table=bestfit_nsmerger_table, mgal=bestfit_mgal) for binaryfrac in loa_binaryfrac_vals ] #omega-instances with new parameter loa_omega_names = [ "$f_{bin}$=%1.2e" % binaryfrac for binaryfrac in loa_binaryfrac_vals ] #visualize masses and sfr with 'visualize'
from visualize import visualize #get some bestfit-value from 'Eris' from bestfit_param_omega.find_bestfit_param_v0.set_param_from_eris import * #get values and filename from document from read_parameter_space import read_param #Set cmdline arguments try: from sys import argv n = int(argv[1]) except IndexError: n = 10 #Run calculations of GCE with 'Omega' timesteps = n loa_coaltime_vals, save_name = read_param( "t_nsm_coal") #get values from 'parameter_space.txt' bestfit_mgal = 4.0e+10 loa_omega_inst = [ omega(special_timesteps=timesteps, t_nsm_coal=coaltime, imf_type=bestfit_imf_type, sfh_array=bestfit_sfh_array, ns_merger_on=bestfit_ns_merger_on, nsmerger_table=bestfit_nsmerger_table, mgal=bestfit_mgal) for coaltime in loa_coaltime_vals ] #omega-instances with new parameter loa_omega_names = [ "$t_{coal,nsm}$=%1.2e" % coaltime for coaltime in loa_coaltime_vals ] #visualize masses and sfr with 'visualize'
from visualize import visualize, eris_data #get some bestfit-value from 'Eris' from bestfit_param_omega.find_bestfit_param_v0.set_param_from_eris import * #get values and filename from document from read_parameter_space import read_param #Set cmdline arguments try: from sys import argv n = int(argv[1]) except IndexError: n = 10 #Run calculations of GCE with 'Omega' timesteps = n loa_norm_vals, save_name = read_param( "sfh_array_norm") #get values from 'parameter_space.txt' bestfit_mgal = 4.0e+10 loa_omega_inst = [ omega(special_timesteps=timesteps, sfh_array_norm=norm, mgal=bestfit_mgal, imf_type=bestfit_imf_type, sfh_array=bestfit_sfh_array, ns_merger_on=bestfit_ns_merger_on, nsmerger_table=bestfit_nsmerger_table) for norm in loa_norm_vals ] #omega-instances with new m_gal loa_omega_names = ["$f_{norm,sfr}$=%1.2f" % norm for norm in loa_norm_vals] norm_eris = eris_data().sfr["m_growth"][-1] print norm_eris loa_omega_inst += [ omega(special_timesteps=timesteps,
from visualize import visualize #get some bestfit-value from 'Eris' from bestfit_param_omega.find_bestfit_param_v0.set_param_from_eris import * #get values and filename from document from read_parameter_space import read_param #Set cmdline arguments try: from sys import argv n = int(argv[1]) except IndexError: n = 10 #Run calculations of GCE with 'Omega' timesteps = n loa_imfboundary_vals, save_name = read_param( "imf_bdys") #get values from 'parameter_space.txt' bestfit_mgal = 4.0e+10 loa_omega_inst = [ omega(special_timesteps=timesteps, imf_bdys=imfboundary, imf_type=bestfit_imf_type, sfh_array=bestfit_sfh_array, ns_merger_on=bestfit_ns_merger_on, nsmerger_table=bestfit_nsmerger_table, mgal=bestfit_mgal) for imfboundary in loa_imfboundary_vals ] #omega-instances with new outflow loa_omega_names = [ "$IMF\in$%s" % imf_bound for imf_bound in loa_imfboundary_vals ] #visualize masses and sfr with 'visualize'
from visualize import visualize #get some bestfit-value from 'Eris' from bestfit_param_omega.find_bestfit_param_v0.set_param_from_eris import * #get values and filename from document from read_parameter_space import read_param #Set cmdline arguments try: from sys import argv n = int(argv[1]) except IndexError: n = 10 #Run calculations of GCE with 'Omega' timesteps = n loa_pop3yields_vals, save_name = read_param( "imf_yields_range_pop3") #get values from 'parameter_space.txt' bestfit_mgal = 4.0e+10 loa_omega_inst = [ omega(special_timesteps=timesteps, imf_yields_range_pop3=pop3yields, imf_type=bestfit_imf_type, sfh_array=bestfit_sfh_array, ns_merger_on=bestfit_ns_merger_on, nsmerger_table=bestfit_nsmerger_table, mgal=bestfit_mgal) for pop3yields in loa_pop3yields_vals ] #omega-instances with new parameter loa_omega_names = [ "$Y_{pop3}\in$%s" % pop3yields for pop3yields in loa_pop3yields_vals ] #visualize masses and sfr with 'visualize'
from visualize import visualize #get some bestfit-value from 'Eris' from bestfit_param_omega.find_bestfit_param_v0.set_param_from_eris import * #get values and filename from document from read_parameter_space import read_param #Set cmdline arguments try: from sys import argv n = int(argv[1]) except IndexError: n = 10 #Run calculations of GCE with 'Omega' timesteps = n loa_mgal_vals, save_name = read_param("mgal") #get values from 'parameter_space.txt' loa_omega_inst = [omega(special_timesteps=timesteps, mgal=mgal, imf_type=bestfit_imf_type, sfh_array=bestfit_sfh_array, ns_merger_on=bestfit_ns_merger_on, nsmerger_table=bestfit_nsmerger_table) for mgal in loa_mgal_vals] #omega-instances with new m_gal loa_omega_names = ["$M_{gal}$=%1.2e"%mgal for mgal in loa_mgal_vals] #visualize masses and sfr with 'visualize' title = "Vary initial galactic mass $M_{gal}(t_0)$" #plot sfr, ism-mass, locked_mass, total_mass plot_obj = visualize(loa_omega_inst, loa_omega_names, num_yaxes=4) plot_obj.add_time_mass("total", index_yaxis=0) plot_obj.add_time_mass("ism", index_yaxis=1) plot_obj.add_time_mass("locked", index_yaxis=2) plot_obj.add_time_rate("sf", index_yaxis=3) plot_obj.finalize(show=True, title=title, save=save_name)
from visualize import visualize #get some bestfit-value from 'Eris' from bestfit_param_omega.find_bestfit_param_v0.set_param_from_eris import * #get values and filename from document from read_parameter_space import read_param #Set cmdline arguments try: from sys import argv n = int(argv[1]) except IndexError: n = 10 #Run calculations of GCE with 'Omega' timesteps = n loa_pop3bdys_vals, save_name = read_param( "imf_bdys_pop3") #get values from 'parameter_space.txt' bestfit_mgal = 4.0e+10 loa_omega_inst = [ omega(special_timesteps=timesteps, imf_bdys_pop3=pop3bdys, imf_type=bestfit_imf_type, sfh_array=bestfit_sfh_array, ns_merger_on=bestfit_ns_merger_on, nsmerger_table=bestfit_nsmerger_table, mgal=bestfit_mgal) for pop3bdys in loa_pop3bdys_vals ] #omega-instances with new parameter loa_omega_names = [ "$IMF_{pop3}\in$%s" % pop3bdys for pop3bdys in loa_pop3bdys_vals ] #visualize masses and sfr with 'visualize'
from visualize import visualize #get some bestfit-value from 'Eris' from bestfit_param_omega.find_bestfit_param_v0.set_param_from_eris import * #get values and filename from document from read_parameter_space import read_param #Set cmdline arguments try: from sys import argv n = int(argv[1]) except IndexError: n = 10 #Run calculations of GCE with 'Omega' timesteps = n loa_eject_vals, save_name = read_param("m_ej_nsm") #get values from 'parameter_space.txt' bestfit_mgal = 4.0e+10 loa_omega_inst = [omega(special_timesteps=timesteps, m_ej_nsm=eject, imf_type=bestfit_imf_type, sfh_array=bestfit_sfh_array, ns_merger_on=bestfit_ns_merger_on, nsmerger_table=bestfit_nsmerger_table, mgal=bestfit_mgal) for eject in loa_eject_vals] #omega-instances with new parameter loa_omega_names = ["$\dot{M}_{NSM}$=%1.2e"%eject for eject in loa_eject_vals] #visualize masses and sfr with 'visualize' title = "Vary ejecta mass from NSM event" #plot sfr, ism-mass, locked_mass, total_mass plot_obj = visualize(loa_omega_inst, loa_omega_names, num_yaxes=4, yields=True) plot_obj.add_time_relabu("[Eu/H]", index_yaxis=0) plot_obj.add_yields("Eu", index_yaxis=1, time="sum")
from visualize import visualize #get some bestfit-value from 'Eris' from bestfit_param_omega.find_bestfit_param_v0.set_param_from_eris import * #get values and filename from document from read_parameter_space import read_param #Set cmdline arguments try: from sys import argv n = int(argv[1]) except IndexError: n = 10 #Run calculations of GCE with 'Omega' timesteps = n loa_time_vals, save_name = read_param( "t_merger_max") #get values from 'parameter_space.txt' bestfit_mgal = 4.0e+10 loa_omega_inst = [ omega(special_timesteps=timesteps, t_merger_max=time, imf_type=bestfit_imf_type, sfh_array=bestfit_sfh_array, ns_merger_on=bestfit_ns_merger_on, nsmerger_table=bestfit_nsmerger_table, mgal=bestfit_mgal) for time in loa_time_vals ] #omega-instances with new parameter loa_omega_names = ["Max($t_{merge}$)=%1.2e" % time for time in loa_time_vals] #visualize masses and sfr with 'visualize' title = "Vary maximum merging time" #plot sfr, ism-mass, locked_mass, total_mass
from visualize import visualize #get some bestfit-value from 'Eris' from bestfit_param_omega.find_bestfit_param_v0.set_param_from_eris import * #get values and filename from document from read_parameter_space import read_param #Set cmdline arguments try: from sys import argv n = int(argv[1]) except IndexError: n = 10 #Run calculations of GCE with 'Omega' timesteps = n loa_nbnsm_vals, save_name = read_param( "nb_nsm_per_m") #get values from 'parameter_space.txt' bestfit_mgal = 3.0e+10 bestfit_t_nsm_coal = 1e+9 loa_omega_inst = [ omega(special_timesteps=timesteps, nb_nsm_per_m=nbnsm, imf_type=bestfit_imf_type, sfh_array=bestfit_sfh_array, ns_merger_on=bestfit_ns_merger_on, nsmerger_table=bestfit_nsmerger_table, mgal=bestfit_mgal) for nbnsm in loa_nbnsm_vals ] #omega-instances with new outflow loa_omega_names = ["#NSM/M=%1.2f" % nbnsm for nbnsm in loa_nbnsm_vals] #visualize masses and sfr with 'visualize' title = "Vary number of NSM per mass"
from visualize import visualize #get some bestfit-value from 'Eris' from bestfit_param_omega.find_bestfit_param_v0.set_param_from_eris import * #get values and filename from document from read_parameter_space import read_param #Set cmdline arguments try: from sys import argv n = int(argv[1]) except IndexError: n = 10 #Run calculations of GCE with 'Omega' timesteps = n loa_fraction_vals, save_name = read_param( "f_arfo") #get values from 'parameter_space.txt' bestfit_mgal = 3.5e+10 loa_omega_inst = [ omega(special_timesteps=timesteps, f_arfo=fraction, imf_type=bestfit_imf_type, sfh_array=bestfit_sfh_array, ns_merger_on=bestfit_ns_merger_on, nsmerger_table=bestfit_nsmerger_table, mgal=bestfit_mgal) for fraction in loa_fraction_vals ] #omega-instances with new m_gal loa_omega_names = ["f=%1.2e" % fraction for fraction in loa_fraction_vals] #visualize masses and sfr with 'visualize' title = "Vary fractional yield of massive stars" #plot sfr, ism-mass, locked_mass, total_mass
from visualize import visualize #get some bestfit-value from 'Eris' from bestfit_param_omega.find_bestfit_param_v0.set_param_from_eris import * #get values and filename from document from read_parameter_space import read_param #Set cmdline arguments try: from sys import argv n = int(argv[1]) except IndexError: n = 10 #Run calculations of GCE with 'Omega' timesteps = n loa_massload_vals, save_name = read_param( "mass_loading") #get values from 'parameter_space.txt' bestfit_mgal = 4.0e+10 loa_omega_inst = [ omega(special_timesteps=timesteps, in_out_control=True, mass_loading=massload, out_follows_E_rate=False, in_out_ratio=0.0, imf_type=bestfit_imf_type, sfh_array=bestfit_sfh_array, ns_merger_on=bestfit_ns_merger_on, nsmerger_table=bestfit_nsmerger_table, mgal=bestfit_mgal) for massload in loa_massload_vals ] #omega-instances with new outflow loa_omega_names = [ "$\dot{M}_{out}$=%1.2e*SFR" % massload for massload in loa_massload_vals