all_results = []
galaxies = []
for item in inits:
ch = ChemModel(**item)
'''
call modules to run the model:
snrate: SN rate at each time step - this also sets time array
so ch.supernova_rate() must be called first to set
time array for the entire code
all results: t, mg, m*, mz, Z, md, md/mz, sfr,
dust_source(all), dust_source(stars),
dust_source(ism), destruction_time, graingrowth_time
'''
snrate = ch.supernova_rate()
all_results = ch.gas_metal_dust_mass(snrate)
# write all the parameters to a dictionary for each init set
params = {'time' : all_results[:,0],
'mgas' : all_results[:,1],
'mstars' : all_results[:,2],
'metalmass' : all_results[:,3],
'metallicity' : all_results[:,4],
'dustmass' : all_results[:,5],
'dust_metals_ratio' : all_results[:,6],
'sfr' : all_results[:,7],
'dust_all' : all_results[:,8],
'dust_stars' : all_results[:,9],
'dust_ism' : all_results[:,10],
'time_destroy' : all_results[:,11],
'time_gg' : all_results[:,12]}
