I'm trying to load in two Gizmo cosmo simulations that are exactly identical except that one has a 25 Mpc/h size (call it sim) and the other has 100 Mpc/h (sim2). Weirdly, the 100 Mpc/h loads fine, but the 25 Mpc/h doesn't load any of its particle data, it seems. Below are some outputs from an ipython session to show what I mean.
Any thoughts/pointers on what could be going wrong with the 25 box would be much appreciated.
thanks,
romeel
In [4]: sim2 = yt.load("/cosma/home/dc-dave2/data/m25n1024/s43/snap_m25n1024_042.hdf5")
yt : [INFO ] 2018-09-11 19:22:52,554 Calculating time from 1.660e-01 to be 3.729e+16 seconds
yt : [INFO ] 2018-09-11 19:22:52,554 Assuming length units are in kpc/h (comoving)
yt : [INFO ] 2018-09-11 19:22:52,565 Parameters: current_time = 3.72864881334e+16 s
yt : [INFO ] 2018-09-11 19:22:52,565 Parameters: domain_dimensions = [2 2 2]
yt : [INFO ] 2018-09-11 19:22:52,565 Parameters: domain_left_edge = [ 0. 0. 0.]
yt : [INFO ] 2018-09-11 19:22:52,566 Parameters: domain_right_edge = [ 25000. 25000. 25000.]
yt : [INFO ] 2018-09-11 19:22:52,566 Parameters: cosmological_simulation = 1
yt : [INFO ] 2018-09-11 19:22:52,566 Parameters: current_redshift = 5.02438675155
yt : [INFO ] 2018-09-11 19:22:52,566 Parameters: omega_lambda = 0.7
yt : [INFO ] 2018-09-11 19:22:52,566 Parameters: omega_matter = 0.3
yt : [INFO ] 2018-09-11 19:22:52,566 Parameters: hubble_constant = 0.68
In [6]: sim2 = yt.load("/cosma/home/dc-dave2/data/m100n1024/s43/snap_m100n1024_042.hdf5")
yt : [INFO ] 2018-09-11 19:23:33,838 Calculating time from 1.660e-01 to be 3.729e+16 seconds
yt : [INFO ] 2018-09-11 19:23:33,838 Assuming length units are in kpc/h (comoving)
yt : [INFO ] 2018-09-11 19:23:33,847 Parameters: current_time = 3.72864881334e+16 s
yt : [INFO ] 2018-09-11 19:23:33,847 Parameters: domain_dimensions = [2 2 2]
yt : [INFO ] 2018-09-11 19:23:33,847 Parameters: domain_left_edge = [ 0. 0. 0.]
yt : [INFO ] 2018-09-11 19:23:33,848 Parameters: domain_right_edge = [ 100000. 100000. 100000.]
yt : [INFO ] 2018-09-11 19:23:33,848 Parameters: cosmological_simulation = 1
yt : [INFO ] 2018-09-11 19:23:33,848 Parameters: current_redshift = 5.02438675155
yt : [INFO ] 2018-09-11 19:23:33,848 Parameters: omega_lambda = 0.7
yt : [INFO ] 2018-09-11 19:23:33,848 Parameters: omega_matter = 0.3
yt : [INFO ] 2018-09-11 19:23:33,848 Parameters: hubble_constant = 0.68
In [23]: sim.dataset_type
Out[23]: 'gadget_hdf5'
In [24]: sim2.dataset_type
Out[24]: 'gadget_hdf5'
*********** all looks good so far
In [17]: sim.parameters
Out[17]:
{'BoxSize': 25000.0,
'Flag_Cooling': 1,
'Flag_DoublePrecision': 0,
'Flag_Feedback': 1,
'Flag_IC_Info': 3,
'Flag_Metals': 11,
'Flag_Sfr': 1,
'Flag_StellarAge': 1,
'HubbleParam': 0.68000000000000005,
'MassTable': array([ 0., 0., 0., 0., 0., 0.]),
'Massarr': array([ 0., 0., 0., 0., 0., 0.]),
'NumFiles': 1,
'NumFilesPerSnapshot': 1,
'NumPart_ThisFile': array([1072810110, 1073741824, 0, 0, 931728,
5], dtype=int32),
'NumPart_Total': array([1072810110, 1073741824, 0, 0, 931728,
5], dtype=uint32),
'NumPart_Total_HighWord': array([0, 0, 0, 0, 0, 0], dtype=uint32),
'Omega0': 0.29999999999999999,
'OmegaLambda': 0.69999999999999996,
'Redshift': 5.0243867515502965,
'Time': 0.16599199905993139}
In [18]: sim2.parameters
Out[18]:
{'BoxSize': 100000.0,
'Flag_Cooling': 1,
'Flag_DoublePrecision': 0,
'Flag_Feedback': 1,
'Flag_IC_Info': 3,
'Flag_Metals': 11,
'Flag_Sfr': 1,
'Flag_StellarAge': 1,
'HubbleParam': 0.68000000000000005,
'MassTable': array([ 0., 0., 0., 0., 0., 0.]),
'Massarr': array([ 0., 0., 0., 0., 0., 0.]),
'NumFiles': 1,
'NumFilesPerSnapshot': 1,
'NumPart_ThisFile': array([1071819062, 1073741824, 0, 0, 1922044,
482], dtype=int32),
'NumPart_Total': array([1071819062, 1073741824, 0, 0, 1922044,
482], dtype=uint32),
'NumPart_Total_HighWord': array([0, 0, 0, 0, 0, 0], dtype=uint32),
'Omega0': 0.29999999999999999,
'OmegaLambda': 0.69999999999999996,
'Redshift': 5.0243867515502965,
'Time': 0.16599199905993139}
************** both look reasonable. Now we look at the loaded data ....
In [31]: sim.particle_types
Out[31]: ()
In [32]: sim2.particle_types
Out[32]: ('PartType4', 'PartType5', 'PartType0', 'PartType1', 'all')
********** sim has no particle_types! Why?
In [19]: sim.field_info
Out[19]:
{('gas',
'cell_volume'): Alias Field for "('index', 'cell_volume')" (gas, cell_volume): (units: code_length**3),
('gas',
'dx'): Alias Field for "('index', 'dx')" (gas, dx): (units: code_length),
('gas',
'dy'): Alias Field for "('index', 'dy')" (gas, dy): (units: code_length),
('gas',
'dz'): Alias Field for "('index', 'dz')" (gas, dz): (units: code_length),
('gas',
'path_element_x'): Alias Field for "('index', 'path_element_x')" (gas, path_element_x): (units: code_length),
('gas',
'path_element_y'): Alias Field for "('index', 'path_element_y')" (gas, path_element_y): (units: code_length),
('gas',
'path_element_z'): Alias Field for "('index', 'path_element_z')" (gas, path_element_z): (units: code_length),
('gas',
'vertex_x'): Alias Field for "('index', 'vertex_x')" (gas, vertex_x): (units: code_length),
('gas',
'vertex_y'): Alias Field for "('index', 'vertex_y')" (gas, vertex_y): (units: code_length),
('gas',
'vertex_z'): Alias Field for "('index', 'vertex_z')" (gas, vertex_z): (units: code_length),
('gas', 'x'): Alias Field for "('index', 'x')" (gas, x): (units: code_length),
('gas', 'y'): Alias Field for "('index', 'y')" (gas, y): (units: code_length),
('gas', 'z'): Alias Field for "('index', 'z')" (gas, z): (units: code_length),
('index',
'cell_volume'): Derived Field (index, cell_volume): (units: code_length**3),
('index', 'cylindrical_r'): Derived Field (index, cylindrical_r): (units: cm),
('index',
'cylindrical_radius'): Derived Field (index, cylindrical_radius): (units: cm),
('index',
'cylindrical_theta'): Derived Field (index, cylindrical_theta): (units: ),
('index', 'cylindrical_z'): Derived Field (index, cylindrical_z): (units: cm),
('index', 'disk_angle'): Derived Field (index, disk_angle): (units: ),
('index', 'dx'): Derived Field (index, dx): (units: code_length),
('index', 'dy'): Derived Field (index, dy): (units: code_length),
('index', 'dz'): Derived Field (index, dz): (units: code_length),
('index', 'grid_indices'): Derived Field (index, grid_indices): (units: ),
('index', 'grid_level'): Derived Field (index, grid_level): (units: ),
('index', 'height'): Derived Field (index, height): (units: cm),
('index', 'morton_index'): Derived Field (index, morton_index): (units: ),
('index', 'ones'): Derived Field (index, ones): (units: ),
('index', 'ones_over_dx'): Derived Field (index, ones_over_dx): (units: 1/cm),
('index',
'path_element_x'): Derived Field (index, path_element_x): (units: code_length),
('index',
'path_element_y'): Derived Field (index, path_element_y): (units: code_length),
('index',
'path_element_z'): Derived Field (index, path_element_z): (units: code_length),
('index', 'radius'): Derived Field (index, radius): (units: cm),
('index', 'spherical_phi'): Derived Field (index, spherical_phi): (units: ),
('index', 'spherical_r'): Derived Field (index, spherical_r): (units: cm),
('index',
'spherical_radius'): Derived Field (index, spherical_radius): (units: cm),
('index',
'spherical_theta'): Derived Field (index, spherical_theta): (units: ),
('index', 'vertex_x'): Derived Field (index, vertex_x): (units: code_length),
('index', 'vertex_y'): Derived Field (index, vertex_y): (units: code_length),
('index', 'vertex_z'): Derived Field (index, vertex_z): (units: code_length),
('index',
'virial_radius_fraction'): Derived Field (index, virial_radius_fraction): (units: ),
('index', 'x'): Derived Field (index, x): (units: code_length),
('index', 'y'): Derived Field (index, y): (units: code_length),
('index', 'z'): Derived Field (index, z): (units: code_length),
('index', 'zeros'): Derived Field (index, zeros): (units: )}
In [20]: sim2.field_info
---------------------------------------------------------------------------
AttributeError Traceback (most recent call last)
<ipython-input-20-7ba53e5d35de> in <module>()
----> 1 sim2.field_info
AttributeError: 'GizmoDataset' object has no attribute 'field_info'
********* So sim2 (which is behaving properly) appears to be a GizmoDataset, as expected. But sim isn't? Then what is it?
In [21]: sim.field_list
Out[21]: []
In [22]: sim2.field_list
yt : [INFO ] 2018-09-11 19:33:26,796 Allocating for 2.147e+09 particles (index particle type 'all')
yt : [INFO ] 2018-09-11 19:39:02,984 Identified 1.516e+08 octs
Out[22]:
[('PartType0', 'AGS-Softening'),
('PartType0', 'Coordinates'),
('PartType0', 'DelayTime'),
('PartType0', 'Density'),
('PartType0', 'Dust_Masses'),
('PartType0', 'Dust_Metallicity'),
('PartType0', 'ElectronAbundance'),
('PartType0', 'FractionH2'),
('PartType0', 'GrackleHI'),
('PartType0', 'GrackleHII'),
('PartType0', 'GrackleHM'),
('PartType0', 'GrackleHeI'),
('PartType0', 'GrackleHeII'),
('PartType0', 'GrackleHeIII'),
('PartType0', 'HaloID'),
('PartType0', 'InternalEnergy'),
('PartType0', 'Masses'),
('PartType0', 'Metallicity_00'),
('PartType0', 'Metallicity_01'),
('PartType0', 'Metallicity_02'),
('PartType0', 'Metallicity_03'),
('PartType0', 'Metallicity_04'),
('PartType0', 'Metallicity_05'),
('PartType0', 'Metallicity_06'),
('PartType0', 'Metallicity_07'),
('PartType0', 'Metallicity_08'),
('PartType0', 'Metallicity_09'),
('PartType0', 'Metallicity_10'),
('PartType0', 'NWindLaunches'),
('PartType0', 'NeutralHydrogenAbundance'),
('PartType0', 'ParticleIDs'),
('PartType0', 'Potential'),
('PartType0', 'Sigma'),
('PartType0', 'SmoothingLength'),
('PartType0', 'StarFormationRate'),
('PartType0', 'Velocities'),
('PartType1', 'AGS-Softening'),
('PartType1', 'Coordinates'),
('PartType1', 'HaloID'),
('PartType1', 'Masses'),
('PartType1', 'ParticleIDs'),
('PartType1', 'Potential'),
('PartType1', 'Velocities'),
('PartType4', 'AGS-Softening'),
('PartType4', 'Coordinates'),
('PartType4', 'Dust_Masses'),
('PartType4', 'Dust_Metallicity'),
('PartType4', 'HaloID'),
('PartType4', 'Masses'),
('PartType4', 'Metallicity_00'),
('PartType4', 'Metallicity_01'),
('PartType4', 'Metallicity_02'),
('PartType4', 'Metallicity_03'),
('PartType4', 'Metallicity_04'),
('PartType4', 'Metallicity_05'),
('PartType4', 'Metallicity_06'),
('PartType4', 'Metallicity_07'),
('PartType4', 'Metallicity_08'),
('PartType4', 'Metallicity_09'),
('PartType4', 'Metallicity_10'),
('PartType4', 'ParticleIDs'),
('PartType4', 'Potential'),
('PartType4', 'StellarFormationTime'),
('PartType4', 'Velocities'),
('PartType5', 'AGS-Softening'),
('PartType5', 'BH_AccretionLength'),
('PartType5', 'BH_Mass'),
('PartType5', 'BH_Mass_AlphaDisk'),
('PartType5', 'BH_Mdot'),
('PartType5', 'BH_NProgs'),
('PartType5', 'Coordinates'),
('PartType5', 'HaloID'),
('PartType5', 'Masses'),
('PartType5', 'ParticleIDs'),
('PartType5', 'Potential'),
('PartType5', 'StellarFormationTime'),
('PartType5', 'Velocities'),
('all', 'AGS-Softening'),
('all', 'Coordinates'),
('all', 'HaloID'),
('all', 'Masses'),
('all', 'ParticleIDs'),
('all', 'Potential'),
('all', 'Velocities')]
*********** sim2 has all the right fields. sim has no fields?!?