yt for Cosmology

  • Accessing Cosmological Parameters
  • Cosmological Coordinates
  • Halo Finding
  • Halo Profiling
In [1]:
# This will make sure all plots get showed in the notebook.
%matplotlib inline  
import yt
In [13]:
ds = yt.load("../sample_data/Enzo_64/DD0043/data0043")
print ds
print ds.current_redshift
data0043
0.00139308806408

In [3]:
for k,v in ds.parameters.iteritems():
    print "%s: %s" % (k, v)
ClusterSMBHJetAngleTheta: 0.0
DataLabel[5]: GasEnergy
H2FormationOnDust: 0
HistoryDumpName: history
Coordinate: 0
H2StarMakerAssumeColdWarmPressureBalance: 0
CourantSafetyNumber: 0.5
OutputCoolingTime: 0
UnigridTranspose: 0
SinkMergeDistance: 1e+16
CycleLastHistoryDump: 0
ClusterSMBHJetOpenAngleRadius: 0.0
MetaDataDatasetUUID: ad0ae2b4-3428-432f-9428-b220d4e9a6ce
MBHFeedbackMassEjectionFraction: 0.1
ExternalGravityOrientation: [0 0 0]
EOSType: 0
MHDCTSlopeLimiter: 1
ClusterSMBHFeedbackSwitch: 0
CoolDataIpiht: 1
SmallRho: 1e-30
ResetMagneticField: 0
MinimumSlopeForRefinement: [ 0.3  0.3  0.3  0.3  0.3  0.3  0.3]
VelAnyl: 0
TimeLastInterpolatedDataDump: -99999
DrivingEfficiency: 1.0
DomainRightEdge: [1 1 1]
ClusterSMBHJetLaunchOffset: 10.0
StarEnergyToThermalFeedback: 1e-05
SubcycleNumber: 1374
SubcycleLastDataDump: 0
PointSourceGravity: 0
RiemannSolver: 5
MHD_WriteElectric: 1
RadiationXRayComptonHeating: 0
ResubmitOn: 0
H2StarMakerSigmaOverR: 0.0333333
MustRefineRegionRightEdge: [1 1 1]
DivergenceCleaningBoundaryBuffer: 0
LoadBalancingCycleSkip: 10
RadiationRedshiftDropOff: 0.0
CMBTemperatureFloor: 1
StarEnergyToQuasarUV: 5e-06
NumberOfGhostZones: 3
#DataCGSConversionFactor[1]: 1448940000.0
StellarWindTurnOnMass: 0.1
MBHMinDynamicalTime: 10000000.0
UseAmbipolarDiffusion: 0
RadiativeTransferHydrogenOnly: 0
RestartDumpDir: RS
ClusterSMBHAccretionTime: 5.0
CorrectParentBoundaryFlux: 0
ClusterSMBHJetPrecessionPeriod: 0.0
RadiationRedshiftOff: 0.0
InlineHaloFinder: 0
SimpleRampTime: 0.1
DataDumpDir: DD
ShockwaveRefinementMaxLevel: 0
BigStarFormationDone: 0
Mu: 0.6
DivBDampingLength: 1
PopIIIBlackHoles: 0
EOSSoundSpeed: 26500
PointSourceGravityPosition: [0 0 0]
StarMakerMinimumMass: 1000000.0
CosmologyCurrentRedshift: 2.22044604925e-16
RestartDumpNumber: 0
CycleSkipGlobalDataDump: 0
dtDataDump: 5
TopGridRank: 3
MBHAccretionRadius: 50
InitialCPUTime: 749.883
RadiationXRaySecondaryIon: 0
GadgetEquilibriumCooling: 0
PopIIIColorDensityThreshold: 1000000.0
MSStellarWindTurnOnMass: 10.0
MHDPowellSource: 0
GreensFunctionMaxSize: 1
MinimumEnergyRatioForRefinement: 0.1
BigStarSeparation: 0.25
RestartDumpName: restart
RadiativeTransferSourceRadius: 0
MultiRefineRegionMinimumOuterLevel: -99999
NumberOfDustTemperatureBins: 250
CoolDataCompXray: 0.0
DataLabel[1]: x-velocity
DeuteriumToHydrogenRatio: 6.8e-05
dtPhoton: -99999
OutputParticleTypeGrouping: 0
MetalCoolingTable: metal_cool.dat
ZEUSLinearArtificialViscosity: 0
PopIIIInitialMassFunction: 0
ShockTemperatureFloor: 1.0
MaximumAlvenSpeed: 1e+30
GasDragCoefficient: 0.0
TimeLastRestartDump: 0
dtInterpolatedDataDump: 0
RadiativeTransferTraceSpectrum: 0
RadiationFieldType: 1
ClusterSMBHEnoughColdGas: 10000000.0
DualEnergyFormalism: 1
CycleSkipDataDump: 0
#DataCGSConversionFactor[11]: 2.49552e-30
ClusterSMBHStartTime: 1.0
NumberOfParticles: 32768
H2StarMakerStochastic: 1
RefineByJeansLengthUnits: 0
PopIIISupernovaRadius: 1
ReconstructionMethod: 1
TimeLastHistoryDump: 0.816510360156
MBHFeedbackThermalRadius: 50
RadiativeTransferHIIRestrictedTimestep: 0
MustRefineParticlesRefineToLevel: 0
MHDCTUseSpecificEnergy: 1
tiny_pressure: 1e-20
MoveParticlesBetweenSiblings: 1
ClusterSMBHJetRadius: 6.0
StopFirstTimeAtLevel: 0
RefineRegionLeftEdge: [0 0 0]
MBHFeedbackRadiativeEfficiency: 0.1
StarClusterMinDynamicalTime: 10000000.0
StarClusterUseMetalField: 0
ShearingBoxProblemType: 0
SelfGravity: 1
CosmologyOmegaLambdaNow: 0.732
CycleSkipHistoryDump: 0
AnisotropicConduction: 0
MHD_CT_Method: 2
MultiRefineRegionMaximumOuterLevel: -99999
StarClusterSNEnergy: 6.8e+48
MinimumShearForRefinement: 1.0
PoissonBoundaryType: 0
MBHCombineRadius: 50
StarEnergyToStellarUV: 3e-06
ClusterSMBHJetMdot: 3.0
CoolingPowerCutOffDensity1: 0
CoolingPowerCutOffDensity2: 10000000000.0
SmallP: 1e-35
MetallicityRefinementMinLevel: -1
SetUVBAmplitude: 1
FastSiblingLocatorEntireDomain: 1
S2ParticleSize: 3
MovieTimestepCounter: 1374
DualEnergyFormalismEta2: 0.1
DualEnergyFormalismEta1: 0.001
SmoothedDarkMatterNeighbors: 32
PopIIISupernovaMustRefineResolution: 32
DivergenceCleaningThreshold: 0.001
CIECooling: 1
AdjustUVBackgroundHighRedshift: 0
RadiationSpectrumNormalization: 1e-21
CoolingCutOffTemperature: 0
MaximumRefinementLevel: 5
StarMakerMinimumDynamicalTime: 1000000.0
UseConstantAcceleration: 0
MBHFeedbackEnergyCoupling: 0.05
ProblemType: 30
ResetLoadBalancing: 0
WriteBoundary: 0
GlobalDir: .
StorePreShockFields: 0
ClusterSMBHJetDim: 2
ResubmitCommand: (null)
RadiationRedshiftOn: 7.0
MustRefineRegionLeftEdge: [0 0 0]
CosmologyComovingBoxSize: 32
PhotoelectricHeatingRate: 8.5e-26
MetallicityRefinementMinMetallicity: 1e-05
InitialTime: 229.708469913
NumberOfRootGridTilesPerDimensionPerProcessor: 1
RadiativeTransferPeriodicBoundary: 0
SubgridSizeAutoAdjust: 1
InterpolationMethod: 1
PPMSteepeningParameter: 0
UseFloor: 0
HaloFinderCycleSkip: 3
ParticleSplitterIterations: 0
#TemperatureUnits: 25434000000.0
ReadGhostZones: 0
ClusterSMBHJetEdot: 1.0
SecondDerivativeFlaggingFields: []
RadiativeTransferTimestepVelocityLimit: 100.0
LoadBalancingMinLevel: 0
RefineByResistiveLengthSafetyFactor: 2.0
MaximumGravityRefinementLevel: 5
PressureFree: 0
DataLabel[4]: TotalEnergy
RadiativeTransferOpticallyThinH2: 0
IncrementDensityOutput: 999
WritePotential: 0
AngularVelocity: 0.001
TracerParticleDumpNumber: 0
StartDensityOutputs: 999
PopIIIBHLuminosityEfficiency: 0.1
MustRefineRegionMinRefinementLevel: -1
VersionNumber: 2.0
FixedTimestep: -1
StarFormationOncePerRootGridTimeStep: 0
HydrogenFractionByMass: 0.76
MustRefineParticlesMinimumMass: 0.0
RadiativeTransferPhotonMergeRadius: 10.0
ParticleSubgridDepositMode: 1
TemperatureEnd: 100000000.0
SolarMetalFractionByMass: 0.02041
MBHMinimumMass: 1000
tiny_number: 1e-20
NewMovieParticleOn: 0
ParallelParticleIO: 1
ViscosityCoefficient: 0
NewMovieDumpNumber: 0
CosmologyHubbleConstantNow: 0.704
MBHFeedbackMetalYield: 0.02
MinimumMassForRefinement: [  4.01738864e-05   2.03966739e-04   7.53260371e-06   7.53260371e-06
   7.53260371e-06   7.53260371e-06   7.53260371e-06]
DataLabel[9]: HeI_Density
MinimumPressureSupportParameter: 100.0
PopIIIMassRange: [   1.  300.]
LoadBalancingMaxLevel: 50
CurrentDensityOutput: 999
ShearingVelocityDirection: -1
PopIIIInitialMassFunctionSlope: -1.3
BigStarFormation: 0
StarClusterMinimumMass: 1000
RadiativeTransferRadiationPressure: 0
HaloFinderLastTime: 0.0
#DataCGSConversionFactor[2]: 1448940000.0
TracerParticleDumpName: TracerOutput
SecondDerivativeEpsilon: 0.01
HistoryDumpDir: HD
MultiMetals: 0
Unigrid: 0
HaloFinderRunAfterOutput: 1
DataLabel[7]: HI_Density
UseViscosity: 0
DataLabel[8]: HII_Density
ExternalGravityRadius: 0.0
dtHistoryDump: 0
H2StarMakerColdGasTemperature: 10000
MHDCTDualEnergyMethod: 2
NumberOfBufferZones: 1
ClusterSMBHKineticFraction: 1.0
StarClusterRegionRightEdge: [ 1.  1.  1.]
EquationOfState: 0
CloudyCoolingGridFile: solar_2008_3D_metals.h5
NewMovieName: MoviePack
GreensFunctionMaxNumber: 40
MBHAccretionFixedTemperature: 300000
RadiativeTransferSplitPhotonRadius: -99999.0
#DataCGSConversionFactor[10]: 2.49552e-30
PointSourceGravityCoreRadius: 0
TracerParticleOn: 0
TimeLastTracerParticleDump: 0.816510360156
PythonTopGridSkip: 0
MaximumTopGridTimeStep: 1e+20
TimeLastDataDump: 225.81646036
HierarchyFileInputFormat: 1
PythonReloadScript: 0
RadiativeTransferSourceBeamAngle: 30.0
CosmologyInitialRedshift: 50
MinimumEfficiency: 0.35
ResetMagneticFieldAmplitude: [0 0 0]
WriteGhostZones: 0
ShockwaveRefinementMinMach: 1.3
DataLabel[0]: Density
PopIIISupernovaMustRefine: 0
HaloFinderMinimumSize: 50
RadiativeTransferInitialHEALPixLevel: 3
PPMDiffusionParameter: 0
H2StarMakerMinimumH2FractionForStarFormation: 1e-05
TracerParticleDumpDir: TD
MBHParticleIOFilename: mbh_particle_io.dat
RefineRegionAutoAdjust: 0
UsePhysicalUnit: 0
RefineBy: 2
MetallicityRefinementMinDensity: -99999
DataLabel[10]: HeII_Density
NumberOfTemperatureBins: 600
RootGridCourantSafetyNumber: 1.0
RadiationFieldRedshift: 0.0
MBHAccretingMassRatio: 1
UseGasDrag: 0
SetHeIIHeatingScale: 1.8
HaloFinderSubfind: 0
CheckpointRestart: 0
PotentialIterations: 4
MinimumSecondDerivativeForRefinement: [ 0.3  0.3  0.3  0.3  0.3  0.3  0.3]
ParallelRootGridIO: 1
Debug1: 0
Debug2: 0
PopIIIOverDensityThreshold: 1000000.0
SlopeFlaggingFields: []
#DataCGSConversionFactor[6]: 2.49552e-30
MBHAccretion: 0
MemoryLimit: 4000000000
ClusterSMBHFastJetVelocity: 10000.0
IsotropicConduction: 0
PopIIIInitialMassFunctionCalls: 0
ParticleBoundaryType: 3
ThreeBodyRate: 0
ClusterSMBHTramp: 0.1
StellarWindFeedback: 0
AnisotropicConductionSpitzerFraction: 1.0
DataLabel[3]: z-velocity
ShockwaveRefinementMinVelocity: 10000000.0
IsotropicConductionSpitzerFraction: 1.0
CycleSkipRestartDump: 0
RadiationFieldLevelRecompute: 0
RadiativeCoolingModel: 1
dtTracerParticleDump: 0
MBHTurnOffStarFormation: 0
LoadBalancing: 1
H2StarMakerEfficiency: 0.01
CosmologyOmegaMatterNow: 0.268
StarMakerTypeIaSNe: 0
ImplicitProblem: 0
RefineByJeansLengthSafetyFactor: 4.0
StarClusterRegionLeftEdge: [ 0.  0.  0.]
#DataCGSConversionFactor[8]: 2.49552e-30
ConstantAcceleration: [0 0 0]
StarMakerOverDensityThreshold: 100
RadiativeTransferLoadBalance: 0
MaximumParticleRefinementLevel: 5
RadiativeTransferRaysPerCell: 5.1
ClusterSMBHFastJetRadius: 0.1
StarClusterFormEfficiency: 0.1
StarClusterIonizingLuminosity: 1e+47
HistoryDumpNumber: 0
OutputTemperature: 1
GravityResolution: 1
PointSourceGravityConstant: 1
CosmologyFinalRedshift: 2.22044604925e-16
BoundaryConditionName: ./DD0046/DD0046.boundary
RefineRegionRightEdge: [1 1 1]
SelfGravityGasOff: 0
RightFaceBoundaryCondition: [3 3 3]
AccretionKernal: 0
UseHydro: 1
Theta_Limiter: 1.5
MetalCooling: 3
RefineRegionTimeType: -1
AddParticleAttributes: 1
CoolingCutOffDensity2: 10000000000.0
CoolingCutOffDensity1: 0
StarFeedbackDistCellStep: 3
ConductionCourantSafetyNumber: 0.5
ClusterSMBHJetVelocity: 10000.0
CurrentTimeIdentifier: 1376404298
UniformGravityDirection: 0
StarFeedbackDistRadius: 1
UseMinimumPressureSupport: 0
RadiativeTransferPropagationDistance: 0.1
H2StarMakerH2FloorInColdGas: 0
JeansRefinementColdTemperature: -1.0
StarClusterSNRadius: 10
TabulatedLWBackground: 0
MultiSpecies: 1
GravitationalConstant: 1.0
ComputePotential: 0
StarMakerMassEfficiency: 0.1
CosmologyMaxExpansionRate: 0.02
HaloFinderTimestep: -99999.0
PopIIIMetalCriticalFraction: 0.0001
OutputFirstTimeAtLevel: 0
PhotoelectricHeating: 0
PopIIISupernovaUseColour: 0
PopIIIInitialMassFunctionSeed: -99999
AdjustUVBackground: 1
StopCPUTime: 2592000.0
ExternalGravityDensity: 0.0
PoissonApproximationThreshold: 0.001
ExternalGravityConstant: 0.0
DustTemperatureStart: 1.0
StarParticleFeedback: 1
H2StarMakerH2DissociationFlux_MW: 1
TopGridDimensions: [32 32 32]
DataDumpName: DD
PopIIIColorMass: 1000000.0
RadiativeTransferTraceSpectrumTable: spectrum_table.dat
NumberOfOutputsBeforeExit: 0
MBHFeedback: 0
CloudyElectronFractionFactor: 0.009154
CoolDataf0to3: 0.1
DataDumpNumber: 46
RadiationRedshiftFullOn: 6.0
FluxCorrection: 1
RadiativeTransferFluxBackgroundLimit: 0.01
MinimumMassForRefinementLevelExponent: [ 0.  0.  0.  0.  0.  0.  0.]
PPMFlatteningParameter: 0
#TimeUnits: 1.89806e+15
dtRestartDump: -99999
DustTemperatureEnd: 1500.0
Gamma: 1.6667
OptimalSubgridsPerProcessor: 16
ConservativeReconstruction: 0
SubcycleSkipDataDump: 0
UniformGravityConstant: 1
MovieDataField: []
FileDirectedOutput: 1
EOSCriticalDensity: 1e-13
HydroMethod: 0
RandomForcingEdot: -1
H2StarMakerMinimumMass: 0
CellFlaggingMethod: [2 4]
RadiativeTransferSourceClustering: 0
Movie3DVolumes: 0
PythonSubcycleSkip: 1
StarClusterHeliumIonization: 0
MHD_ProjectE: 1
StaticHierarchy: 0
ConservativeInterpolation: 1
#DataCGSConversionFactor[7]: 2.49552e-30
UniformGravity: 0
MHD_ProjectB: 0
TemperatureStart: 1.0
CycleLastRestartDump: 0
HaloFinderLinkingLength: 0.1
DataLabel[6]: Electron_Density
VelocityGradient: 1.0
UseResistivity: 0
CoolDataTempXray: 0.0
MovieSkipTimestep: -99999
RadiativeTransferAdaptiveTimestep: 0
StopTime: 229.708469913
H2StarMakerUseSobolevColumn: 0
SinkMergeMass: 0.1
ZEUSQuadraticArtificialViscosity: 2
RadiationSpectrumSlope: 1.5
RadiativeTransferCoupledRateSolver: 1
BaryonSelfGravityApproximation: 1
GridVelocity: [0 0 0]
TopGridGravityBoundary: 0
NumberOfParticleAttributes: 3
#DataCGSConversionFactor[9]: 2.49552e-30
CycleLastDataDump: 0
OutputDustTemperature: 0
CT_AthenaDissipation: 0.1
RadiativeTransferInterpolateField: 0
H2OpticalDepthApproximation: 1
UseDivergenceCleaning: 0
OutputSmoothedDarkMatter: 0
SimpleQ: 1e+50
H2StarMakerNumberDensityThreshold: 0
ClusterSMBHJetAnglePhi: 0.0
TimingCycleSkip: 1
PopIIIHeliumIonization: 0
RadiativeTransferFLD: 0
#DataCGSConversionFactor[12]: 2.49552e-30
StarMakerPlanetaryNebulae: 0
ProcessorTopology: []
IncludeCloudyHeating: 0
ClusterSMBHAccretionEpsilon: 0.001
HaloFinderOutputParticleList: 0
RiemannSolverFallback: 1
RedshiftDumpName: RD
StarClusterUnresolvedModel: 0
ClusterSMBHCalculateGasMass: 1
#DataCGSConversionFactor[0]: 2.49552e-30
StopSteps: 10000
StarMakerSHDensityThreshold: 7e-26
DataLabel[11]: HeIII_Density
#DataCGSConversionFactor[3]: 1448940000.0
MBHAccretionFixedRate: 0.001
BAnyl: 0
ParticleSplitterChildrenParticleSeparation: 1.0
PositiveReconstruction: 0
StarMetalYield: 0.1
PopIIIH2CriticalFraction: 0.0005
MovieVertexCentered: 0
SmallT: 1e-10
MBHInsertLocationFilename: mbh_insert_location.in
CubeDumpEnabled: 0
CoolDataIh2co: 1
FindShocksOnlyOnOutput: 0
UseDrivingField: 0
MinimumPressureJumpForRefinement: 0.33
PopIIIStarMass: 100
InitialCycleNumber: 235
NewMovieRightEdge: [1 1 1]
ParticleTypeInFile: 1
MetaDataSimulationUUID: 8268b492-df69-4853-b27e-4ea4665d9279
ExternalGravity: 0
DomainLeftEdge: [0 0 0]
huge_number: 1e+20
MinimumSubgridEdge: 4
ComovingCoordinates: 1
MBHParticleIO: 0
NewMovieLeftEdge: [0 0 0]
RadiativeTransferRadiationPressureScale: 1.0
DataLabel[2]: y-velocity
LeftFaceBoundaryCondition: [3 3 3]
StarClusterCombineRadius: 10
ExternalGravityPosition: [0 0 0]
StarMassEjectionFraction: 0.25
ShockMethod: 0
ExtractFieldsOnly: 1
RadiativeTransferPropagationSpeedFraction: 1
HierarchyFileOutputFormat: 2
RadiationShield: 0
MinimumOverDensityForRefinement: [ 8.   8.   1.5  1.5  1.5  1.5  1.5]
RadiativeTransferPhotonEscapeRadius: 0.0
ParticleCourantSafetyNumber: 0.8
RandomForcing: 0
DataLabel[12]: Metal_Density
EOSGamma: 1.667
RedshiftDumpDir: RD
RadiativeCooling: 1
RadiativeTransfer: 0
StopCycle: 100000
RadiativeTransferFLDCallOnLevel: 0
RKOrder: 2
IsothermalSoundSpeed: 1
OutputGriddedStarParticle: 0
MBHFeedbackJetsThresholdMass: 10
MaximumSubgridSize: 2000
MustRefineParticlesRefineToLevelAutoAdjust: 0
StarParticleCreation: 1
PartitionNestedGrids: 0
OutputOnDensity: 0
ClusterSMBHFeedback: 0

In [14]:
ds2 = yt.load("../sample_data/Enzo_64/DD0033/data0033")
print ds2.current_redshift
0.270690474165

In [15]:
## At z=0, the comoving Mpc/h should be equal to the proper Mpc/h:
print "ds:"
print ds.domain_width.in_units("Mpccm/h")
print ds.domain_width.in_units("Mpc/h")

## Here, at z-0.27, the comoving and proper box sizes should be different.
print "ds2:"
print ds2.domain_width.in_units("Mpccm/h")
print ds2.domain_width.in_units("Mpc/h")
ds:
[ 128.  128.  128.] Mpccm/h
[ 127.82193279  127.82193279  127.82193279] Mpc/h
ds2:
[ 128.  128.  128.] Mpccm/h
[ 100.73263521  100.73263521  100.73263521] Mpc/h

In [16]:
p = yt.ProjectionPlot(ds, 'x', 'density', weight_field='density')
p.show()

In [17]:
import yt
from yt.analysis_modules.halo_analysis.api import HaloCatalog
hc = HaloCatalog(data_ds=ds, finder_method='hop')
hc.create()
In [18]:
ad = hc.halos_ds.all_data()
masses = ad['particle_mass'][:]
print 'Found %i halos, ranging from %e to %e' % (masses.size, masses.min(), masses.max())
print 'Found %i halos, ranging from %e to %e' % (masses.size, masses.min().in_units("Msun/h"), masses.max().in_units("Msun/h"))
Found 138 halos, ranging from 6.064129e+46 to 2.696205e+48
Found 138 halos, ranging from 2.165307e+13 to 9.627290e+14

In [19]:
p.annotate_halos(hc, annotate_field='particle_mass')
Out[19]:

In [21]:
hc = HaloCatalog(data_ds=ds, finder_method='hop')

# Filter out less massive halos
hc.add_filter("quantity_value", "particle_mass", ">", 1e14, "Msun")

# attach a sphere object to each halo whose radius extends
#   to twice the radius of the halo
hc.add_callback("sphere", factor=2.0)

# use the sphere to calculate radial profiles of gas density
# weighted by cell volume in terms of the virial radius
hc.add_callback("profile", ["radius"],
                [("gas", "overdensity")],
                weight_field="cell_volume",
                accumulation=False,
                storage="virial_quantities_profiles")


hc.add_callback("virial_quantities", ["radius"],
                profile_storage="virial_quantities_profiles")
hc.add_callback('delete_attribute', 'virial_quantities_profiles')

field_params = dict(virial_radius=('quantity', 'radius_200'))
hc.add_callback('sphere', radius_field='radius_200', factor=5,
                field_parameters=field_params)
hc.add_callback('profile', ['virial_radius_fraction'], 
                [('gas', 'temperature')],
                storage='virial_profiles',
                weight_field='cell_mass',
                accumulation=False, output_dir='profiles')

# Save the profiles
hc.add_callback("save_profiles", storage="virial_profiles",
                output_dir="profiles")

