In [1]:
%matplotlib inline

In [2]:
import numpy as np
import matplotlib.pyplot as plt
import yt

In [3]:
ts = yt.load("/srv/nbgrader/data/enzo_cosmology_plus/*/*.hierarchy")

In [4]:
len(ts)


Out[4]:
57

In [10]:
ts[0].current_time.in_units("Myr")


yt : [INFO     ] 2017-04-14 16:41:47,595 Parameters: current_time              = 0.81651036015641
yt : [INFO     ] 2017-04-14 16:41:47,597 Parameters: domain_dimensions         = [32 32 32]
yt : [INFO     ] 2017-04-14 16:41:47,600 Parameters: domain_left_edge          = [ 0.  0.  0.]
yt : [INFO     ] 2017-04-14 16:41:47,603 Parameters: domain_right_edge         = [ 1.  1.  1.]
yt : [INFO     ] 2017-04-14 16:41:47,606 Parameters: cosmological_simulation   = 1
yt : [INFO     ] 2017-04-14 16:41:47,607 Parameters: current_redshift          = 50
yt : [INFO     ] 2017-04-14 16:41:47,609 Parameters: omega_lambda              = 0.732
yt : [INFO     ] 2017-04-14 16:41:47,611 Parameters: omega_matter              = 0.268
yt : [INFO     ] 2017-04-14 16:41:47,613 Parameters: hubble_constant           = 0.704
Out[10]:
49.1205196929592 Myr

In [11]:
ts[10].current_time.in_units("Myr")


yt : [INFO     ] 2017-04-14 16:41:55,763 Parameters: current_time              = 50.816545265311
yt : [INFO     ] 2017-04-14 16:41:55,764 Parameters: domain_dimensions         = [32 32 32]
yt : [INFO     ] 2017-04-14 16:41:55,765 Parameters: domain_left_edge          = [ 0.  0.  0.]
yt : [INFO     ] 2017-04-14 16:41:55,767 Parameters: domain_right_edge         = [ 1.  1.  1.]
yt : [INFO     ] 2017-04-14 16:41:55,768 Parameters: cosmological_simulation   = 1
yt : [INFO     ] 2017-04-14 16:41:55,770 Parameters: current_redshift          = 2.2188733677898
yt : [INFO     ] 2017-04-14 16:41:55,772 Parameters: omega_lambda              = 0.732
yt : [INFO     ] 2017-04-14 16:41:55,774 Parameters: omega_matter              = 0.268
yt : [INFO     ] 2017-04-14 16:41:55,776 Parameters: hubble_constant           = 0.704
Out[11]:
3057.077085898462 Myr

In [14]:
ds = ts[56]


yt : [INFO     ] 2017-04-14 16:43:58,124 Parameters: current_time              = 228.1000985764
yt : [INFO     ] 2017-04-14 16:43:58,126 Parameters: domain_dimensions         = [32 32 32]
yt : [INFO     ] 2017-04-14 16:43:58,127 Parameters: domain_left_edge          = [ 0.  0.  0.]
yt : [INFO     ] 2017-04-14 16:43:58,128 Parameters: domain_right_edge         = [ 1.  1.  1.]
yt : [INFO     ] 2017-04-14 16:43:58,130 Parameters: cosmological_simulation   = 1
yt : [INFO     ] 2017-04-14 16:43:58,130 Parameters: current_redshift          = 0.0069990069523933
yt : [INFO     ] 2017-04-14 16:43:58,131 Parameters: omega_lambda              = 0.732
yt : [INFO     ] 2017-04-14 16:43:58,132 Parameters: omega_matter              = 0.268
yt : [INFO     ] 2017-04-14 16:43:58,133 Parameters: hubble_constant           = 0.704

In [15]:
p = yt.ProjectionPlot(ds, "x", "density")
p.show()


Parsing Hierarchy : 100%|██████████| 135/135 [00:00<00:00, 9883.77it/s]
yt : [INFO     ] 2017-04-14 16:43:59,122 Gathering a field list (this may take a moment.)
yt : [INFO     ] 2017-04-14 16:44:00,523 Projection completed
yt : [INFO     ] 2017-04-14 16:44:00,525 xlim = 0.000000 1.000000
yt : [INFO     ] 2017-04-14 16:44:00,527 ylim = 0.000000 1.000000
yt : [INFO     ] 2017-04-14 16:44:00,530 xlim = 0.000000 1.000000
yt : [INFO     ] 2017-04-14 16:44:00,531 ylim = 0.000000 1.000000
yt : [INFO     ] 2017-04-14 16:44:00,534 Making a fixed resolution buffer of (('gas', 'density')) 800 by 800


In [16]:
ds.field_list


Out[16]:
[('all', 'creation_time'),
 ('all', 'dynamical_time'),
 ('all', 'metallicity_fraction'),
 ('all', 'particle_index'),
 ('all', 'particle_mass'),
 ('all', 'particle_position_x'),
 ('all', 'particle_position_y'),
 ('all', 'particle_position_z'),
 ('all', 'particle_type'),
 ('all', 'particle_velocity_x'),
 ('all', 'particle_velocity_y'),
 ('all', 'particle_velocity_z'),
 ('enzo', 'Dark_Matter_Density'),
 ('enzo', 'Density'),
 ('enzo', 'Electron_Density'),
 ('enzo', 'GasEnergy'),
 ('enzo', 'HII_Density'),
 ('enzo', 'HI_Density'),
 ('enzo', 'HeIII_Density'),
 ('enzo', 'HeII_Density'),
 ('enzo', 'HeI_Density'),
 ('enzo', 'Metal_Density'),
 ('enzo', 'Temperature'),
 ('enzo', 'TotalEnergy'),
 ('enzo', 'x-velocity'),
 ('enzo', 'y-velocity'),
 ('enzo', 'z-velocity'),
 ('io', 'creation_time'),
 ('io', 'dynamical_time'),
 ('io', 'metallicity_fraction'),
 ('io', 'particle_index'),
 ('io', 'particle_mass'),
 ('io', 'particle_position_x'),
 ('io', 'particle_position_y'),
 ('io', 'particle_position_z'),
 ('io', 'particle_type'),
 ('io', 'particle_velocity_x'),
 ('io', 'particle_velocity_y'),
 ('io', 'particle_velocity_z')]

In [17]:
ds.derived_field_list


Out[17]:
[('all', 'age'),
 ('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'),
 ('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'),
 ('enzo', 'Dark_Matter_Density'),
 ('enzo', 'Density'),
 ('enzo', 'Electron_Density'),
 ('enzo', 'GasEnergy'),
 ('enzo', 'HII_Density'),
 ('enzo', 'HI_Density'),
 ('enzo', 'HeIII_Density'),
 ('enzo', 'HeII_Density'),
 ('enzo', 'HeI_Density'),
 ('enzo', 'Metal_Density'),
 ('enzo', 'Temperature'),
 ('enzo', 'TotalEnergy'),
 ('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', 'El_density'),
 ('gas', 'El_fraction'),
 ('gas', 'El_mass'),
 ('gas', 'El_number_density'),
 ('gas', 'H_density'),
 ('gas', 'H_fraction'),
 ('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', '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', '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', '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', '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', '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', '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_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', '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', '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', 'creation_time'),
 ('io', 'dynamical_time'),
 ('io', 'mesh_id'),
 ('io', 'metallicity_fraction'),
 ('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_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')]

In [32]:
sp = ds.sphere([0.5, 0.5, 0.5], 0.1)

In [33]:
sp.max("density")


Out[33]:
1.8989874577524162e-30 g/cm**3

In [34]:
sp.max("density", axis="x").to_pw()


yt : [INFO     ] 2017-04-14 16:53:57,724 Projection completed
yt : [INFO     ] 2017-04-14 16:53:57,726 xlim = 0.000000 1.000000
yt : [INFO     ] 2017-04-14 16:53:57,727 ylim = 0.000000 1.000000
yt : [INFO     ] 2017-04-14 16:53:57,729 Making a fixed resolution buffer of (('gas', 'density')) 800 by 800
Out[34]:


In [35]:
df = sp.to_dataframe(["density", "temperature", "velocity_magnitude"])

In [38]:
pplot = yt.ProfilePlot(sp, "density", "temperature", weight_field="cell_mass")
pplot.show()




In [40]:
source = ds.r[:, :, :]

In [41]:
pplot = yt.ProfilePlot(source, "density", "temperature", weight_field="cell_mass")

In [42]:
pplot.show()




In [45]:
phase_plot = yt.PhasePlot(source, "density", "temperature", "cell_mass", weight_field=None)
phase_plot.show()




In [46]:
phase_avg = yt.PhasePlot(source, "density", "temperature", "velocity_magnitude", weight_field="cell_mass")
phase_avg.show()




In [48]:
phase_avg.plots["velocity_magnitude"].figure


Out[48]:

In [49]:
source


Out[49]:
YTRegion (RD0009): , center=[  6.96416138e+25   6.96416138e+25   6.96416138e+25] cm, left_edge=[ 0.  0.  0.] cm, right_edge=[  1.39283228e+26   1.39283228e+26   1.39283228e+26] cm

In [51]:
source.argmax("density", axis="temperature")


Out[51]:
13489700.909542166 K

In [52]:
source.argmax("density")


Out[52]:
(0.46240234375 code_length,
 0.21240234375 code_length,
 0.57568359375 code_length)

In [54]:
p = yt.SlicePlot(ds, "x", "density", center = source.argmax("density"))


yt : [INFO     ] 2017-04-14 17:04:46,399 xlim = -0.287598 0.712402
yt : [INFO     ] 2017-04-14 17:04:46,400 ylim = 0.075684 1.075684
yt : [INFO     ] 2017-04-14 17:04:46,402 xlim = -0.287598 0.712402
yt : [INFO     ] 2017-04-14 17:04:46,403 ylim = 0.075684 1.075684
yt : [INFO     ] 2017-04-14 17:04:46,405 Making a fixed resolution buffer of (('gas', 'density')) 800 by 800

In [55]:
p.show()




In [56]:
p.zoom(5)


yt : [INFO     ] 2017-04-14 17:04:59,187 Making a fixed resolution buffer of (('gas', 'density')) 800 by 800
Out[56]:


In [58]:
p.annotate_velocity()


Out[58]:


In [62]:
p = yt.SlicePlot(ds, "x", "density", center = source.argmax("density"))


yt : [INFO     ] 2017-04-14 17:06:33,241 xlim = -0.287598 0.712402
yt : [INFO     ] 2017-04-14 17:06:33,241 ylim = 0.075684 1.075684
yt : [INFO     ] 2017-04-14 17:06:33,243 xlim = -0.287598 0.712402
yt : [INFO     ] 2017-04-14 17:06:33,244 ylim = 0.075684 1.075684
yt : [INFO     ] 2017-04-14 17:06:33,246 Making a fixed resolution buffer of (('gas', 'density')) 800 by 800

In [63]:
p.annotate_line_integral_convolution("velocity_y", "velocity_z")
p.show()




In [65]:
p.zoom(5)


yt : [INFO     ] 2017-04-14 17:07:29,305 Making a fixed resolution buffer of (('gas', 'density')) 800 by 800
Out[65]:


In [66]:
p.annotate_clear()


Out[66]:


In [67]:
p.annotate_contour("temperature")


/opt/conda/lib/python3.5/site-packages/numpy/ma/core.py:867: RuntimeWarning: invalid value encountered in less_equal
  return umath.less_equal(x, self.critical_value)
Out[67]:


In [68]:
ds.current_redshift


Out[68]:
0.0069990069523933

In [70]:
ds.current_time.in_units("Myr")


Out[70]:
13722.294205724085 Myr

In [72]:
rho = []
temp = []
z = []
t = []
for ds in ts:
    source = ds.r[:, :, :]
    rho.append(source.max("density"))
    temp.append(source.argmax("density", axis="temperature"))
    z.append(ds.current_redshift)
    t.append(ds.current_time.in_units("Myr"))


yt : [INFO     ] 2017-04-14 17:12:17,950 Parameters: current_time              = 0.81651036015641
yt : [INFO     ] 2017-04-14 17:12:17,951 Parameters: domain_dimensions         = [32 32 32]
yt : [INFO     ] 2017-04-14 17:12:17,952 Parameters: domain_left_edge          = [ 0.  0.  0.]
yt : [INFO     ] 2017-04-14 17:12:17,953 Parameters: domain_right_edge         = [ 1.  1.  1.]
yt : [INFO     ] 2017-04-14 17:12:17,955 Parameters: cosmological_simulation   = 1
yt : [INFO     ] 2017-04-14 17:12:17,955 Parameters: current_redshift          = 50
yt : [INFO     ] 2017-04-14 17:12:17,956 Parameters: omega_lambda              = 0.732
yt : [INFO     ] 2017-04-14 17:12:17,957 Parameters: omega_matter              = 0.268
yt : [INFO     ] 2017-04-14 17:12:17,958 Parameters: hubble_constant           = 0.704
Parsing Hierarchy : 100%|██████████| 2/2 [00:00<00:00, 15768.06it/s]
yt : [INFO     ] 2017-04-14 17:12:17,963 Gathering a field list (this may take a moment.)
yt : [INFO     ] 2017-04-14 17:12:19,079 Parameters: current_time              = 5.8164677145709
yt : [INFO     ] 2017-04-14 17:12:19,080 Parameters: domain_dimensions         = [32 32 32]
yt : [INFO     ] 2017-04-14 17:12:19,081 Parameters: domain_left_edge          = [ 0.  0.  0.]
yt : [INFO     ] 2017-04-14 17:12:19,082 Parameters: domain_right_edge         = [ 1.  1.  1.]
yt : [INFO     ] 2017-04-14 17:12:19,083 Parameters: cosmological_simulation   = 1
yt : [INFO     ] 2017-04-14 17:12:19,084 Parameters: current_redshift          = 12.773860186605
yt : [INFO     ] 2017-04-14 17:12:19,085 Parameters: omega_lambda              = 0.732
yt : [INFO     ] 2017-04-14 17:12:19,085 Parameters: omega_matter              = 0.268
yt : [INFO     ] 2017-04-14 17:12:19,086 Parameters: hubble_constant           = 0.704
Parsing Hierarchy : 100%|██████████| 2/2 [00:00<00:00, 16946.68it/s]
yt : [INFO     ] 2017-04-14 17:12:19,091 Gathering a field list (this may take a moment.)
yt : [INFO     ] 2017-04-14 17:12:20,013 Parameters: current_time              = 10.816460428131
yt : [INFO     ] 2017-04-14 17:12:20,014 Parameters: domain_dimensions         = [32 32 32]
yt : [INFO     ] 2017-04-14 17:12:20,015 Parameters: domain_left_edge          = [ 0.  0.  0.]
yt : [INFO     ] 2017-04-14 17:12:20,016 Parameters: domain_right_edge         = [ 1.  1.  1.]
yt : [INFO     ] 2017-04-14 17:12:20,018 Parameters: cosmological_simulation   = 1
yt : [INFO     ] 2017-04-14 17:12:20,019 Parameters: current_redshift          = 8.1057247460185
yt : [INFO     ] 2017-04-14 17:12:20,019 Parameters: omega_lambda              = 0.732
yt : [INFO     ] 2017-04-14 17:12:20,020 Parameters: omega_matter              = 0.268
yt : [INFO     ] 2017-04-14 17:12:20,021 Parameters: hubble_constant           = 0.704
Parsing Hierarchy : 100%|██████████| 2/2 [00:00<00:00, 16100.98it/s]
yt : [INFO     ] 2017-04-14 17:12:20,026 Gathering a field list (this may take a moment.)
yt : [INFO     ] 2017-04-14 17:12:21,013 Parameters: current_time              = 15.816500763285
yt : [INFO     ] 2017-04-14 17:12:21,014 Parameters: domain_dimensions         = [32 32 32]
yt : [INFO     ] 2017-04-14 17:12:21,015 Parameters: domain_left_edge          = [ 0.  0.  0.]
yt : [INFO     ] 2017-04-14 17:12:21,017 Parameters: domain_right_edge         = [ 1.  1.  1.]
yt : [INFO     ] 2017-04-14 17:12:21,018 Parameters: cosmological_simulation   = 1
yt : [INFO     ] 2017-04-14 17:12:21,018 Parameters: current_redshift          = 6.0647907838631
yt : [INFO     ] 2017-04-14 17:12:21,019 Parameters: omega_lambda              = 0.732
yt : [INFO     ] 2017-04-14 17:12:21,020 Parameters: omega_matter              = 0.268
yt : [INFO     ] 2017-04-14 17:12:21,021 Parameters: hubble_constant           = 0.704
Parsing Hierarchy : 100%|██████████| 2/2 [00:00<00:00, 12018.06it/s]
yt : [INFO     ] 2017-04-14 17:12:21,026 Gathering a field list (this may take a moment.)
yt : [INFO     ] 2017-04-14 17:12:22,209 Parameters: current_time              = 20.816480098769
yt : [INFO     ] 2017-04-14 17:12:22,210 Parameters: domain_dimensions         = [32 32 32]
yt : [INFO     ] 2017-04-14 17:12:22,211 Parameters: domain_left_edge          = [ 0.  0.  0.]
yt : [INFO     ] 2017-04-14 17:12:22,213 Parameters: domain_right_edge         = [ 1.  1.  1.]
yt : [INFO     ] 2017-04-14 17:12:22,214 Parameters: cosmological_simulation   = 1
yt : [INFO     ] 2017-04-14 17:12:22,215 Parameters: current_redshift          = 4.8788890436873
yt : [INFO     ] 2017-04-14 17:12:22,216 Parameters: omega_lambda              = 0.732
yt : [INFO     ] 2017-04-14 17:12:22,218 Parameters: omega_matter              = 0.268
yt : [INFO     ] 2017-04-14 17:12:22,219 Parameters: hubble_constant           = 0.704
Parsing Hierarchy : 100%|██████████| 3/3 [00:00<00:00, 14631.29it/s]
yt : [INFO     ] 2017-04-14 17:12:22,225 Gathering a field list (this may take a moment.)
yt : [INFO     ] 2017-04-14 17:12:23,183 Parameters: current_time              = 25.816484630092
yt : [INFO     ] 2017-04-14 17:12:23,184 Parameters: domain_dimensions         = [32 32 32]
yt : [INFO     ] 2017-04-14 17:12:23,186 Parameters: domain_left_edge          = [ 0.  0.  0.]
yt : [INFO     ] 2017-04-14 17:12:23,188 Parameters: domain_right_edge         = [ 1.  1.  1.]
yt : [INFO     ] 2017-04-14 17:12:23,190 Parameters: cosmological_simulation   = 1
yt : [INFO     ] 2017-04-14 17:12:23,191 Parameters: current_redshift          = 4.0888731438237
yt : [INFO     ] 2017-04-14 17:12:23,192 Parameters: omega_lambda              = 0.732
yt : [INFO     ] 2017-04-14 17:12:23,194 Parameters: omega_matter              = 0.268
yt : [INFO     ] 2017-04-14 17:12:23,195 Parameters: hubble_constant           = 0.704
Parsing Hierarchy : 100%|██████████| 6/6 [00:00<00:00, 24818.37it/s]
yt : [INFO     ] 2017-04-14 17:12:23,202 Gathering a field list (this may take a moment.)
yt : [INFO     ] 2017-04-14 17:12:24,153 Parameters: current_time              = 30.816463652018
yt : [INFO     ] 2017-04-14 17:12:24,154 Parameters: domain_dimensions         = [32 32 32]
yt : [INFO     ] 2017-04-14 17:12:24,157 Parameters: domain_left_edge          = [ 0.  0.  0.]
yt : [INFO     ] 2017-04-14 17:12:24,158 Parameters: domain_right_edge         = [ 1.  1.  1.]
yt : [INFO     ] 2017-04-14 17:12:24,160 Parameters: cosmological_simulation   = 1
yt : [INFO     ] 2017-04-14 17:12:24,161 Parameters: current_redshift          = 3.5179684066031
yt : [INFO     ] 2017-04-14 17:12:24,162 Parameters: omega_lambda              = 0.732
yt : [INFO     ] 2017-04-14 17:12:24,164 Parameters: omega_matter              = 0.268
yt : [INFO     ] 2017-04-14 17:12:24,165 Parameters: hubble_constant           = 0.704
Parsing Hierarchy : 100%|██████████| 13/13 [00:00<00:00, 24266.11it/s]
yt : [INFO     ] 2017-04-14 17:12:24,173 Gathering a field list (this may take a moment.)
yt : [INFO     ] 2017-04-14 17:12:25,413 Parameters: current_time              = 35.816475202361
yt : [INFO     ] 2017-04-14 17:12:25,414 Parameters: domain_dimensions         = [32 32 32]
yt : [INFO     ] 2017-04-14 17:12:25,416 Parameters: domain_left_edge          = [ 0.  0.  0.]
yt : [INFO     ] 2017-04-14 17:12:25,417 Parameters: domain_right_edge         = [ 1.  1.  1.]
yt : [INFO     ] 2017-04-14 17:12:25,419 Parameters: cosmological_simulation   = 1
yt : [INFO     ] 2017-04-14 17:12:25,420 Parameters: current_redshift          = 3.0823895109346
yt : [INFO     ] 2017-04-14 17:12:25,422 Parameters: omega_lambda              = 0.732
yt : [INFO     ] 2017-04-14 17:12:25,423 Parameters: omega_matter              = 0.268
yt : [INFO     ] 2017-04-14 17:12:25,424 Parameters: hubble_constant           = 0.704
Parsing Hierarchy : 100%|██████████| 15/15 [00:00<00:00, 26468.05it/s]
yt : [INFO     ] 2017-04-14 17:12:25,431 Gathering a field list (this may take a moment.)
yt : [INFO     ] 2017-04-14 17:12:26,473 Parameters: current_time              = 40.816517503487
yt : [INFO     ] 2017-04-14 17:12:26,474 Parameters: domain_dimensions         = [32 32 32]
yt : [INFO     ] 2017-04-14 17:12:26,476 Parameters: domain_left_edge          = [ 0.  0.  0.]
yt : [INFO     ] 2017-04-14 17:12:26,478 Parameters: domain_right_edge         = [ 1.  1.  1.]
yt : [INFO     ] 2017-04-14 17:12:26,479 Parameters: cosmological_simulation   = 1
yt : [INFO     ] 2017-04-14 17:12:26,481 Parameters: current_redshift          = 2.7368848656654
yt : [INFO     ] 2017-04-14 17:12:26,482 Parameters: omega_lambda              = 0.732
yt : [INFO     ] 2017-04-14 17:12:26,483 Parameters: omega_matter              = 0.268
yt : [INFO     ] 2017-04-14 17:12:26,484 Parameters: hubble_constant           = 0.704
Parsing Hierarchy : 100%|██████████| 18/18 [00:00<00:00, 28739.05it/s]
yt : [INFO     ] 2017-04-14 17:12:26,492 Gathering a field list (this may take a moment.)
yt : [INFO     ] 2017-04-14 17:12:27,482 Parameters: current_time              = 45.816588819493
yt : [INFO     ] 2017-04-14 17:12:27,483 Parameters: domain_dimensions         = [32 32 32]
yt : [INFO     ] 2017-04-14 17:12:27,485 Parameters: domain_left_edge          = [ 0.  0.  0.]
yt : [INFO     ] 2017-04-14 17:12:27,486 Parameters: domain_right_edge         = [ 1.  1.  1.]
yt : [INFO     ] 2017-04-14 17:12:27,487 Parameters: cosmological_simulation   = 1
yt : [INFO     ] 2017-04-14 17:12:27,489 Parameters: current_redshift          = 2.4546875080899
yt : [INFO     ] 2017-04-14 17:12:27,490 Parameters: omega_lambda              = 0.732
yt : [INFO     ] 2017-04-14 17:12:27,491 Parameters: omega_matter              = 0.268
yt : [INFO     ] 2017-04-14 17:12:27,494 Parameters: hubble_constant           = 0.704
Parsing Hierarchy : 100%|██████████| 25/25 [00:00<00:00, 23327.61it/s]
yt : [INFO     ] 2017-04-14 17:12:27,502 Gathering a field list (this may take a moment.)
yt : [INFO     ] 2017-04-14 17:12:28,713 Parameters: current_time              = 50.816545265311
yt : [INFO     ] 2017-04-14 17:12:28,714 Parameters: domain_dimensions         = [32 32 32]
yt : [INFO     ] 2017-04-14 17:12:28,716 Parameters: domain_left_edge          = [ 0.  0.  0.]
yt : [INFO     ] 2017-04-14 17:12:28,717 Parameters: domain_right_edge         = [ 1.  1.  1.]
yt : [INFO     ] 2017-04-14 17:12:28,719 Parameters: cosmological_simulation   = 1
yt : [INFO     ] 2017-04-14 17:12:28,720 Parameters: current_redshift          = 2.2188733677898
yt : [INFO     ] 2017-04-14 17:12:28,721 Parameters: omega_lambda              = 0.732
yt : [INFO     ] 2017-04-14 17:12:28,722 Parameters: omega_matter              = 0.268
yt : [INFO     ] 2017-04-14 17:12:28,723 Parameters: hubble_constant           = 0.704
Parsing Hierarchy : 100%|██████████| 32/32 [00:00<00:00, 24447.67it/s]
yt : [INFO     ] 2017-04-14 17:12:28,731 Gathering a field list (this may take a moment.)
yt : [INFO     ] 2017-04-14 17:12:29,857 Parameters: current_time              = 55.81650237435
yt : [INFO     ] 2017-04-14 17:12:29,858 Parameters: domain_dimensions         = [32 32 32]
yt : [INFO     ] 2017-04-14 17:12:29,860 Parameters: domain_left_edge          = [ 0.  0.  0.]
yt : [INFO     ] 2017-04-14 17:12:29,861 Parameters: domain_right_edge         = [ 1.  1.  1.]
yt : [INFO     ] 2017-04-14 17:12:29,862 Parameters: cosmological_simulation   = 1
yt : [INFO     ] 2017-04-14 17:12:29,864 Parameters: current_redshift          = 2.0181601157575
yt : [INFO     ] 2017-04-14 17:12:29,865 Parameters: omega_lambda              = 0.732
yt : [INFO     ] 2017-04-14 17:12:29,866 Parameters: omega_matter              = 0.268
yt : [INFO     ] 2017-04-14 17:12:29,867 Parameters: hubble_constant           = 0.704
Parsing Hierarchy : 100%|██████████| 42/42 [00:00<00:00, 25664.45it/s]
yt : [INFO     ] 2017-04-14 17:12:29,876 Gathering a field list (this may take a moment.)
yt : [INFO     ] 2017-04-14 17:12:30,887 Parameters: current_time              = 60.816631390806
yt : [INFO     ] 2017-04-14 17:12:30,888 Parameters: domain_dimensions         = [32 32 32]
yt : [INFO     ] 2017-04-14 17:12:30,889 Parameters: domain_left_edge          = [ 0.  0.  0.]
yt : [INFO     ] 2017-04-14 17:12:30,890 Parameters: domain_right_edge         = [ 1.  1.  1.]
yt : [INFO     ] 2017-04-14 17:12:30,892 Parameters: cosmological_simulation   = 1
yt : [INFO     ] 2017-04-14 17:12:30,893 Parameters: current_redshift          = 1.8447188773741
yt : [INFO     ] 2017-04-14 17:12:30,895 Parameters: omega_lambda              = 0.732
yt : [INFO     ] 2017-04-14 17:12:30,896 Parameters: omega_matter              = 0.268
yt : [INFO     ] 2017-04-14 17:12:30,897 Parameters: hubble_constant           = 0.704
Parsing Hierarchy : 100%|██████████| 49/49 [00:00<00:00, 26632.23it/s]
yt : [INFO     ] 2017-04-14 17:12:30,906 Gathering a field list (this may take a moment.)
yt : [INFO     ] 2017-04-14 17:12:32,120 Parameters: current_time              = 65.816604038397
yt : [INFO     ] 2017-04-14 17:12:32,121 Parameters: domain_dimensions         = [32 32 32]
yt : [INFO     ] 2017-04-14 17:12:32,123 Parameters: domain_left_edge          = [ 0.  0.  0.]
yt : [INFO     ] 2017-04-14 17:12:32,125 Parameters: domain_right_edge         = [ 1.  1.  1.]
yt : [INFO     ] 2017-04-14 17:12:32,126 Parameters: cosmological_simulation   = 1
yt : [INFO     ] 2017-04-14 17:12:32,127 Parameters: current_redshift          = 1.6929484501255
yt : [INFO     ] 2017-04-14 17:12:32,129 Parameters: omega_lambda              = 0.732
yt : [INFO     ] 2017-04-14 17:12:32,130 Parameters: omega_matter              = 0.268
yt : [INFO     ] 2017-04-14 17:12:32,131 Parameters: hubble_constant           = 0.704
Parsing Hierarchy : 100%|██████████| 60/60 [00:00<00:00, 27676.04it/s]
yt : [INFO     ] 2017-04-14 17:12:32,142 Gathering a field list (this may take a moment.)
yt : [INFO     ] 2017-04-14 17:12:33,183 Parameters: current_time              = 70.816577212426
yt : [INFO     ] 2017-04-14 17:12:33,184 Parameters: domain_dimensions         = [32 32 32]
yt : [INFO     ] 2017-04-14 17:12:33,187 Parameters: domain_left_edge          = [ 0.  0.  0.]
yt : [INFO     ] 2017-04-14 17:12:33,188 Parameters: domain_right_edge         = [ 1.  1.  1.]
yt : [INFO     ] 2017-04-14 17:12:33,190 Parameters: cosmological_simulation   = 1
yt : [INFO     ] 2017-04-14 17:12:33,191 Parameters: current_redshift          = 1.5587032944109
yt : [INFO     ] 2017-04-14 17:12:33,192 Parameters: omega_lambda              = 0.732
yt : [INFO     ] 2017-04-14 17:12:33,193 Parameters: omega_matter              = 0.268
yt : [INFO     ] 2017-04-14 17:12:33,195 Parameters: hubble_constant           = 0.704
Parsing Hierarchy : 100%|██████████| 61/61 [00:00<00:00, 28196.22it/s]
yt : [INFO     ] 2017-04-14 17:12:33,205 Gathering a field list (this may take a moment.)
yt : [INFO     ] 2017-04-14 17:12:34,343 Parameters: current_time              = 75.81655092386
yt : [INFO     ] 2017-04-14 17:12:34,344 Parameters: domain_dimensions         = [32 32 32]
yt : [INFO     ] 2017-04-14 17:12:34,347 Parameters: domain_left_edge          = [ 0.  0.  0.]
yt : [INFO     ] 2017-04-14 17:12:34,349 Parameters: domain_right_edge         = [ 1.  1.  1.]
yt : [INFO     ] 2017-04-14 17:12:34,350 Parameters: cosmological_simulation   = 1
yt : [INFO     ] 2017-04-14 17:12:34,352 Parameters: current_redshift          = 1.4388577442766
yt : [INFO     ] 2017-04-14 17:12:34,353 Parameters: omega_lambda              = 0.732
yt : [INFO     ] 2017-04-14 17:12:34,354 Parameters: omega_matter              = 0.268
yt : [INFO     ] 2017-04-14 17:12:34,356 Parameters: hubble_constant           = 0.704
Parsing Hierarchy : 100%|██████████| 60/60 [00:00<00:00, 12779.72it/s]
yt : [INFO     ] 2017-04-14 17:12:34,369 Gathering a field list (this may take a moment.)
yt : [INFO     ] 2017-04-14 17:12:35,772 Parameters: current_time              = 80.816525197391
yt : [INFO     ] 2017-04-14 17:12:35,773 Parameters: domain_dimensions         = [32 32 32]
yt : [INFO     ] 2017-04-14 17:12:35,775 Parameters: domain_left_edge          = [ 0.  0.  0.]
yt : [INFO     ] 2017-04-14 17:12:35,777 Parameters: domain_right_edge         = [ 1.  1.  1.]
yt : [INFO     ] 2017-04-14 17:12:35,779 Parameters: cosmological_simulation   = 1
yt : [INFO     ] 2017-04-14 17:12:35,780 Parameters: current_redshift          = 1.3310040593995
yt : [INFO     ] 2017-04-14 17:12:35,782 Parameters: omega_lambda              = 0.732
yt : [INFO     ] 2017-04-14 17:12:35,783 Parameters: omega_matter              = 0.268
yt : [INFO     ] 2017-04-14 17:12:35,785 Parameters: hubble_constant           = 0.704
Parsing Hierarchy : 100%|██████████| 61/61 [00:00<00:00, 26225.15it/s]
yt : [INFO     ] 2017-04-14 17:12:35,795 Gathering a field list (this may take a moment.)
yt : [INFO     ] 2017-04-14 17:12:36,914 Parameters: current_time              = 85.816500054852
yt : [INFO     ] 2017-04-14 17:12:36,915 Parameters: domain_dimensions         = [32 32 32]
yt : [INFO     ] 2017-04-14 17:12:36,917 Parameters: domain_left_edge          = [ 0.  0.  0.]
yt : [INFO     ] 2017-04-14 17:12:36,918 Parameters: domain_right_edge         = [ 1.  1.  1.]
yt : [INFO     ] 2017-04-14 17:12:36,920 Parameters: cosmological_simulation   = 1
yt : [INFO     ] 2017-04-14 17:12:36,921 Parameters: current_redshift          = 1.2332551244551
yt : [INFO     ] 2017-04-14 17:12:36,922 Parameters: omega_lambda              = 0.732
yt : [INFO     ] 2017-04-14 17:12:36,924 Parameters: omega_matter              = 0.268
yt : [INFO     ] 2017-04-14 17:12:36,925 Parameters: hubble_constant           = 0.704
Parsing Hierarchy : 100%|██████████| 70/70 [00:00<00:00, 26235.48it/s]
yt : [INFO     ] 2017-04-14 17:12:36,936 Gathering a field list (this may take a moment.)
yt : [INFO     ] 2017-04-14 17:12:37,972 Parameters: current_time              = 90.816475515226
yt : [INFO     ] 2017-04-14 17:12:37,973 Parameters: domain_dimensions         = [32 32 32]
yt : [INFO     ] 2017-04-14 17:12:37,976 Parameters: domain_left_edge          = [ 0.  0.  0.]
yt : [INFO     ] 2017-04-14 17:12:37,977 Parameters: domain_right_edge         = [ 1.  1.  1.]
yt : [INFO     ] 2017-04-14 17:12:37,979 Parameters: cosmological_simulation   = 1
yt : [INFO     ] 2017-04-14 17:12:37,981 Parameters: current_redshift          = 1.1441093230359
yt : [INFO     ] 2017-04-14 17:12:37,982 Parameters: omega_lambda              = 0.732
yt : [INFO     ] 2017-04-14 17:12:37,983 Parameters: omega_matter              = 0.268
yt : [INFO     ] 2017-04-14 17:12:37,985 Parameters: hubble_constant           = 0.704
Parsing Hierarchy : 100%|██████████| 76/76 [00:00<00:00, 27656.35it/s]
yt : [INFO     ] 2017-04-14 17:12:37,996 Gathering a field list (this may take a moment.)
yt : [INFO     ] 2017-04-14 17:12:39,267 Parameters: current_time              = 95.816708910656
yt : [INFO     ] 2017-04-14 17:12:39,268 Parameters: domain_dimensions         = [32 32 32]
yt : [INFO     ] 2017-04-14 17:12:39,270 Parameters: domain_left_edge          = [ 0.  0.  0.]
yt : [INFO     ] 2017-04-14 17:12:39,271 Parameters: domain_right_edge         = [ 1.  1.  1.]
yt : [INFO     ] 2017-04-14 17:12:39,273 Parameters: cosmological_simulation   = 1
yt : [INFO     ] 2017-04-14 17:12:39,275 Parameters: current_redshift          = 1.0623517145997
yt : [INFO     ] 2017-04-14 17:12:39,276 Parameters: omega_lambda              = 0.732
yt : [INFO     ] 2017-04-14 17:12:39,277 Parameters: omega_matter              = 0.268
yt : [INFO     ] 2017-04-14 17:12:39,279 Parameters: hubble_constant           = 0.704
Parsing Hierarchy : 100%|██████████| 84/84 [00:00<00:00, 24849.88it/s]
yt : [INFO     ] 2017-04-14 17:12:39,291 Gathering a field list (this may take a moment.)
yt : [INFO     ] 2017-04-14 17:12:40,392 Parameters: current_time              = 100.81669729468
yt : [INFO     ] 2017-04-14 17:12:40,393 Parameters: domain_dimensions         = [32 32 32]
yt : [INFO     ] 2017-04-14 17:12:40,395 Parameters: domain_left_edge          = [ 0.  0.  0.]
yt : [INFO     ] 2017-04-14 17:12:40,397 Parameters: domain_right_edge         = [ 1.  1.  1.]
yt : [INFO     ] 2017-04-14 17:12:40,398 Parameters: cosmological_simulation   = 1
yt : [INFO     ] 2017-04-14 17:12:40,399 Parameters: current_redshift          = 0.98700240373113
yt : [INFO     ] 2017-04-14 17:12:40,400 Parameters: omega_lambda              = 0.732
yt : [INFO     ] 2017-04-14 17:12:40,401 Parameters: omega_matter              = 0.268
yt : [INFO     ] 2017-04-14 17:12:40,403 Parameters: hubble_constant           = 0.704
Parsing Hierarchy : 100%|██████████| 88/88 [00:00<00:00, 28914.90it/s]
yt : [INFO     ] 2017-04-14 17:12:40,414 Gathering a field list (this may take a moment.)
yt : [INFO     ] 2017-04-14 17:12:41,672 Parameters: current_time              = 105.81668601342
yt : [INFO     ] 2017-04-14 17:12:41,674 Parameters: domain_dimensions         = [32 32 32]
yt : [INFO     ] 2017-04-14 17:12:41,676 Parameters: domain_left_edge          = [ 0.  0.  0.]
yt : [INFO     ] 2017-04-14 17:12:41,678 Parameters: domain_right_edge         = [ 1.  1.  1.]
yt : [INFO     ] 2017-04-14 17:12:41,680 Parameters: cosmological_simulation   = 1
yt : [INFO     ] 2017-04-14 17:12:41,681 Parameters: current_redshift          = 0.91724233768264
yt : [INFO     ] 2017-04-14 17:12:41,682 Parameters: omega_lambda              = 0.732
yt : [INFO     ] 2017-04-14 17:12:41,684 Parameters: omega_matter              = 0.268
yt : [INFO     ] 2017-04-14 17:12:41,685 Parameters: hubble_constant           = 0.704
Parsing Hierarchy : 100%|██████████| 91/91 [00:00<00:00, 26128.26it/s]
yt : [INFO     ] 2017-04-14 17:12:41,699 Gathering a field list (this may take a moment.)
yt : [INFO     ] 2017-04-14 17:12:42,789 Parameters: current_time              = 110.81667504624
yt : [INFO     ] 2017-04-14 17:12:42,790 Parameters: domain_dimensions         = [32 32 32]
yt : [INFO     ] 2017-04-14 17:12:42,791 Parameters: domain_left_edge          = [ 0.  0.  0.]
yt : [INFO     ] 2017-04-14 17:12:42,793 Parameters: domain_right_edge         = [ 1.  1.  1.]
yt : [INFO     ] 2017-04-14 17:12:42,795 Parameters: cosmological_simulation   = 1
yt : [INFO     ] 2017-04-14 17:12:42,796 Parameters: current_redshift          = 0.85239308960285
yt : [INFO     ] 2017-04-14 17:12:42,798 Parameters: omega_lambda              = 0.732
yt : [INFO     ] 2017-04-14 17:12:42,799 Parameters: omega_matter              = 0.268
yt : [INFO     ] 2017-04-14 17:12:42,800 Parameters: hubble_constant           = 0.704
Parsing Hierarchy : 100%|██████████| 97/97 [00:00<00:00, 18451.13it/s]
yt : [INFO     ] 2017-04-14 17:12:42,814 Gathering a field list (this may take a moment.)
yt : [INFO     ] 2017-04-14 17:12:43,937 Parameters: current_time              = 115.81666439643
yt : [INFO     ] 2017-04-14 17:12:43,938 Parameters: domain_dimensions         = [32 32 32]
yt : [INFO     ] 2017-04-14 17:12:43,940 Parameters: domain_left_edge          = [ 0.  0.  0.]
yt : [INFO     ] 2017-04-14 17:12:43,942 Parameters: domain_right_edge         = [ 1.  1.  1.]
yt : [INFO     ] 2017-04-14 17:12:43,944 Parameters: cosmological_simulation   = 1
yt : [INFO     ] 2017-04-14 17:12:43,945 Parameters: current_redshift          = 0.79188495366179
yt : [INFO     ] 2017-04-14 17:12:43,946 Parameters: omega_lambda              = 0.732
yt : [INFO     ] 2017-04-14 17:12:43,947 Parameters: omega_matter              = 0.268
yt : [INFO     ] 2017-04-14 17:12:43,948 Parameters: hubble_constant           = 0.704
Parsing Hierarchy : 100%|██████████| 101/101 [00:00<00:00, 28804.29it/s]
yt : [INFO     ] 2017-04-14 17:12:43,963 Gathering a field list (this may take a moment.)
yt : [INFO     ] 2017-04-14 17:12:45,315 Parameters: current_time              = 120.81665406615
yt : [INFO     ] 2017-04-14 17:12:45,316 Parameters: domain_dimensions         = [32 32 32]
yt : [INFO     ] 2017-04-14 17:12:45,319 Parameters: domain_left_edge          = [ 0.  0.  0.]
yt : [INFO     ] 2017-04-14 17:12:45,320 Parameters: domain_right_edge         = [ 1.  1.  1.]
yt : [INFO     ] 2017-04-14 17:12:45,322 Parameters: cosmological_simulation   = 1
yt : [INFO     ] 2017-04-14 17:12:45,323 Parameters: current_redshift          = 0.73523577485216
yt : [INFO     ] 2017-04-14 17:12:45,325 Parameters: omega_lambda              = 0.732
yt : [INFO     ] 2017-04-14 17:12:45,326 Parameters: omega_matter              = 0.268
yt : [INFO     ] 2017-04-14 17:12:45,328 Parameters: hubble_constant           = 0.704
Parsing Hierarchy : 100%|██████████| 110/110 [00:00<00:00, 26175.73it/s]
yt : [INFO     ] 2017-04-14 17:12:45,343 Gathering a field list (this may take a moment.)
yt : [INFO     ] 2017-04-14 17:12:46,494 Parameters: current_time              = 125.81664405645
yt : [INFO     ] 2017-04-14 17:12:46,496 Parameters: domain_dimensions         = [32 32 32]
yt : [INFO     ] 2017-04-14 17:12:46,498 Parameters: domain_left_edge          = [ 0.  0.  0.]
yt : [INFO     ] 2017-04-14 17:12:46,500 Parameters: domain_right_edge         = [ 1.  1.  1.]
yt : [INFO     ] 2017-04-14 17:12:46,501 Parameters: cosmological_simulation   = 1
yt : [INFO     ] 2017-04-14 17:12:46,503 Parameters: current_redshift          = 0.68203458028522
yt : [INFO     ] 2017-04-14 17:12:46,505 Parameters: omega_lambda              = 0.732
yt : [INFO     ] 2017-04-14 17:12:46,507 Parameters: omega_matter              = 0.268
yt : [INFO     ] 2017-04-14 17:12:46,509 Parameters: hubble_constant           = 0.704
Parsing Hierarchy : 100%|██████████| 111/111 [00:00<00:00, 23013.73it/s]
yt : [INFO     ] 2017-04-14 17:12:46,527 Gathering a field list (this may take a moment.)
yt : [INFO     ] 2017-04-14 17:12:47,635 Parameters: current_time              = 130.81663436738
yt : [INFO     ] 2017-04-14 17:12:47,636 Parameters: domain_dimensions         = [32 32 32]
yt : [INFO     ] 2017-04-14 17:12:47,638 Parameters: domain_left_edge          = [ 0.  0.  0.]
yt : [INFO     ] 2017-04-14 17:12:47,639 Parameters: domain_right_edge         = [ 1.  1.  1.]
yt : [INFO     ] 2017-04-14 17:12:47,641 Parameters: cosmological_simulation   = 1
yt : [INFO     ] 2017-04-14 17:12:47,642 Parameters: current_redshift          = 0.63192878032802
yt : [INFO     ] 2017-04-14 17:12:47,643 Parameters: omega_lambda              = 0.732
yt : [INFO     ] 2017-04-14 17:12:47,644 Parameters: omega_matter              = 0.268
yt : [INFO     ] 2017-04-14 17:12:47,646 Parameters: hubble_constant           = 0.704
Parsing Hierarchy : 100%|██████████| 122/122 [00:00<00:00, 28485.03it/s]
yt : [INFO     ] 2017-04-14 17:12:47,659 Gathering a field list (this may take a moment.)
yt : [INFO     ] 2017-04-14 17:12:48,905 Parameters: current_time              = 135.81662904138
yt : [INFO     ] 2017-04-14 17:12:48,906 Parameters: domain_dimensions         = [32 32 32]
yt : [INFO     ] 2017-04-14 17:12:48,909 Parameters: domain_left_edge          = [ 0.  0.  0.]
yt : [INFO     ] 2017-04-14 17:12:48,910 Parameters: domain_right_edge         = [ 1.  1.  1.]
yt : [INFO     ] 2017-04-14 17:12:48,912 Parameters: cosmological_simulation   = 1
yt : [INFO     ] 2017-04-14 17:12:48,913 Parameters: current_redshift          = 0.58461402038788
yt : [INFO     ] 2017-04-14 17:12:48,915 Parameters: omega_lambda              = 0.732
yt : [INFO     ] 2017-04-14 17:12:48,916 Parameters: omega_matter              = 0.268
yt : [INFO     ] 2017-04-14 17:12:48,918 Parameters: hubble_constant           = 0.704
Parsing Hierarchy : 100%|██████████| 126/126 [00:00<00:00, 23610.88it/s]
yt : [INFO     ] 2017-04-14 17:12:48,935 Gathering a field list (this may take a moment.)
yt : [INFO     ] 2017-04-14 17:12:49,987 Parameters: current_time              = 140.8166291087
yt : [INFO     ] 2017-04-14 17:12:49,988 Parameters: domain_dimensions         = [32 32 32]
yt : [INFO     ] 2017-04-14 17:12:49,991 Parameters: domain_left_edge          = [ 0.  0.  0.]
yt : [INFO     ] 2017-04-14 17:12:49,992 Parameters: domain_right_edge         = [ 1.  1.  1.]
yt : [INFO     ] 2017-04-14 17:12:49,993 Parameters: cosmological_simulation   = 1
yt : [INFO     ] 2017-04-14 17:12:49,995 Parameters: current_redshift          = 0.53982618866656
yt : [INFO     ] 2017-04-14 17:12:49,996 Parameters: omega_lambda              = 0.732
yt : [INFO     ] 2017-04-14 17:12:49,997 Parameters: omega_matter              = 0.268
yt : [INFO     ] 2017-04-14 17:12:49,999 Parameters: hubble_constant           = 0.704
Parsing Hierarchy : 100%|██████████| 118/118 [00:00<00:00, 22716.66it/s]
yt : [INFO     ] 2017-04-14 17:12:50,013 Gathering a field list (this may take a moment.)
yt : [INFO     ] 2017-04-14 17:12:51,169 Parameters: current_time              = 145.81663725255
yt : [INFO     ] 2017-04-14 17:12:51,170 Parameters: domain_dimensions         = [32 32 32]
yt : [INFO     ] 2017-04-14 17:12:51,171 Parameters: domain_left_edge          = [ 0.  0.  0.]
yt : [INFO     ] 2017-04-14 17:12:51,173 Parameters: domain_right_edge         = [ 1.  1.  1.]
yt : [INFO     ] 2017-04-14 17:12:51,175 Parameters: cosmological_simulation   = 1
yt : [INFO     ] 2017-04-14 17:12:51,176 Parameters: current_redshift          = 0.49733488223364
yt : [INFO     ] 2017-04-14 17:12:51,177 Parameters: omega_lambda              = 0.732
yt : [INFO     ] 2017-04-14 17:12:51,179 Parameters: omega_matter              = 0.268
yt : [INFO     ] 2017-04-14 17:12:51,180 Parameters: hubble_constant           = 0.704
Parsing Hierarchy : 100%|██████████| 122/122 [00:00<00:00, 27598.57it/s]
yt : [INFO     ] 2017-04-14 17:12:51,196 Gathering a field list (this may take a moment.)
yt : [INFO     ] 2017-04-14 17:12:52,529 Parameters: current_time              = 150.81666762016
yt : [INFO     ] 2017-04-14 17:12:52,530 Parameters: domain_dimensions         = [32 32 32]
yt : [INFO     ] 2017-04-14 17:12:52,531 Parameters: domain_left_edge          = [ 0.  0.  0.]
yt : [INFO     ] 2017-04-14 17:12:52,533 Parameters: domain_right_edge         = [ 1.  1.  1.]
yt : [INFO     ] 2017-04-14 17:12:52,535 Parameters: cosmological_simulation   = 1
yt : [INFO     ] 2017-04-14 17:12:52,536 Parameters: current_redshift          = 0.4569380718956
yt : [INFO     ] 2017-04-14 17:12:52,538 Parameters: omega_lambda              = 0.732
yt : [INFO     ] 2017-04-14 17:12:52,539 Parameters: omega_matter              = 0.268
yt : [INFO     ] 2017-04-14 17:12:52,540 Parameters: hubble_constant           = 0.704
Parsing Hierarchy : 100%|██████████| 124/124 [00:00<00:00, 19911.70it/s]
yt : [INFO     ] 2017-04-14 17:12:52,558 Gathering a field list (this may take a moment.)
yt : [INFO     ] 2017-04-14 17:12:53,735 Parameters: current_time              = 155.81669000627
yt : [INFO     ] 2017-04-14 17:12:53,736 Parameters: domain_dimensions         = [32 32 32]
yt : [INFO     ] 2017-04-14 17:12:53,739 Parameters: domain_left_edge          = [ 0.  0.  0.]
yt : [INFO     ] 2017-04-14 17:12:53,742 Parameters: domain_right_edge         = [ 1.  1.  1.]
yt : [INFO     ] 2017-04-14 17:12:53,745 Parameters: cosmological_simulation   = 1
yt : [INFO     ] 2017-04-14 17:12:53,747 Parameters: current_redshift          = 0.418458243302
yt : [INFO     ] 2017-04-14 17:12:53,749 Parameters: omega_lambda              = 0.732
yt : [INFO     ] 2017-04-14 17:12:53,751 Parameters: omega_matter              = 0.268
yt : [INFO     ] 2017-04-14 17:12:53,753 Parameters: hubble_constant           = 0.704
Parsing Hierarchy : 100%|██████████| 121/121 [00:00<00:00, 18773.75it/s]
yt : [INFO     ] 2017-04-14 17:12:53,771 Gathering a field list (this may take a moment.)
yt : [INFO     ] 2017-04-14 17:12:54,920 Parameters: current_time              = 160.81646253715
yt : [INFO     ] 2017-04-14 17:12:54,922 Parameters: domain_dimensions         = [32 32 32]
yt : [INFO     ] 2017-04-14 17:12:54,924 Parameters: domain_left_edge          = [ 0.  0.  0.]
yt : [INFO     ] 2017-04-14 17:12:54,925 Parameters: domain_right_edge         = [ 1.  1.  1.]
yt : [INFO     ] 2017-04-14 17:12:54,927 Parameters: cosmological_simulation   = 1
yt : [INFO     ] 2017-04-14 17:12:54,928 Parameters: current_redshift          = 0.38173998206753
yt : [INFO     ] 2017-04-14 17:12:54,929 Parameters: omega_lambda              = 0.732
yt : [INFO     ] 2017-04-14 17:12:54,930 Parameters: omega_matter              = 0.268
yt : [INFO     ] 2017-04-14 17:12:54,931 Parameters: hubble_constant           = 0.704
Parsing Hierarchy : 100%|██████████| 121/121 [00:00<00:00, 24032.14it/s]
yt : [INFO     ] 2017-04-14 17:12:54,948 Gathering a field list (this may take a moment.)
yt : [INFO     ] 2017-04-14 17:12:56,336 Parameters: current_time              = 165.81648811368
yt : [INFO     ] 2017-04-14 17:12:56,337 Parameters: domain_dimensions         = [32 32 32]
yt : [INFO     ] 2017-04-14 17:12:56,339 Parameters: domain_left_edge          = [ 0.  0.  0.]
yt : [INFO     ] 2017-04-14 17:12:56,341 Parameters: domain_right_edge         = [ 1.  1.  1.]
yt : [INFO     ] 2017-04-14 17:12:56,343 Parameters: cosmological_simulation   = 1
yt : [INFO     ] 2017-04-14 17:12:56,344 Parameters: current_redshift          = 0.34664015640906
yt : [INFO     ] 2017-04-14 17:12:56,346 Parameters: omega_lambda              = 0.732
yt : [INFO     ] 2017-04-14 17:12:56,347 Parameters: omega_matter              = 0.268
yt : [INFO     ] 2017-04-14 17:12:56,348 Parameters: hubble_constant           = 0.704
Parsing Hierarchy : 100%|██████████| 115/115 [00:00<00:00, 18918.46it/s]
yt : [INFO     ] 2017-04-14 17:12:56,364 Gathering a field list (this may take a moment.)
yt : [INFO     ] 2017-04-14 17:12:57,522 Parameters: current_time              = 170.81647985764
yt : [INFO     ] 2017-04-14 17:12:57,523 Parameters: domain_dimensions         = [32 32 32]
yt : [INFO     ] 2017-04-14 17:12:57,525 Parameters: domain_left_edge          = [ 0.  0.  0.]
yt : [INFO     ] 2017-04-14 17:12:57,526 Parameters: domain_right_edge         = [ 1.  1.  1.]
yt : [INFO     ] 2017-04-14 17:12:57,528 Parameters: cosmological_simulation   = 1
yt : [INFO     ] 2017-04-14 17:12:57,529 Parameters: current_redshift          = 0.31303669951134
yt : [INFO     ] 2017-04-14 17:12:57,531 Parameters: omega_lambda              = 0.732
yt : [INFO     ] 2017-04-14 17:12:57,532 Parameters: omega_matter              = 0.268
yt : [INFO     ] 2017-04-14 17:12:57,534 Parameters: hubble_constant           = 0.704
Parsing Hierarchy : 100%|██████████| 128/128 [00:00<00:00, 21936.38it/s]
yt : [INFO     ] 2017-04-14 17:12:57,551 Gathering a field list (this may take a moment.)
yt : [INFO     ] 2017-04-14 17:12:58,954 Parameters: current_time              = 175.81671540261
yt : [INFO     ] 2017-04-14 17:12:58,955 Parameters: domain_dimensions         = [32 32 32]
yt : [INFO     ] 2017-04-14 17:12:58,958 Parameters: domain_left_edge          = [ 0.  0.  0.]
yt : [INFO     ] 2017-04-14 17:12:58,959 Parameters: domain_right_edge         = [ 1.  1.  1.]
yt : [INFO     ] 2017-04-14 17:12:58,961 Parameters: cosmological_simulation   = 1
yt : [INFO     ] 2017-04-14 17:12:58,962 Parameters: current_redshift          = 0.28081685085325
yt : [INFO     ] 2017-04-14 17:12:58,964 Parameters: omega_lambda              = 0.732
yt : [INFO     ] 2017-04-14 17:12:58,965 Parameters: omega_matter              = 0.268
yt : [INFO     ] 2017-04-14 17:12:58,967 Parameters: hubble_constant           = 0.704
Parsing Hierarchy : 100%|██████████| 127/127 [00:00<00:00, 24336.47it/s]
yt : [INFO     ] 2017-04-14 17:12:58,983 Gathering a field list (this may take a moment.)
yt : [INFO     ] 2017-04-14 17:13:00,236 Parameters: current_time              = 180.81669075462
yt : [INFO     ] 2017-04-14 17:13:00,237 Parameters: domain_dimensions         = [32 32 32]
yt : [INFO     ] 2017-04-14 17:13:00,239 Parameters: domain_left_edge          = [ 0.  0.  0.]
yt : [INFO     ] 2017-04-14 17:13:00,241 Parameters: domain_right_edge         = [ 1.  1.  1.]
yt : [INFO     ] 2017-04-14 17:13:00,242 Parameters: cosmological_simulation   = 1
yt : [INFO     ] 2017-04-14 17:13:00,244 Parameters: current_redshift          = 0.24988435558784
yt : [INFO     ] 2017-04-14 17:13:00,245 Parameters: omega_lambda              = 0.732
yt : [INFO     ] 2017-04-14 17:13:00,246 Parameters: omega_matter              = 0.268
yt : [INFO     ] 2017-04-14 17:13:00,248 Parameters: hubble_constant           = 0.704
Parsing Hierarchy : 100%|██████████| 132/132 [00:00<00:00, 28224.31it/s]
yt : [INFO     ] 2017-04-14 17:13:00,267 Gathering a field list (this may take a moment.)
yt : [INFO     ] 2017-04-14 17:13:01,483 Parameters: current_time              = 185.81668581456
yt : [INFO     ] 2017-04-14 17:13:01,484 Parameters: domain_dimensions         = [32 32 32]
yt : [INFO     ] 2017-04-14 17:13:01,486 Parameters: domain_left_edge          = [ 0.  0.  0.]
yt : [INFO     ] 2017-04-14 17:13:01,488 Parameters: domain_right_edge         = [ 1.  1.  1.]
yt : [INFO     ] 2017-04-14 17:13:01,489 Parameters: cosmological_simulation   = 1
yt : [INFO     ] 2017-04-14 17:13:01,490 Parameters: current_redshift          = 0.22014792325267
yt : [INFO     ] 2017-04-14 17:13:01,492 Parameters: omega_lambda              = 0.732
yt : [INFO     ] 2017-04-14 17:13:01,493 Parameters: omega_matter              = 0.268
yt : [INFO     ] 2017-04-14 17:13:01,494 Parameters: hubble_constant           = 0.704
Parsing Hierarchy : 100%|██████████| 127/127 [00:00<00:00, 25341.42it/s]
yt : [INFO     ] 2017-04-14 17:13:01,511 Gathering a field list (this may take a moment.)
yt : [INFO     ] 2017-04-14 17:13:02,966 Parameters: current_time              = 190.81666319251
yt : [INFO     ] 2017-04-14 17:13:02,967 Parameters: domain_dimensions         = [32 32 32]
yt : [INFO     ] 2017-04-14 17:13:02,970 Parameters: domain_left_edge          = [ 0.  0.  0.]
yt : [INFO     ] 2017-04-14 17:13:02,971 Parameters: domain_right_edge         = [ 1.  1.  1.]
yt : [INFO     ] 2017-04-14 17:13:02,972 Parameters: cosmological_simulation   = 1
yt : [INFO     ] 2017-04-14 17:13:02,974 Parameters: current_redshift          = 0.19152709865797
yt : [INFO     ] 2017-04-14 17:13:02,975 Parameters: omega_lambda              = 0.732
yt : [INFO     ] 2017-04-14 17:13:02,976 Parameters: omega_matter              = 0.268
yt : [INFO     ] 2017-04-14 17:13:02,978 Parameters: hubble_constant           = 0.704
Parsing Hierarchy : 100%|██████████| 132/132 [00:00<00:00, 22684.92it/s]
yt : [INFO     ] 2017-04-14 17:13:02,996 Gathering a field list (this may take a moment.)
yt : [INFO     ] 2017-04-14 17:13:04,177 Parameters: current_time              = 195.81664708888
yt : [INFO     ] 2017-04-14 17:13:04,178 Parameters: domain_dimensions         = [32 32 32]
yt : [INFO     ] 2017-04-14 17:13:04,180 Parameters: domain_left_edge          = [ 0.  0.  0.]
yt : [INFO     ] 2017-04-14 17:13:04,182 Parameters: domain_right_edge         = [ 1.  1.  1.]
yt : [INFO     ] 2017-04-14 17:13:04,184 Parameters: cosmological_simulation   = 1
yt : [INFO     ] 2017-04-14 17:13:04,185 Parameters: current_redshift          = 0.16394878625904
yt : [INFO     ] 2017-04-14 17:13:04,187 Parameters: omega_lambda              = 0.732
yt : [INFO     ] 2017-04-14 17:13:04,188 Parameters: omega_matter              = 0.268
yt : [INFO     ] 2017-04-14 17:13:04,189 Parameters: hubble_constant           = 0.704
Parsing Hierarchy : 100%|██████████| 134/134 [00:00<00:00, 26866.00it/s]
yt : [INFO     ] 2017-04-14 17:13:04,208 Gathering a field list (this may take a moment.)
yt : [INFO     ] 2017-04-14 17:13:05,400 Parameters: current_time              = 200.81662398535
yt : [INFO     ] 2017-04-14 17:13:05,402 Parameters: domain_dimensions         = [32 32 32]
yt : [INFO     ] 2017-04-14 17:13:05,404 Parameters: domain_left_edge          = [ 0.  0.  0.]
yt : [INFO     ] 2017-04-14 17:13:05,405 Parameters: domain_right_edge         = [ 1.  1.  1.]
yt : [INFO     ] 2017-04-14 17:13:05,406 Parameters: cosmological_simulation   = 1
yt : [INFO     ] 2017-04-14 17:13:05,408 Parameters: current_redshift          = 0.13734692208715
yt : [INFO     ] 2017-04-14 17:13:05,409 Parameters: omega_lambda              = 0.732
yt : [INFO     ] 2017-04-14 17:13:05,410 Parameters: omega_matter              = 0.268
yt : [INFO     ] 2017-04-14 17:13:05,411 Parameters: hubble_constant           = 0.704
Parsing Hierarchy : 100%|██████████| 134/134 [00:00<00:00, 28404.34it/s]
yt : [INFO     ] 2017-04-14 17:13:05,426 Gathering a field list (this may take a moment.)
yt : [INFO     ] 2017-04-14 17:13:06,825 Parameters: current_time              = 205.8165768785
yt : [INFO     ] 2017-04-14 17:13:06,827 Parameters: domain_dimensions         = [32 32 32]
yt : [INFO     ] 2017-04-14 17:13:06,828 Parameters: domain_left_edge          = [ 0.  0.  0.]
yt : [INFO     ] 2017-04-14 17:13:06,830 Parameters: domain_right_edge         = [ 1.  1.  1.]
yt : [INFO     ] 2017-04-14 17:13:06,831 Parameters: cosmological_simulation   = 1
yt : [INFO     ] 2017-04-14 17:13:06,832 Parameters: current_redshift          = 0.11166148005496
yt : [INFO     ] 2017-04-14 17:13:06,834 Parameters: omega_lambda              = 0.732
yt : [INFO     ] 2017-04-14 17:13:06,835 Parameters: omega_matter              = 0.268
yt : [INFO     ] 2017-04-14 17:13:06,836 Parameters: hubble_constant           = 0.704
Parsing Hierarchy : 100%|██████████| 129/129 [00:00<00:00, 18946.85it/s]
yt : [INFO     ] 2017-04-14 17:13:06,854 Gathering a field list (this may take a moment.)
yt : [INFO     ] 2017-04-14 17:13:08,077 Parameters: current_time              = 210.81650455891
yt : [INFO     ] 2017-04-14 17:13:08,078 Parameters: domain_dimensions         = [32 32 32]
yt : [INFO     ] 2017-04-14 17:13:08,080 Parameters: domain_left_edge          = [ 0.  0.  0.]
yt : [INFO     ] 2017-04-14 17:13:08,082 Parameters: domain_right_edge         = [ 1.  1.  1.]
yt : [INFO     ] 2017-04-14 17:13:08,083 Parameters: cosmological_simulation   = 1
yt : [INFO     ] 2017-04-14 17:13:08,085 Parameters: current_redshift          = 0.086837698033633
yt : [INFO     ] 2017-04-14 17:13:08,086 Parameters: omega_lambda              = 0.732
yt : [INFO     ] 2017-04-14 17:13:08,087 Parameters: omega_matter              = 0.268
yt : [INFO     ] 2017-04-14 17:13:08,089 Parameters: hubble_constant           = 0.704
Parsing Hierarchy : 100%|██████████| 131/131 [00:00<00:00, 26423.67it/s]
yt : [INFO     ] 2017-04-14 17:13:08,105 Gathering a field list (this may take a moment.)
yt : [INFO     ] 2017-04-14 17:13:09,309 Parameters: current_time              = 215.81654844752
yt : [INFO     ] 2017-04-14 17:13:09,310 Parameters: domain_dimensions         = [32 32 32]
yt : [INFO     ] 2017-04-14 17:13:09,313 Parameters: domain_left_edge          = [ 0.  0.  0.]
yt : [INFO     ] 2017-04-14 17:13:09,314 Parameters: domain_right_edge         = [ 1.  1.  1.]
yt : [INFO     ] 2017-04-14 17:13:09,316 Parameters: cosmological_simulation   = 1
yt : [INFO     ] 2017-04-14 17:13:09,318 Parameters: current_redshift          = 0.062824903812594
yt : [INFO     ] 2017-04-14 17:13:09,319 Parameters: omega_lambda              = 0.732
yt : [INFO     ] 2017-04-14 17:13:09,320 Parameters: omega_matter              = 0.268
yt : [INFO     ] 2017-04-14 17:13:09,322 Parameters: hubble_constant           = 0.704
Parsing Hierarchy : 100%|██████████| 122/122 [00:00<00:00, 26227.84it/s]
yt : [INFO     ] 2017-04-14 17:13:09,338 Gathering a field list (this may take a moment.)
yt : [INFO     ] 2017-04-14 17:13:10,774 Parameters: current_time              = 220.81662170548
yt : [INFO     ] 2017-04-14 17:13:10,775 Parameters: domain_dimensions         = [32 32 32]
yt : [INFO     ] 2017-04-14 17:13:10,778 Parameters: domain_left_edge          = [ 0.  0.  0.]
yt : [INFO     ] 2017-04-14 17:13:10,780 Parameters: domain_right_edge         = [ 1.  1.  1.]
yt : [INFO     ] 2017-04-14 17:13:10,781 Parameters: cosmological_simulation   = 1
yt : [INFO     ] 2017-04-14 17:13:10,783 Parameters: current_redshift          = 0.039577805405602
yt : [INFO     ] 2017-04-14 17:13:10,784 Parameters: omega_lambda              = 0.732
yt : [INFO     ] 2017-04-14 17:13:10,785 Parameters: omega_matter              = 0.268
yt : [INFO     ] 2017-04-14 17:13:10,787 Parameters: hubble_constant           = 0.704
Parsing Hierarchy : 100%|██████████| 124/124 [00:00<00:00, 18801.06it/s]
yt : [INFO     ] 2017-04-14 17:13:10,807 Gathering a field list (this may take a moment.)
yt : [INFO     ] 2017-04-14 17:13:12,037 Parameters: current_time              = 225.81668302743
yt : [INFO     ] 2017-04-14 17:13:12,039 Parameters: domain_dimensions         = [32 32 32]
yt : [INFO     ] 2017-04-14 17:13:12,041 Parameters: domain_left_edge          = [ 0.  0.  0.]
yt : [INFO     ] 2017-04-14 17:13:12,042 Parameters: domain_right_edge         = [ 1.  1.  1.]
yt : [INFO     ] 2017-04-14 17:13:12,044 Parameters: cosmological_simulation   = 1
yt : [INFO     ] 2017-04-14 17:13:12,045 Parameters: current_redshift          = 0.017054656507229
yt : [INFO     ] 2017-04-14 17:13:12,046 Parameters: omega_lambda              = 0.732
yt : [INFO     ] 2017-04-14 17:13:12,048 Parameters: omega_matter              = 0.268
yt : [INFO     ] 2017-04-14 17:13:12,049 Parameters: hubble_constant           = 0.704
Parsing Hierarchy : 100%|██████████| 126/126 [00:00<00:00, 25288.65it/s]
yt : [INFO     ] 2017-04-14 17:13:12,065 Gathering a field list (this may take a moment.)
yt : [INFO     ] 2017-04-14 17:13:13,281 Parameters: current_time              = 229.70846991283
yt : [INFO     ] 2017-04-14 17:13:13,282 Parameters: domain_dimensions         = [32 32 32]
yt : [INFO     ] 2017-04-14 17:13:13,285 Parameters: domain_left_edge          = [ 0.  0.  0.]
yt : [INFO     ] 2017-04-14 17:13:13,287 Parameters: domain_right_edge         = [ 1.  1.  1.]
yt : [INFO     ] 2017-04-14 17:13:13,288 Parameters: cosmological_simulation   = 1
yt : [INFO     ] 2017-04-14 17:13:13,289 Parameters: current_redshift          = 2.2204460492503e-16
yt : [INFO     ] 2017-04-14 17:13:13,291 Parameters: omega_lambda              = 0.732
yt : [INFO     ] 2017-04-14 17:13:13,292 Parameters: omega_matter              = 0.268
yt : [INFO     ] 2017-04-14 17:13:13,293 Parameters: hubble_constant           = 0.704
Parsing Hierarchy : 100%|██████████| 124/124 [00:00<00:00, 25578.80it/s]
yt : [INFO     ] 2017-04-14 17:13:13,309 Gathering a field list (this may take a moment.)
yt : [INFO     ] 2017-04-14 17:13:14,753 Parameters: current_time              = 208.14467318643
yt : [INFO     ] 2017-04-14 17:13:14,754 Parameters: domain_dimensions         = [32 32 32]
yt : [INFO     ] 2017-04-14 17:13:14,756 Parameters: domain_left_edge          = [ 0.  0.  0.]
yt : [INFO     ] 2017-04-14 17:13:14,758 Parameters: domain_right_edge         = [ 1.  1.  1.]
yt : [INFO     ] 2017-04-14 17:13:14,760 Parameters: cosmological_simulation   = 1
yt : [INFO     ] 2017-04-14 17:13:14,761 Parameters: current_redshift          = 0.099998843013062
yt : [INFO     ] 2017-04-14 17:13:14,762 Parameters: omega_lambda              = 0.732
yt : [INFO     ] 2017-04-14 17:13:14,764 Parameters: omega_matter              = 0.268
yt : [INFO     ] 2017-04-14 17:13:14,765 Parameters: hubble_constant           = 0.704
Parsing Hierarchy : 100%|██████████| 133/133 [00:00<00:00, 16984.09it/s]
yt : [INFO     ] 2017-04-14 17:13:14,783 Gathering a field list (this may take a moment.)
yt : [INFO     ] 2017-04-14 17:13:15,961 Parameters: current_time              = 210.37447283232
yt : [INFO     ] 2017-04-14 17:13:15,963 Parameters: domain_dimensions         = [32 32 32]
yt : [INFO     ] 2017-04-14 17:13:15,966 Parameters: domain_left_edge          = [ 0.  0.  0.]
yt : [INFO     ] 2017-04-14 17:13:15,968 Parameters: domain_right_edge         = [ 1.  1.  1.]
yt : [INFO     ] 2017-04-14 17:13:15,970 Parameters: cosmological_simulation   = 1
yt : [INFO     ] 2017-04-14 17:13:15,972 Parameters: current_redshift          = 0.088998908295534
yt : [INFO     ] 2017-04-14 17:13:15,973 Parameters: omega_lambda              = 0.732
yt : [INFO     ] 2017-04-14 17:13:15,975 Parameters: omega_matter              = 0.268
yt : [INFO     ] 2017-04-14 17:13:15,977 Parameters: hubble_constant           = 0.704
Parsing Hierarchy : 100%|██████████| 134/134 [00:00<00:00, 16506.22it/s]
yt : [INFO     ] 2017-04-14 17:13:15,996 Gathering a field list (this may take a moment.)
yt : [INFO     ] 2017-04-14 17:13:17,185 Parameters: current_time              = 212.63786229367
yt : [INFO     ] 2017-04-14 17:13:17,186 Parameters: domain_dimensions         = [32 32 32]
yt : [INFO     ] 2017-04-14 17:13:17,189 Parameters: domain_left_edge          = [ 0.  0.  0.]
yt : [INFO     ] 2017-04-14 17:13:17,190 Parameters: domain_right_edge         = [ 1.  1.  1.]
yt : [INFO     ] 2017-04-14 17:13:17,191 Parameters: cosmological_simulation   = 1
yt : [INFO     ] 2017-04-14 17:13:17,192 Parameters: current_redshift          = 0.077999121535763
yt : [INFO     ] 2017-04-14 17:13:17,194 Parameters: omega_lambda              = 0.732
yt : [INFO     ] 2017-04-14 17:13:17,196 Parameters: omega_matter              = 0.268
yt : [INFO     ] 2017-04-14 17:13:17,197 Parameters: hubble_constant           = 0.704
Parsing Hierarchy : 100%|██████████| 136/136 [00:00<00:00, 26891.63it/s]
yt : [INFO     ] 2017-04-14 17:13:17,213 Gathering a field list (this may take a moment.)
yt : [INFO     ] 2017-04-14 17:13:18,711 Parameters: current_time              = 214.93558667652
yt : [INFO     ] 2017-04-14 17:13:18,712 Parameters: domain_dimensions         = [32 32 32]
yt : [INFO     ] 2017-04-14 17:13:18,714 Parameters: domain_left_edge          = [ 0.  0.  0.]
yt : [INFO     ] 2017-04-14 17:13:18,716 Parameters: domain_right_edge         = [ 1.  1.  1.]
yt : [INFO     ] 2017-04-14 17:13:18,717 Parameters: cosmological_simulation   = 1
yt : [INFO     ] 2017-04-14 17:13:18,719 Parameters: current_redshift          = 0.066998908324564
yt : [INFO     ] 2017-04-14 17:13:18,720 Parameters: omega_lambda              = 0.732
yt : [INFO     ] 2017-04-14 17:13:18,721 Parameters: omega_matter              = 0.268
yt : [INFO     ] 2017-04-14 17:13:18,723 Parameters: hubble_constant           = 0.704
Parsing Hierarchy : 100%|██████████| 125/125 [00:00<00:00, 24980.37it/s]
yt : [INFO     ] 2017-04-14 17:13:18,743 Gathering a field list (this may take a moment.)
yt : [INFO     ] 2017-04-14 17:13:19,945 Parameters: current_time              = 217.05453517381
yt : [INFO     ] 2017-04-14 17:13:19,946 Parameters: domain_dimensions         = [32 32 32]
yt : [INFO     ] 2017-04-14 17:13:19,948 Parameters: domain_left_edge          = [ 0.  0.  0.]
yt : [INFO     ] 2017-04-14 17:13:19,949 Parameters: domain_right_edge         = [ 1.  1.  1.]
yt : [INFO     ] 2017-04-14 17:13:19,951 Parameters: cosmological_simulation   = 1
yt : [INFO     ] 2017-04-14 17:13:19,952 Parameters: current_redshift          = 0.056999419833864
yt : [INFO     ] 2017-04-14 17:13:19,954 Parameters: omega_lambda              = 0.732
yt : [INFO     ] 2017-04-14 17:13:19,956 Parameters: omega_matter              = 0.268
yt : [INFO     ] 2017-04-14 17:13:19,957 Parameters: hubble_constant           = 0.704
Parsing Hierarchy : 100%|██████████| 124/124 [00:00<00:00, 26293.92it/s]
yt : [INFO     ] 2017-04-14 17:13:19,974 Gathering a field list (this may take a moment.)
yt : [INFO     ] 2017-04-14 17:13:21,406 Parameters: current_time              = 219.20300709899
yt : [INFO     ] 2017-04-14 17:13:21,407 Parameters: domain_dimensions         = [32 32 32]
yt : [INFO     ] 2017-04-14 17:13:21,409 Parameters: domain_left_edge          = [ 0.  0.  0.]
yt : [INFO     ] 2017-04-14 17:13:21,410 Parameters: domain_right_edge         = [ 1.  1.  1.]
yt : [INFO     ] 2017-04-14 17:13:21,412 Parameters: cosmological_simulation   = 1
yt : [INFO     ] 2017-04-14 17:13:21,413 Parameters: current_redshift          = 0.046999006913759
yt : [INFO     ] 2017-04-14 17:13:21,415 Parameters: omega_lambda              = 0.732
yt : [INFO     ] 2017-04-14 17:13:21,416 Parameters: omega_matter              = 0.268
yt : [INFO     ] 2017-04-14 17:13:21,418 Parameters: hubble_constant           = 0.704
Parsing Hierarchy : 100%|██████████| 135/135 [00:00<00:00, 25752.99it/s]
yt : [INFO     ] 2017-04-14 17:13:21,434 Gathering a field list (this may take a moment.)
yt : [INFO     ] 2017-04-14 17:13:22,660 Parameters: current_time              = 221.38106136178
yt : [INFO     ] 2017-04-14 17:13:22,661 Parameters: domain_dimensions         = [32 32 32]
yt : [INFO     ] 2017-04-14 17:13:22,663 Parameters: domain_left_edge          = [ 0.  0.  0.]
yt : [INFO     ] 2017-04-14 17:13:22,664 Parameters: domain_right_edge         = [ 1.  1.  1.]
yt : [INFO     ] 2017-04-14 17:13:22,666 Parameters: cosmological_simulation   = 1
yt : [INFO     ] 2017-04-14 17:13:22,668 Parameters: current_redshift          = 0.036999742678557
yt : [INFO     ] 2017-04-14 17:13:22,669 Parameters: omega_lambda              = 0.732
yt : [INFO     ] 2017-04-14 17:13:22,670 Parameters: omega_matter              = 0.268
yt : [INFO     ] 2017-04-14 17:13:22,671 Parameters: hubble_constant           = 0.704
Parsing Hierarchy : 100%|██████████| 126/126 [00:00<00:00, 25217.46it/s]
yt : [INFO     ] 2017-04-14 17:13:22,687 Gathering a field list (this may take a moment.)
yt : [INFO     ] 2017-04-14 17:13:23,902 Parameters: current_time              = 223.58978764375
yt : [INFO     ] 2017-04-14 17:13:23,903 Parameters: domain_dimensions         = [32 32 32]
yt : [INFO     ] 2017-04-14 17:13:23,906 Parameters: domain_left_edge          = [ 0.  0.  0.]
yt : [INFO     ] 2017-04-14 17:13:23,908 Parameters: domain_right_edge         = [ 1.  1.  1.]
yt : [INFO     ] 2017-04-14 17:13:23,910 Parameters: cosmological_simulation   = 1
yt : [INFO     ] 2017-04-14 17:13:23,911 Parameters: current_redshift          = 0.026999006928412
yt : [INFO     ] 2017-04-14 17:13:23,912 Parameters: omega_lambda              = 0.732
yt : [INFO     ] 2017-04-14 17:13:23,913 Parameters: omega_matter              = 0.268
yt : [INFO     ] 2017-04-14 17:13:23,915 Parameters: hubble_constant           = 0.704
Parsing Hierarchy : 100%|██████████| 132/132 [00:00<00:00, 18057.08it/s]
yt : [INFO     ] 2017-04-14 17:13:23,933 Gathering a field list (this may take a moment.)
yt : [INFO     ] 2017-04-14 17:13:25,337 Parameters: current_time              = 225.82900916478
yt : [INFO     ] 2017-04-14 17:13:25,338 Parameters: domain_dimensions         = [32 32 32]
yt : [INFO     ] 2017-04-14 17:13:25,340 Parameters: domain_left_edge          = [ 0.  0.  0.]
yt : [INFO     ] 2017-04-14 17:13:25,341 Parameters: domain_right_edge         = [ 1.  1.  1.]
yt : [INFO     ] 2017-04-14 17:13:25,342 Parameters: cosmological_simulation   = 1
yt : [INFO     ] 2017-04-14 17:13:25,344 Parameters: current_redshift          = 0.016999994534557
yt : [INFO     ] 2017-04-14 17:13:25,345 Parameters: omega_lambda              = 0.732
yt : [INFO     ] 2017-04-14 17:13:25,346 Parameters: omega_matter              = 0.268
yt : [INFO     ] 2017-04-14 17:13:25,348 Parameters: hubble_constant           = 0.704
Parsing Hierarchy : 100%|██████████| 126/126 [00:00<00:00, 20295.80it/s]
yt : [INFO     ] 2017-04-14 17:13:25,367 Gathering a field list (this may take a moment.)
yt : [INFO     ] 2017-04-14 17:13:26,598 Parameters: current_time              = 228.1000985764
yt : [INFO     ] 2017-04-14 17:13:26,599 Parameters: domain_dimensions         = [32 32 32]
yt : [INFO     ] 2017-04-14 17:13:26,602 Parameters: domain_left_edge          = [ 0.  0.  0.]
yt : [INFO     ] 2017-04-14 17:13:26,603 Parameters: domain_right_edge         = [ 1.  1.  1.]
yt : [INFO     ] 2017-04-14 17:13:26,604 Parameters: cosmological_simulation   = 1
yt : [INFO     ] 2017-04-14 17:13:26,605 Parameters: current_redshift          = 0.0069990069523933
yt : [INFO     ] 2017-04-14 17:13:26,607 Parameters: omega_lambda              = 0.732
yt : [INFO     ] 2017-04-14 17:13:26,613 Parameters: omega_matter              = 0.268
yt : [INFO     ] 2017-04-14 17:13:26,614 Parameters: hubble_constant           = 0.704
Parsing Hierarchy : 100%|██████████| 135/135 [00:00<00:00, 27102.77it/s]
yt : [INFO     ] 2017-04-14 17:13:26,631 Gathering a field list (this may take a moment.)

In [73]:
print(rho)


[7.182167246767479e-26 g/cm**3, 2.9059611911967146e-27 g/cm**3, 1.3005399720818743e-27 g/cm**3, 8.665855780016177e-28 g/cm**3, 6.891936512431531e-28 g/cm**3, 5.909674099074075e-28 g/cm**3, 6.437978652578952e-28 g/cm**3, 7.781435079625696e-28 g/cm**3, 9.521005081655007e-28 g/cm**3, 1.2057546479919688e-27 g/cm**3, 2.0161676691464046e-27 g/cm**3, 3.380572129927601e-27 g/cm**3, 3.295075322999272e-27 g/cm**3, 9.666878720323255e-27 g/cm**3, 1.825961802282628e-26 g/cm**3, 1.1726898858891321e-26 g/cm**3, 6.813775120263505e-27 g/cm**3, 1.4230549555324727e-26 g/cm**3, 1.078782931550446e-26 g/cm**3, 6.745072238857048e-27 g/cm**3, 6.808817705910852e-27 g/cm**3, 7.561628519554906e-27 g/cm**3, 8.434088433480284e-27 g/cm**3, 8.517243269487163e-27 g/cm**3, 1.1286984884893572e-26 g/cm**3, 1.3194951024096913e-26 g/cm**3, 9.119251590395433e-27 g/cm**3, 9.965776408071898e-27 g/cm**3, 8.333596630921615e-27 g/cm**3, 8.151365986808684e-27 g/cm**3, 7.787286493148639e-27 g/cm**3, 8.300723955850137e-27 g/cm**3, 7.801276432330143e-27 g/cm**3, 9.373997930365866e-27 g/cm**3, 1.515729240904511e-26 g/cm**3, 1.1629496442549857e-26 g/cm**3, 1.6189246512340015e-26 g/cm**3, 1.894826495302461e-26 g/cm**3, 1.6602405287644688e-26 g/cm**3, 9.617606745970016e-27 g/cm**3, 7.144667095938952e-27 g/cm**3, 1.3164275978638868e-26 g/cm**3, 8.037914530585359e-27 g/cm**3, 1.331117721827776e-26 g/cm**3, 8.20353531245526e-27 g/cm**3, 1.192355294538153e-26 g/cm**3, 1.1263226910520681e-26 g/cm**3, 9.65861944978284e-27 g/cm**3, 8.694137708035828e-27 g/cm**3, 1.0150175873599733e-26 g/cm**3, 1.4849862845110662e-26 g/cm**3, 1.2333428016349913e-26 g/cm**3, 1.1793796248186496e-26 g/cm**3, 7.578845626351999e-27 g/cm**3, 6.3709097451064186e-27 g/cm**3, 1.1923266136800247e-26 g/cm**3, 1.0779097150543768e-26 g/cm**3]

In [89]:
plt.style.use("ggplot")
plt.plot(t, rho, '-x')
plt.xscale("log")
plt.yscale("log")



In [90]:
import h5py

In [94]:
f = h5py.File("/srv/nbgrader/data/zebrafish/ViBE-Z_96hpf_v1.h5", "r")
list(f.keys())


Out[94]:
['anatomy']

In [95]:
list(f["/anatomy"].keys())


Out[95]:
['average_brain', 'reference_larva', 'segment_lines', 'segment_regions']

In [100]:
data = {}
for k in f["/anatomy"]:
    data[k] = f["/anatomy"][k][:].astype("f8")
    print(data[k].shape)


(500, 500, 1000)
(500, 500, 1000)
(500, 500, 1000)
(500, 500, 1000)

In [104]:
ds = yt.load_uniform_grid(data, [500, 500, 1000], nprocs=20)


yt : [INFO     ] 2017-04-14 17:25:19,522 Parameters: current_time              = 0.0
yt : [INFO     ] 2017-04-14 17:25:19,524 Parameters: domain_dimensions         = [ 500  500 1000]
yt : [INFO     ] 2017-04-14 17:25:19,525 Parameters: domain_left_edge          = [ 0.  0.  0.]
yt : [INFO     ] 2017-04-14 17:25:19,527 Parameters: domain_right_edge         = [ 1.  1.  1.]
yt : [INFO     ] 2017-04-14 17:25:19,529 Parameters: cosmological_simulation   = 0.0

In [105]:
ds


Out[105]:
UniformGridData

In [106]:
p = yt.ProjectionPlot(ds, "z", "average_brain", method="mip")
p.show()


yt : [INFO     ] 2017-04-14 17:26:10,428 Projection completed
yt : [INFO     ] 2017-04-14 17:26:10,430 xlim = 0.000000 1.000000
yt : [INFO     ] 2017-04-14 17:26:10,431 ylim = 0.000000 1.000000
yt : [INFO     ] 2017-04-14 17:26:10,433 xlim = 0.000000 1.000000
yt : [INFO     ] 2017-04-14 17:26:10,434 ylim = 0.000000 1.000000
yt : [INFO     ] 2017-04-14 17:26:10,436 Making a fixed resolution buffer of (('stream', 'average_brain')) 800 by 800


In [107]:
p.set_log("all", False)


Out[107]:


In [108]:
ds.field_list


Out[108]:
[('stream', 'average_brain'),
 ('stream', 'reference_larva'),
 ('stream', 'segment_lines'),
 ('stream', 'segment_regions')]

In [110]:
p2 = yt.ProjectionPlot(ds, "z", "reference_larva", method="mip")


yt : [INFO     ] 2017-04-14 17:27:56,299 Projection completed
yt : [INFO     ] 2017-04-14 17:27:56,301 xlim = 0.000000 1.000000
yt : [INFO     ] 2017-04-14 17:27:56,302 ylim = 0.000000 1.000000
yt : [INFO     ] 2017-04-14 17:27:56,304 xlim = 0.000000 1.000000
yt : [INFO     ] 2017-04-14 17:27:56,305 ylim = 0.000000 1.000000
yt : [INFO     ] 2017-04-14 17:27:56,307 Making a fixed resolution buffer of (('stream', 'reference_larva')) 800 by 800

In [111]:
p2.set_log("all", False)


Out[111]:


In [113]:
phase = yt.PhasePlot(ds.r[:,:,:], "average_brain", "reference_larva", "cell_volume", weight_field=None)
phase.show()


---------------------------------------------------------------------------
MemoryError                               Traceback (most recent call last)
<ipython-input-113-9db7637d0a52> in <module>()
----> 1 phase = yt.PhasePlot(ds.r[:,:,:], "average_brain", "reference_larva", "cell_volume", weight_field=None)
      2 phase.show()

/opt/conda/lib/python3.5/site-packages/yt/visualization/profile_plotter.py in __init__(self, data_source, x_field, y_field, z_fields, weight_field, x_bins, y_bins, accumulation, fractional, fontsize, figure_size)
    748                 weight_field=weight_field,
    749                 accumulation=accumulation,
--> 750                 fractional=fractional)
    751 
    752         type(self)._initialize_instance(self, data_source, profile, fontsize,

/opt/conda/lib/python3.5/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)
   1040         setattr(obj, "fractional", fractional)
   1041     if fields is not None:
-> 1042         obj.add_fields([field for field in fields])
   1043     for field in fields:
   1044         if fractional:

/opt/conda/lib/python3.5/site-packages/yt/data_objects/profiles.py in add_fields(self, fields)
    113         citer = self.data_source.chunks([], "io")
    114         for chunk in parallel_objects(citer):
