In [1]:
%matplotlib inline
import yt
from galaxy_analysis.analysis import Galaxy
from galaxy_analysis.yt_fields import field_generators as fg
import matplotlib.pyplot as plt
import numpy as np

In [2]:
#
#
#

dsname = 'DD0088'
wdir   = '~/work/enzo_runs/Wise12/Nested-Test-L2'

#gal = Galaxy(dsname, wdir = wdir)

In [3]:
ds = yt.load(wdir + '/' + dsname + '/' + dsname)
data = ds.all_data()
print(ds.domain_width.to('kpc'))


fg.generate_derived_fields(ds)

ds = yt.load(wdir + '/' + dsname + '/' + dsname)
data = ds.all_data()

fg.generate_particle_filters(ds)


[102.51652059 102.51652059 102.51652059] kpc
Parsing Hierarchy : 100%|██████████| 1638/1638 [00:00<00:00, 11002.98it/s]
/home/aemerick/.local/lib/python3.7/site-packages/yt/fields/local_fields.py:46: UserWarning: Because 'sampling_type' not specified, yt will assume a cell 'sampling_type'
  warnings.warn("Because 'sampling_type' not specified, yt will "
tracer species present:  ['C', 'N', 'O', 'Na', 'Mg', 'Si', 'S', 'Ca', 'Mn', 'Fe', 'Ni', 'As', 'Sr', 'Y', 'Ba']
16 mass fields defined
16 mass fraction fields defined
15 number density fields defined
141 abundance ratio fields defined
5 additional helper fields defined
Parsing Hierarchy : 100%|██████████| 1638/1638 [00:00<00:00, 10476.65it/s]
/home/aemerick/code/onezone/onezone/data_tables.py:98: FutureWarning: elementwise comparison failed; returning scalar instead, but in the future will perform elementwise comparison
  if  c == flag or id == flag:

In [78]:
ds.print_stats()
ds.current_time.to('Myr')


level	# grids	       # cells	     # cells^3
----------------------------------------------
  0	    64	        262144	            64
  1	    16	        195112	            58
  2	    41	        512000	            80
  3	   199	       1319752	           110
  4	   373	       1085016	           103
  5	   324	        503704	            80
  6	   174	        153352	            54
  7	    33	         19896	            28
  8	     1	           216	             6
----------------------------------------------
   	  1225	       4051192


t = 1.95089178e+01 = 1.32238181e+16 s = 4.19037509e+08 years

Smallest Cell:
	Width: 5.014e-06 Mpc
	Width: 5.014e+00 pc
	Width: 1.034e+06 AU
	Width: 1.547e+19 cm
Out[78]:
419.03750890206084 Myr

In [50]:
ds.derived_field_list


Out[50]:
[('all', 'age'),
 ('all', 'birth_mass'),
 ('all', 'creation_time'),
 ('all', 'dynamical_time'),
 ('all', 'mesh_id'),
 ('all', 'metallicity_fraction'),
 ('all', 'particle_angular_momentum'),
 ('all', 'particle_angular_momentum_magnitude'),
 ('all', 'particle_angular_momentum_x'),
 ('all', 'particle_angular_momentum_y'),
 ('all', 'particle_angular_momentum_z'),
 ('all', 'particle_cylindrical_velocity_theta'),
 ('all', 'particle_cylindrical_velocity_z'),
 ('all', 'particle_index'),
 ('all', 'particle_mass'),
 ('all', 'particle_ones'),
 ('all', 'particle_position'),
 ('all', 'particle_position_cylindrical_radius'),
 ('all', 'particle_position_cylindrical_theta'),
 ('all', 'particle_position_cylindrical_z'),
 ('all', 'particle_position_relative'),
 ('all', 'particle_position_relative_x'),
 ('all', 'particle_position_relative_y'),
 ('all', 'particle_position_relative_z'),
 ('all', 'particle_position_spherical_phi'),
 ('all', 'particle_position_spherical_radius'),
 ('all', 'particle_position_spherical_theta'),
 ('all', 'particle_position_x'),
 ('all', 'particle_position_y'),
 ('all', 'particle_position_z'),
 ('all', 'particle_radial_velocity'),
 ('all', 'particle_radius'),
 ('all', 'particle_specific_angular_momentum'),
 ('all', 'particle_specific_angular_momentum_x'),
 ('all', 'particle_specific_angular_momentum_y'),
 ('all', 'particle_specific_angular_momentum_z'),
 ('all', 'particle_spherical_position_phi'),
 ('all', 'particle_spherical_position_radius'),
 ('all', 'particle_spherical_position_theta'),
 ('all', 'particle_spherical_velocity_phi'),
 ('all', 'particle_spherical_velocity_radius'),
 ('all', 'particle_spherical_velocity_theta'),
 ('all', 'particle_type'),
 ('all', 'particle_velocity'),
 ('all', 'particle_velocity_cylindrical_radius'),
 ('all', 'particle_velocity_cylindrical_theta'),
 ('all', 'particle_velocity_cylindrical_z'),
 ('all', 'particle_velocity_magnitude'),
 ('all', 'particle_velocity_relative'),
 ('all', 'particle_velocity_relative_x'),
 ('all', 'particle_velocity_relative_y'),
 ('all', 'particle_velocity_relative_z'),
 ('all', 'particle_velocity_spherical_phi'),
 ('all', 'particle_velocity_spherical_radius'),
 ('all', 'particle_velocity_spherical_theta'),
 ('all', 'particle_velocity_x'),
 ('all', 'particle_velocity_y'),
 ('all', 'particle_velocity_z'),
 ('all', 'relative_particle_position'),
 ('all', 'relative_particle_position_x'),
 ('all', 'relative_particle_position_y'),
 ('all', 'relative_particle_position_z'),
 ('all', 'relative_particle_velocity'),
 ('all', 'relative_particle_velocity_x'),
 ('all', 'relative_particle_velocity_y'),
 ('all', 'relative_particle_velocity_z'),
 ('all', 'sn_mass_ejected'),
 ('all', 'wind_mass_ejected'),
 ('deposit', 'all_cic'),
 ('deposit', 'all_cic_age'),
 ('deposit', 'all_cic_velocity_x'),
 ('deposit', 'all_cic_velocity_y'),
 ('deposit', 'all_cic_velocity_z'),
 ('deposit', 'all_count'),
 ('deposit', 'all_density'),
 ('deposit', 'all_mass'),
 ('deposit', 'all_nn_age'),
 ('deposit', 'all_nn_velocity_x'),
 ('deposit', 'all_nn_velocity_y'),
 ('deposit', 'all_nn_velocity_z'),
 ('deposit', 'io_cic'),
 ('deposit', 'io_cic_age'),
 ('deposit', 'io_cic_velocity_x'),
 ('deposit', 'io_cic_velocity_y'),
 ('deposit', 'io_cic_velocity_z'),
 ('deposit', 'io_count'),
 ('deposit', 'io_density'),
 ('deposit', 'io_mass'),
 ('deposit', 'io_nn_age'),
 ('deposit', 'io_nn_velocity_x'),
 ('deposit', 'io_nn_velocity_y'),
 ('deposit', 'io_nn_velocity_z'),
 ('deposit', 'low_mass_stars_cic'),
 ('deposit', 'low_mass_stars_cic_age'),
 ('deposit', 'low_mass_stars_cic_velocity_x'),
 ('deposit', 'low_mass_stars_cic_velocity_y'),
 ('deposit', 'low_mass_stars_cic_velocity_z'),
 ('deposit', 'low_mass_stars_count'),
 ('deposit', 'low_mass_stars_density'),
 ('deposit', 'low_mass_stars_mass'),
 ('deposit', 'low_mass_stars_nn_age'),
 ('deposit', 'low_mass_stars_nn_velocity_x'),
 ('deposit', 'low_mass_stars_nn_velocity_y'),
 ('deposit', 'low_mass_stars_nn_velocity_z'),
 ('deposit', 'low_mass_unresolved_stars_cic'),
 ('deposit', 'low_mass_unresolved_stars_cic_age'),
 ('deposit', 'low_mass_unresolved_stars_cic_velocity_x'),
 ('deposit', 'low_mass_unresolved_stars_cic_velocity_y'),
 ('deposit', 'low_mass_unresolved_stars_cic_velocity_z'),
 ('deposit', 'low_mass_unresolved_stars_count'),
 ('deposit', 'low_mass_unresolved_stars_density'),
 ('deposit', 'low_mass_unresolved_stars_mass'),
 ('deposit', 'low_mass_unresolved_stars_nn_age'),
 ('deposit', 'low_mass_unresolved_stars_nn_velocity_x'),
 ('deposit', 'low_mass_unresolved_stars_nn_velocity_y'),
 ('deposit', 'low_mass_unresolved_stars_nn_velocity_z'),
 ('deposit', 'main_sequence_popIII_stars_cic'),
 ('deposit', 'main_sequence_popIII_stars_cic_age'),
 ('deposit', 'main_sequence_popIII_stars_cic_velocity_x'),
 ('deposit', 'main_sequence_popIII_stars_cic_velocity_y'),
 ('deposit', 'main_sequence_popIII_stars_cic_velocity_z'),
 ('deposit', 'main_sequence_popIII_stars_count'),
 ('deposit', 'main_sequence_popIII_stars_density'),
 ('deposit', 'main_sequence_popIII_stars_mass'),
 ('deposit', 'main_sequence_popIII_stars_nn_age'),
 ('deposit', 'main_sequence_popIII_stars_nn_velocity_x'),
 ('deposit', 'main_sequence_popIII_stars_nn_velocity_y'),
 ('deposit', 'main_sequence_popIII_stars_nn_velocity_z'),
 ('deposit', 'main_sequence_stars_cic'),
 ('deposit', 'main_sequence_stars_cic_age'),
 ('deposit', 'main_sequence_stars_cic_velocity_x'),
 ('deposit', 'main_sequence_stars_cic_velocity_y'),
 ('deposit', 'main_sequence_stars_cic_velocity_z'),
 ('deposit', 'main_sequence_stars_count'),
 ('deposit', 'main_sequence_stars_density'),
 ('deposit', 'main_sequence_stars_mass'),
 ('deposit', 'main_sequence_stars_nn_age'),
 ('deposit', 'main_sequence_stars_nn_velocity_x'),
 ('deposit', 'main_sequence_stars_nn_velocity_y'),
 ('deposit', 'main_sequence_stars_nn_velocity_z'),
 ('deposit', 'remnant_stars_cic'),
 ('deposit', 'remnant_stars_cic_age'),
 ('deposit', 'remnant_stars_cic_velocity_x'),
 ('deposit', 'remnant_stars_cic_velocity_y'),
 ('deposit', 'remnant_stars_cic_velocity_z'),
 ('deposit', 'remnant_stars_count'),
 ('deposit', 'remnant_stars_density'),
 ('deposit', 'remnant_stars_mass'),
 ('deposit', 'remnant_stars_nn_age'),
 ('deposit', 'remnant_stars_nn_velocity_x'),
 ('deposit', 'remnant_stars_nn_velocity_y'),
 ('deposit', 'remnant_stars_nn_velocity_z'),
 ('enzo', 'AGB_Metal_Density'),
 ('enzo', 'As_Density'),
 ('enzo', 'Ba_Density'),
 ('enzo', 'C_Density'),
 ('enzo', 'Ca_Density'),
 ('enzo', 'Dark_Matter_Density'),
 ('enzo', 'Density'),
 ('enzo', 'Electron_Density'),
 ('enzo', 'FUV_FluxDensity'),
 ('enzo', 'Fe_Density'),
 ('enzo', 'GasEnergy'),
 ('enzo', 'H2II_Density'),
 ('enzo', 'H2II_kdiss'),
 ('enzo', 'H2I_Density'),
 ('enzo', 'H2I_kdiss'),
 ('enzo', 'HII_Density'),
 ('enzo', 'HI_Density'),
 ('enzo', 'HI_kph'),
 ('enzo', 'HM_Density'),
 ('enzo', 'HM_kph'),
 ('enzo', 'HeIII_Density'),
 ('enzo', 'HeII_Density'),
 ('enzo', 'HeII_kph'),
 ('enzo', 'HeI_Density'),
 ('enzo', 'HeI_kph'),
 ('enzo', 'Metal_Density'),
 ('enzo', 'Mg_Density'),
 ('enzo', 'Mn_Density'),
 ('enzo', 'N_Density'),
 ('enzo', 'Na_Density'),
 ('enzo', 'Ni_Density'),
 ('enzo', 'O_Density'),
 ('enzo', 'Pe_heating_rate'),
 ('enzo', 'PhotoGamma'),
 ('enzo', 'PopIII_Metal_Density'),
 ('enzo', 'RadAccel1'),
 ('enzo', 'RadAccel2'),
 ('enzo', 'RadAccel3'),
 ('enzo', 'Ray_Segments'),
 ('enzo', 'SNII_Metal_Density'),
 ('enzo', 'SNIa_Metal_Density'),
 ('enzo', 'S_Density'),
 ('enzo', 'Si_Density'),
 ('enzo', 'Sr_Density'),
 ('enzo', 'Temperature'),
 ('enzo', 'TotalEnergy'),
 ('enzo', 'Y_Density'),
 ('enzo', 'cell_volume'),
 ('enzo', 'dx'),
 ('enzo', 'dy'),
 ('enzo', 'dz'),
 ('enzo', 'path_element_x'),
 ('enzo', 'path_element_y'),
 ('enzo', 'path_element_z'),
 ('enzo', 'vertex_x'),
 ('enzo', 'vertex_y'),
 ('enzo', 'vertex_z'),
 ('enzo', 'x'),
 ('enzo', 'x-velocity'),
 ('enzo', 'y'),
 ('enzo', 'y-velocity'),
 ('enzo', 'z'),
 ('enzo', 'z-velocity'),
 ('gas', 'As_Abundance'),
 ('gas', 'As_Density'),
 ('gas', 'As_Fraction'),
 ('gas', 'As_Mass'),
 ('gas', 'As_Number_Density'),
 ('gas', 'As_over_Ba'),
 ('gas', 'As_over_C'),
 ('gas', 'As_over_Fe'),
 ('gas', 'As_over_H'),
 ('gas', 'As_over_Mg'),
 ('gas', 'As_over_N'),
 ('gas', 'As_over_Na'),
 ('gas', 'As_over_Ni'),
 ('gas', 'As_over_O'),
 ('gas', 'Ba_Abundance'),
 ('gas', 'Ba_Density'),
 ('gas', 'Ba_Fraction'),
 ('gas', 'Ba_Mass'),
 ('gas', 'Ba_Number_Density'),
 ('gas', 'Ba_over_As'),
 ('gas', 'Ba_over_C'),
 ('gas', 'Ba_over_Ca'),
 ('gas', 'Ba_over_Fe'),
 ('gas', 'Ba_over_H'),
 ('gas', 'Ba_over_Mg'),
 ('gas', 'Ba_over_Mn'),
 ('gas', 'Ba_over_N'),
 ('gas', 'Ba_over_Na'),
 ('gas', 'Ba_over_Ni'),
 ('gas', 'Ba_over_O'),
 ('gas', 'Ba_over_S'),
 ('gas', 'Ba_over_Si'),
 ('gas', 'Ba_over_Sr'),
 ('gas', 'Ba_over_Y'),
 ('gas', 'C_Abundance'),
 ('gas', 'C_Density'),
 ('gas', 'C_Fraction'),
 ('gas', 'C_Mass'),
 ('gas', 'C_Number_Density'),
 ('gas', 'C_over_Ba'),
 ('gas', 'C_over_Fe'),
 ('gas', 'C_over_H'),
 ('gas', 'C_over_Mg'),
 ('gas', 'C_over_N'),
 ('gas', 'C_over_Na'),
 ('gas', 'C_over_Ni'),
 ('gas', 'C_over_O'),
 ('gas', 'Ca_Abundance'),
 ('gas', 'Ca_Density'),
 ('gas', 'Ca_Fraction'),
 ('gas', 'Ca_Mass'),
 ('gas', 'Ca_Number_Density'),
 ('gas', 'Ca_over_Ba'),
 ('gas', 'Ca_over_C'),
 ('gas', 'Ca_over_Fe'),
 ('gas', 'Ca_over_H'),
 ('gas', 'Ca_over_Mg'),
 ('gas', 'Ca_over_N'),
 ('gas', 'Ca_over_Na'),
 ('gas', 'Ca_over_Ni'),
 ('gas', 'Ca_over_O'),
 ('gas', 'El_density'),
 ('gas', 'El_fraction'),
 ('gas', 'El_mass'),
 ('gas', 'El_number_density'),
 ('gas', 'FUV_flux'),
 ('gas', 'Fe_Abundance'),
 ('gas', 'Fe_Density'),
 ('gas', 'Fe_Fraction'),
 ('gas', 'Fe_Mass'),
 ('gas', 'Fe_Number_Density'),
 ('gas', 'Fe_over_As'),
 ('gas', 'Fe_over_Ba'),
 ('gas', 'Fe_over_C'),
 ('gas', 'Fe_over_Ca'),
 ('gas', 'Fe_over_H'),
 ('gas', 'Fe_over_Mg'),
 ('gas', 'Fe_over_Mn'),
 ('gas', 'Fe_over_N'),
 ('gas', 'Fe_over_Na'),
 ('gas', 'Fe_over_Ni'),
 ('gas', 'Fe_over_O'),
 ('gas', 'Fe_over_S'),
 ('gas', 'Fe_over_Si'),
 ('gas', 'Fe_over_Sr'),
 ('gas', 'Fe_over_Y'),
 ('gas', 'G_eff'),
 ('gas', 'G_o'),
 ('gas', 'H2_density'),
 ('gas', 'H2_fraction'),
 ('gas', 'H2_mass'),
 ('gas', 'H2_number_density'),
 ('gas', 'H2_p0_density'),
 ('gas', 'H2_p0_fraction'),
 ('gas', 'H2_p0_mass'),
 ('gas', 'H2_p0_number_density'),
 ('gas', 'H2_p1_density'),
 ('gas', 'H2_p1_fraction'),
 ('gas', 'H2_p1_mass'),
 ('gas', 'H2_p1_number_density'),
 ('gas', 'H_Abundance'),
 ('gas', 'H_Mass'),
 ('gas', 'H_density'),
 ('gas', 'H_fraction'),
 ('gas', 'H_m1_density'),
 ('gas', 'H_m1_fraction'),
 ('gas', 'H_m1_mass'),
 ('gas', 'H_m1_number_density'),
 ('gas', 'H_mass'),
 ('gas', 'H_nuclei_density'),
 ('gas', 'H_number_density'),
 ('gas', 'H_p0_density'),
 ('gas', 'H_p0_fraction'),
 ('gas', 'H_p0_mass'),
 ('gas', 'H_p0_number_density'),
 ('gas', 'H_p1_density'),
 ('gas', 'H_p1_fraction'),
 ('gas', 'H_p1_mass'),
 ('gas', 'H_p1_number_density'),
 ('gas', 'H_total_mass'),
 ('gas', 'H_total_number_density'),
 ('gas', 'He_density'),
 ('gas', 'He_fraction'),
 ('gas', 'He_mass'),
 ('gas', 'He_nuclei_density'),
 ('gas', 'He_number_density'),
 ('gas', 'He_p0_density'),
 ('gas', 'He_p0_fraction'),
 ('gas', 'He_p0_mass'),
 ('gas', 'He_p0_number_density'),
 ('gas', 'He_p1_density'),
 ('gas', 'He_p1_fraction'),
 ('gas', 'He_p1_mass'),
 ('gas', 'He_p1_number_density'),
 ('gas', 'He_p2_density'),
 ('gas', 'He_p2_fraction'),
 ('gas', 'He_p2_mass'),
 ('gas', 'He_p2_number_density'),
 ('gas', 'He_total_mass'),
 ('gas', 'LW_flux'),
 ('gas', 'Metal_Density'),
 ('gas', 'Metal_Number_Density'),
 ('gas', 'Mg_Abundance'),
 ('gas', 'Mg_Density'),
 ('gas', 'Mg_Fraction'),
 ('gas', 'Mg_Mass'),
 ('gas', 'Mg_Number_Density'),
 ('gas', 'Mg_over_Ba'),
 ('gas', 'Mg_over_C'),
 ('gas', 'Mg_over_Fe'),
 ('gas', 'Mg_over_H'),
 ('gas', 'Mg_over_N'),
 ('gas', 'Mg_over_Na'),
 ('gas', 'Mg_over_Ni'),
 ('gas', 'Mg_over_O'),
 ('gas', 'Mn_Abundance'),
 ('gas', 'Mn_Density'),
 ('gas', 'Mn_Fraction'),
 ('gas', 'Mn_Mass'),
 ('gas', 'Mn_Number_Density'),
 ('gas', 'Mn_over_Ba'),
 ('gas', 'Mn_over_C'),
 ('gas', 'Mn_over_Fe'),
 ('gas', 'Mn_over_H'),
 ('gas', 'Mn_over_Mg'),
 ('gas', 'Mn_over_N'),
 ('gas', 'Mn_over_Na'),
 ('gas', 'Mn_over_Ni'),
 ('gas', 'Mn_over_O'),
 ('gas', 'N_Abundance'),
 ('gas', 'N_Density'),
 ('gas', 'N_Fraction'),
 ('gas', 'N_Mass'),
 ('gas', 'N_Number_Density'),
 ('gas', 'N_over_Ba'),
 ('gas', 'N_over_C'),
 ('gas', 'N_over_Fe'),
 ('gas', 'N_over_H'),
 ('gas', 'N_over_Mg'),
 ('gas', 'N_over_Na'),
 ('gas', 'N_over_Ni'),
 ('gas', 'N_over_O'),
 ('gas', 'Na_Abundance'),
 ('gas', 'Na_Density'),
 ('gas', 'Na_Fraction'),
 ('gas', 'Na_Mass'),
 ('gas', 'Na_Number_Density'),
 ('gas', 'Na_over_Ba'),
 ('gas', 'Na_over_C'),
 ('gas', 'Na_over_Fe'),
 ('gas', 'Na_over_H'),
 ('gas', 'Na_over_Mg'),
 ('gas', 'Na_over_N'),
 ('gas', 'Na_over_Ni'),
 ('gas', 'Na_over_O'),
 ('gas', 'Ni_Abundance'),
 ('gas', 'Ni_Density'),
 ('gas', 'Ni_Fraction'),
 ('gas', 'Ni_Mass'),
 ('gas', 'Ni_Number_Density'),
 ('gas', 'Ni_over_Ba'),
 ('gas', 'Ni_over_C'),
 ('gas', 'Ni_over_Fe'),
 ('gas', 'Ni_over_H'),
 ('gas', 'Ni_over_Mg'),
 ('gas', 'Ni_over_N'),
 ('gas', 'Ni_over_Na'),
 ('gas', 'Ni_over_O'),
 ('gas', 'OTLW_kdissH2I'),
 ('gas', 'O_Abundance'),
 ('gas', 'O_Density'),
 ('gas', 'O_Fraction'),
 ('gas', 'O_Mass'),
 ('gas', 'O_Number_Density'),
 ('gas', 'O_over_Ba'),
 ('gas', 'O_over_C'),
 ('gas', 'O_over_Fe'),
 ('gas', 'O_over_H'),
 ('gas', 'O_over_Mg'),
 ('gas', 'O_over_N'),
 ('gas', 'O_over_Na'),
 ('gas', 'O_over_Ni'),
 ('gas', 'Pe_heating_rate'),
 ('gas', 'Pe_heating_rate_masked'),
 ('gas', 'Q0_flux'),
 ('gas', 'Q1_flux'),
 ('gas', 'S_Abundance'),
 ('gas', 'S_Density'),
 ('gas', 'S_Fraction'),
 ('gas', 'S_Mass'),
 ('gas', 'S_Number_Density'),
 ('gas', 'S_over_Ba'),
 ('gas', 'S_over_C'),
 ('gas', 'S_over_Fe'),
 ('gas', 'S_over_H'),
 ('gas', 'S_over_Mg'),
 ('gas', 'S_over_N'),
 ('gas', 'S_over_Na'),
 ('gas', 'S_over_Ni'),
 ('gas', 'S_over_O'),
 ('gas', 'Si_Abundance'),
 ('gas', 'Si_Density'),
 ('gas', 'Si_Fraction'),
 ('gas', 'Si_Mass'),
 ('gas', 'Si_Number_Density'),
 ('gas', 'Si_over_Ba'),
 ('gas', 'Si_over_C'),
 ('gas', 'Si_over_Fe'),
 ('gas', 'Si_over_H'),
 ('gas', 'Si_over_Mg'),
 ('gas', 'Si_over_N'),
 ('gas', 'Si_over_Na'),
 ('gas', 'Si_over_Ni'),
 ('gas', 'Si_over_O'),
 ('gas', 'Sr_Abundance'),
 ('gas', 'Sr_Density'),
 ('gas', 'Sr_Fraction'),
 ('gas', 'Sr_Mass'),
 ('gas', 'Sr_Number_Density'),
 ('gas', 'Sr_over_Ba'),
 ('gas', 'Sr_over_C'),
 ('gas', 'Sr_over_Fe'),
 ('gas', 'Sr_over_H'),
 ('gas', 'Sr_over_Mg'),
 ('gas', 'Sr_over_N'),
 ('gas', 'Sr_over_Na'),
 ('gas', 'Sr_over_Ni'),
 ('gas', 'Sr_over_O'),
 ('gas', 'Y_Abundance'),
 ('gas', 'Y_Density'),
 ('gas', 'Y_Fraction'),
 ('gas', 'Y_Mass'),
 ('gas', 'Y_Number_Density'),
 ('gas', 'Y_over_Ba'),
 ('gas', 'Y_over_C'),
 ('gas', 'Y_over_Fe'),
 ('gas', 'Y_over_H'),
 ('gas', 'Y_over_Mg'),
 ('gas', 'Y_over_N'),
 ('gas', 'Y_over_Na'),
 ('gas', 'Y_over_Ni'),
 ('gas', 'Y_over_O'),
 ('gas', 'a_rad'),
 ('gas', 'alpha_Abundance'),
 ('gas', 'alpha_Mass'),
 ('gas', 'alpha_over_As'),
 ('gas', 'alpha_over_Ba'),
 ('gas', 'alpha_over_C'),
 ('gas', 'alpha_over_Ca'),
 ('gas', 'alpha_over_Fe'),
 ('gas', 'alpha_over_H'),
 ('gas', 'alpha_over_Mg'),
 ('gas', 'alpha_over_Mn'),
 ('gas', 'alpha_over_N'),
 ('gas', 'alpha_over_Na'),
 ('gas', 'alpha_over_Ni'),
 ('gas', 'alpha_over_O'),
 ('gas', 'alpha_over_S'),
 ('gas', 'alpha_over_Si'),
 ('gas', 'alpha_over_Sr'),
 ('gas', 'alpha_over_Y'),
 ('gas', 'angular_momentum_magnitude'),
 ('gas', 'angular_momentum_x'),
 ('gas', 'angular_momentum_y'),
 ('gas', 'angular_momentum_z'),
 ('gas', 'averaged_density'),
 ('gas', 'baroclinic_vorticity_magnitude'),
 ('gas', 'baroclinic_vorticity_x'),
 ('gas', 'baroclinic_vorticity_y'),
 ('gas', 'baroclinic_vorticity_z'),
 ('gas', 'baryon_overdensity'),
 ('gas', 'cell_mass'),
 ('gas', 'cell_volume'),
 ('gas', 'chandra_emissivity'),
 ('gas', 'courant_time_step'),
 ('gas', 'cutting_plane_velocity_x'),
 ('gas', 'cutting_plane_velocity_y'),
 ('gas', 'cutting_plane_velocity_z'),
 ('gas', 'cylindrical_radial_velocity'),
 ('gas', 'cylindrical_radial_velocity_absolute'),
 ('gas', 'cylindrical_tangential_velocity'),
 ('gas', 'cylindrical_tangential_velocity_absolute'),
 ('gas', 'dark_matter_density'),
 ('gas', 'density'),
 ('gas', 'density_gradient_magnitude'),
 ('gas', 'density_gradient_x'),
 ('gas', 'density_gradient_y'),
 ('gas', 'density_gradient_z'),
 ('gas', 'dx'),
 ('gas', 'dy'),
 ('gas', 'dynamical_time'),
 ('gas', 'dz'),
 ('gas', 'emission_measure'),
 ('gas', 'entropy'),
 ('gas', 'is_star_forming'),
 ('gas', 'jeans_mass'),
 ('gas', 'kT'),
 ('gas', 'kinetic_energy'),
 ('gas', 'mach_number'),
 ('gas', 'matter_density'),
 ('gas', 'matter_mass'),
 ('gas', 'matter_overdensity'),
 ('gas', 'mazzotta_weighting'),
 ('gas', 'mean_molecular_weight'),
 ('gas', 'metal_density'),
 ('gas', 'metal_mass'),
 ('gas', 'metallicity'),
 ('gas', 'number_density'),
 ('gas', 'overdensity'),
 ('gas', 'path_element_x'),
 ('gas', 'path_element_y'),
 ('gas', 'path_element_z'),
 ('gas', 'photo_gamma'),
 ('gas', 'pressure'),
 ('gas', 'pressure_gradient_magnitude'),
 ('gas', 'pressure_gradient_x'),
 ('gas', 'pressure_gradient_y'),
 ('gas', 'pressure_gradient_z'),
 ('gas', 'radial_mach_number'),
 ('gas', 'radial_velocity'),
 ('gas', 'radial_velocity_absolute'),
 ('gas', 'radiation_acceleration_x'),
 ('gas', 'radiation_acceleration_y'),
 ('gas', 'radiation_acceleration_z'),
 ('gas', 'relative_velocity_x'),
 ('gas', 'relative_velocity_y'),
 ('gas', 'relative_velocity_z'),
 ('gas', 'shear'),
 ('gas', 'shear_criterion'),
 ('gas', 'shear_mach'),
 ('gas', 'sound_speed'),
 ('gas', 'specific_angular_momentum_magnitude'),
 ('gas', 'specific_angular_momentum_x'),
 ('gas', 'specific_angular_momentum_y'),
 ('gas', 'specific_angular_momentum_z'),
 ('gas', 'sz_kinetic'),
 ('gas', 'szy'),
 ('gas', 'tangential_over_velocity_magnitude'),
 ('gas', 'tangential_velocity'),
 ('gas', 'temperature'),
 ('gas', 'thermal_energy'),
 ('gas', 'total_energy'),
 ('gas', 'velocity_cylindrical_radius'),
 ('gas', 'velocity_cylindrical_theta'),
 ('gas', 'velocity_cylindrical_z'),
 ('gas', 'velocity_divergence'),
 ('gas', 'velocity_divergence_absolute'),
 ('gas', 'velocity_magnitude'),
 ('gas', 'velocity_spherical_phi'),
 ('gas', 'velocity_spherical_radius'),
 ('gas', 'velocity_spherical_theta'),
 ('gas', 'velocity_x'),
 ('gas', 'velocity_y'),
 ('gas', 'velocity_z'),
 ('gas', 'vertex_x'),
 ('gas', 'vertex_y'),
 ('gas', 'vertex_z'),
 ('gas', 'vorticity_growth_magnitude'),
 ('gas', 'vorticity_growth_magnitude_absolute'),
 ('gas', 'vorticity_growth_timescale'),
 ('gas', 'vorticity_growth_x'),
 ('gas', 'vorticity_growth_y'),
 ('gas', 'vorticity_growth_z'),
 ('gas', 'vorticity_magnitude'),
 ('gas', 'vorticity_radiation_pressure_growth_magnitude'),
 ('gas', 'vorticity_radiation_pressure_growth_magnitude_absolute'),
 ('gas', 'vorticity_radiation_pressure_growth_timescale'),
 ('gas', 'vorticity_radiation_pressure_growth_x'),
 ('gas', 'vorticity_radiation_pressure_growth_y'),
 ('gas', 'vorticity_radiation_pressure_growth_z'),
 ('gas', 'vorticity_radiation_pressure_magnitude'),
 ('gas', 'vorticity_radiation_pressure_x'),
 ('gas', 'vorticity_radiation_pressure_y'),
 ('gas', 'vorticity_radiation_pressure_z'),
 ('gas', 'vorticity_squared'),
 ('gas', 'vorticity_stretching_magnitude'),
 ('gas', 'vorticity_stretching_x'),
 ('gas', 'vorticity_stretching_y'),
 ('gas', 'vorticity_stretching_z'),
 ('gas', 'vorticity_x'),
 ('gas', 'vorticity_y'),
 ('gas', 'vorticity_z'),
 ('gas', 'weak_lensing_convergence'),
 ('gas', 'x'),
 ('gas', 'xray_emissivity'),
 ('gas', 'y'),
 ('gas', 'z'),
 ('index', 'DM_background_density'),
 ('index', 'DM_background_potential'),
 ('index', 'cell_volume'),
 ('index', 'cylindrical_r'),
 ('index', 'cylindrical_radius'),
 ('index', 'cylindrical_theta'),
 ('index', 'cylindrical_z'),
 ('index', 'disk_angle'),
 ('index', 'dx'),
 ('index', 'dy'),
 ('index', 'dz'),
 ('index', 'grid_indices'),
 ('index', 'grid_level'),
 ('index', 'height'),
 ('index', 'magnitude_cylindrical_radius'),
 ('index', 'magnitude_cylindrical_z'),
 ('index', 'morton_index'),
 ('index', 'ones'),
 ('index', 'ones_over_dx'),
 ('index', 'path_element_x'),
 ('index', 'path_element_y'),
 ('index', 'path_element_z'),
 ('index', 'radius'),
 ('index', 'spherical_phi'),
 ('index', 'spherical_r'),
 ('index', 'spherical_radius'),
 ('index', 'spherical_theta'),
 ('index', 'vertex_x'),
 ('index', 'vertex_y'),
 ('index', 'vertex_z'),
 ('index', 'virial_radius_fraction'),
 ('index', 'x'),
 ('index', 'y'),
 ('index', 'z'),
 ('index', 'zeros'),
 ('io', 'age'),
 ('io', 'birth_mass'),
 ('io', 'creation_time'),
 ('io', 'dynamical_time'),
 ('io', 'mesh_id'),
 ('io', 'metallicity_fraction'),
 ('io', 'particle_age'),
 ('io', 'particle_angular_momentum'),
 ('io', 'particle_angular_momentum_magnitude'),
 ('io', 'particle_angular_momentum_x'),
 ('io', 'particle_angular_momentum_y'),
 ('io', 'particle_angular_momentum_z'),
 ('io', 'particle_cylindrical_velocity_theta'),
 ('io', 'particle_cylindrical_velocity_z'),
 ('io', 'particle_index'),
 ('io', 'particle_mass'),
 ('io', 'particle_model_E0'),
 ('io', 'particle_model_E1'),
 ('io', 'particle_model_L0'),
 ('io', 'particle_model_L1'),
 ('io', 'particle_model_L_1_3eV'),
 ('io', 'particle_model_L_FUV'),
 ('io', 'particle_model_L_LW'),
 ('io', 'particle_model_Q0'),
 ('io', 'particle_model_Q1'),
 ('io', 'particle_model_R'),
 ('io', 'particle_model_Teff'),
 ('io', 'particle_model_lifetime'),
 ('io', 'particle_model_luminosity'),
 ('io', 'particle_ones'),
 ('io', 'particle_position'),
 ('io', 'particle_position_cylindrical_radius'),
 ('io', 'particle_position_cylindrical_theta'),
 ('io', 'particle_position_cylindrical_z'),
 ('io', 'particle_position_relative'),
 ('io', 'particle_position_relative_x'),
 ('io', 'particle_position_relative_y'),
 ('io', 'particle_position_relative_z'),
 ('io', 'particle_position_spherical_phi'),
 ('io', 'particle_position_spherical_radius'),
 ('io', 'particle_position_spherical_theta'),
 ('io', 'particle_position_x'),
 ('io', 'particle_position_y'),
 ('io', 'particle_position_z'),
 ('io', 'particle_radial_velocity'),
 ('io', 'particle_radius'),
 ('io', 'particle_specific_angular_momentum'),
 ('io', 'particle_specific_angular_momentum_x'),
 ('io', 'particle_specific_angular_momentum_y'),
 ('io', 'particle_specific_angular_momentum_z'),
 ('io', 'particle_spherical_position_phi'),
 ('io', 'particle_spherical_position_radius'),
 ('io', 'particle_spherical_position_theta'),
 ('io', 'particle_spherical_velocity_phi'),
 ('io', 'particle_spherical_velocity_radius'),
 ('io', 'particle_spherical_velocity_theta'),
 ('io', 'particle_type'),
 ('io', 'particle_velocity'),
 ('io', 'particle_velocity_cylindrical_radius'),
 ('io', 'particle_velocity_cylindrical_theta'),
 ('io', 'particle_velocity_cylindrical_z'),
 ('io', 'particle_velocity_magnitude'),
 ('io', 'particle_velocity_relative'),
 ('io', 'particle_velocity_relative_x'),
 ('io', 'particle_velocity_relative_y'),
 ('io', 'particle_velocity_relative_z'),
 ('io', 'particle_velocity_spherical_phi'),
 ('io', 'particle_velocity_spherical_radius'),
 ('io', 'particle_velocity_spherical_theta'),
 ('io', 'particle_velocity_x'),
 ('io', 'particle_velocity_y'),
 ('io', 'particle_velocity_z'),
 ('io', 'relative_particle_position'),
 ('io', 'relative_particle_position_x'),
 ('io', 'relative_particle_position_y'),
 ('io', 'relative_particle_position_z'),
 ('io', 'relative_particle_velocity'),
 ('io', 'relative_particle_velocity_x'),
 ('io', 'relative_particle_velocity_y'),
 ('io', 'relative_particle_velocity_z'),
 ('io', 'sn_mass_ejected'),
 ('io', 'wind_mass_ejected'),
 ('low_mass_stars', 'age'),
 ('low_mass_stars', 'birth_mass'),
 ('low_mass_stars', 'creation_time'),
 ('low_mass_stars', 'dynamical_time'),
 ('low_mass_stars', 'mesh_id'),
 ('low_mass_stars', 'metallicity_fraction'),
 ('low_mass_stars', 'particle_angular_momentum'),
 ('low_mass_stars', 'particle_angular_momentum_magnitude'),
 ('low_mass_stars', 'particle_angular_momentum_x'),
 ('low_mass_stars', 'particle_angular_momentum_y'),
 ('low_mass_stars', 'particle_angular_momentum_z'),
 ('low_mass_stars', 'particle_cylindrical_velocity_theta'),
 ('low_mass_stars', 'particle_cylindrical_velocity_z'),
 ('low_mass_stars', 'particle_index'),
 ('low_mass_stars', 'particle_mass'),
 ('low_mass_stars', 'particle_ones'),
 ('low_mass_stars', 'particle_position'),
 ('low_mass_stars', 'particle_position_cylindrical_radius'),
 ('low_mass_stars', 'particle_position_cylindrical_theta'),
 ('low_mass_stars', 'particle_position_cylindrical_z'),
 ('low_mass_stars', 'particle_position_relative'),
 ('low_mass_stars', 'particle_position_relative_x'),
 ('low_mass_stars', 'particle_position_relative_y'),
 ('low_mass_stars', 'particle_position_relative_z'),
 ('low_mass_stars', 'particle_position_spherical_phi'),
 ('low_mass_stars', 'particle_position_spherical_radius'),
 ('low_mass_stars', 'particle_position_spherical_theta'),
 ('low_mass_stars', 'particle_position_x'),
 ('low_mass_stars', 'particle_position_y'),
 ('low_mass_stars', 'particle_position_z'),
 ('low_mass_stars', 'particle_radial_velocity'),
 ('low_mass_stars', 'particle_radius'),
 ('low_mass_stars', 'particle_specific_angular_momentum'),
 ('low_mass_stars', 'particle_specific_angular_momentum_x'),
 ('low_mass_stars', 'particle_specific_angular_momentum_y'),
 ('low_mass_stars', 'particle_specific_angular_momentum_z'),
 ('low_mass_stars', 'particle_spherical_position_phi'),
 ('low_mass_stars', 'particle_spherical_position_radius'),
 ('low_mass_stars', 'particle_spherical_position_theta'),
 ('low_mass_stars', 'particle_spherical_velocity_phi'),
 ('low_mass_stars', 'particle_spherical_velocity_radius'),
 ('low_mass_stars', 'particle_spherical_velocity_theta'),
 ('low_mass_stars', 'particle_type'),
 ('low_mass_stars', 'particle_velocity'),
 ('low_mass_stars', 'particle_velocity_cylindrical_radius'),
 ('low_mass_stars', 'particle_velocity_cylindrical_theta'),
 ('low_mass_stars', 'particle_velocity_cylindrical_z'),
 ('low_mass_stars', 'particle_velocity_magnitude'),
 ('low_mass_stars', 'particle_velocity_relative'),
 ('low_mass_stars', 'particle_velocity_relative_x'),
 ('low_mass_stars', 'particle_velocity_relative_y'),
 ('low_mass_stars', 'particle_velocity_relative_z'),
 ('low_mass_stars', 'particle_velocity_spherical_phi'),
 ('low_mass_stars', 'particle_velocity_spherical_radius'),
 ('low_mass_stars', 'particle_velocity_spherical_theta'),
 ('low_mass_stars', 'particle_velocity_x'),
 ('low_mass_stars', 'particle_velocity_y'),
 ('low_mass_stars', 'particle_velocity_z'),
 ('low_mass_stars', 'relative_particle_position'),
 ('low_mass_stars', 'relative_particle_position_x'),
 ('low_mass_stars', 'relative_particle_position_y'),
 ('low_mass_stars', 'relative_particle_position_z'),
 ('low_mass_stars', 'relative_particle_velocity'),
 ('low_mass_stars', 'relative_particle_velocity_x'),
 ('low_mass_stars', 'relative_particle_velocity_y'),
 ('low_mass_stars', 'relative_particle_velocity_z'),
 ('low_mass_stars', 'sn_mass_ejected'),
 ('low_mass_stars', 'wind_mass_ejected'),
 ('low_mass_unresolved_stars', 'age'),
 ('low_mass_unresolved_stars', 'birth_mass'),
 ('low_mass_unresolved_stars', 'creation_time'),
 ('low_mass_unresolved_stars', 'dynamical_time'),
 ('low_mass_unresolved_stars', 'mesh_id'),
 ('low_mass_unresolved_stars', 'metallicity_fraction'),
 ('low_mass_unresolved_stars', 'particle_angular_momentum'),
 ('low_mass_unresolved_stars', 'particle_angular_momentum_magnitude'),
 ('low_mass_unresolved_stars', 'particle_angular_momentum_x'),
 ('low_mass_unresolved_stars', 'particle_angular_momentum_y'),
 ('low_mass_unresolved_stars', 'particle_angular_momentum_z'),
 ('low_mass_unresolved_stars', 'particle_cylindrical_velocity_theta'),
 ('low_mass_unresolved_stars', 'particle_cylindrical_velocity_z'),
 ('low_mass_unresolved_stars', 'particle_index'),
 ('low_mass_unresolved_stars', 'particle_mass'),
 ('low_mass_unresolved_stars', 'particle_ones'),
 ('low_mass_unresolved_stars', 'particle_position'),
 ('low_mass_unresolved_stars', 'particle_position_cylindrical_radius'),
 ('low_mass_unresolved_stars', 'particle_position_cylindrical_theta'),
 ('low_mass_unresolved_stars', 'particle_position_cylindrical_z'),
 ('low_mass_unresolved_stars', 'particle_position_relative'),
 ('low_mass_unresolved_stars', 'particle_position_relative_x'),
 ('low_mass_unresolved_stars', 'particle_position_relative_y'),
 ('low_mass_unresolved_stars', 'particle_position_relative_z'),
 ('low_mass_unresolved_stars', 'particle_position_spherical_phi'),
 ('low_mass_unresolved_stars', 'particle_position_spherical_radius'),
 ('low_mass_unresolved_stars', 'particle_position_spherical_theta'),
 ('low_mass_unresolved_stars', 'particle_position_x'),
 ('low_mass_unresolved_stars', 'particle_position_y'),
 ('low_mass_unresolved_stars', 'particle_position_z'),
 ('low_mass_unresolved_stars', 'particle_radial_velocity'),
 ('low_mass_unresolved_stars', 'particle_radius'),
 ('low_mass_unresolved_stars', 'particle_specific_angular_momentum'),
 ('low_mass_unresolved_stars', 'particle_specific_angular_momentum_x'),
 ('low_mass_unresolved_stars', 'particle_specific_angular_momentum_y'),
 ('low_mass_unresolved_stars', 'particle_specific_angular_momentum_z'),
 ('low_mass_unresolved_stars', 'particle_spherical_position_phi'),
 ('low_mass_unresolved_stars', 'particle_spherical_position_radius'),
 ('low_mass_unresolved_stars', 'particle_spherical_position_theta'),
 ('low_mass_unresolved_stars', 'particle_spherical_velocity_phi'),
 ('low_mass_unresolved_stars', 'particle_spherical_velocity_radius'),
 ('low_mass_unresolved_stars', 'particle_spherical_velocity_theta'),
 ('low_mass_unresolved_stars', 'particle_type'),
 ('low_mass_unresolved_stars', 'particle_velocity'),
 ('low_mass_unresolved_stars', 'particle_velocity_cylindrical_radius'),
 ('low_mass_unresolved_stars', 'particle_velocity_cylindrical_theta'),
 ('low_mass_unresolved_stars', 'particle_velocity_cylindrical_z'),
 ('low_mass_unresolved_stars', 'particle_velocity_magnitude'),
 ('low_mass_unresolved_stars', 'particle_velocity_relative'),
 ('low_mass_unresolved_stars', 'particle_velocity_relative_x'),
 ('low_mass_unresolved_stars', 'particle_velocity_relative_y'),
 ('low_mass_unresolved_stars', 'particle_velocity_relative_z'),
 ('low_mass_unresolved_stars', 'particle_velocity_spherical_phi'),
 ('low_mass_unresolved_stars', 'particle_velocity_spherical_radius'),
 ('low_mass_unresolved_stars', 'particle_velocity_spherical_theta'),
 ('low_mass_unresolved_stars', 'particle_velocity_x'),
 ('low_mass_unresolved_stars', 'particle_velocity_y'),
 ('low_mass_unresolved_stars', 'particle_velocity_z'),
 ('low_mass_unresolved_stars', 'relative_particle_position'),
 ('low_mass_unresolved_stars', 'relative_particle_position_x'),
 ('low_mass_unresolved_stars', 'relative_particle_position_y'),
 ('low_mass_unresolved_stars', 'relative_particle_position_z'),
 ('low_mass_unresolved_stars', 'relative_particle_velocity'),
 ('low_mass_unresolved_stars', 'relative_particle_velocity_x'),
 ('low_mass_unresolved_stars', 'relative_particle_velocity_y'),
 ('low_mass_unresolved_stars', 'relative_particle_velocity_z'),
 ('low_mass_unresolved_stars', 'sn_mass_ejected'),
 ('low_mass_unresolved_stars', 'wind_mass_ejected'),
 ('main_sequence_popIII_stars', 'age'),
 ('main_sequence_popIII_stars', 'birth_mass'),
 ('main_sequence_popIII_stars', 'creation_time'),
 ('main_sequence_popIII_stars', 'dynamical_time'),
 ('main_sequence_popIII_stars', 'mesh_id'),
 ('main_sequence_popIII_stars', 'metallicity_fraction'),
 ('main_sequence_popIII_stars', 'particle_angular_momentum'),
 ('main_sequence_popIII_stars', 'particle_angular_momentum_magnitude'),
 ('main_sequence_popIII_stars', 'particle_angular_momentum_x'),
 ('main_sequence_popIII_stars', 'particle_angular_momentum_y'),
 ('main_sequence_popIII_stars', 'particle_angular_momentum_z'),
 ('main_sequence_popIII_stars', 'particle_cylindrical_velocity_theta'),
 ('main_sequence_popIII_stars', 'particle_cylindrical_velocity_z'),
 ('main_sequence_popIII_stars', 'particle_index'),
 ('main_sequence_popIII_stars', 'particle_mass'),
 ('main_sequence_popIII_stars', 'particle_ones'),
 ('main_sequence_popIII_stars', 'particle_position'),
 ('main_sequence_popIII_stars', 'particle_position_cylindrical_radius'),
 ('main_sequence_popIII_stars', 'particle_position_cylindrical_theta'),
 ('main_sequence_popIII_stars', 'particle_position_cylindrical_z'),
 ('main_sequence_popIII_stars', 'particle_position_relative'),
 ('main_sequence_popIII_stars', 'particle_position_relative_x'),
 ('main_sequence_popIII_stars', 'particle_position_relative_y'),
 ('main_sequence_popIII_stars', 'particle_position_relative_z'),
 ('main_sequence_popIII_stars', 'particle_position_spherical_phi'),
 ('main_sequence_popIII_stars', 'particle_position_spherical_radius'),
 ('main_sequence_popIII_stars', 'particle_position_spherical_theta'),
 ('main_sequence_popIII_stars', 'particle_position_x'),
 ('main_sequence_popIII_stars', 'particle_position_y'),
 ('main_sequence_popIII_stars', 'particle_position_z'),
 ('main_sequence_popIII_stars', 'particle_radial_velocity'),
 ('main_sequence_popIII_stars', 'particle_radius'),
 ('main_sequence_popIII_stars', 'particle_specific_angular_momentum'),
 ('main_sequence_popIII_stars', 'particle_specific_angular_momentum_x'),
 ('main_sequence_popIII_stars', 'particle_specific_angular_momentum_y'),
 ('main_sequence_popIII_stars', 'particle_specific_angular_momentum_z'),
 ('main_sequence_popIII_stars', 'particle_spherical_position_phi'),
 ('main_sequence_popIII_stars', 'particle_spherical_position_radius'),
 ('main_sequence_popIII_stars', 'particle_spherical_position_theta'),
 ('main_sequence_popIII_stars', 'particle_spherical_velocity_phi'),
 ('main_sequence_popIII_stars', 'particle_spherical_velocity_radius'),
 ('main_sequence_popIII_stars', 'particle_spherical_velocity_theta'),
 ('main_sequence_popIII_stars', 'particle_type'),
 ('main_sequence_popIII_stars', 'particle_velocity'),
 ('main_sequence_popIII_stars', 'particle_velocity_cylindrical_radius'),
 ('main_sequence_popIII_stars', 'particle_velocity_cylindrical_theta'),
 ('main_sequence_popIII_stars', 'particle_velocity_cylindrical_z'),
 ('main_sequence_popIII_stars', 'particle_velocity_magnitude'),
 ('main_sequence_popIII_stars', 'particle_velocity_relative'),
 ('main_sequence_popIII_stars', 'particle_velocity_relative_x'),
 ('main_sequence_popIII_stars', 'particle_velocity_relative_y'),
 ('main_sequence_popIII_stars', 'particle_velocity_relative_z'),
 ('main_sequence_popIII_stars', 'particle_velocity_spherical_phi'),
 ('main_sequence_popIII_stars', 'particle_velocity_spherical_radius'),
 ('main_sequence_popIII_stars', 'particle_velocity_spherical_theta'),
 ('main_sequence_popIII_stars', 'particle_velocity_x'),
 ('main_sequence_popIII_stars', 'particle_velocity_y'),
 ('main_sequence_popIII_stars', 'particle_velocity_z'),
 ('main_sequence_popIII_stars', 'relative_particle_position'),
 ('main_sequence_popIII_stars', 'relative_particle_position_x'),
 ('main_sequence_popIII_stars', 'relative_particle_position_y'),
 ('main_sequence_popIII_stars', 'relative_particle_position_z'),
 ('main_sequence_popIII_stars', 'relative_particle_velocity'),
 ('main_sequence_popIII_stars', 'relative_particle_velocity_x'),
 ('main_sequence_popIII_stars', 'relative_particle_velocity_y'),
 ('main_sequence_popIII_stars', 'relative_particle_velocity_z'),
 ('main_sequence_popIII_stars', 'sn_mass_ejected'),
 ('main_sequence_popIII_stars', 'wind_mass_ejected'),
 ('main_sequence_stars', 'age'),
 ('main_sequence_stars', 'birth_mass'),
 ('main_sequence_stars', 'creation_time'),
 ('main_sequence_stars', 'dynamical_time'),
 ('main_sequence_stars', 'mesh_id'),
 ('main_sequence_stars', 'metallicity_fraction'),
 ('main_sequence_stars', 'particle_angular_momentum'),
 ('main_sequence_stars', 'particle_angular_momentum_magnitude'),
 ('main_sequence_stars', 'particle_angular_momentum_x'),
 ('main_sequence_stars', 'particle_angular_momentum_y'),
 ('main_sequence_stars', 'particle_angular_momentum_z'),
 ('main_sequence_stars', 'particle_cylindrical_velocity_theta'),
 ('main_sequence_stars', 'particle_cylindrical_velocity_z'),
 ('main_sequence_stars', 'particle_index'),
 ('main_sequence_stars', 'particle_mass'),
 ('main_sequence_stars', 'particle_ones'),
 ('main_sequence_stars', 'particle_position'),
 ('main_sequence_stars', 'particle_position_cylindrical_radius'),
 ('main_sequence_stars', 'particle_position_cylindrical_theta'),
 ('main_sequence_stars', 'particle_position_cylindrical_z'),
 ('main_sequence_stars', 'particle_position_relative'),
 ('main_sequence_stars', 'particle_position_relative_x'),
 ('main_sequence_stars', 'particle_position_relative_y'),
 ('main_sequence_stars', 'particle_position_relative_z'),
 ('main_sequence_stars', 'particle_position_spherical_phi'),
 ('main_sequence_stars', 'particle_position_spherical_radius'),
 ('main_sequence_stars', 'particle_position_spherical_theta'),
 ('main_sequence_stars', 'particle_position_x'),
 ('main_sequence_stars', 'particle_position_y'),
 ('main_sequence_stars', 'particle_position_z'),
 ('main_sequence_stars', 'particle_radial_velocity'),
 ('main_sequence_stars', 'particle_radius'),
 ('main_sequence_stars', 'particle_specific_angular_momentum'),
 ('main_sequence_stars', 'particle_specific_angular_momentum_x'),
 ('main_sequence_stars', 'particle_specific_angular_momentum_y'),
 ('main_sequence_stars', 'particle_specific_angular_momentum_z'),
 ('main_sequence_stars', 'particle_spherical_position_phi'),
 ...]

In [ ]:
#if np.size(px)==0:: # main sequnce stars
select = data['particle_type'] == 11
px = data['particle_position_x'][select]
py = data['particle_position_y'][select]
pz = data['particle_position_z'][select]

print("number of main sequnce stars ", np.size(px))

if np.size(px)==0: # unresolved stars
    select = data['particle_type'] == 14
    px = data['particle_position_x'][select]
    py = data['particle_position_y'][select]
    pz = data['particle_position_z'][select]
    
    if np.size(px) > 0:
        print("using MS popIII stars")        
else:
    print("using main sequnce stars")
    
print("number of MS popIII stars ", np.size(data['birth_mass'][ data['particle_type'] == 14]))


if np.size(px)==0: # unresolved stars
    select = data['particle_type'] == 15
    px = data['particle_position_x'][select]
    py = data['particle_position_y'][select]
    pz = data['particle_position_z'][select]
    
    if np.size(px) > 0:
        print("using unresolved stars")
        

if np.size(px)==0: # remnants
    select = data['particle_type'] == 13
    px = data['particle_position_x'][select]
    py = data['particle_position_y'][select]
    pz = data['particle_position_z'][select]
    
    if np.size(px) > 0:
    
        print("using remnants")

if np.size(px) == 0: # else, max density
    center = 'max'
    print("using max density")
else:
    center = [px[0],py[0],pz[0]]    


center = ds.domain_center.value



#
#
#
select = (data['birth_mass'] > 3.0) * (data['birth_mass'] < 8.0) * (data['particle_type'] == 13)

print("Number of AGB events: ",np.size(data['birth_mass'][select]))
select = (data['birth_mass'] > 2.0) * (data['birth_mass'] < 8.0) * (data['particle_type'] == 11)

print("Number of stars that will have AGB events : ",np.size(data['birth_mass'][select]))
t_life = data['dynamical_time'].to('Myr')[select]
t_birth = data['creation_time'].to('Myr')[select]

death = t_life + t_birth

time_to_death = (death - ds.current_time)
print("First will occur in : ", np.min(time_to_death), np.min(time_to_death[time_to_death>0]),"Myr")

In [18]:
l_kpc = 15.0 * yt.units.kpc

ds.domain_width[0].to('kpc') * 0.5 # 10.0*yt.units.kpc
l = (l_kpc / ds.length_unit.to('kpc')).value

box = ds.region( center, center - l, center + l)


all_fields = ['number_density','Temperature','Fe_over_H','alpha_over_Fe',
              #,'FUV_FluxDensity','HI_kph',
              'PopIII_Metal_Density']#,'H2I_kdiss','Ba_over_Fe']

sp = yt.ProjectionPlot(ds, 'z', all_fields, center = center,
                       data_source = box, weight_field = ('gas','metal_density'))
sp.set_width(l_kpc)
sp.annotate_particles(0.95, ptype = 'main_sequence_stars', marker='*', col='C1', p_size=100)
sp.annotate_particles(0.95, ptype = 'main_sequence_popIII_stars', marker='*', col='black', p_size=100)
sp.annotate_particles(0.95, ptype = 'remnant_stars', marker='D', col='black', p_size = 10)
sp.annotate_particles(0.95, ptype = 'low_mass_unresolved_stars', marker='o', col = 'black', p_size=10)



if 'number_density' in sp.fields:
    field = 'number_density'
    sp.set_zlim(field, 1.0E-4, 1.0E1)
    sp.set_cmap(field,'viridis')
    
if 'Temperature' in sp.fields:
    field = 'Temperature'
    sp.set_zlim(field, 100.0, 1.0E5)
    sp.set_cmap(field,'RdYlBu_r')    
    
if 'Fe_over_H' in sp.fields:
    field = 'Fe_over_H'
    sp.set_log(field,False)
    sp.set_zlim(field, -10, -1)
    sp.set_cmap(field,'cubehelix') 
    
if 'alpha_over_Fe' in sp.fields:
    field = 'alpha_over_Fe'
    sp.set_log(field,False)
    sp.set_zlim(field, -1, 2)
    sp.set_cmap(field,'cubehelix') 

        
if 'Ba_over_Fe' in sp.fields:
    field = 'Ba_over_Fe'
    sp.set_log(field,False)
    sp.set_zlim(field, -3, 3)
    sp.set_cmap(field,'cubehelix')

    
sp.show()
sp.set_buff_size(1024)
sp.save()


/home/aemerick/.local/lib/python3.7/site-packages/yt/data_objects/construction_data_containers.py:379: RuntimeWarning: divide by zero encountered in true_divide
  np.divide(nvals, nwvals[:,None], nvals)





Out[18]:
['DD0088_Projection_z_Temperature_metal_density.png',
 'DD0088_Projection_z_Fe_over_H_metal_density.png',
 'DD0088_Projection_z_alpha_over_Fe_metal_density.png',
 'DD0088_Projection_z_PopIII_Metal_Density_metal_density.png',
 'DD0088_Projection_z_number_density_metal_density.png']

In [15]:
np.max(data['HI_kph'].to('1/s'))


Out[15]:
3.9943893670630954e-11 1/s

In [81]:
#
#
#
#
#
pp = yt.PhasePlot( ds, 'number_density', 'Temperature','cell_mass')
pp.set_xlim(1.0E-6,1.0E3)
pp.set_ylim(1.0,1.0E7)
pp.set_unit('cell_mass','Msun')
pp.set_cmap('cell_mass','viridis')

pp.show()




In [59]:



Out[59]:
42.40710538525811 dimensionless

In [22]:
#
#
#
#
#
reg = ds.region(ds.domain_center, ds.domain_center - 20.0*yt.units.kpc, ds.domain_center+20.0*yt.units.kpc)
pp = yt.PhasePlot( reg, 'number_density', 'O_over_Fe','cell_mass')#, ybins
pp.set_log("O_over_Fe",False)
pp.set_ylim(-3,3)

pp.set_xlim(1.0E-6,1.0E3)
pp.set_unit('cell_mass','Msun')
pp.set_cmap('cell_mass','viridis')

pp.show()




In [30]:
#
#
#
#
#
reg = ds.region(ds.domain_center, ds.domain_center - 20.0*yt.units.kpc, ds.domain_center+20.0*yt.units.kpc)
pp = yt.PhasePlot( reg, 'Fe_over_H','O_over_Fe', 'cell_mass')

pp.set_ylim(1,3)
pp.set_log("O_over_Fe",False)    
pp.set_ylim(-3,3)

pp.set_xlim(0.1,1)
pp.set_log("Fe_over_H",False)
pp.set_xlim(-12,1)


pp.set_unit('cell_mass','Msun')
pp.set_cmap('cell_mass','viridis')

pp.show()


---------------------------------------------------------------------------
YTIllDefinedBounds                        Traceback (most recent call last)
<ipython-input-30-66604e55ec00> in <module>
      8 
      9 pp.set_ylim(1,3)
---> 10 pp.set_log("O_over_Fe",False)
     11 pp.set_ylim(-3,3)
     12 

~/.local/lib/python3.7/site-packages/yt/visualization/plot_container.py in newfunc(*args, **kwargs)
     73     @wraps(f)
     74     def newfunc(*args, **kwargs):
---> 75         rv = f(*args, **kwargs)
     76         args[0]._plot_valid = False
     77         return rv

~/.local/lib/python3.7/site-packages/yt/visualization/profile_plotter.py in set_log(self, field, log)
   1402             self._profile_valid = False
   1403         else:
-> 1404             field, = self.profile.data_source._determine_fields([field])
   1405             if field == p.x_field:
   1406                 self.x_log = log

~/.local/lib/python3.7/site-packages/yt/visualization/profile_plotter.py in profile(self)
   1017     def profile(self):
   1018         if not self._profile_valid:
-> 1019             self._recreate_profile()
   1020         return self._profile
   1021 

~/.local/lib/python3.7/site-packages/yt/visualization/profile_plotter.py in _recreate_profile(self)
   1542             extrema=extrema,
   1543             logs=logs,
-> 1544             **additional_kwargs)
   1545         for field in zunits:
   1546             self._profile.set_field_unit(field, zunits[field])

~/.local/lib/python3.7/site-packages/yt/data_objects/profiles.py in create_profile(data_source, bin_fields, fields, n_bins, extrema, logs, units, weight_field, accumulation, fractional, deposition, override_bins)
   1193     for f, n, (mi, ma), l in zip(bin_fields, n_bins, ex, logs):
   1194         if mi <= 0 and l:
-> 1195             raise YTIllDefinedBounds(mi, ma)
   1196         args += [f, n, mi, ma, l]
   1197     kwargs = dict(weight_field=weight_field)

YTIllDefinedBounds: The bounds -9.414e+00 and 5.496e-01 are ill-defined. Typically this happens when a log binning is specified and zero or negative values are given for the bounds.

In [28]:
yt.PhasePlot?

In [19]:
#
#
#
#
#
reg = ds.region(ds.domain_center, ds.domain_center - 20.0*yt.units.kpc, ds.domain_center+20.0*yt.units.kpc)
pp = yt.PhasePlot( reg, 'number_density', 'Fe_over_H','cell_mass')

                  
pp.set_log("Fe_over_H",False)
pp.set_ylim(-12,1)

pp.set_xlim(1.0E-6,1.0E3)
pp.set_unit('cell_mass','Msun')
pp.set_cmap('cell_mass','viridis')

pp.show()




In [28]:
#
#
#

x  =   data[('io','particle_Fe_over_H')]

bins, hist = np.histogram(x, bins = np.arange(-10,1,0.25))

fig, ax = plt.subplots()
fig.set_size_inches(6,6)

plt.step(bins, hist, where = 'post')


---------------------------------------------------------------------------
YTFieldNotFound                           Traceback (most recent call last)
<ipython-input-28-4aebe1463991> in <module>
      3 #
      4 
----> 5 x  =   data[('io','particle_Fe_over_H')]
      6 
      7 bins, hist = np.histogram(x, bins = np.arange(-10,1,0.25))

~/.local/lib/python3.7/site-packages/yt/data_objects/data_containers.py in __getitem__(self, key)
    247         Returns a single field.  Will add if necessary.
    248         """
--> 249         f = self._determine_fields([key])[0]
    250         if f not in self.field_data and key not in self.field_data:
    251             if f in self._container_fields:

~/.local/lib/python3.7/site-packages/yt/data_objects/data_containers.py in _determine_fields(self, fields)
   1333                     raise YTFieldNotParseable(field)
   1334                 ftype, fname = field
-> 1335                 finfo = self.ds._get_field_info(ftype, fname)
   1336             elif isinstance(field, DerivedField):
   1337                 ftype, fname = field.name

~/.local/lib/python3.7/site-packages/yt/data_objects/static_output.py in _get_field_info(self, ftype, fname)
    737                     self._last_finfo = self.field_info[(ftype, fname)]
    738                     return self._last_finfo
--> 739         raise YTFieldNotFound((ftype, fname), self)
    740 
    741     def _setup_classes(self):

YTFieldNotFound: Could not find field '('io', 'particle_Fe_over_H')' in DD0037.

In [ ]:
f= data['PopIII_Metal_Density'] / data['Density'].value
print(np.min(f.value), np.max(f.value), np.sum(f*data['cell_mass']).value/np.sum(data['cell_mass']).value)
f= data['AGB_Metal_Density'] / data['Density'].value
print(np.min(f.value), np.max(f.value), np.sum(f*data['cell_mass']).value/np.sum(data['cell_mass']).value)
f= data['SNII_Metal_Density'] / data['Density'].value
print(np.min(f.value), np.max(f.value), np.sum(f*data['cell_mass']).value/np.sum(data['cell_mass']).value)
f= data['SNIa_Metal_Density'] / data['Density'].value
print(np.min(f.value), np.max(f.value), np.sum(f*data['cell_mass']).value/np.sum(data['cell_mass']).value)

In [88]:
select = data['particle_type']>=11
print(np.sum(data['birth_mass'][select]))

print(np.sum(data['birth_mass'][data['birth_mass']<100]))


15547.292579297628 dimensionless
-88085809382.84601 dimensionless

In [ ]: