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