--> 115             self._bin_chunk(chunk, fields, temp_storage)
    116         self._finalize_storage(fields, temp_storage)
    117 

/opt/conda/lib/python3.5/site-packages/yt/data_objects/profiles.py in _bin_chunk(self, chunk, fields, storage)
    543 
    544     def _bin_chunk(self, chunk, fields, storage):
--> 545         rv = self._get_data(chunk, fields)
    546         if rv is None: return
    547         fdata, wdata, (bf_x, bf_y) = rv

/opt/conda/lib/python3.5/site-packages/yt/data_objects/profiles.py in _get_data(self, chunk, fields)
    241         for i, field in enumerate(fields):
    242             units = chunk.ds.field_info[field].output_units
--> 243             arr[:,i] = chunk[field][filter].in_units(units)
    244         if self.weight_field is not None:
    245             units = chunk.ds.field_info[self.weight_field].output_units

/opt/conda/lib/python3.5/site-packages/yt/data_objects/data_containers.py in __getitem__(self, key)
    270                 return self.field_data[f]
    271             else:
--> 272                 self.get_data(f)
    273         # fi.units is the unit expression string. We depend on the registry
    274         # hanging off the dataset to define this unit object.

/opt/conda/lib/python3.5/site-packages/yt/data_objects/data_containers.py in get_data(self, fields)
   1190 
   1191         fields_to_generate += gen_fluids + gen_particles
-> 1192         self._generate_fields(fields_to_generate)
   1193         for field in list(self.field_data.keys()):
   1194             if field not in ofields:

/opt/conda/lib/python3.5/site-packages/yt/data_objects/data_containers.py in _generate_fields(self, fields_to_generate)
   1210                 fi = self.ds._get_field_info(*field)
   1211                 try:
-> 1212                     fd = self._generate_field(field)
   1213                     if fd is None:
   1214                         raise RuntimeError

/opt/conda/lib/python3.5/site-packages/yt/data_objects/data_containers.py in _generate_field(self, field)
    307                 tr = self._generate_particle_field(field)
    308             else:
--> 309                 tr = self._generate_fluid_field(field)
    310             if tr is None:
    311                 raise YTCouldNotGenerateField(field, self.ds)

/opt/conda/lib/python3.5/site-packages/yt/data_objects/data_containers.py in _generate_fluid_field(self, field)
    327             rv = self._generate_spatial_fluid(field, ngt_exception.ghost_zones)
    328         else:
--> 329             rv = finfo(gen_obj)
    330         return rv
    331 

/opt/conda/lib/python3.5/site-packages/yt/fields/derived_field.py in __call__(self, data)
    202                 "for %s" % (self.name,))
    203         with self.unit_registry(data):
--> 204             dd = self._function(self, data)
    205         for field_name in data.keys():
    206             if field_name not in original_fields:

/opt/conda/lib/python3.5/site-packages/yt/geometry/coordinates/cartesian_coordinates.py in _cell_volume(field, data)
     52                                units = "code_length")
     53         def _cell_volume(field, data):
---> 54             rv  = data["index", "dx"].copy(order='K')
     55             rv *= data["index", "dy"]
     56             rv *= data["index", "dz"]

MemoryError: 

In [ ]: