{
"metadata": {
"name": ""
},
"nbformat": 3,
"nbformat_minor": 0,
"worksheets": [
{
"cells": [
{
"cell_type": "code",
"collapsed": false,
"input": [
"from yt.imods import *"
],
"language": "python",
"metadata": {},
"outputs": [
{
"output_type": "stream",
"stream": "stderr",
"text": [
"/home/mturk/yt-x86_64/lib/python2.7/site-packages/matplotlib-1.3.0-py2.7-linux-x86_64.egg/matplotlib/__init__.py:727: UserWarning: Found matplotlib configuration in ~/.matplotlib/. To conform with the XDG base directory standard, this configuration location has been deprecated on Linux, and the new location is now '/home/mturk/.config'/matplotlib/. Please move your configuration there to ensure that matplotlib will continue to find it in the future.\n",
" _get_xdg_config_dir())\n"
]
}
],
"prompt_number": 2
},
{
"cell_type": "code",
"collapsed": false,
"input": [
"pf = load(\"DD0046/DD0046\")"
],
"language": "python",
"metadata": {},
"outputs": [],
"prompt_number": 9
},
{
"cell_type": "code",
"collapsed": false,
"input": [
"pf.current_redshift"
],
"language": "python",
"metadata": {},
"outputs": [
{
"metadata": {},
"output_type": "pyout",
"prompt_number": 10,
"text": [
"-2.7810863612032e-09"
]
}
],
"prompt_number": 10
},
{
"cell_type": "code",
"collapsed": false,
"input": [
"pf2 = load(\"DD0027/DD0027\")"
],
"language": "python",
"metadata": {},
"outputs": [],
"prompt_number": 11
},
{
"cell_type": "code",
"collapsed": false,
"input": [
"pf2.current_redshift"
],
"language": "python",
"metadata": {},
"outputs": [
{
"metadata": {},
"output_type": "pyout",
"prompt_number": 12,
"text": [
"0.58703627083642"
]
}
],
"prompt_number": 12
},
{
"cell_type": "code",
"collapsed": false,
"input": [
"pf.domain_left_edge"
],
"language": "python",
"metadata": {},
"outputs": [
{
"metadata": {},
"output_type": "pyout",
"prompt_number": 13,
"text": [
"array([ 0., 0., 0.])"
]
}
],
"prompt_number": 13
},
{
"cell_type": "code",
"collapsed": false,
"input": [
"pf.domain_right_edge"
],
"language": "python",
"metadata": {},
"outputs": [
{
"metadata": {},
"output_type": "pyout",
"prompt_number": 14,
"text": [
"array([ 1., 1., 1.])"
]
}
],
"prompt_number": 14
},
{
"cell_type": "code",
"collapsed": false,
"input": [
"pf.parameters"
],
"language": "python",
"metadata": {},
"outputs": [
{
"metadata": {},
"output_type": "pyout",
"prompt_number": 15,
"text": [
"{'#DataCGSConversionFactor[0]': 2.52765e-30,\n",
" '#DataCGSConversionFactor[1]': 1462390000.0,\n",
" '#DataCGSConversionFactor[2]': 1462390000.0,\n",
" '#DataCGSConversionFactor[3]': 1462390000.0,\n",
" '#DataCGSConversionFactor[5]': 2.52765e-30,\n",
" '#TemperatureUnits': 26803400000.0,\n",
" '#TimeUnits': 1885950000000000.0,\n",
" 'AccretionKernal': 0,\n",
" 'AddParticleAttributes': 0,\n",
" 'AdjustUVBackground': 1,\n",
" 'AdjustUVBackgroundHighRedshift': 0,\n",
" 'AngularVelocity': 0.001,\n",
" 'AnisotropicConduction': 0,\n",
" 'AnisotropicConductionSpitzerFraction': 1.0,\n",
" 'BAnyl': 0,\n",
" 'BaryonSelfGravityApproximation': 1,\n",
" 'BigStarFormation': 0,\n",
" 'BigStarFormationDone': 0,\n",
" 'BigStarSeparation': 0.25,\n",
" 'BoundaryConditionName': './DD0046/DD0046.boundary',\n",
" 'CIECooling': 1,\n",
" 'CMBTemperatureFloor': 1,\n",
" 'CellFlaggingMethod': array([2, 4]),\n",
" 'CheckpointRestart': 0,\n",
" 'CloudyCoolingGridFile': '',\n",
" 'CloudyElectronFractionFactor': 0.009154,\n",
" 'ComovingCoordinates': 1,\n",
" 'ComputePotential': 0,\n",
" 'ConductionCourantSafetyNumber': 0.5,\n",
" 'ConservativeInterpolation': 0,\n",
" 'ConservativeReconstruction': 0,\n",
" 'ConstantAcceleration': array([0, 0, 0]),\n",
" 'CoolDataCompXray': 0.0,\n",
" 'CoolDataIh2co': 1,\n",
" 'CoolDataIpiht': 1,\n",
" 'CoolDataTempXray': 0.0,\n",
" 'CoolDataf0to3': 0.1,\n",
" 'CoolingCutOffDensity1': 0,\n",
" 'CoolingCutOffDensity2': 10000000000.0,\n",
" 'CoolingCutOffTemperature': 0,\n",
" 'CoolingPowerCutOffDensity1': 0,\n",
" 'CoolingPowerCutOffDensity2': 10000000000.0,\n",
" 'Coordinate': 0,\n",
" 'CosmologyComovingBoxSize': 32,\n",
" 'CosmologyCurrentRedshift': -2.7810863612032e-09,\n",
" 'CosmologyFinalRedshift': -2.6598269142752e-09,\n",
" 'CosmologyHubbleConstantNow': 0.702,\n",
" 'CosmologyInitialRedshift': 50,\n",
" 'CosmologyMaxExpansionRate': 0.02,\n",
" 'CosmologyOmegaLambdaNow': 0.727,\n",
" 'CosmologyOmegaMatterNow': 0.273,\n",
" 'CourantSafetyNumber': 0.3,\n",
" 'CubeDumpEnabled': 0,\n",
" 'CurrentDensityOutput': 999,\n",
" 'CurrentTimeIdentifier': 1326853192,\n",
" 'CycleLastDataDump': 0,\n",
" 'CycleLastHistoryDump': 0,\n",
" 'CycleLastRestartDump': 0,\n",
" 'CycleSkipDataDump': 0,\n",
" 'CycleSkipGlobalDataDump': 0,\n",
" 'CycleSkipHistoryDump': 0,\n",
" 'CycleSkipRestartDump': 0,\n",
" 'DataDumpDir': 'DD',\n",
" 'DataDumpName': 'DD',\n",
" 'DataDumpNumber': 46,\n",
" 'DataLabel[0]': 'Density',\n",
" 'DataLabel[1]': 'x-velocity',\n",
" 'DataLabel[2]': 'y-velocity',\n",
" 'DataLabel[3]': 'z-velocity',\n",
" 'DataLabel[4]': 'TotalEnergy',\n",
" 'DataLabel[5]': 'Metal_Density',\n",
" 'Debug1': 0,\n",
" 'Debug2': 0,\n",
" 'DeuteriumToHydrogenRatio': 6.8e-05,\n",
" 'DivBDampingLength': 1,\n",
" 'DivergenceCleaningBoundaryBuffer': 0,\n",
" 'DivergenceCleaningThreshold': 0.001,\n",
" 'DomainLeftEdge': array([0, 0, 0]),\n",
" 'DomainRightEdge': array([1, 1, 1]),\n",
" 'DrivingEfficiency': 1.0,\n",
" 'DualEnergyFormalism': 0,\n",
" 'DualEnergyFormalismEta1': 0.001,\n",
" 'DualEnergyFormalismEta2': 0.1,\n",
" 'DustTemperatureEnd': 1500.0,\n",
" 'DustTemperatureStart': 1.0,\n",
" 'EOSCriticalDensity': 1e-13,\n",
" 'EOSGamma': 1.667,\n",
" 'EOSSoundSpeed': 26500,\n",
" 'EOSType': 0,\n",
" 'ExternalGravity': 0,\n",
" 'ExternalGravityConstant': 0.0,\n",
" 'ExternalGravityDensity': 0.0,\n",
" 'ExternalGravityOrientation': array([0, 0, 0]),\n",
" 'ExternalGravityPosition': array([0, 0, 0]),\n",
" 'ExternalGravityRadius': 0.0,\n",
" 'ExtractFieldsOnly': 1,\n",
" 'FastSiblingLocatorEntireDomain': 1,\n",
" 'FileDirectedOutput': 1,\n",
" 'FluxCorrection': 1,\n",
" 'GadgetEquilibriumCooling': 0,\n",
" 'Gamma': 1.6667,\n",
" 'GasDragCoefficient': 0.0,\n",
" 'GlobalDir': '.',\n",
" 'GravitationalConstant': 1.0,\n",
" 'GravityResolution': 1,\n",
" 'GreensFunctionMaxNumber': 40,\n",
" 'GreensFunctionMaxSize': 1,\n",
" 'GridVelocity': array([0, 0, 0]),\n",
" 'H2FormationOnDust': 0,\n",
" 'H2OpticalDepthApproximation': 1,\n",
" 'HaloFinderCycleSkip': 3,\n",
" 'HaloFinderLastTime': 0.0,\n",
" 'HaloFinderLinkingLength': 0.1,\n",
" 'HaloFinderMinimumSize': 50,\n",
" 'HaloFinderOutputParticleList': 0,\n",
" 'HaloFinderRunAfterOutput': 1,\n",
" 'HaloFinderSubfind': 0,\n",
" 'HaloFinderTimestep': -99999.0,\n",
" 'HierarchyFileInputFormat': 1,\n",
" 'HierarchyFileOutputFormat': 2,\n",
" 'HistoryDumpDir': 'HD',\n",
" 'HistoryDumpName': 'history',\n",
" 'HistoryDumpNumber': 0,\n",
" 'HydroMethod': 2,\n",
" 'HydrogenFractionByMass': 0.76,\n",
" 'ImplicitProblem': 0,\n",
" 'IncludeCloudyHeating': 0,\n",
" 'IncrementDensityOutput': 999,\n",
" 'InitialCPUTime': 1845.52,\n",
" 'InitialCycleNumber': 232,\n",
" 'InitialTime': 230.66527489234,\n",
" 'InlineHaloFinder': 0,\n",
" 'InterpolationMethod': 1,\n",
" 'IsothermalSoundSpeed': 1,\n",
" 'IsotropicConduction': 0,\n",
" 'IsotropicConductionSpitzerFraction': 1.0,\n",
" 'JeansRefinementColdTemperature': -1.0,\n",
" 'LeftFaceBoundaryCondition': array([3, 3, 3]),\n",
" 'LoadBalancing': 1,\n",
" 'LoadBalancingCycleSkip': 10,\n",
" 'LoadBalancingMaxLevel': 50,\n",
" 'LoadBalancingMinLevel': 0,\n",
" 'MBHAccretingMassRatio': 1,\n",
" 'MBHAccretion': 0,\n",
" 'MBHAccretionFixedRate': 0.001,\n",
" 'MBHAccretionFixedTemperature': 300000,\n",
" 'MBHAccretionRadius': 50,\n",
" 'MBHCombineRadius': 50,\n",
" 'MBHFeedback': 0,\n",
" 'MBHFeedbackEnergyCoupling': 0.05,\n",
" 'MBHFeedbackJetsThresholdMass': 10,\n",
" 'MBHFeedbackMassEjectionFraction': 0.1,\n",
" 'MBHFeedbackMetalYield': 0.02,\n",
" 'MBHFeedbackRadiativeEfficiency': 0.1,\n",
" 'MBHFeedbackThermalRadius': 50,\n",
" 'MBHInsertLocationFilename': 'mbh_insert_location.in',\n",
" 'MBHMinDynamicalTime': 10000000.0,\n",
" 'MBHMinimumMass': 1000,\n",
" 'MBHParticleIO': 0,\n",
" 'MBHParticleIOFilename': 'mbh_particle_io.dat',\n",
" 'MBHTurnOffStarFormation': 0,\n",
" 'MSStellarWindTurnOnMass': 10.0,\n",
" 'MaximumAlvenSpeed': 1e+30,\n",
" 'MaximumGravityRefinementLevel': 5,\n",
" 'MaximumParticleRefinementLevel': 5,\n",
" 'MaximumRefinementLevel': 5,\n",
" 'MaximumSubgridSize': 64,\n",
" 'MaximumTopGridTimeStep': 1e+20,\n",
" 'MemoryLimit': 4000000000,\n",
" 'MetaDataDatasetUUID': 'f953c44f-c6da-4999-8e1b-d41fb7538cbe',\n",
" 'MetaDataSimulationUUID': '8d3f46fb-5f61-4dd9-b486-7f3e7da124c8',\n",
" 'MetalCooling': 0,\n",
" 'MetalCoolingTable': 'metal_cool.dat',\n",
" 'MetallicityRefinementMinDensity': -99999,\n",
" 'MetallicityRefinementMinLevel': -1,\n",
" 'MetallicityRefinementMinMetallicity': 1e-05,\n",
" 'MinimumEfficiency': 0.35,\n",
" 'MinimumEnergyRatioForRefinement': 0.1,\n",
" 'MinimumMassForRefinement': array([ 2.03450509e-05, 2.03450531e-04, 7.62939408e-06,\n",
" 7.62939408e-06, 7.62939408e-06, 7.62939408e-06,\n",
" 7.62939408e-06]),\n",
" 'MinimumMassForRefinementLevelExponent': array([ 0., 0., 0., 0., 0., 0., 0.]),\n",
" 'MinimumOverDensityForRefinement': array([ 4. , 8. , 1.5, 1.5, 1.5, 1.5, 1.5]),\n",
" 'MinimumPressureJumpForRefinement': 0.33,\n",
" 'MinimumPressureSupportParameter': 100.0,\n",
" 'MinimumShearForRefinement': 1.0,\n",
" 'MinimumSlopeForRefinement': array([ 0.3, 0.3, 0.3, 0.3, 0.3, 0.3, 0.3]),\n",
" 'MinimumSubgridEdge': 2,\n",
" 'MoveParticlesBetweenSiblings': 1,\n",
" 'Movie3DVolumes': 0,\n",
" 'MovieDataField': array([], dtype=float64),\n",
" 'MovieSkipTimestep': -99999,\n",
" 'MovieTimestepCounter': 1340,\n",
" 'MovieVertexCentered': 0,\n",
" 'Mu': 0.6,\n",
" 'MultiMetals': 0,\n",
" 'MultiSpecies': 0,\n",
" 'MustRefineParticlesMinimumMass': 0.0,\n",
" 'MustRefineParticlesRefineToLevel': 0,\n",
" 'MustRefineParticlesRefineToLevelAutoAdjust': 0,\n",
" 'MustRefineRegionLeftEdge': array([0, 0, 0]),\n",
" 'MustRefineRegionMinRefinementLevel': -1,\n",
" 'MustRefineRegionRightEdge': array([1, 1, 1]),\n",
" 'NewMovieDumpNumber': 0,\n",
" 'NewMovieLeftEdge': array([0, 0, 0]),\n",
" 'NewMovieName': 'MoviePack',\n",
" 'NewMovieParticleOn': 0,\n",
" 'NewMovieRightEdge': array([1, 1, 1]),\n",
" 'NumberOfBufferZones': 1,\n",
" 'NumberOfDustTemperatureBins': 250,\n",
" 'NumberOfOutputsBeforeExit': 0,\n",
" 'NumberOfParticleAttributes': 3,\n",
" 'NumberOfParticles': 113027,\n",
" 'NumberOfRootGridTilesPerDimensionPerProcessor': 1,\n",
" 'NumberOfTemperatureBins': 502,\n",
" 'OptimalSubgridsPerProcessor': 16,\n",
" 'OutputCoolingTime': 0,\n",
" 'OutputDustTemperature': 0,\n",
" 'OutputFirstTimeAtLevel': 0,\n",
" 'OutputGriddedStarParticle': 0,\n",
" 'OutputOnDensity': 0,\n",
" 'OutputParticleTypeGrouping': 0,\n",
" 'OutputSmoothedDarkMatter': 0,\n",
" 'OutputTemperature': 1,\n",
" 'PPMDiffusionParameter': 0,\n",
" 'PPMFlatteningParameter': 0,\n",
" 'PPMSteepeningParameter': 0,\n",
" 'ParallelParticleIO': 1,\n",
" 'ParallelRootGridIO': 1,\n",
" 'ParticleBoundaryType': 3,\n",
" 'ParticleCourantSafetyNumber': 0.5,\n",
" 'ParticleSplitterChildrenParticleSeparation': 1.0,\n",
" 'ParticleSplitterIterations': 0,\n",
" 'ParticleTypeInFile': 1,\n",
" 'PartitionNestedGrids': 0,\n",
" 'PhotoelectricHeating': 0,\n",
" 'PhotoelectricHeatingRate': 8.5e-26,\n",
" 'PointSourceGravity': 0,\n",
" 'PointSourceGravityConstant': 1,\n",
" 'PointSourceGravityCoreRadius': 0,\n",
" 'PointSourceGravityPosition': array([0, 0, 0]),\n",
" 'PoissonApproximationThreshold': 0.001,\n",
" 'PoissonBoundaryType': 0,\n",
" 'PopIIIBHLuminosityEfficiency': 0.1,\n",
" 'PopIIIBlackHoles': 0,\n",
" 'PopIIIColorDensityThreshold': 1000000.0,\n",
" 'PopIIIColorMass': 1000000.0,\n",
" 'PopIIIH2CriticalFraction': 0.0005,\n",
" 'PopIIIHeliumIonization': 0,\n",
" 'PopIIIInitialMassFunction': 0,\n",
" 'PopIIIInitialMassFunctionCalls': 0,\n",
" 'PopIIIInitialMassFunctionSeed': -99999,\n",
" 'PopIIIInitialMassFunctionSlope': -1.3,\n",
" 'PopIIIMassRange': array([ 1., 300.]),\n",
" 'PopIIIMetalCriticalFraction': 0.0001,\n",
" 'PopIIIOverDensityThreshold': 1000000.0,\n",
" 'PopIIIStarMass': 100,\n",
" 'PopIIISupernovaMustRefine': 0,\n",
" 'PopIIISupernovaMustRefineResolution': 32,\n",
" 'PopIIISupernovaRadius': 1,\n",
" 'PopIIISupernovaUseColour': 0,\n",
" 'PositiveReconstruction': 0,\n",
" 'PotentialIterations': 4,\n",
" 'PressureFree': 0,\n",
" 'PrimordialChemistrySolver': 0,\n",
" 'ProblemType': 30,\n",
" 'PythonSubcycleSkip': 1,\n",
" 'PythonTopGridSkip': 0,\n",
" 'RKOrder': 2,\n",
" 'RadiationFieldLevelRecompute': 0,\n",
" 'RadiationFieldRedshift': 0.0,\n",
" 'RadiationFieldType': 0,\n",
" 'RadiationRedshiftDropOff': 0.0,\n",
" 'RadiationRedshiftFullOn': 6.0,\n",
" 'RadiationRedshiftOff': 0.0,\n",
" 'RadiationRedshiftOn': 7.0,\n",
" 'RadiationShield': 0,\n",
" 'RadiationSpectrumNormalization': 1e-21,\n",
" 'RadiationSpectrumSlope': 1.5,\n",
" 'RadiationXRayComptonHeating': 0,\n",
" 'RadiationXRaySecondaryIon': 0,\n",
" 'RadiativeCooling': 1,\n",
" 'RadiativeCoolingModel': 1,\n",
" 'RadiativeTransfer': 0,\n",
" 'RadiativeTransferAdaptiveTimestep': 0,\n",
" 'RadiativeTransferCoupledRateSolver': 1,\n",
" 'RadiativeTransferFLD': 0,\n",
" 'RadiativeTransferFLDCallOnLevel': 0,\n",
" 'RadiativeTransferFluxBackgroundLimit': 0.01,\n",
" 'RadiativeTransferHIIRestrictedTimestep': 0,\n",
" 'RadiativeTransferHydrogenOnly': 0,\n",
" 'RadiativeTransferInitialHEALPixLevel': 3,\n",
" 'RadiativeTransferInterpolateField': 0,\n",
" 'RadiativeTransferLoadBalance': 0,\n",
" 'RadiativeTransferOpticallyThinH2': 0,\n",
" 'RadiativeTransferPeriodicBoundary': 0,\n",
" 'RadiativeTransferPhotonEscapeRadius': 0.0,\n",
" 'RadiativeTransferPhotonMergeRadius': 10.0,\n",
" 'RadiativeTransferPropagationDistance': 0.1,\n",
" 'RadiativeTransferPropagationSpeedFraction': 1,\n",
" 'RadiativeTransferRadiationPressure': 0,\n",
" 'RadiativeTransferRadiationPressureScale': 1.0,\n",
" 'RadiativeTransferRaysPerCell': 5.1,\n",
" 'RadiativeTransferSourceBeamAngle': 30.0,\n",
" 'RadiativeTransferSourceClustering': 0,\n",
" 'RadiativeTransferSourceRadius': 0,\n",
" 'RadiativeTransferSplitPhotonRadius': -99999.0,\n",
" 'RadiativeTransferTimestepVelocityLimit': 100.0,\n",
" 'RadiativeTransferTraceSpectrum': 0,\n",
" 'RadiativeTransferTraceSpectrumTable': 'spectrum_table.dat',\n",
" 'RandomForcing': 0,\n",
" 'RandomForcingEdot': -1,\n",
" 'ReadGhostZones': 0,\n",
" 'ReconstructionMethod': -99999,\n",
" 'RedshiftDumpDir': 'RD',\n",
" 'RedshiftDumpName': 'RD',\n",
" 'RefineBy': 2,\n",
" 'RefineByJeansLengthSafetyFactor': 4.0,\n",
" 'RefineByJeansLengthUnits': 0,\n",
" 'RefineByResistiveLengthSafetyFactor': 2.0,\n",
" 'RefineRegionAutoAdjust': 0,\n",
" 'RefineRegionLeftEdge': array([0, 0, 0]),\n",
" 'RefineRegionRightEdge': array([1, 1, 1]),\n",
" 'RefineRegionTimeType': -1,\n",
" 'ResetLoadBalancing': 0,\n",
" 'ResetMagneticField': 0,\n",
" 'ResetMagneticFieldAmplitude': array([0, 0, 0]),\n",
" 'RestartDumpDir': 'RS',\n",
" 'RestartDumpName': 'restart',\n",
" 'RestartDumpNumber': 0,\n",
" 'ResubmitCommand': '(null)',\n",
" 'ResubmitOn': 0,\n",
" 'RiemannSolver': -99999,\n",
" 'RiemannSolverFallback': 1,\n",
" 'RightFaceBoundaryCondition': array([3, 3, 3]),\n",
" 'RootGridCourantSafetyNumber': 1.0,\n",
" 'S2ParticleSize': 3,\n",
" 'SelfGravity': 1,\n",
" 'SelfGravityGasOff': 0,\n",
" 'SetHeIIHeatingScale': 1.8,\n",
" 'SetUVBAmplitude': 1,\n",
" 'ShearingBoxProblemType': 0,\n",
" 'ShearingVelocityDirection': -1,\n",
" 'ShockMethod': 0,\n",
" 'ShockTemperatureFloor': 1.0,\n",
" 'ShockwaveRefinementMaxLevel': 0,\n",
" 'ShockwaveRefinementMinMach': 1.3,\n",
" 'ShockwaveRefinementMinVelocity': 10000000.0,\n",
" 'SimpleQ': 1e+50,\n",
" 'SimpleRampTime': 0.1,\n",
" 'SinkMergeDistance': 1.0000000272564224e+16,\n",
" 'SinkMergeMass': 0.1,\n",
" 'SlopeFlaggingFields': array([], dtype=float64),\n",
" 'SmallP': 1e-35,\n",
" 'SmallRho': 1e-30,\n",
" 'SmallT': 1e-10,\n",
" 'SmoothedDarkMatterNeighbors': 32,\n",
" 'SolarMetalFractionByMass': 0.02041,\n",
" 'StarClusterCombineRadius': 10,\n",
" 'StarClusterFormEfficiency': 0.1,\n",
" 'StarClusterHeliumIonization': 0,\n",
" 'StarClusterIonizingLuminosity': 1e+47,\n",
" 'StarClusterMinDynamicalTime': 10000000.0,\n",
" 'StarClusterMinimumMass': 1000,\n",
" 'StarClusterRegionLeftEdge': array([ 0., 0., 0.]),\n",
" 'StarClusterRegionRightEdge': array([ 1., 1., 1.]),\n",
" 'StarClusterSNEnergy': 6.8e+48,\n",
" 'StarClusterSNRadius': 10,\n",
" 'StarClusterUnresolvedModel': 0,\n",
" 'StarClusterUseMetalField': 0,\n",
" 'StarEnergyToQuasarUV': 5e-06,\n",
" 'StarEnergyToStellarUV': 3e-06,\n",
" 'StarEnergyToThermalFeedback': 1e-05,\n",
" 'StarFeedbackDistCellStep': 0,\n",
" 'StarFeedbackDistRadius': 0,\n",
" 'StarMakerMassEfficiency': 0.1,\n",
" 'StarMakerMinimumDynamicalTime': 1000000.0,\n",
" 'StarMakerMinimumMass': 1000000.0,\n",
" 'StarMakerOverDensityThreshold': 100,\n",
" 'StarMakerPlanetaryNebulae': 0,\n",
" 'StarMakerSHDensityThreshold': 7e-26,\n",
" 'StarMakerTypeIaSNe': 0,\n",
" 'StarMassEjectionFraction': 0.25,\n",
" 'StarMetalYield': 0.1,\n",
" 'StarParticleCreation': 1,\n",
" 'StarParticleFeedback': 1,\n",
" 'StartDensityOutputs': 999,\n",
" 'StaticHierarchy': 0,\n",
" 'StellarWindFeedback': 0,\n",
" 'StellarWindTurnOnMass': 0.1,\n",
" 'StopCPUTime': 2592000.0,\n",
" 'StopCycle': 100000,\n",
" 'StopFirstTimeAtLevel': 0,\n",
" 'StopSteps': 10000,\n",
" 'StopTime': 230.66527486408,\n",
" 'StorePreShockFields': 0,\n",
" 'SubcycleLastDataDump': 0,\n",
" 'SubcycleNumber': 1340,\n",
" 'SubcycleSkipDataDump': 0,\n",
" 'SubgridSizeAutoAdjust': 1,\n",
" 'TabulatedLWBackground': 0,\n",
" 'TemperatureEnd': 977236736.0,\n",
" 'TemperatureStart': 9549.926758,\n",
" 'Theta_Limiter': 1.5,\n",
" 'ThreeBodyRate': 0,\n",
" 'TimeLastDataDump': 225.81646043367,\n",
" 'TimeLastHistoryDump': 0.8165104336731,\n",
" 'TimeLastInterpolatedDataDump': -99999,\n",
" 'TimeLastRestartDump': 0,\n",
" 'TimeLastTracerParticleDump': 0.8165104336731,\n",
" 'TopGridDimensions': array([32, 32, 32]),\n",
" 'TopGridGravityBoundary': 0,\n",
" 'TopGridRank': 3,\n",
" 'TracerParticleDumpDir': 'TD',\n",
" 'TracerParticleDumpName': 'TracerOutput',\n",
" 'TracerParticleDumpNumber': 0,\n",
" 'TracerParticleOn': 0,\n",
" 'UniformGravity': 0,\n",
" 'UniformGravityConstant': 1,\n",
" 'UniformGravityDirection': 0,\n",
" 'Unigrid': 0,\n",
" 'UnigridTranspose': 0,\n",
" 'UseAmbipolarDiffusion': 0,\n",
" 'UseConstantAcceleration': 0,\n",
" 'UseDivergenceCleaning': 0,\n",
" 'UseDrivingField': 0,\n",
" 'UseFloor': 0,\n",
" 'UseGasDrag': 0,\n",
" 'UseHydro': 1,\n",
" 'UseMinimumPressureSupport': 0,\n",
" 'UsePhysicalUnit': 0,\n",
" 'UseResistivity': 0,\n",
" 'UseViscosity': 0,\n",
" 'VelAnyl': 0,\n",
" 'VelocityGradient': 1.0,\n",
" 'VersionNumber': 2.0,\n",
" 'ViscosityCoefficient': 0,\n",
" 'WriteGhostZones': 0,\n",
" 'WritePotential': 0,\n",
" 'ZEUSLinearArtificialViscosity': 0,\n",
" 'ZEUSQuadraticArtificialViscosity': 2,\n",
" 'dtDataDump': 5,\n",
" 'dtHistoryDump': 0,\n",
" 'dtInterpolatedDataDump': 0,\n",
" 'dtPhoton': -99999,\n",
" 'dtRestartDump': -99999,\n",
" 'dtTracerParticleDump': 0,\n",
" 'huge_number': 1e+20,\n",
" 'tiny_number': 1e-20}"
]
}
],
"prompt_number": 15
},
{
"cell_type": "code",
"collapsed": false,
"input": [
"pf.parameters[\"CellFlaggingMethod\"]"
],
"language": "python",
"metadata": {},
"outputs": [
{
"metadata": {},
"output_type": "pyout",
"prompt_number": 16,
"text": [
"array([2, 4])"
]
}
],
"prompt_number": 16
},
{
"cell_type": "code",
"collapsed": false,
"input": [
"pf.domain_width * pf[\"mpch\"]"
],
"language": "python",
"metadata": {},
"outputs": [
{
"metadata": {},
"output_type": "pyout",
"prompt_number": 23,
"text": [
"array([ 32.00000009, 32.00000009, 32.00000009])"
]
}
],
"prompt_number": 23
},
{
"cell_type": "code",
"collapsed": false,
"input": [
"pf.h.print_stats()"
],
"language": "python",
"metadata": {},
"outputs": [
{
"output_type": "stream",
"stream": "stdout",
"text": [
"level\t# grids\t # cells\t # cells^3\n",
"----------------------------------------------\n",
" 0\t 2\t 32768\t 31\n",
" 1\t 54\t 103992\t 47\n",
" 2\t 53\t 69336\t 41\n",
" 3\t 48\t 29816\t 31\n",
" 4\t 35\t 17304\t 25\n",
" 5\t 19\t 7072\t 19\n",
"----------------------------------------------\n",
" \t 211\t 260288\n",
"\n",
"\n",
"t = 2.30665275e+02 = 4.35120015e+17 s = 1.37881212e+10 years\n",
"\n",
"Smallest Cell:\n",
"\tWidth: 9.766e-04 1\n",
"\tWidth: 9.766e-04 unitary\n",
"\tWidth: 3.125e-02 Mpchcm\n",
"\tWidth: 3.125e-02 mpchcm\n",
"\tWidth: 3.125e-02 Mpch\n",
"\tWidth: 3.125e-02 mpch\n",
"\tWidth: 4.452e-02 Mpccm\n",
"\tWidth: 4.452e-02 mpccm\n",
"\tWidth: 4.452e-02 Mpc\n",
"\tWidth: 4.452e-02 mpc\n",
"\tWidth: 4.980e-02 aye\n",
"\tWidth: 3.125e+01 kpchcm\n",
"\tWidth: 3.125e+01 kpch\n",
"\tWidth: 4.452e+01 kpccm\n",
"\tWidth: 4.452e+01 kpc\n",
"\tWidth: 3.125e+04 pchcm\n",
"\tWidth: 3.125e+04 pch\n",
"\tWidth: 4.452e+04 pccm\n",
"\tWidth: 4.452e+04 pc\n",
"\tWidth: 6.446e+09 auhcm\n",
"\tWidth: 6.446e+09 auh\n",
"\tWidth: 9.182e+09 aucm\n",
"\tWidth: 9.182e+09 au\n",
"\tWidth: 1.386e+12 rsunhcm\n",
"\tWidth: 1.386e+12 rsunh\n",
"\tWidth: 1.975e+12 rsuncm\n",
"\tWidth: 1.975e+12 rsun\n",
"\tWidth: 5.992e+17 mileshcm\n",
"\tWidth: 5.992e+17 milesh\n",
"\tWidth: 8.535e+17 milescm\n",
"\tWidth: 8.535e+17 miles\n",
"\tWidth: 9.643e+17 kmhcm\n",
"\tWidth: 9.643e+17 kmh\n",
"\tWidth: 1.374e+18 kmcm\n",
"\tWidth: 1.374e+18 km\n",
"\tWidth: 9.643e+22 cmhcm\n",
"\tWidth: 9.643e+22 cmh\n",
"\tWidth: 1.374e+23 cmcm\n",
"\tWidth: 1.374e+23 cm\n"
]
}
],
"prompt_number": 24
},
{
"cell_type": "code",
"collapsed": false,
"input": [
"10.0 / pf[\"mpc\"]"
],
"language": "python",
"metadata": {},
"outputs": [
{
"metadata": {},
"output_type": "pyout",
"prompt_number": 25,
"text": [
"0.21937499938989916"
]
}
],
"prompt_number": 25
},
{
"cell_type": "code",
"collapsed": false,
"input": [
"2000 / pf[\"kpc\"]"
],
"language": "python",
"metadata": {},
"outputs": [
{
"metadata": {},
"output_type": "pyout",
"prompt_number": 26,
"text": [
"0.04387499987797983"
]
}
],
"prompt_number": 26
},
{
"cell_type": "code",
"collapsed": false,
"input": [
"p1 = ProjectionPlot(pf, \"y\", \"Density\")"
],
"language": "python",
"metadata": {},
"outputs": [],
"prompt_number": 31
},
{
"cell_type": "code",
"collapsed": false,
"input": [
"p1.show()"
],
"language": "python",
"metadata": {},
"outputs": [
{
"html": [
"
"
],
"metadata": {},
"output_type": "display_data",
"text": [
""
]
}
],
"prompt_number": 32
},
{
"cell_type": "code",
"collapsed": false,
"input": [
"p2 = ProjectionPlot(pf, \"y\", \"Density\", weight_field=\"Density\")"
],
"language": "python",
"metadata": {},
"outputs": [],
"prompt_number": 33
},
{
"cell_type": "code",
"collapsed": false,
"input": [
"p2.show()"
],
"language": "python",
"metadata": {},
"outputs": [
{
"html": [
"
"
],
"metadata": {},
"output_type": "display_data",
"text": [
""
]
}
],
"prompt_number": 34
},
{
"cell_type": "code",
"collapsed": false,
"input": [
"halos = HaloFinder(pf)"
],
"language": "python",
"metadata": {},
"outputs": [],
"prompt_number": 35
},
{
"cell_type": "code",
"collapsed": false,
"input": [
"len(halos)"
],
"language": "python",
"metadata": {},
"outputs": [
{
"metadata": {},
"output_type": "pyout",
"prompt_number": 36,
"text": [
"30"
]
}
],
"prompt_number": 36
},
{
"cell_type": "code",
"collapsed": false,
"input": [
"for halo in halos:\n",
" print halo.total_mass()"
],
"language": "python",
"metadata": {},
"outputs": [
{
"output_type": "stream",
"stream": "stdout",
"text": [
"3.93029e+13\n",
"3.66947e+13\n",
"3.64249e+13\n",
"3.63349e+13\n",
"2.04158e+13\n",
"1.88869e+13\n",
"1.76277e+13\n",
"1.70881e+13\n",
"1.40302e+13\n",
"1.39403e+13\n",
"1.32208e+13\n",
"1.25912e+13\n",
"1.22315e+13\n",
"1.12422e+13\n",
"1.09724e+13\n",
"1.04328e+13\n",
"1.02529e+13\n",
"8.634e+12\n",
"8.54406e+12\n",
"8.27425e+12\n",
"8.09438e+12\n",
"8.09438e+12\n",
"7.10507e+12\n",
"6.20569e+12\n",
"5.84594e+12\n",
"5.57613e+12\n",
"5.03651e+12\n",
"4.49688e+12\n",
"4.49688e+12\n",
"4.13713e+12\n"
]
}
],
"prompt_number": 37
},
{
"cell_type": "code",
"collapsed": false,
"input": [
"for halo in halos:\n",
" print halo.center_of_mass()"
],
"language": "python",
"metadata": {},
"outputs": [
{
"output_type": "stream",
"stream": "stdout",
"text": [
"[ 0.58943546 0.9174071 0.81971047]\n",
"[ 0.52445403 0.02436666 0.77793467]\n",
"[ 0.46436154 0.20728192 0.581331 ]\n",
"[ 0.53894666 0.72162902 0.47548559]\n",
"[ 0.554628 0.66257589 0.85877835]\n",
"[ 0.55891204 0.75135646 0.67284367]\n",
"[ 0.65487748 0.37812087 0.17422425]\n",
"[ 0.8800615 0.59401279 0.93420305]\n",
"[ 0.20572017 0.21106693 0.00481071]\n",
"[ 0.49429624 0.06680771 0.80341858]\n",
"[ 0.65042379 0.53749597 0.2839634 ]\n",
"[ 0.79709084 0.87278154 0.00839798]\n",
"[ 0.60951889 0.84400326 0.1132075 ]\n",
"[ 0.51798954 0.38783484 0.8616772 ]\n",
"[ 0.52271743 0.76397941 0.29785954]\n",
"[ 0.43688137 0.67239363 0.0975858 ]\n",
"[ 0.05467522 0.19381474 0.04019758]\n",
"[ 0.64962371 0.65253555 0.40207628]\n",
"[ 0.53175531 0.40996185 0.82469474]\n",
"[ 0.48761615 0.96540439 0.59194927]\n",
"[ 0.23046686 0.20758259 0.99864345]\n",
"[ 0.85845639 0.87086715 0.03634816]\n",
"[ 0.54048269 0.87113514 0.68829863]\n",
"[ 0.08484187 0.21368094 0.09348651]\n",
"[ 0.97079234 0.83473201 0.58975742]\n",
"[ 0.44184821 0.12351559 0.06352914]\n",
"[ 0.52165652 0.29968925 0.8076959 ]\n",
"[ 0.97708911 0.62956642 0.9168672 ]\n",
"[ 0.98735729 0.2374026 0.03843924]\n",
"[ 0.28647738 0.5686128 0.12947163]\n"
]
}
],
"prompt_number": 38
},
{
"cell_type": "code",
"collapsed": false,
"input": [
"for halo in halos:\n",
" print halo.maximum_radius() * pf[\"kpc\"]"
],
"language": "python",
"metadata": {},
"outputs": [
{
"output_type": "stream",
"stream": "stdout",
"text": [
"1253.27620754\n",
"1060.15445463\n",
"1152.50099889\n",
"852.809341693\n",
"649.811016763\n",
"747.383554145\n",
"915.310091066\n",
"705.965246075\n",
"751.415936696\n",
"628.671854139\n",
"824.972502769\n",
"699.167793007\n",
"759.904797014\n",
"1290.2863654\n",
"628.080568798\n",
"536.680863178\n",
"625.742475612\n",
"527.243215075\n",
"643.293092082\n",
"652.333377183\n",
"639.181599501\n",
"569.552561232\n",
"525.532295951\n",
"441.647555839\n",
"479.552815762\n",
"415.160397889\n",
"502.249842686\n",
"407.178605949\n",
"475.999165244\n",
"376.943006496\n"
]
}
],
"prompt_number": 40
},
{
"cell_type": "code",
"collapsed": false,
"input": [
"p1.annotate_hop_circles(halos)\n",
"p1.show()"
],
"language": "python",
"metadata": {},
"outputs": [
{
"html": [
"
"
],
"metadata": {},
"output_type": "pyout",
"prompt_number": 41,
"text": [
""
]
}
],
"prompt_number": 41
},
{
"cell_type": "code",
"collapsed": false,
"input": [
"most_massive = halos[0]"
],
"language": "python",
"metadata": {},
"outputs": [],
"prompt_number": 42
},
{
"cell_type": "code",
"collapsed": false,
"input": [
"new_center = most_massive.center_of_mass()"
],
"language": "python",
"metadata": {},
"outputs": [],
"prompt_number": 61
},
{
"cell_type": "code",
"collapsed": false,
"input": [
"p1.set_center((new_center[0], new_center[2]))\n",
"p1.set_width(10, 'mpc')\n",
"p1.show()"
],
"language": "python",
"metadata": {},
"outputs": [
{
"html": [
"
"
],
"metadata": {},
"output_type": "display_data",
"text": [
""
]
}
],
"prompt_number": 63
},
{
"cell_type": "code",
"collapsed": false,
"input": [
"p2.show()"
],
"language": "python",
"metadata": {},
"outputs": [
{
"html": [
"
"
],
"metadata": {},
"output_type": "display_data",
"text": [
""
]
}
],
"prompt_number": 64
},
{
"cell_type": "code",
"collapsed": false,
"input": [
"p2.annotate_contour(\"Temperature\")\n",
"p2.show()"
],
"language": "python",
"metadata": {},
"outputs": [
{
"html": [
"
"
],
"metadata": {},
"output_type": "display_data",
"text": [
""
]
}
],
"prompt_number": 65
},
{
"cell_type": "code",
"collapsed": false,
"input": [
"pf.h.field_list"
],
"language": "python",
"metadata": {},
"outputs": [
{
"metadata": {},
"output_type": "pyout",
"prompt_number": 66,
"text": [
"['particle_type',\n",
" 'Temperature',\n",
" 'Density',\n",
" 'creation_time',\n",
" 'TotalEnergy',\n",
" 'Dark_Matter_Density',\n",
" 'Metal_Density',\n",
" 'x-velocity',\n",
" 'dynamical_time',\n",
" 'particle_index',\n",
" 'particle_mass',\n",
" 'z-velocity',\n",
" 'metallicity_fraction',\n",
" 'particle_position_z',\n",
" 'particle_position_x',\n",
" 'particle_position_y',\n",
" 'y-velocity',\n",
" 'particle_velocity_z',\n",
" 'particle_velocity_x',\n",
" 'particle_velocity_y']"
]
}
],
"prompt_number": 66
},
{
"cell_type": "code",
"collapsed": false,
"input": [
"s = SlicePlot(pf, \"x\", \"x-velocity\", center=\"max\")\n",
"s.show()"
],
"language": "python",
"metadata": {},
"outputs": [
{
"html": [
"
"
],
"metadata": {},
"output_type": "display_data",
"text": [
""
]
}
],
"prompt_number": 67
},
{
"cell_type": "code",
"collapsed": false,
"input": [
"s.set_cmap(\"x-velocity\", \"kamae\")\n",
"s.show()"
],
"language": "python",
"metadata": {},
"outputs": [
{
"html": [
"
"
],
"metadata": {},
"output_type": "display_data",
"text": [
""
]
}
],
"prompt_number": 69
},
{
"cell_type": "code",
"collapsed": false,
"input": [
"s.zoom(2)\n",
"s.show()"
],
"language": "python",
"metadata": {},
"outputs": [
{
"html": [
"
"
],
"metadata": {},
"output_type": "display_data",
"text": [
""
]
}
],
"prompt_number": 70
},
{
"cell_type": "code",
"collapsed": false,
"input": [
"pf.h.derived_field_list"
],
"language": "python",
"metadata": {},
"outputs": [
{
"metadata": {},
"output_type": "pyout",
"prompt_number": 71,
"text": [
"['particle_type',\n",
" 'z-velocity',\n",
" 'Temperature',\n",
" 'metallicity_fraction',\n",
" 'Density',\n",
" 'particle_velocity_z',\n",
" 'creation_time',\n",
" 'TotalEnergy',\n",
" 'Dark_Matter_Density',\n",
" 'particle_position_z',\n",
" 'Metal_Density',\n",
" 'x-velocity',\n",
" 'particle_position_y',\n",
" 'dynamical_time',\n",
" 'y-velocity',\n",
" 'particle_index',\n",
" 'particle_position_x',\n",
" 'particle_mass',\n",
" 'particle_velocity_x',\n",
" 'particle_velocity_y',\n",
" 'StarMetallicity',\n",
" 'ParticleAge',\n",
" 'Dark_Matter_Mass',\n",
" 'star_creation_time',\n",
" 'star_metallicity_fraction',\n",
" 'Metal_MassMsun',\n",
" 'star_density',\n",
" 'ParticleMass',\n",
" 'StarDynamicalTimeYears',\n",
" 'Dark_Matter_MassMsun',\n",
" 'StarAgeYears',\n",
" 'tracer_number_density',\n",
" 'Total_Energy',\n",
" 'ParticleMassMsun',\n",
" 'TotalEnergy',\n",
" 'cic_particle_velocity_x',\n",
" 'cic_particle_velocity_y',\n",
" 'dm_density',\n",
" 'cic_particle_velocity_z',\n",
" 'IsStarParticle',\n",
" 'KineticEnergy',\n",
" 'ThermalEnergy',\n",
" 'Metal_Mass',\n",
" 'Metal_Fraction',\n",
" 'star_dynamical_time',\n",
" 'H_NumberDensity',\n",
" 'Metallicity',\n",
" 'NumberDensity',\n",
" 'particle_mass',\n",
" 'StarCreationTimeYears',\n",
" 'RadialVelocityABS',\n",
" 'SZKinetic',\n",
" 'RadialVelocity',\n",
" 'VorticityZ',\n",
" 'CellVolume',\n",
" 'Zeros',\n",
" 'CellMassCode',\n",
" 'ParticleRadiuskpc',\n",
" 'ParticleRadiusAU',\n",
" 'Pressure',\n",
" 'Ones',\n",
" 'CellMassMsun',\n",
" 'DynamicalTime',\n",
" 'VorticityStretchingY',\n",
" 'dz',\n",
" 'Baryon_Overdensity',\n",
" 'dx',\n",
" 'dy',\n",
" 'StarMassMsun',\n",
" 'XRayEmissivity',\n",
" 'x',\n",
" 'GridLevel',\n",
" 'TangentialVelocity',\n",
" 'AngularMomentumX',\n",
" 'AngularMomentumY',\n",
" 'AngularMomentumZ',\n",
" 'sph_r',\n",
" 'ParticleMass',\n",
" 'Contours',\n",
" 'Radiuspc',\n",
" 'DensityPerturbation',\n",
" 'CellVolumeCode',\n",
" 'VorticityY',\n",
" 'particle_position_z',\n",
" 'sph_phi',\n",
" 'MeanMolecularWeight',\n",
" 'CellsPerBin',\n",
" 'AveragedDensity',\n",
" 'SZY',\n",
" 'TempkeV',\n",
" 'VorticityStretchingX',\n",
" 'cyl_RCode',\n",
" 'VorticityStretchingZ',\n",
" 'VorticityMagnitude',\n",
" 'gradDensityMagnitude',\n",
" 'ParticleSpecificAngularMomentumXKMSMPC',\n",
" 'ParticleRadiuskpch',\n",
" 'cyl_RadialVelocityKMSABS',\n",
" 'VorticitySquared',\n",
" 'cyl_z',\n",
" 'ParticleAngularMomentumZ',\n",
" 'ParticleAngularMomentumX',\n",
" 'ParticleAngularMomentumY',\n",
" 'AbsDivV',\n",
" 'RadialMachNumber',\n",
" 'Radius',\n",
" 'cyl_R',\n",
" 'tempContours',\n",
" 'ParticleSpecificAngularMomentumZKMSMPC',\n",
" 'CourantTimeStep',\n",
" 'RadiusMpc',\n",
" 'RadialVelocityKMS',\n",
" 'MachNumber',\n",
" 'sph_theta',\n",
" 'DivV',\n",
" 'CellVolumeMpc',\n",
" 'JeansMassMsun',\n",
" 'gradDensityZ',\n",
" 'ParticleRadiusMpc',\n",
" 'Radiuskpch',\n",
" 'Matter_Density',\n",
" 'ParticleVelocityMagnitude',\n",
" 'TotalMass',\n",
" 'particle_position_y',\n",
" 'ParticleRadiusCode',\n",
" 'SpecificAngularMomentumX',\n",
" 'SpecificAngularMomentumY',\n",
" 'SpecificAngularMomentumZ',\n",
" 'OnesOverDx',\n",
" 'TotalMassMsun',\n",
" 'RadiusAU',\n",
" 'gradDensityX',\n",
" 'gradDensityY',\n",
" 'cyl_TangentialVelocityABS',\n",
" 'Height',\n",
" 'cyl_TangentialVelocityKMSABS',\n",
" 'Radiuskpc',\n",
" 'particle_position_x',\n",
" 'HeightAU',\n",
" 'VelocityMagnitude',\n",
" 'RadialVelocityKMSABS',\n",
" 'z',\n",
" 'RadiusCode',\n",
" 'SoundSpeed',\n",
" 'cyl_RadialVelocity',\n",
" 'cyl_RadialVelocityKMS',\n",
" 'cyl_RadialVelocityABS',\n",
" 'ParticleSpecificAngularMomentumX',\n",
" 'ParticleSpecificAngularMomentumY',\n",
" 'ParticleRadius',\n",
" 'DiskAngle',\n",
" 'TangentialOverVelocityMagnitude',\n",
" 'Entropy',\n",
" 'ParticleSpecificAngularMomentumZ',\n",
" 'GridIndices',\n",
" 'ChandraEmissivity',\n",
" 'VorticityStretchingMagnitude',\n",
" 'cyl_TangentialVelocity',\n",
" 'CellMass',\n",
" 'ComovingDensity',\n",
" 'particle_density',\n",
" 'ParticleSpecificAngularMomentumYKMSMPC',\n",
" 'ParticleRadiuspc',\n",
" 'cyl_TangentialVelocityKMS',\n",
" 'Overdensity',\n",
" 'y',\n",
" 'cyl_theta',\n",
" 'VorticityX']"
]
}
],
"prompt_number": 71
},
{
"cell_type": "code",
"collapsed": false,
"input": [
"s = SlicePlot(pf, \"x\", \"VelocityMagnitude\", center=\"max\")\n",
"s.show()"
],
"language": "python",
"metadata": {},
"outputs": [
{
"html": [
"
"
],
"metadata": {},
"output_type": "display_data",
"text": [
""
]
}
],
"prompt_number": 72
},
{
"cell_type": "code",
"collapsed": false,
"input": [
"s.zoom(2)\n",
"s.set_axes_unit(\"kpc\")\n",
"s.set_log(\"VelocityMagnitude\", True)\n",
"s.set_cmap(\"VelocityMagnitude\", \"cubehelix\")\n",
"s.annotate_velocity()\n",
"s.show()"
],
"language": "python",
"metadata": {},
"outputs": [
{
"html": [
"
"
],
"metadata": {},
"output_type": "display_data",
"text": [
""
]
}
],
"prompt_number": 74
},
{
"cell_type": "code",
"collapsed": false,
"input": [
"s.zoom(2)\n",
"s.show()"
],
"language": "python",
"metadata": {},
"outputs": [
{
"html": [
"
"
],
"metadata": {},
"output_type": "display_data",
"text": [
""
]
}
],
"prompt_number": 75
},
{
"cell_type": "code",
"collapsed": false,
"input": [
"s.annotate_title(\"Il Bruno, Il Brutto, Il Cattivo\")\n",
"s.show()"
],
"language": "python",
"metadata": {},
"outputs": [
{
"html": [
"
"
],
"metadata": {},
"output_type": "display_data",
"text": [
""
]
}
],
"prompt_number": 76
},
{
"cell_type": "code",
"collapsed": false,
"input": [
"p3 = ProjectionPlot(pf, \"y\", [\"Density\", \"Temperature\", \"AbsDivV\", \"MachNumber\"],\n",
" weight_field = \"Density\")\n",
"p3.show()"
],
"language": "python",
"metadata": {},
"outputs": [
{
"html": [
"
"
],
"metadata": {},
"output_type": "display_data",
"text": [
""
]
}
],
"prompt_number": 83
},
{
"cell_type": "code",
"collapsed": false,
"input": [
"p3.set_log(\"AbsDivV\", True)\n",
"p3.set_cmap(\"Temperature\", \"hot\")\n",
"p3.set_cmap(\"AbsDivV\", \"kamae\")\n",
"p3.set_cmap(\"MachNumber\", \"cubehelix\")\n",
"p3.show()"
],
"language": "python",
"metadata": {},
"outputs": [
{
"html": [
"
"
],
"metadata": {},
"output_type": "display_data",
"text": [
""
]
}
],
"prompt_number": 84
},
{
"cell_type": "code",
"collapsed": false,
"input": [
"p3.zoom(2)\n",
"p3.show()"
],
"language": "python",
"metadata": {},
"outputs": [
{
"html": [
"
"
],
"metadata": {},
"output_type": "display_data",
"text": [
""
]
}
],
"prompt_number": 85
},
{
"cell_type": "code",
"collapsed": false,
"input": [
"p3.save()"
],
"language": "python",
"metadata": {},
"outputs": [
{
"metadata": {},
"output_type": "pyout",
"prompt_number": 86,
"text": [
"['DD0046/DD0046_Projection_y_AbsDivV_Density.png',\n",
" 'DD0046/DD0046_Projection_y_Temperature_Density.png',\n",
" 'DD0046/DD0046_Projection_y_MachNumber_Density.png',\n",
" 'DD0046/DD0046_Projection_y_Density_Density.png']"
]
}
],
"prompt_number": 86
},
{
"cell_type": "code",
"collapsed": false,
"input": [
"halos"
],
"language": "python",
"metadata": {},
"outputs": [
{
"metadata": {},
"output_type": "pyout",
"prompt_number": 87,
"text": [
""
]
}
],
"prompt_number": 87
},
{
"cell_type": "code",
"collapsed": false,
"input": [
"most_massive"
],
"language": "python",
"metadata": {},
"outputs": [
{
"metadata": {},
"output_type": "pyout",
"prompt_number": 88,
"text": [
""
]
}
],
"prompt_number": 88
},
{
"cell_type": "code",
"collapsed": false,
"input": [
"my_sphere = most_massive.get_sphere()"
],
"language": "python",
"metadata": {},
"outputs": [],
"prompt_number": 89
},
{
"cell_type": "code",
"collapsed": false,
"input": [
"my_sphere.radius * pf[\"mpc\"]"
],
"language": "python",
"metadata": {},
"outputs": [
{
"metadata": {},
"output_type": "pyout",
"prompt_number": 91,
"text": [
"1.2532762075390431"
]
}
],
"prompt_number": 91
},
{
"cell_type": "code",
"collapsed": false,
"input": [
"my_sphere[\"Density\"]"
],
"language": "python",
"metadata": {},
"outputs": [
{
"metadata": {},
"output_type": "pyout",
"prompt_number": 92,
"text": [
"array([ 2.49070207e-30, 2.98411297e-30, 5.97460307e-30, ...,\n",
" 1.27477000e-27, 1.04645690e-27, 8.89098805e-28])"
]
}
],
"prompt_number": 92
},
{
"cell_type": "code",
"collapsed": false,
"input": [
"my_sphere[\"Temperature\"]"
],
"language": "python",
"metadata": {},
"outputs": [
{
"metadata": {},
"output_type": "pyout",
"prompt_number": 93,
"text": [
"array([ 5115659.5 , 2222620.25 , 1687170.875, ..., 6484474.5 ,\n",
" 5702430.5 , 5362007. ])"
]
}
],
"prompt_number": 93
},
{
"cell_type": "code",
"collapsed": false,
"input": [
"my_sphere[\"Density\"].size"
],
"language": "python",
"metadata": {},
"outputs": [
{
"metadata": {},
"output_type": "pyout",
"prompt_number": 94,
"text": [
"3652"
]
}
],
"prompt_number": 94
},
{
"cell_type": "code",
"collapsed": false,
"input": [
"my_sphere.quantities[\"AngularMomentumVector\"]()"
],
"language": "python",
"metadata": {},
"outputs": [
{
"metadata": {},
"output_type": "pyout",
"prompt_number": 95,
"text": [
"array([-0.47974874, -0.87363042, -0.08130824])"
]
}
],
"prompt_number": 95
},
{
"cell_type": "code",
"collapsed": false,
"input": [
"my_sphere.quantities[\"Extrema\"](\"Density\")"
],
"language": "python",
"metadata": {},
"outputs": [
{
"metadata": {},
"output_type": "pyout",
"prompt_number": 96,
"text": [
"[(8.386588622152806e-31, 2.1837533438671876e-26)]"
]
}
],
"prompt_number": 96
},
{
"cell_type": "code",
"collapsed": false,
"input": [
"my_sphere.quantities[\"WeightedAverageQuantity\"](\"Temperature\", \"CellMassMsun\")"
],
"language": "python",
"metadata": {},
"outputs": [
{
"metadata": {},
"output_type": "pyout",
"prompt_number": 97,
"text": [
"5215186.0081834383"
]
}
],
"prompt_number": 97
},
{
"cell_type": "code",
"collapsed": false,
"input": [
"\"%0.3e\" % _"
],
"language": "python",
"metadata": {},
"outputs": [
{
"metadata": {},
"output_type": "pyout",
"prompt_number": 98,
"text": [
"'5.215e+06'"
]
}
],
"prompt_number": 98
},
{
"cell_type": "code",
"collapsed": false,
"input": [
"masses = []\n",
"temperatures = []\n",
"\n",
"for halo in halos:\n",
" sp = halo.get_sphere()\n",
" total_mass = sp.quantities[\"TotalMass\"]()\n",
" avg_temp = sp.quantities[\"WeightedAverageQuantity\"](\"Temperature\", \"CellMassMsun\")\n",
" masses.append(total_mass)\n",
" temperatures.append(avg_temp)"
],
"language": "python",
"metadata": {},
"outputs": [],
"prompt_number": 99
},
{
"cell_type": "code",
"collapsed": false,
"input": [
"print masses"
],
"language": "python",
"metadata": {},
"outputs": [
{
"output_type": "stream",
"stream": "stdout",
"text": [
"[56720465549339.719, 51715498622456.859, 52913970018595.75, 46411942429289.195, 25165745747348.477, 25157848649524.148, 25371615424521.797, 22088425117487.18, 18121377062920.234, 17276803105900.252, 19042307282655.898, 17268860517603.289, 17131043706439.941, 22100966965195.246, 13930128498809.535, 12825467934881.279, 13289385249261.043, 10053954890574.197, 11883248433277.059, 10240374488563.861, 10935265277880.361, 10655707196405.16, 8766743352949.1523, 6819245124060.0752, 7012482418455.0039, 6217686381087.4873, 6475544709904.2734, 5251568025607.2451, 5397269614598.1016, 4789926416747.5547]\n"
]
}
],
"prompt_number": 100
},
{
"cell_type": "code",
"collapsed": false,
"input": [
"print temperatures"
],
"language": "python",
"metadata": {},
"outputs": [
{
"output_type": "stream",
"stream": "stdout",
"text": [
"[5215186.0081834383, 6584243.9661670364, 7517465.6976599572, 7815751.5203913832, 5990948.2334709335, 5038441.5903957412, 3913081.3255644245, 5545732.4161329363, 2457276.8501548967, 4421172.6048297463, 3237673.2658461425, 5287235.9263523361, 3135118.0320181036, 1947864.6866669746, 3438434.4668884226, 3800040.3361177375, 3248708.3538134927, 2563261.214621475, 1839855.3541096128, 1648128.6719817019, 2603173.321045273, 3265123.0486687133, 954723.63230793586, 1877172.9734493182, 1977533.9141875124, 2310657.8547804183, 2067395.6738823995, 550105.42744388885, 970254.9868437025, 573102.356046036]\n"
]
}
],
"prompt_number": 101
},
{
"cell_type": "code",
"collapsed": false,
"input": [
"pylab.clf()\n",
"pylab.loglog(masses, temperatures, 'o')\n",
"pylab.xlabel(r\"$\\mathrm{Mass}\\/[M_\\odot]$\")\n",
"pylab.ylabel(r\"$\\mathrm{T}\\/[\\mathrm{K}]$\")\n",
"pylab.show()"
],
"language": "python",
"metadata": {},
"outputs": [
{
"metadata": {},
"output_type": "display_data",
"png": "iVBORw0KGgoAAAANSUhEUgAAAZMAAAEcCAYAAAAC+llsAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAF4NJREFUeJzt3X9s1Hfhx/FXB4waKVodQ0fvK9uVCRU2lB86ZseZCAWb\nsYRko0wDAxHhm17JnIkRSnrdWDLUudAS0IgygkIqstmUcyfd4rUyge4HjSCdjEKhVg0/5g+20W6t\n9/1jX2qvvbZ3fd/d58c9H8kl3qf3+dz78LN73ft3ViQSiQgAAAM3WV0AAIDzESYAAGOECQDAGGEC\nADBGmAAAjBEmAABjhAkAwBhhAgAwNtrqAsTjyJEj+sUvfqHu7m6dPn1aL7/8stVFAgD0keWkGfC1\ntbW6dOmSvvGNb1hdFABAH5Y1c61evVoTJ07UjBkzoo6HQiFNnTpVU6ZM0datW6P+tm/fPj388MPp\nLCYAIA6WhcmqVasUCoWijvX09Ki0tFShUEinT5/W/v371dLSIkm6ePGiPvKRj+jDH/6wFcUFAAzB\nsjApLCxUbm5u1LGmpibl5+dr8uTJGjNmjEpKSlRbWytJ+tnPfqbVq1dbUVQAwDBs1QHf0dEhj8fT\n+zwvL0/Hjx+XJAUCgSHPzc/PV2trayqLBwCu4/V6dfbsWePr2GpocFZW1ojPbW1tVSQSccWjoqLC\nFe9pes2Rnp/IefG+drjXmf7dKQ+rPodb7k8r7s3hXpOsH+G2CpNJkyapvb2993l7e7vy8vIsLJE1\nfD6fK97T9JojPT+R8+J97XCvs+L/MytY9Tndcn9acW8m+r4jZenQ4La2Nt1///06efKkJKm7u1uf\n/vSn9dJLL+m2227T3LlztX//fk2bNm3Ya2VlZcnCjwIMKRAIDNtUC1ghWd+dltVMli9frnnz5unM\nmTPyeDzavXu3Ro8ere3bt6uoqEgFBQVatmxZXEEC2F2m1FyQuRw1aXEoWVlZqqiokM/n4z9cABhG\nOBxWOBxWZWVlUmomrgoTl3wUAEgbxzdzAQDcgzABABiz1aRFALCLYLBRVVWH1dU1WmPHdqusbKGK\ni++zuli2RZgAQD/BYKM2bPitWluf7D3W2rpJkgiUQbiqmSsQCCgcDltdDAAOV1V1OCpIJKm19UlV\nV9cPeG0w2KiionL5fAEVFZUrGGxMVzGNhMPhpM59clXNhElhAJKhqyv2V2Nn56io506uwdyYRlFZ\nWZmU67mqZgIAyTB2bHfM49nZPVHPE6nBuB1hAgD9lJUtlNe7KeqY17tRfv+CqGPx1mAygauauQAg\nGW40UVVXb1Zn5yhlZ/fI7180oOkq3hpMJmAGPACMUKw+E693o7ZtGxg8dpWs705X1UwCgQBrcwFI\nixvzULKzr+jjH1+mT37yk5o0KSdmDcaObqzNlSzUTAAgQbFrJJu0bVuRI4KkL9bmAgCLMIprIMIE\nABLEKK6BCBMASBCjuAYiTACkhVOXHYkl3nkomcRVo7kA2JOTlx2JJd55KJnEVaO52LYXsKeionId\nPrwlxvHNCoWesKBEYNveQTA0GLBGPPt++HwBNTQEBpw7f35A4fDA40gfJi0CsFy8zVd0WLsfHfAA\nRize+RZ0WLsfNRMgQyVjW9p451tkWod1Jm75S5gAGShZo6sSab4qLr7P9V+okvtGrsWLZi4gAyVr\nOZBEmq/cNM9kKJm61Ao1EyADJWs5kHibrzLp13qmLrVCmAAZKJmjq+Jpvhr81/pm14VJpo5cc1Uz\nVyAQSOr6/IBbpXt0VSb9WnfKyLVwOKxAIJC06zFpEchQwWCjqqvr+zRPLUhZLSHTZsCn89/WVLK+\nOwkTACnnhu1t3Yow6YcwARKXzvkQTvq1nkkIk34IEyAxbtp6FiPHtr0AjGTqfAikBmECZKhMGmGF\n1CNMgAyVqfMhkBqECZChnDIfAs5ABzyQwdI5wioTV9J1AkZz9UOYAPbFyDH7YjRXDCynArdy+oq7\njByzn2Qvp+KqhR6T+Q8D2IUbVtxl5Jj9+Hw++Xw+VVZWJuV6rqqZAG7khl/1jBxzP8IEsDk3/Kpn\n5Jj7uaqZC3CjoX7VO2WEVKbtAZ+JGM0F2NxgK+5+7Wt5+vnPOxghBSMMDe6HMIGbxZoPUlV1OKP2\nCEFqECb9ECZwAtNmqb7n//GPb+gf//hfSdHnz58fUDgcsLyscIZkfXfSZwKkiekQ31jnSzc6tf97\nfjJGSLlhODLSi9FcQJqYDvGNdb70pKT/np/ICKmhJkK6YTgy0ouaCZAmpkN8Bzs/N/ei7rorkNAI\nqeFqHm4Yjoz0IkyANDGduDfY+XPn/o9CoUBCZRm85rFZxcX3MckQCaOZC0gT04l7yZz4N1zNg0mG\nSBQ1EyBNTCfuJXPi33A1DyYZIlGuGhpcUVHRu3gZgMENNhFy2zYCI1OEw2GFw2FVVlYyz6Qv5pkA\niUnnxliwLyYt9kOYAEDi2BwLAGAbdMADFmCpErgNYQKkUKzQkMRSJXAd+kyAFIk9YmqTxo//h06c\n2DHg9az2Cyuw0CNgc4PNMs/NXRnz9VYvVULTG0wQJkCKDDbLXOqKedTKpUpYJRimGM0FpMhgs8wn\nTx5nu6VKWCUYpqiZwHXs0lxTVrZQra2bBswyf+KJFZLstVQJqwTDFGECV7FTc81w61vZqfmIVYJh\nitFccJWionL2RR8B1urKXIzmAmKguWZkWCUYpggTuArNNSNXXHwf4YERYzQXXGXgpk6N+tCHlqmj\n49qAfc4BJA81E7hK3+aav/zlks6dy9L16zU6dUo6dYq5E0Cq0AEP1xqqM97vX2CL4cOA1eiAB4Yx\nWGd8R8c12wwfdgO7zOuBtQgTuNZgnfF/+9vfdPVqTdSxD2Z7b+ZLMEF2mtcDa7mqAz4QCCgcDltd\nDNjEwM546ROfeFTXr3fGfD3DhxPHMizOFQ6HFQgEknY9V9VMkvkPA+frP3fi2rXL+utfO/XuuzNi\nvt4Ow4ed1mTEvB7n8vl88vl8qqysTMr1XBUmQH99504UFZXr9de3SGqUtElS9Gxvv3+RJWW8wYlN\nRszrwQ2uauYChvLfX9H3SSqStFlSQLm5y22xbIgTm4xiNSVavQIyrEHNBBkj+lf0ff//kObOtUfH\nuxObjFiGBTcQJsgYgy0Jb3Xz1g1ObTJiGRZIhAkyiFW/ouPtVLd72AFDYQY8kEKxl3bfpG3bimIG\nSjDYqOrq+j5ht4Bf/UipZH13EiZACrG/CuwuWd+djOYCUsiJnerASBAmQAo5tVMdSBRhAqQQ8zCQ\nKegzAVKMTnXYGR3w/RAmAJA4OuABALZBmAAAjBEmAABjhAkAwBhhAgAwRpgAAIwRJgAAY4QJAMAY\nYQIAMEaYAACMDbvT4rPPPqusrKy4LhaJRPTII4+YlgkA4DDDhklubq4eeOCBuC5WW1trXCAAgPPE\ntdDjlStXdMstt8T827lz53THHXckvWCJYqFHAEhcWhd6LC8vj3n8nXfe0WOPPWZcCACAs8UVJi0t\nLaqrq4s6tnfvXs2cOVMnTpxIScEAAM4R934mFy5c0B/+8AfdfvvtevTRR3X77bdr69atkiSPx5PS\nQkYiEZWXl+vatWuaPXu2VqxYMeA1NHMBQOKS9d05bAe8JHV0dOhTn/qUcnJytHz5cn3ve99TYWFh\n799S7de//rU6Ojp0yy23KC8vL+XvBwBITFzNXIFAQI2NjTp58qS+9a1v6cCBA2psbFRDQ4MqKytH\n9MarV6/WxIkTNWPGjKjjoVBIU6dO1ZQpU3prPmfOnNG9996rH/zgB9q5c+eI3g8AkDpx1UyOHTum\n7u7uqGO7d++WJL366qsjeuNVq1bJ7/dHNVn19PSotLRUL774oiZNmqQ5c+ZoyZIlysvL08033yxJ\nuukm5llmqmCwUVVVh9XVNVpjx3arrGwhe6kDNhFXmPzwhz/UggULYv6tvr5+RG9cWFiotra2qGNN\nTU3Kz8/X5MmTJUklJSWqra3Vhg0b5Pf79fvf/14+n2/QawYCgd7/7fP5hnwtnCUYbNSGDb9Va+uT\nvcdaWzdJkiWBQrDBqcLhsMLhcPIvHBnGq6++OtxLRvTaSCQSOX/+fGT69Om9zw8cOBBZs2ZN7/O9\ne/dGSktL47pWHB8FDrZw4aaIFBnwKCoqT3tZDh1qiHi9G6PK4fVujBw61JD2sgCmkvXdOWzN5IUX\nXtCpU6fiCqb29nbNmjVrxMEW77ItyDxdXbFv1c7OUWkuiVRVdTiqhiRJra1Pqrp6M7UTZKxhw2Sw\nCYupMGnSJLW3t/c+b29vZ/QWJEljx3bHPJ6d3ZPmktgr2AC7sFVv9uzZs/Xmm2+qra1N7733nmpq\narRkyRKriwUbKCtbKK93U9Qxr3ej/P7YfXmpZKdgA+wirg74VFi+fLkaGhp09epVeTwePf7441q1\napW2b9+uoqIi9fT06Otf/7qmTZsW9zUDgQAd7y51o/mounqzOjtHKTu7R37/IkualcrKFqq1dVNU\nU9cHwbYo7WUBRirZHfFxz4C3O2bAI52CwUZVV9f3CbYF9JfAkZL13RlXmKxfv1733HOPPvOZzxh1\nsKcSYQIAiUvrciq33nprzPWwJOnkyZMDZrEDADJLXB3wEydOHPRv/VcTBgBknriauXJycgbdHOvK\nlSu6du1a0guWKJq5ACBxaW3mevrpp7V27dqYf/v+979vXIhkYTQXAMTHktFcjz32mJ5++umYf2tr\na+tdS8tK1EwAIHFp3bZ3/Pjx2rdvX8wVgu0QJAAAazHPBAAyWFprJgAADIUwAQAYI0wAAMZcFSaB\nQCA1O4gBgMuEw+Go3WlN0QEPABmMDngAgG0QJgAAY4QJAMAYYQIAMGbZtr1AX8Fgo6qqDqura7TG\nju1WWdlCdi4EHMRVYcKqwc4UDDZqw4bfRu2p3tq6SZIIFCBF2AN+EAwNdq6ionIdPrwlxvHNCoWe\noNYCpFBa9zMBUqmrK/Zt2Nk5iloL4BB0wMNyY8d2xzyend2jqqrDUUEiSa2tT6q6uj4dRQMQJ8IE\nlisrWyivd1PUMa93o/z+BUPWWgDYB81csNyN5qrq6s3q7Byl7Owe+f2LVFx8n6qqDsc8Jzu7J51F\nBDAMOuBha7H6TLzejdq2bRF9JkASJOu7kzCB7QWDjaquru9Ta1lAkABJQpj0Q5gAQOJYNTgG9jMB\ngPiwn8kgqJkAQOKomQAAbIMwAQAYI0wAAMYIEwCAMcIEAGCMMAEAGCNMAADGCBMAgDHCBABgzFVh\nwnIqABAfllMZBMupAEDiWE4FAGAbhAkAwBhhAgAwRpgAAIwRJgAAY4QJAMAYYQIAMEaYAACMESYA\nAGOECQDAGGECADBGmAAAjLkqTFg1GADiw6rBg2DVYABIHKsGAwBsgzABABgjTAAAxggTAIAxwgQA\nYIwwAQAYI0wAAMYIEwCAMcIEAGCMMAEAGCNMAADGCBMAgDHCBABgjDABABgjTAAAxggTAIAxwgQA\nYMxVYcK2vQAQH7btHQTb9gJA4ti2FwBgG4QJAMAYYQIAMEaYAACMESYAAGOECQDAGGECADBGmAAA\njBEmAABjhAkAwBhhAgAwRpgAAIwRJgAAY4QJAMAYYQIAMEaYAACMESYAAGOECQDAGGECADBGmAAA\njBEmAABjhAkAwBhhAgAwRpgAAIwRJgAAY4QJAMAYYQIAMOaIMAmHwyosLNT69evV0NBgdXEAAP04\nIkxuuukm5eTkqKurS3l5eVYXB0hYOBy2ughASlkWJqtXr9bEiRM1Y8aMqOOhUEhTp07VlClTtHXr\nVklSYWGhfvOb3+ipp55SRUWFFcUFjBAmcDvLwmTVqlUKhUJRx3p6elRaWqpQKKTTp09r//79amlp\nUVZWliTpox/9qLq6uqwoblpZ8cWTivc0veZIz0/kvHhfO9zrMiUsrPqcbrk/rbg3E33fkbIsTAoL\nC5Wbmxt1rKmpSfn5+Zo8ebLGjBmjkpIS1dbW6vnnn9e6deu0YsUK+f1+i0qcPoSJ2fmESeoQJmbn\nuzlMFLHQ+fPnI9OnT+99fuDAgciaNWt6n+/duzdSWloa17W8Xm9EEg8ePHjwSODh9XqT8n0+WjZy\nozlrJM6ePZvEkgAAEmGr0VyTJk1Se3t77/P29nZGbwGAA9gqTGbPnq0333xTbW1teu+991RTU6Ml\nS5ZYXSwAwDAsC5Ply5dr3rx5OnPmjDwej3bv3q3Ro0dr+/btKioqUkFBgZYtW6Zp06ZZVUQAQJyy\nIpFIxOpCAACczVbNXAAAZ3JdmJw/f15r1qzRgw8+2HustrZWa9euVUlJierr6y0sHTJZrHvzjTfe\n0Pr16/XQQw/ppz/9qYWlQ6aLdX9K0jvvvKM5c+YoGAwOeb5rm7kefPBBHThwIOrYP//5T33729/W\nrl27LCoVEPve/M9//qOSkhL98pe/tKhUwAf6358VFRXKycnRtGnTVFxcPOh5jqiZJLKO11C2bNmi\n0tLSVBUTGSgZ92ZdXZ2Ki4tVUlKSyqIiA5nen/X19SooKNCECROGfS9HhEki63jFEolE9J3vfEeL\nFy/WzJkz01FkZAjTe1OS7r//fr3wwgvas2dPqouLDGN6fzY0NOjYsWPat2+ffvKTn2iohixHhEki\n63i99dZbWrdunU6cONGbuNXV1XrppZf0q1/9Sj/+8Y+t+AhwKdN7s6GhQRs2bNA3v/lNfelLX7Li\nI8DFRnJ/Njc3996fW7Zs0TPPPKOHH35Ya9euHXKVElstp5KIjo4OeTye3ud5eXk6fvy4Pvaxj+lH\nP/pR1GvLyspUVlaW7iIiQyVyb86fP1/z589PdxGRwRK5P29YuXLlsNd1RM0kFpN1vIBU4t6EnaXq\n/nRsmLCOF+yKexN2lqr707FhwjpesCvuTdhZqu5PR4QJ63jBrrg3YWfpvD9dO2kRAJA+jqiZAADs\njTABABgjTAAAxggTAIAxwgQAYIwwAQAYI0wAAMYIEwCAMcIESKMtW7botddeS8q1amtr9eyzzybl\nWoApwgTQBzvK3XPPPQoEAlHHn3/+eXk8Hj333HNJeR+Px6NZs2b1Pt+xY4fGjx+vq1evRr3uoYce\n0iOPPDLkploPPPAAKxTDNggTQNKCBQu0fv16/e53v4s6fv36deXn52vp0qUped+5c+dq8eLFamtr\n6z32+uuv6+2339bmzZtZ0wuOQZgA/2/UqFG68847e2sD//rXvwbsUpdsFy5cUGFhoS5evNh77O23\n39alS5fk9XpT+t5AMjl2p0UgFZYuXaqDBw+qvLxcf/rTn1RQUND7tz//+c/auXOnPvvZz+r69eta\nt26djh49qrfeekudnZ3q7OyU1+vV1atXe59/9atfHfL9IpGIPB5Pb83k6NGjuuOOO3Trrbem8mMC\nSUfNBOjjy1/+sl588UVJH3zR93XnnXfqqaee0sqVK7Vnzx5J0qFDh3TzzTdryZIlmj59uurq6qKe\nx8Pj8ai9vV3vv/++srKy1NzcrLlz5yb3gwEpRpgAfYwZM0Yej0fNzc0aN25c1N+ysrJ05MgRPffc\nc3r33XclSX6/XzU1Nfrc5z6nzs7OAc+H8u9//1u5ubnyeDy6ePGijh07pi984QtqamrS5z//+ZR9\nRiAVCBOgn6VLlyoQCOjuu++Oqp1s3bpVLS0tWrp0qXJycnTx4kXV1dVp165dOnr0qILB4IDnQ3nl\nlVc0a9YsTZgwQefOnVNOTo4kqampSXPmzIl6bUtLi5qbm/XKK6/o/fffT/6HBgwRJoCkUCikZ555\nRgcPHtSiRYt011136dq1a9q5c6fOnj2rgwcPKj8/XxcuXFAoFFJBQYHq6ur097//XTU1NTp06JAW\nL14c9fwrX/nKoO935MgRffe739WhQ4ckSV/84hc1c+ZM7dixQ83NzTpy5Ejva5ubmzVhwgTNnDlT\nc+bMUWNjY8r/PYBEsdMikEZ79uzRypUrEzrn5Zdf1r333qs1a9bo8ccfV1tbm+bNmzfi6wGpQM0E\nsLlx48bp8uXL2rVrl2677TZ1dXVZXSRgAMIESKP29vaEl1O5++67dfnyZb322ms6fvy4CgsLJX2w\nnAoNC7ALmrkAAMaomQAAjBEmAABjhAkAwBhhAgAwRpgAAIwRJgAAY4QJAMAYYQIAMPZ/tdII0gRP\n36EAAAAASUVORK5CYII=\n",
"text": [
""
]
}
],
"prompt_number": 104
},
{
"cell_type": "code",
"collapsed": false,
"input": [],
"language": "python",
"metadata": {},
"outputs": []
}
],
"metadata": {}
}
]
}