hc.create()
---------------------------------------------------------------------------
YTFieldNotFound                           Traceback (most recent call last)
<ipython-input-21-457de4bd14a7> in <module>()
     34                 output_dir="profiles")
     35 
---> 36 hc.create()

/home/skillman/local/yt-x86_64/src/yt-3.0/yt/analysis_modules/halo_analysis/halo_catalog.pyc in create(self, save_halos, save_catalog, njobs, dynamic)
    295 
    296         """
--> 297         self._run(save_halos, save_catalog, njobs=njobs, dynamic=dynamic)
    298 
    299     def load(self, save_halos=True, save_catalog=False, njobs=-1, dynamic=False):

/home/skillman/local/yt-x86_64/src/yt-3.0/yt/utilities/parallel_tools/parallel_analysis_interface.pyc in barrierize(*args, **kwargs)
    265     def barrierize(*args, **kwargs):
    266         if not parallel_capable:
--> 267             return func(*args, **kwargs)
    268         mylog.debug("Entering barrier before %s", func.func_name)
    269         comm = _get_comm(args)

/home/skillman/local/yt-x86_64/src/yt-3.0/yt/analysis_modules/halo_analysis/halo_catalog.pyc in _run(self, save_halos, save_catalog, njobs, dynamic)
    387             for action_type, action in self.actions:
    388                 if action_type == "callback":
--> 389                     action(new_halo)
    390                 elif action_type == "filter":
    391                     halo_filter = action(new_halo)

/home/skillman/local/yt-x86_64/src/yt-3.0/yt/analysis_modules/halo_analysis/halo_callbacks.pyc in __call__(self, halo)
     54 
     55     def __call__(self, halo):
---> 56         self.function(halo, *self.args, **self.kwargs)
     57         return True
     58 

/home/skillman/local/yt-x86_64/src/yt-3.0/yt/analysis_modules/halo_analysis/halo_callbacks.pyc in profile(halo, bin_fields, profile_fields, n_bins, extrema, logs, weight_field, accumulation, fractional, storage, output_dir)
    224     my_profile = create_profile(halo.data_object, bin_fields, profile_fields, n_bins=n_bins,
    225                                 extrema=extrema, logs=logs, weight_field=weight_field,
--> 226                                 accumulation=accumulation, fractional=fractional)
    227 
    228     prof_store = dict([(field, my_profile[field]) \

/home/skillman/local/yt-x86_64/src/yt-3.0/yt/data_objects/profiles.pyc in create_profile(data_source, bin_fields, fields, n_bins, extrema, logs, units, weight_field, accumulation, fractional)
   1304     else:
   1305         raise NotImplementedError
-> 1306     bin_fields = data_source._determine_fields(bin_fields)
   1307     fields = data_source._determine_fields(fields)
   1308     if units is not None:

/home/skillman/local/yt-x86_64/src/yt-3.0/yt/data_objects/data_containers.pyc in _determine_fields(self, fields)
    487             else:
    488                 fname = field
--> 489                 finfo = self.ds._get_field_info("unknown", fname)
    490                 if finfo.particle_type:
    491                     ftype = self._current_particle_type

/home/skillman/local/yt-x86_64/src/yt-3.0/yt/data_objects/static_output.pyc in _get_field_info(self, ftype, fname)
    480                     self._last_finfo = self.field_info[(ftype, fname)]
    481                     return self._last_finfo
--> 482         raise YTFieldNotFound((ftype, fname), self)
    483 
    484     def _setup_classes(self):

YTFieldNotFound: Could not find field '('all', 'virial_radius_fraction')' in data0043.
In []:
!h5ls halo_catalogs/catalog/catalog.0.h5
In []:
hc
In []:
hc.halos_ds.field_list
In []: