notebook.community

Edit and run 





In [24]:



    


%matplotlib inline
import yt
import numpy as np
import matplotlib.pyplot as plt
from galaxy_analysis import Galaxy



    


















    






/home/aemerick/code/galaxy_analysis/particle_analysis/IMF.py:2: UserWarning: 
This call to matplotlib.use() has no effect because the backend has already
been chosen; matplotlib.use() must be called *before* pylab, matplotlib.pyplot,
or matplotlib.backends is imported for the first time.

The backend was *originally* set to 'module://ipykernel.pylab.backend_inline' by the following code:
  File "/home/aemerick/anaconda2/lib/python2.7/runpy.py", line 174, in _run_module_as_main
    "__main__", fname, loader, pkg_name)
  File "/home/aemerick/anaconda2/lib/python2.7/runpy.py", line 72, in _run_code
    exec code in run_globals
  File "/home/aemerick/anaconda2/lib/python2.7/site-packages/ipykernel_launcher.py", line 16, in <module>
    app.launch_new_instance()
  File "/home/aemerick/anaconda2/lib/python2.7/site-packages/traitlets/config/application.py", line 658, in launch_instance
    app.start()
  File "/home/aemerick/anaconda2/lib/python2.7/site-packages/ipykernel/kernelapp.py", line 477, in start
    ioloop.IOLoop.instance().start()
  File "/home/aemerick/anaconda2/lib/python2.7/site-packages/zmq/eventloop/ioloop.py", line 177, in start
    super(ZMQIOLoop, self).start()
  File "/home/aemerick/anaconda2/lib/python2.7/site-packages/tornado/ioloop.py", line 888, in start
    handler_func(fd_obj, events)
  File "/home/aemerick/anaconda2/lib/python2.7/site-packages/tornado/stack_context.py", line 277, in null_wrapper
    return fn(*args, **kwargs)
  File "/home/aemerick/anaconda2/lib/python2.7/site-packages/zmq/eventloop/zmqstream.py", line 440, in _handle_events
    self._handle_recv()
  File "/home/aemerick/anaconda2/lib/python2.7/site-packages/zmq/eventloop/zmqstream.py", line 472, in _handle_recv
    self._run_callback(callback, msg)
  File "/home/aemerick/anaconda2/lib/python2.7/site-packages/zmq/eventloop/zmqstream.py", line 414, in _run_callback
    callback(*args, **kwargs)
  File "/home/aemerick/anaconda2/lib/python2.7/site-packages/tornado/stack_context.py", line 277, in null_wrapper
    return fn(*args, **kwargs)
  File "/home/aemerick/anaconda2/lib/python2.7/site-packages/ipykernel/kernelbase.py", line 283, in dispatcher
    return self.dispatch_shell(stream, msg)
  File "/home/aemerick/anaconda2/lib/python2.7/site-packages/ipykernel/kernelbase.py", line 235, in dispatch_shell
    handler(stream, idents, msg)
  File "/home/aemerick/anaconda2/lib/python2.7/site-packages/ipykernel/kernelbase.py", line 399, in execute_request
    user_expressions, allow_stdin)
  File "/home/aemerick/anaconda2/lib/python2.7/site-packages/ipykernel/ipkernel.py", line 196, in do_execute
    res = shell.run_cell(code, store_history=store_history, silent=silent)
  File "/home/aemerick/anaconda2/lib/python2.7/site-packages/ipykernel/zmqshell.py", line 533, in run_cell
    return super(ZMQInteractiveShell, self).run_cell(*args, **kwargs)
  File "/home/aemerick/anaconda2/lib/python2.7/site-packages/IPython/core/interactiveshell.py", line 2718, in run_cell
    interactivity=interactivity, compiler=compiler, result=result)
  File "/home/aemerick/anaconda2/lib/python2.7/site-packages/IPython/core/interactiveshell.py", line 2822, in run_ast_nodes
    if self.run_code(code, result):
  File "/home/aemerick/anaconda2/lib/python2.7/site-packages/IPython/core/interactiveshell.py", line 2882, in run_code
    exec(code_obj, self.user_global_ns, self.user_ns)
  File "<ipython-input-24-f61cb47bd4fe>", line 1, in <module>
    get_ipython().magic(u'matplotlib inline')
  File "/home/aemerick/anaconda2/lib/python2.7/site-packages/IPython/core/interactiveshell.py", line 2160, in magic
    return self.run_line_magic(magic_name, magic_arg_s)
  File "/home/aemerick/anaconda2/lib/python2.7/site-packages/IPython/core/interactiveshell.py", line 2081, in run_line_magic
    result = fn(*args,**kwargs)
  File "<decorator-gen-105>", line 2, in matplotlib
  File "/home/aemerick/anaconda2/lib/python2.7/site-packages/IPython/core/magic.py", line 188, in <lambda>
    call = lambda f, *a, **k: f(*a, **k)
  File "/home/aemerick/anaconda2/lib/python2.7/site-packages/IPython/core/magics/pylab.py", line 100, in matplotlib
    gui, backend = self.shell.enable_matplotlib(args.gui)
  File "/home/aemerick/anaconda2/lib/python2.7/site-packages/IPython/core/interactiveshell.py", line 2950, in enable_matplotlib
    pt.activate_matplotlib(backend)
  File "/home/aemerick/anaconda2/lib/python2.7/site-packages/IPython/core/pylabtools.py", line 308, in activate_matplotlib
    matplotlib.pyplot.switch_backend(backend)
  File "/home/aemerick/anaconda2/lib/python2.7/site-packages/matplotlib/pyplot.py", line 229, in switch_backend
    matplotlib.use(newbackend, warn=False, force=True)
  File "/home/aemerick/anaconda2/lib/python2.7/site-packages/matplotlib/__init__.py", line 1305, in use
    reload(sys.modules['matplotlib.backends'])
  File "/home/aemerick/anaconda2/lib/python2.7/site-packages/matplotlib/backends/__init__.py", line 14, in <module>
    line for line in traceback.format_stack()


  mpl.use('Agg')
/home/aemerick/code/galaxy_analysis/particle_analysis/sfrFromParticles.py:4: UserWarning: 
This call to matplotlib.use() has no effect because the backend has already
been chosen; matplotlib.use() must be called *before* pylab, matplotlib.pyplot,
or matplotlib.backends is imported for the first time.

The backend was *originally* set to 'module://ipykernel.pylab.backend_inline' by the following code:
  File "/home/aemerick/anaconda2/lib/python2.7/runpy.py", line 174, in _run_module_as_main
    "__main__", fname, loader, pkg_name)
  File "/home/aemerick/anaconda2/lib/python2.7/runpy.py", line 72, in _run_code
    exec code in run_globals
  File "/home/aemerick/anaconda2/lib/python2.7/site-packages/ipykernel_launcher.py", line 16, in <module>
    app.launch_new_instance()
  File "/home/aemerick/anaconda2/lib/python2.7/site-packages/traitlets/config/application.py", line 658, in launch_instance
    app.start()
  File "/home/aemerick/anaconda2/lib/python2.7/site-packages/ipykernel/kernelapp.py", line 477, in start
    ioloop.IOLoop.instance().start()
  File "/home/aemerick/anaconda2/lib/python2.7/site-packages/zmq/eventloop/ioloop.py", line 177, in start
    super(ZMQIOLoop, self).start()
  File "/home/aemerick/anaconda2/lib/python2.7/site-packages/tornado/ioloop.py", line 888, in start
    handler_func(fd_obj, events)
  File "/home/aemerick/anaconda2/lib/python2.7/site-packages/tornado/stack_context.py", line 277, in null_wrapper
    return fn(*args, **kwargs)
  File "/home/aemerick/anaconda2/lib/python2.7/site-packages/zmq/eventloop/zmqstream.py", line 440, in _handle_events
    self._handle_recv()
  File "/home/aemerick/anaconda2/lib/python2.7/site-packages/zmq/eventloop/zmqstream.py", line 472, in _handle_recv
    self._run_callback(callback, msg)
  File "/home/aemerick/anaconda2/lib/python2.7/site-packages/zmq/eventloop/zmqstream.py", line 414, in _run_callback
    callback(*args, **kwargs)
  File "/home/aemerick/anaconda2/lib/python2.7/site-packages/tornado/stack_context.py", line 277, in null_wrapper
    return fn(*args, **kwargs)
  File "/home/aemerick/anaconda2/lib/python2.7/site-packages/ipykernel/kernelbase.py", line 283, in dispatcher
    return self.dispatch_shell(stream, msg)
  File "/home/aemerick/anaconda2/lib/python2.7/site-packages/ipykernel/kernelbase.py", line 235, in dispatch_shell
    handler(stream, idents, msg)
  File "/home/aemerick/anaconda2/lib/python2.7/site-packages/ipykernel/kernelbase.py", line 399, in execute_request
    user_expressions, allow_stdin)
  File "/home/aemerick/anaconda2/lib/python2.7/site-packages/ipykernel/ipkernel.py", line 196, in do_execute
    res = shell.run_cell(code, store_history=store_history, silent=silent)
  File "/home/aemerick/anaconda2/lib/python2.7/site-packages/ipykernel/zmqshell.py", line 533, in run_cell
    return super(ZMQInteractiveShell, self).run_cell(*args, **kwargs)
  File "/home/aemerick/anaconda2/lib/python2.7/site-packages/IPython/core/interactiveshell.py", line 2718, in run_cell
    interactivity=interactivity, compiler=compiler, result=result)
  File "/home/aemerick/anaconda2/lib/python2.7/site-packages/IPython/core/interactiveshell.py", line 2822, in run_ast_nodes
    if self.run_code(code, result):
  File "/home/aemerick/anaconda2/lib/python2.7/site-packages/IPython/core/interactiveshell.py", line 2882, in run_code
    exec(code_obj, self.user_global_ns, self.user_ns)
  File "<ipython-input-24-f61cb47bd4fe>", line 1, in <module>
    get_ipython().magic(u'matplotlib inline')
  File "/home/aemerick/anaconda2/lib/python2.7/site-packages/IPython/core/interactiveshell.py", line 2160, in magic
    return self.run_line_magic(magic_name, magic_arg_s)
  File "/home/aemerick/anaconda2/lib/python2.7/site-packages/IPython/core/interactiveshell.py", line 2081, in run_line_magic
    result = fn(*args,**kwargs)
  File "<decorator-gen-105>", line 2, in matplotlib
  File "/home/aemerick/anaconda2/lib/python2.7/site-packages/IPython/core/magic.py", line 188, in <lambda>
    call = lambda f, *a, **k: f(*a, **k)
  File "/home/aemerick/anaconda2/lib/python2.7/site-packages/IPython/core/magics/pylab.py", line 100, in matplotlib
    gui, backend = self.shell.enable_matplotlib(args.gui)
  File "/home/aemerick/anaconda2/lib/python2.7/site-packages/IPython/core/interactiveshell.py", line 2950, in enable_matplotlib
    pt.activate_matplotlib(backend)
  File "/home/aemerick/anaconda2/lib/python2.7/site-packages/IPython/core/pylabtools.py", line 308, in activate_matplotlib
    matplotlib.pyplot.switch_backend(backend)
  File "/home/aemerick/anaconda2/lib/python2.7/site-packages/matplotlib/pyplot.py", line 229, in switch_backend
    matplotlib.use(newbackend, warn=False, force=True)
  File "/home/aemerick/anaconda2/lib/python2.7/site-packages/matplotlib/__init__.py", line 1305, in use
    reload(sys.modules['matplotlib.backends'])
  File "/home/aemerick/anaconda2/lib/python2.7/site-packages/matplotlib/backends/__init__.py", line 14, in <module>
    line for line in traceback.format_stack()


  mpl.use('Agg')
/home/aemerick/code/galaxy_analysis/particle_analysis/sn_rate.py:4: UserWarning: 
This call to matplotlib.use() has no effect because the backend has already
been chosen; matplotlib.use() must be called *before* pylab, matplotlib.pyplot,
or matplotlib.backends is imported for the first time.

The backend was *originally* set to 'module://ipykernel.pylab.backend_inline' by the following code:
  File "/home/aemerick/anaconda2/lib/python2.7/runpy.py", line 174, in _run_module_as_main
    "__main__", fname, loader, pkg_name)
  File "/home/aemerick/anaconda2/lib/python2.7/runpy.py", line 72, in _run_code
    exec code in run_globals
  File "/home/aemerick/anaconda2/lib/python2.7/site-packages/ipykernel_launcher.py", line 16, in <module>
    app.launch_new_instance()
  File "/home/aemerick/anaconda2/lib/python2.7/site-packages/traitlets/config/application.py", line 658, in launch_instance
    app.start()
  File "/home/aemerick/anaconda2/lib/python2.7/site-packages/ipykernel/kernelapp.py", line 477, in start
    ioloop.IOLoop.instance().start()
  File "/home/aemerick/anaconda2/lib/python2.7/site-packages/zmq/eventloop/ioloop.py", line 177, in start
    super(ZMQIOLoop, self).start()
  File "/home/aemerick/anaconda2/lib/python2.7/site-packages/tornado/ioloop.py", line 888, in start
    handler_func(fd_obj, events)
  File "/home/aemerick/anaconda2/lib/python2.7/site-packages/tornado/stack_context.py", line 277, in null_wrapper
    return fn(*args, **kwargs)
  File "/home/aemerick/anaconda2/lib/python2.7/site-packages/zmq/eventloop/zmqstream.py", line 440, in _handle_events
    self._handle_recv()
  File "/home/aemerick/anaconda2/lib/python2.7/site-packages/zmq/eventloop/zmqstream.py", line 472, in _handle_recv
    self._run_callback(callback, msg)
  File "/home/aemerick/anaconda2/lib/python2.7/site-packages/zmq/eventloop/zmqstream.py", line 414, in _run_callback
    callback(*args, **kwargs)
  File "/home/aemerick/anaconda2/lib/python2.7/site-packages/tornado/stack_context.py", line 277, in null_wrapper
    return fn(*args, **kwargs)
  File "/home/aemerick/anaconda2/lib/python2.7/site-packages/ipykernel/kernelbase.py", line 283, in dispatcher
    return self.dispatch_shell(stream, msg)
  File "/home/aemerick/anaconda2/lib/python2.7/site-packages/ipykernel/kernelbase.py", line 235, in dispatch_shell
    handler(stream, idents, msg)
  File "/home/aemerick/anaconda2/lib/python2.7/site-packages/ipykernel/kernelbase.py", line 399, in execute_request
    user_expressions, allow_stdin)
  File "/home/aemerick/anaconda2/lib/python2.7/site-packages/ipykernel/ipkernel.py", line 196, in do_execute
    res = shell.run_cell(code, store_history=store_history, silent=silent)
  File "/home/aemerick/anaconda2/lib/python2.7/site-packages/ipykernel/zmqshell.py", line 533, in run_cell
    return super(ZMQInteractiveShell, self).run_cell(*args, **kwargs)
  File "/home/aemerick/anaconda2/lib/python2.7/site-packages/IPython/core/interactiveshell.py", line 2718, in run_cell
    interactivity=interactivity, compiler=compiler, result=result)
  File "/home/aemerick/anaconda2/lib/python2.7/site-packages/IPython/core/interactiveshell.py", line 2822, in run_ast_nodes
    if self.run_code(code, result):
  File "/home/aemerick/anaconda2/lib/python2.7/site-packages/IPython/core/interactiveshell.py", line 2882, in run_code
    exec(code_obj, self.user_global_ns, self.user_ns)
  File "<ipython-input-24-f61cb47bd4fe>", line 1, in <module>
    get_ipython().magic(u'matplotlib inline')
  File "/home/aemerick/anaconda2/lib/python2.7/site-packages/IPython/core/interactiveshell.py", line 2160, in magic
    return self.run_line_magic(magic_name, magic_arg_s)
  File "/home/aemerick/anaconda2/lib/python2.7/site-packages/IPython/core/interactiveshell.py", line 2081, in run_line_magic
    result = fn(*args,**kwargs)
  File "<decorator-gen-105>", line 2, in matplotlib
  File "/home/aemerick/anaconda2/lib/python2.7/site-packages/IPython/core/magic.py", line 188, in <lambda>
    call = lambda f, *a, **k: f(*a, **k)
  File "/home/aemerick/anaconda2/lib/python2.7/site-packages/IPython/core/magics/pylab.py", line 100, in matplotlib
    gui, backend = self.shell.enable_matplotlib(args.gui)
  File "/home/aemerick/anaconda2/lib/python2.7/site-packages/IPython/core/interactiveshell.py", line 2950, in enable_matplotlib
    pt.activate_matplotlib(backend)
  File "/home/aemerick/anaconda2/lib/python2.7/site-packages/IPython/core/pylabtools.py", line 308, in activate_matplotlib
    matplotlib.pyplot.switch_backend(backend)
  File "/home/aemerick/anaconda2/lib/python2.7/site-packages/matplotlib/pyplot.py", line 229, in switch_backend
    matplotlib.use(newbackend, warn=False, force=True)
  File "/home/aemerick/anaconda2/lib/python2.7/site-packages/matplotlib/__init__.py", line 1305, in use
    reload(sys.modules['matplotlib.backends'])
  File "/home/aemerick/anaconda2/lib/python2.7/site-packages/matplotlib/backends/__init__.py", line 14, in <module>
    line for line in traceback.format_stack()


  mpl.use('Agg')
/home/aemerick/code/galaxy_analysis/particle_analysis/sfhFromParticles.py:4: UserWarning: 
This call to matplotlib.use() has no effect because the backend has already
been chosen; matplotlib.use() must be called *before* pylab, matplotlib.pyplot,
or matplotlib.backends is imported for the first time.

The backend was *originally* set to 'module://ipykernel.pylab.backend_inline' by the following code:
  File "/home/aemerick/anaconda2/lib/python2.7/runpy.py", line 174, in _run_module_as_main
    "__main__", fname, loader, pkg_name)
  File "/home/aemerick/anaconda2/lib/python2.7/runpy.py", line 72, in _run_code
    exec code in run_globals
  File "/home/aemerick/anaconda2/lib/python2.7/site-packages/ipykernel_launcher.py", line 16, in <module>
    app.launch_new_instance()
  File "/home/aemerick/anaconda2/lib/python2.7/site-packages/traitlets/config/application.py", line 658, in launch_instance
    app.start()
  File "/home/aemerick/anaconda2/lib/python2.7/site-packages/ipykernel/kernelapp.py", line 477, in start
    ioloop.IOLoop.instance().start()
  File "/home/aemerick/anaconda2/lib/python2.7/site-packages/zmq/eventloop/ioloop.py", line 177, in start
    super(ZMQIOLoop, self).start()
  File "/home/aemerick/anaconda2/lib/python2.7/site-packages/tornado/ioloop.py", line 888, in start
    handler_func(fd_obj, events)
  File "/home/aemerick/anaconda2/lib/python2.7/site-packages/tornado/stack_context.py", line 277, in null_wrapper
    return fn(*args, **kwargs)
  File "/home/aemerick/anaconda2/lib/python2.7/site-packages/zmq/eventloop/zmqstream.py", line 440, in _handle_events
    self._handle_recv()
  File "/home/aemerick/anaconda2/lib/python2.7/site-packages/zmq/eventloop/zmqstream.py", line 472, in _handle_recv
    self._run_callback(callback, msg)
  File "/home/aemerick/anaconda2/lib/python2.7/site-packages/zmq/eventloop/zmqstream.py", line 414, in _run_callback
    callback(*args, **kwargs)
  File "/home/aemerick/anaconda2/lib/python2.7/site-packages/tornado/stack_context.py", line 277, in null_wrapper
    return fn(*args, **kwargs)
  File "/home/aemerick/anaconda2/lib/python2.7/site-packages/ipykernel/kernelbase.py", line 283, in dispatcher
    return self.dispatch_shell(stream, msg)
  File "/home/aemerick/anaconda2/lib/python2.7/site-packages/ipykernel/kernelbase.py", line 235, in dispatch_shell
    handler(stream, idents, msg)
  File "/home/aemerick/anaconda2/lib/python2.7/site-packages/ipykernel/kernelbase.py", line 399, in execute_request
    user_expressions, allow_stdin)
  File "/home/aemerick/anaconda2/lib/python2.7/site-packages/ipykernel/ipkernel.py", line 196, in do_execute
    res = shell.run_cell(code, store_history=store_history, silent=silent)
  File "/home/aemerick/anaconda2/lib/python2.7/site-packages/ipykernel/zmqshell.py", line 533, in run_cell
    return super(ZMQInteractiveShell, self).run_cell(*args, **kwargs)
  File "/home/aemerick/anaconda2/lib/python2.7/site-packages/IPython/core/interactiveshell.py", line 2718, in run_cell
    interactivity=interactivity, compiler=compiler, result=result)
  File "/home/aemerick/anaconda2/lib/python2.7/site-packages/IPython/core/interactiveshell.py", line 2822, in run_ast_nodes
    if self.run_code(code, result):
  File "/home/aemerick/anaconda2/lib/python2.7/site-packages/IPython/core/interactiveshell.py", line 2882, in run_code
    exec(code_obj, self.user_global_ns, self.user_ns)
  File "<ipython-input-24-f61cb47bd4fe>", line 1, in <module>
    get_ipython().magic(u'matplotlib inline')
  File "/home/aemerick/anaconda2/lib/python2.7/site-packages/IPython/core/interactiveshell.py", line 2160, in magic
    return self.run_line_magic(magic_name, magic_arg_s)
  File "/home/aemerick/anaconda2/lib/python2.7/site-packages/IPython/core/interactiveshell.py", line 2081, in run_line_magic
    result = fn(*args,**kwargs)
  File "<decorator-gen-105>", line 2, in matplotlib
  File "/home/aemerick/anaconda2/lib/python2.7/site-packages/IPython/core/magic.py", line 188, in <lambda>
    call = lambda f, *a, **k: f(*a, **k)
  File "/home/aemerick/anaconda2/lib/python2.7/site-packages/IPython/core/magics/pylab.py", line 100, in matplotlib
    gui, backend = self.shell.enable_matplotlib(args.gui)
  File "/home/aemerick/anaconda2/lib/python2.7/site-packages/IPython/core/interactiveshell.py", line 2950, in enable_matplotlib
    pt.activate_matplotlib(backend)
  File "/home/aemerick/anaconda2/lib/python2.7/site-packages/IPython/core/pylabtools.py", line 308, in activate_matplotlib
    matplotlib.pyplot.switch_backend(backend)
  File "/home/aemerick/anaconda2/lib/python2.7/site-packages/matplotlib/pyplot.py", line 229, in switch_backend
    matplotlib.use(newbackend, warn=False, force=True)
  File "/home/aemerick/anaconda2/lib/python2.7/site-packages/matplotlib/__init__.py", line 1305, in use
    reload(sys.modules['matplotlib.backends'])
  File "/home/aemerick/anaconda2/lib/python2.7/site-packages/matplotlib/backends/__init__.py", line 14, in <module>
    line for line in traceback.format_stack()


  mpl.use('Agg')
/home/aemerick/code/galaxy_analysis/plot/plot_styles.py:7: UserWarning: 
This call to matplotlib.use() has no effect because the backend has already
been chosen; matplotlib.use() must be called *before* pylab, matplotlib.pyplot,
or matplotlib.backends is imported for the first time.

The backend was *originally* set to 'module://ipykernel.pylab.backend_inline' by the following code:
  File "/home/aemerick/anaconda2/lib/python2.7/runpy.py", line 174, in _run_module_as_main
    "__main__", fname, loader, pkg_name)
  File "/home/aemerick/anaconda2/lib/python2.7/runpy.py", line 72, in _run_code
    exec code in run_globals
  File "/home/aemerick/anaconda2/lib/python2.7/site-packages/ipykernel_launcher.py", line 16, in <module>
    app.launch_new_instance()
  File "/home/aemerick/anaconda2/lib/python2.7/site-packages/traitlets/config/application.py", line 658, in launch_instance
    app.start()
  File "/home/aemerick/anaconda2/lib/python2.7/site-packages/ipykernel/kernelapp.py", line 477, in start
    ioloop.IOLoop.instance().start()
  File "/home/aemerick/anaconda2/lib/python2.7/site-packages/zmq/eventloop/ioloop.py", line 177, in start
    super(ZMQIOLoop, self).start()
  File "/home/aemerick/anaconda2/lib/python2.7/site-packages/tornado/ioloop.py", line 888, in start
    handler_func(fd_obj, events)
  File "/home/aemerick/anaconda2/lib/python2.7/site-packages/tornado/stack_context.py", line 277, in null_wrapper
    return fn(*args, **kwargs)
  File "/home/aemerick/anaconda2/lib/python2.7/site-packages/zmq/eventloop/zmqstream.py", line 440, in _handle_events
    self._handle_recv()
  File "/home/aemerick/anaconda2/lib/python2.7/site-packages/zmq/eventloop/zmqstream.py", line 472, in _handle_recv
    self._run_callback(callback, msg)
  File "/home/aemerick/anaconda2/lib/python2.7/site-packages/zmq/eventloop/zmqstream.py", line 414, in _run_callback
    callback(*args, **kwargs)
  File "/home/aemerick/anaconda2/lib/python2.7/site-packages/tornado/stack_context.py", line 277, in null_wrapper
    return fn(*args, **kwargs)
  File "/home/aemerick/anaconda2/lib/python2.7/site-packages/ipykernel/kernelbase.py", line 283, in dispatcher
    return self.dispatch_shell(stream, msg)
  File "/home/aemerick/anaconda2/lib/python2.7/site-packages/ipykernel/kernelbase.py", line 235, in dispatch_shell
    handler(stream, idents, msg)
  File "/home/aemerick/anaconda2/lib/python2.7/site-packages/ipykernel/kernelbase.py", line 399, in execute_request
    user_expressions, allow_stdin)
  File "/home/aemerick/anaconda2/lib/python2.7/site-packages/ipykernel/ipkernel.py", line 196, in do_execute
    res = shell.run_cell(code, store_history=store_history, silent=silent)
  File "/home/aemerick/anaconda2/lib/python2.7/site-packages/ipykernel/zmqshell.py", line 533, in run_cell
    return super(ZMQInteractiveShell, self).run_cell(*args, **kwargs)
  File "/home/aemerick/anaconda2/lib/python2.7/site-packages/IPython/core/interactiveshell.py", line 2718, in run_cell
    interactivity=interactivity, compiler=compiler, result=result)
  File "/home/aemerick/anaconda2/lib/python2.7/site-packages/IPython/core/interactiveshell.py", line 2822, in run_ast_nodes
    if self.run_code(code, result):
  File "/home/aemerick/anaconda2/lib/python2.7/site-packages/IPython/core/interactiveshell.py", line 2882, in run_code
    exec(code_obj, self.user_global_ns, self.user_ns)
  File "<ipython-input-24-f61cb47bd4fe>", line 1, in <module>
    get_ipython().magic(u'matplotlib inline')
  File "/home/aemerick/anaconda2/lib/python2.7/site-packages/IPython/core/interactiveshell.py", line 2160, in magic
    return self.run_line_magic(magic_name, magic_arg_s)
  File "/home/aemerick/anaconda2/lib/python2.7/site-packages/IPython/core/interactiveshell.py", line 2081, in run_line_magic
    result = fn(*args,**kwargs)
  File "<decorator-gen-105>", line 2, in matplotlib
  File "/home/aemerick/anaconda2/lib/python2.7/site-packages/IPython/core/magic.py", line 188, in <lambda>
    call = lambda f, *a, **k: f(*a, **k)
  File "/home/aemerick/anaconda2/lib/python2.7/site-packages/IPython/core/magics/pylab.py", line 100, in matplotlib
    gui, backend = self.shell.enable_matplotlib(args.gui)
  File "/home/aemerick/anaconda2/lib/python2.7/site-packages/IPython/core/interactiveshell.py", line 2950, in enable_matplotlib
    pt.activate_matplotlib(backend)
  File "/home/aemerick/anaconda2/lib/python2.7/site-packages/IPython/core/pylabtools.py", line 308, in activate_matplotlib
    matplotlib.pyplot.switch_backend(backend)
  File "/home/aemerick/anaconda2/lib/python2.7/site-packages/matplotlib/pyplot.py", line 229, in switch_backend
    matplotlib.use(newbackend, warn=False, force=True)
  File "/home/aemerick/anaconda2/lib/python2.7/site-packages/matplotlib/__init__.py", line 1305, in use
    reload(sys.modules['matplotlib.backends'])
  File "/home/aemerick/anaconda2/lib/python2.7/site-packages/matplotlib/backends/__init__.py", line 14, in <module>
    line for line in traceback.format_stack()


  mpl.use('Agg')

















In [79]:



    


colors   = {'massive_star_winds' : 'black', 'AGB_winds' : 'C1', 'SN' : 'C4', 'other_stars' : 'black'}
markers  = {'massive_star_winds' : '*',     'AGB_winds' : 'D', 'SN' : '*', 'other_stars' : '.'}
ps       = {'massive_star_winds' : 20, 'AGB_winds' : 100, 'SN' : 40, 'other_stars' : 10}



    











In [26]:



    


#ds = yt.load('./../example_data/DD0401/DD0401')
#data = ds.all_data()

gal = Galaxy('DD0401', wdir = './../example_data/')



    


















    






Parsing Hierarchy : 100%|██████████| 3088/3088 [00:00<00:00, 46370.46it/s]
/home/aemerick/anaconda2/lib/python2.7/site-packages/yt/fields/local_fields.py:46: UserWarning: Because 'sampling_type' not specified, yt will assume a cell 'sampling_type'
  warnings.warn("Because 'sampling_type' not specified, yt will "










    






tracer species present:  ['C', 'N', 'O', 'Na', 'Mg', 'Si', 'S', 'Ca', 'Mn', 'Fe', 'Ni', 'As', 'Sr', 'Y', 'Ba']
16 mass fields defined
16 mass fraction fields defined
15 number density fields defined
141 abundance ratio fields defined
141 particle abundance ratio fields defined
5 additional helper fields defined










    






Parsing Hierarchy : 100%|██████████| 3088/3088 [00:00<00:00, 44787.36it/s]
/home/aemerick/code/onezone/data_tables.py:89: FutureWarning: elementwise comparison failed; returning scalar instead, but in the future will perform elementwise comparison
  if  c == flag or id == flag:

















In [27]:



    


data = gal.df



    











In [73]:



    


gal.ds.derived_field_list



    


















    
Out[73]:








[('all', 'age'),
 ('all', 'birth_mass'),
 ('all', 'creation_time'),
 ('all', 'dynamical_time'),
 ('all', 'mesh_id'),
 ('all', 'metallicity_fraction'),
 ('all', 'particle_As_fraction'),
 ('all', 'particle_Ba_fraction'),
 ('all', 'particle_C_fraction'),
 ('all', 'particle_Ca_fraction'),
 ('all', 'particle_Fe_fraction'),
 ('all', 'particle_H_fraction'),
 ('all', 'particle_He_fraction'),
 ('all', 'particle_Mg_fraction'),
 ('all', 'particle_Mn_fraction'),
 ('all', 'particle_N_fraction'),
 ('all', 'particle_Na_fraction'),
 ('all', 'particle_Ni_fraction'),
 ('all', 'particle_O_fraction'),
 ('all', 'particle_S_fraction'),
 ('all', 'particle_Si_fraction'),
 ('all', 'particle_Sr_fraction'),
 ('all', 'particle_Y_fraction'),
 ('all', 'particle_angular_momentum'),
 ('all', 'particle_angular_momentum_magnitude'),
 ('all', 'particle_angular_momentum_x'),
 ('all', 'particle_angular_momentum_y'),
 ('all', 'particle_angular_momentum_z'),
 ('all', 'particle_cylindrical_velocity_theta'),
 ('all', 'particle_cylindrical_velocity_z'),
 ('all', 'particle_index'),
 ('all', 'particle_mass'),
 ('all', 'particle_ones'),
 ('all', 'particle_position'),
 ('all', 'particle_position_cylindrical_radius'),
 ('all', 'particle_position_cylindrical_theta'),
 ('all', 'particle_position_cylindrical_z'),
 ('all', 'particle_position_relative'),
 ('all', 'particle_position_relative_x'),
 ('all', 'particle_position_relative_y'),
 ('all', 'particle_position_relative_z'),
 ('all', 'particle_position_spherical_phi'),
 ('all', 'particle_position_spherical_radius'),
 ('all', 'particle_position_spherical_theta'),
 ('all', 'particle_position_x'),
 ('all', 'particle_position_y'),
 ('all', 'particle_position_z'),
 ('all', 'particle_radial_velocity'),
 ('all', 'particle_radius'),
 ('all', 'particle_specific_angular_momentum'),
 ('all', 'particle_specific_angular_momentum_x'),
 ('all', 'particle_specific_angular_momentum_y'),
 ('all', 'particle_specific_angular_momentum_z'),
 ('all', 'particle_spherical_position_phi'),
 ('all', 'particle_spherical_position_radius'),
 ('all', 'particle_spherical_position_theta'),
 ('all', 'particle_spherical_velocity_phi'),
 ('all', 'particle_spherical_velocity_radius'),
 ('all', 'particle_spherical_velocity_theta'),
 ('all', 'particle_type'),
 ('all', 'particle_velocity'),
 ('all', 'particle_velocity_cylindrical_radius'),
 ('all', 'particle_velocity_cylindrical_theta'),
 ('all', 'particle_velocity_cylindrical_z'),
 ('all', 'particle_velocity_magnitude'),
 ('all', 'particle_velocity_relative'),
 ('all', 'particle_velocity_relative_x'),
 ('all', 'particle_velocity_relative_y'),
 ('all', 'particle_velocity_relative_z'),
 ('all', 'particle_velocity_spherical_phi'),
 ('all', 'particle_velocity_spherical_radius'),
 ('all', 'particle_velocity_spherical_theta'),
 ('all', 'particle_velocity_x'),
 ('all', 'particle_velocity_y'),
 ('all', 'particle_velocity_z'),
 ('all', 'rad_table_T_pos'),
 ('all', 'rad_table_Z_pos'),
 ('all', 'rad_table_g_pos'),
 ('all', 'se_table_M_pos'),
 ('all', 'se_table_Z_pos'),
 ('all', 'sn_mass_ejected'),
 ('all', 'wind_mass_ejected'),
 ('all', 'yield_table_M_pos'),
 ('all', 'yield_table_Z_pos'),
 ('deposit', 'all_cic'),
 ('deposit', 'all_cic_age'),
 ('deposit', 'all_cic_velocity_x'),
 ('deposit', 'all_cic_velocity_y'),
 ('deposit', 'all_cic_velocity_z'),
 ('deposit', 'all_count'),
 ('deposit', 'all_density'),
 ('deposit', 'all_mass'),
 ('deposit', 'all_nn_age'),
 ('deposit', 'all_nn_velocity_x'),
 ('deposit', 'all_nn_velocity_y'),
 ('deposit', 'all_nn_velocity_z'),
 ('deposit', 'io_cic'),
 ('deposit', 'io_cic_age'),
 ('deposit', 'io_cic_velocity_x'),
 ('deposit', 'io_cic_velocity_y'),
 ('deposit', 'io_cic_velocity_z'),
 ('deposit', 'io_count'),
 ('deposit', 'io_density'),
 ('deposit', 'io_mass'),
 ('deposit', 'io_nn_age'),
 ('deposit', 'io_nn_velocity_x'),
 ('deposit', 'io_nn_velocity_y'),
 ('deposit', 'io_nn_velocity_z'),
 ('deposit', 'main_sequence_stars_cic'),
 ('deposit', 'main_sequence_stars_cic_age'),
 ('deposit', 'main_sequence_stars_cic_velocity_x'),
 ('deposit', 'main_sequence_stars_cic_velocity_y'),
 ('deposit', 'main_sequence_stars_cic_velocity_z'),
 ('deposit', 'main_sequence_stars_count'),
 ('deposit', 'main_sequence_stars_density'),
 ('deposit', 'main_sequence_stars_mass'),
 ('deposit', 'main_sequence_stars_nn_age'),
 ('deposit', 'main_sequence_stars_nn_velocity_x'),
 ('deposit', 'main_sequence_stars_nn_velocity_y'),
 ('deposit', 'main_sequence_stars_nn_velocity_z'),
 ('enzo', 'As_Density'),
 ('enzo', 'Ba_Density'),
 ('enzo', 'C_Density'),
 ('enzo', 'Ca_Density'),
 ('enzo', 'Dark_Matter_Density'),
 ('enzo', 'Density'),
 ('enzo', 'Electron_Density'),
 ('enzo', 'Fe_Density'),
 ('enzo', 'GasEnergy'),
 ('enzo', 'GravPotential'),
 ('enzo', 'H2II_Density'),
 ('enzo', 'H2I_Density'),
 ('enzo', 'H2I_kdiss'),
 ('enzo', 'HII_Density'),
 ('enzo', 'HI_Density'),
 ('enzo', 'HI_kph'),
 ('enzo', 'HM_Density'),
 ('enzo', 'HeIII_Density'),
 ('enzo', 'HeII_Density'),
 ('enzo', 'HeII_kph'),
 ('enzo', 'HeI_Density'),
 ('enzo', 'HeI_kph'),
 ('enzo', 'Metal_Density'),
 ('enzo', 'Mg_Density'),
 ('enzo', 'Mn_Density'),
 ('enzo', 'N_Density'),
 ('enzo', 'Na_Density'),
 ('enzo', 'Ni_Density'),
 ('enzo', 'OTLW_kdissH2I'),
 ('enzo', 'O_Density'),
 ('enzo', 'Pe_heating_rate'),
 ('enzo', 'PhotoGamma'),
 ('enzo', 'RadAccel1'),
 ('enzo', 'RadAccel2'),
 ('enzo', 'RadAccel3'),
 ('enzo', 'Ray_Segments'),
 ('enzo', 'S_Density'),
 ('enzo', 'Si_Density'),
 ('enzo', 'Sr_Density'),
 ('enzo', 'Temperature'),
 ('enzo', 'TotalEnergy'),
 ('enzo', 'Y_Density'),
 ('enzo', 'cell_volume'),
 ('enzo', 'dx'),
 ('enzo', 'dy'),
 ('enzo', 'dz'),
 ('enzo', 'path_element_x'),
 ('enzo', 'path_element_y'),
 ('enzo', 'path_element_z'),
 ('enzo', 'vertex_x'),
 ('enzo', 'vertex_y'),
 ('enzo', 'vertex_z'),
 ('enzo', 'x'),
 ('enzo', 'x-velocity'),
 ('enzo', 'y'),
 ('enzo', 'y-velocity'),
 ('enzo', 'z'),
 ('enzo', 'z-velocity'),
 ('gas', 'As_Abundance'),
 ('gas', 'As_Density'),
 ('gas', 'As_Fraction'),
 ('gas', 'As_Mass'),
 ('gas', 'As_Number_Density'),
 ('gas', 'As_over_Ba'),
 ('gas', 'As_over_C'),
 ('gas', 'As_over_Fe'),
 ('gas', 'As_over_H'),
 ('gas', 'As_over_Mg'),
 ('gas', 'As_over_N'),
 ('gas', 'As_over_Na'),
 ('gas', 'As_over_Ni'),
 ('gas', 'As_over_O'),
 ('gas', 'Ba_Abundance'),
 ('gas', 'Ba_Density'),
 ('gas', 'Ba_Fraction'),
 ('gas', 'Ba_Mass'),
 ('gas', 'Ba_Number_Density'),
 ('gas', 'Ba_over_As'),
 ('gas', 'Ba_over_C'),
 ('gas', 'Ba_over_Ca'),
 ('gas', 'Ba_over_Fe'),
 ('gas', 'Ba_over_H'),
 ('gas', 'Ba_over_Mg'),
 ('gas', 'Ba_over_Mn'),
 ('gas', 'Ba_over_N'),
 ('gas', 'Ba_over_Na'),
 ('gas', 'Ba_over_Ni'),
 ('gas', 'Ba_over_O'),
 ('gas', 'Ba_over_S'),
 ('gas', 'Ba_over_Si'),
 ('gas', 'Ba_over_Sr'),
 ('gas', 'Ba_over_Y'),
 ('gas', 'C_Abundance'),
 ('gas', 'C_Density'),
 ('gas', 'C_Fraction'),
 ('gas', 'C_Mass'),
 ('gas', 'C_Number_Density'),
 ('gas', 'C_over_Ba'),
 ('gas', 'C_over_Fe'),
 ('gas', 'C_over_H'),
 ('gas', 'C_over_Mg'),
 ('gas', 'C_over_N'),
 ('gas', 'C_over_Na'),
 ('gas', 'C_over_Ni'),
 ('gas', 'C_over_O'),
 ('gas', 'Ca_Abundance'),
 ('gas', 'Ca_Density'),
 ('gas', 'Ca_Fraction'),
 ('gas', 'Ca_Mass'),
 ('gas', 'Ca_Number_Density'),
 ('gas', 'Ca_over_Ba'),
 ('gas', 'Ca_over_C'),
 ('gas', 'Ca_over_Fe'),
 ('gas', 'Ca_over_H'),
 ('gas', 'Ca_over_Mg'),
 ('gas', 'Ca_over_N'),
 ('gas', 'Ca_over_Na'),
 ('gas', 'Ca_over_Ni'),
 ('gas', 'Ca_over_O'),
 ('gas', 'El_density'),
 ('gas', 'El_fraction'),
 ('gas', 'El_mass'),
 ('gas', 'El_number_density'),
 ('gas', 'FUV_flux'),
 ('gas', 'Fe_Abundance'),
 ('gas', 'Fe_Density'),
 ('gas', 'Fe_Fraction'),
 ('gas', 'Fe_Mass'),
 ('gas', 'Fe_Number_Density'),
 ('gas', 'Fe_over_As'),
 ('gas', 'Fe_over_Ba'),
 ('gas', 'Fe_over_C'),
 ('gas', 'Fe_over_Ca'),
 ('gas', 'Fe_over_H'),
 ('gas', 'Fe_over_Mg'),
 ('gas', 'Fe_over_Mn'),
 ('gas', 'Fe_over_N'),
 ('gas', 'Fe_over_Na'),
 ('gas', 'Fe_over_Ni'),
 ('gas', 'Fe_over_O'),
 ('gas', 'Fe_over_S'),
 ('gas', 'Fe_over_Si'),
 ('gas', 'Fe_over_Sr'),
 ('gas', 'Fe_over_Y'),
 ('gas', 'G_eff'),
 ('gas', 'G_o'),
 ('gas', 'H2_density'),
 ('gas', 'H2_fraction'),
 ('gas', 'H2_mass'),
 ('gas', 'H2_number_density'),
 ('gas', 'H2_p0_density'),
 ('gas', 'H2_p0_fraction'),
 ('gas', 'H2_p0_mass'),
 ('gas', 'H2_p0_number_density'),
 ('gas', 'H2_p1_density'),
 ('gas', 'H2_p1_fraction'),
 ('gas', 'H2_p1_mass'),
 ('gas', 'H2_p1_number_density'),
 ('gas', 'H_Abundance'),
 ('gas', 'H_Mass'),
 ('gas', 'H_density'),
 ('gas', 'H_fraction'),
 ('gas', 'H_m1_density'),
 ('gas', 'H_m1_fraction'),
 ('gas', 'H_m1_mass'),
 ('gas', 'H_m1_number_density'),
 ('gas', 'H_mass'),
 ('gas', 'H_nuclei_density'),
 ('gas', 'H_number_density'),
 ('gas', 'H_p0_density'),
 ('gas', 'H_p0_fraction'),
 ('gas', 'H_p0_mass'),
 ('gas', 'H_p0_number_density'),
 ('gas', 'H_p1_density'),
 ('gas', 'H_p1_fraction'),
 ('gas', 'H_p1_mass'),
 ('gas', 'H_p1_number_density'),
 ('gas', 'H_total_mass'),
 ('gas', 'H_total_number_density'),
 ('gas', 'He_density'),
 ('gas', 'He_fraction'),
 ('gas', 'He_mass'),
 ('gas', 'He_nuclei_density'),
 ('gas', 'He_number_density'),
 ('gas', 'He_p0_density'),
 ('gas', 'He_p0_fraction'),
 ('gas', 'He_p0_mass'),
 ('gas', 'He_p0_number_density'),
 ('gas', 'He_p1_density'),
 ('gas', 'He_p1_fraction'),
 ('gas', 'He_p1_mass'),
 ('gas', 'He_p1_number_density'),
 ('gas', 'He_p2_density'),
 ('gas', 'He_p2_fraction'),
 ('gas', 'He_p2_mass'),
 ('gas', 'He_p2_number_density'),
 ('gas', 'He_total_mass'),
 ('gas', 'LW_flux'),
 ('gas', 'Metal_Density'),
 ('gas', 'Metal_Number_Density'),
 ('gas', 'Mg_Abundance'),
 ('gas', 'Mg_Density'),
 ('gas', 'Mg_Fraction'),
 ('gas', 'Mg_Mass'),
 ('gas', 'Mg_Number_Density'),
 ('gas', 'Mg_over_Ba'),
 ('gas', 'Mg_over_C'),
 ('gas', 'Mg_over_Fe'),
 ('gas', 'Mg_over_H'),
 ('gas', 'Mg_over_N'),
 ('gas', 'Mg_over_Na'),
 ('gas', 'Mg_over_Ni'),
 ('gas', 'Mg_over_O'),
 ('gas', 'Mn_Abundance'),
 ('gas', 'Mn_Density'),
 ('gas', 'Mn_Fraction'),
 ('gas', 'Mn_Mass'),
 ('gas', 'Mn_Number_Density'),
 ('gas', 'Mn_over_Ba'),
 ('gas', 'Mn_over_C'),
 ('gas', 'Mn_over_Fe'),
 ('gas', 'Mn_over_H'),
 ('gas', 'Mn_over_Mg'),
 ('gas', 'Mn_over_N'),
 ('gas', 'Mn_over_Na'),
 ('gas', 'Mn_over_Ni'),
 ('gas', 'Mn_over_O'),
 ('gas', 'N_Abundance'),
 ('gas', 'N_Density'),
 ('gas', 'N_Fraction'),
 ('gas', 'N_Mass'),
 ('gas', 'N_Number_Density'),
 ('gas', 'N_over_Ba'),
 ('gas', 'N_over_C'),
 ('gas', 'N_over_Fe'),
 ('gas', 'N_over_H'),
 ('gas', 'N_over_Mg'),
 ('gas', 'N_over_Na'),
 ('gas', 'N_over_Ni'),
 ('gas', 'N_over_O'),
 ('gas', 'Na_Abundance'),
 ('gas', 'Na_Density'),
 ('gas', 'Na_Fraction'),
 ('gas', 'Na_Mass'),
 ('gas', 'Na_Number_Density'),
 ('gas', 'Na_over_Ba'),
 ('gas', 'Na_over_C'),
 ('gas', 'Na_over_Fe'),
 ('gas', 'Na_over_H'),
 ('gas', 'Na_over_Mg'),
 ('gas', 'Na_over_N'),
 ('gas', 'Na_over_Ni'),
 ('gas', 'Na_over_O'),
 ('gas', 'Ni_Abundance'),
 ('gas', 'Ni_Density'),
 ('gas', 'Ni_Fraction'),
 ('gas', 'Ni_Mass'),
 ('gas', 'Ni_Number_Density'),
 ('gas', 'Ni_over_Ba'),
 ('gas', 'Ni_over_C'),
 ('gas', 'Ni_over_Fe'),
 ('gas', 'Ni_over_H'),
 ('gas', 'Ni_over_Mg'),
 ('gas', 'Ni_over_N'),
 ('gas', 'Ni_over_Na'),
 ('gas', 'Ni_over_O'),
 ('gas', 'OTLW_kdissH2I'),
 ('gas', 'O_Abundance'),
 ('gas', 'O_Density'),
 ('gas', 'O_Fraction'),
 ('gas', 'O_Mass'),
 ('gas', 'O_Number_Density'),
 ('gas', 'O_over_Ba'),
 ('gas', 'O_over_C'),
 ('gas', 'O_over_Fe'),
 ('gas', 'O_over_H'),
 ('gas', 'O_over_Mg'),
 ('gas', 'O_over_N'),
 ('gas', 'O_over_Na'),
 ('gas', 'O_over_Ni'),
 ('gas', 'Pe_heating_rate'),
 ('gas', 'Pe_heating_rate_masked'),
 ('gas', 'Q0_flux'),
 ('gas', 'Q1_flux'),
 ('gas', 'S_Abundance'),
 ('gas', 'S_Density'),
 ('gas', 'S_Fraction'),
 ('gas', 'S_Mass'),
 ('gas', 'S_Number_Density'),
 ('gas', 'S_over_Ba'),
 ('gas', 'S_over_C'),
 ('gas', 'S_over_Fe'),
 ('gas', 'S_over_H'),
 ('gas', 'S_over_Mg'),
 ('gas', 'S_over_N'),
 ('gas', 'S_over_Na'),
 ('gas', 'S_over_Ni'),
 ('gas', 'S_over_O'),
 ('gas', 'Si_Abundance'),
 ('gas', 'Si_Density'),
 ('gas', 'Si_Fraction'),
 ('gas', 'Si_Mass'),
 ('gas', 'Si_Number_Density'),
 ('gas', 'Si_over_Ba'),
 ('gas', 'Si_over_C'),
 ('gas', 'Si_over_Fe'),
 ('gas', 'Si_over_H'),
 ('gas', 'Si_over_Mg'),
 ('gas', 'Si_over_N'),
 ('gas', 'Si_over_Na'),
 ('gas', 'Si_over_Ni'),
 ('gas', 'Si_over_O'),
 ('gas', 'Sr_Abundance'),
 ('gas', 'Sr_Density'),
 ('gas', 'Sr_Fraction'),
 ('gas', 'Sr_Mass'),
 ('gas', 'Sr_Number_Density'),
 ('gas', 'Sr_over_Ba'),
 ('gas', 'Sr_over_C'),
 ('gas', 'Sr_over_Fe'),
 ('gas', 'Sr_over_H'),
 ('gas', 'Sr_over_Mg'),
 ('gas', 'Sr_over_N'),
 ('gas', 'Sr_over_Na'),
 ('gas', 'Sr_over_Ni'),
 ('gas', 'Sr_over_O'),
 ('gas', 'Y_Abundance'),
 ('gas', 'Y_Density'),
 ('gas', 'Y_Fraction'),
 ('gas', 'Y_Mass'),
 ('gas', 'Y_Number_Density'),
 ('gas', 'Y_over_Ba'),
 ('gas', 'Y_over_C'),
 ('gas', 'Y_over_Fe'),
 ('gas', 'Y_over_H'),
 ('gas', 'Y_over_Mg'),
 ('gas', 'Y_over_N'),
 ('gas', 'Y_over_Na'),
 ('gas', 'Y_over_Ni'),
 ('gas', 'Y_over_O'),
 ('gas', 'a_grav'),
 ('gas', 'a_grav_x'),
 ('gas', 'a_grav_y'),
 ('gas', 'a_grav_z'),
 ('gas', 'a_rad'),
 ('gas', 'a_rad_over_a_grav'),
 ('gas', 'alpha_Abundance'),
 ('gas', 'alpha_Mass'),
 ('gas', 'alpha_over_As'),
 ('gas', 'alpha_over_Ba'),
 ('gas', 'alpha_over_C'),
 ('gas', 'alpha_over_Ca'),
 ('gas', 'alpha_over_Fe'),
 ('gas', 'alpha_over_H'),
 ('gas', 'alpha_over_Mg'),
 ('gas', 'alpha_over_Mn'),
 ('gas', 'alpha_over_N'),
 ('gas', 'alpha_over_Na'),
 ('gas', 'alpha_over_Ni'),
 ('gas', 'alpha_over_O'),
 ('gas', 'alpha_over_S'),
 ('gas', 'alpha_over_Si'),
 ('gas', 'alpha_over_Sr'),
 ('gas', 'alpha_over_Y'),
 ('gas', 'angular_momentum_magnitude'),
 ('gas', 'angular_momentum_x'),
 ('gas', 'angular_momentum_y'),
 ('gas', 'angular_momentum_z'),
 ('gas', 'averaged_density'),
 ('gas', 'baroclinic_vorticity_magnitude'),
 ('gas', 'baroclinic_vorticity_x'),
 ('gas', 'baroclinic_vorticity_y'),
 ('gas', 'baroclinic_vorticity_z'),
 ('gas', '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', 'gas_gravitational_potential'),
 ('gas', 'gas_gravitational_potential_gradient_magnitude'),
 ('gas', 'gas_gravitational_potential_gradient_x'),
 ('gas', 'gas_gravitational_potential_gradient_y'),
 ('gas', 'gas_gravitational_potential_gradient_z'),
 ('gas', 'gravitationally_bound'),
 ('gas', 'is_star_forming'),
 ('gas', 'jeans_mass'),
 ('gas', 'kT'),
 ('gas', 'kinetic_energy'),
 ('gas', 'mach_number'),
 ('gas', 'matter_density'),
 ('gas', 'matter_mass'),
 ('gas', 'mazzotta_weighting'),
 ('gas', 'mean_molecular_weight'),
 ('gas', 'metal_density'),
 ('gas', 'metal_mass'),
 ('gas', 'metallicity'),
 ('gas', 'number_density'),
 ('gas', 'path_element_x'),
 ('gas', 'path_element_y'),
 ('gas', 'path_element_z'),
 ('gas', 'photo_gamma'),
 ('gas', 'pos_gravitational_potential'),
 ('gas', 'pos_total_gravitational_potential'),
 ('gas', 'potential_energy'),
 ('gas', 'pressure'),
 ('gas', 'pressure_gradient_magnitude'),
 ('gas', 'pressure_gradient_x'),
 ('gas', 'pressure_gradient_y'),
 ('gas', 'pressure_gradient_z'),
 ('gas', 'radial_mach_number'),
 ('gas', 'radial_velocity'),
 ('gas', 'radial_velocity_absolute'),
 ('gas', 'radiation_acceleration_x'),
 ('gas', 'radiation_acceleration_y'),
 ('gas', 'radiation_acceleration_z'),
 ('gas', 'shear'),
 ('gas', 'shear_criterion'),
 ('gas', 'shear_mach'),
 ('gas', 'sound_speed'),
 ('gas', 'specific_angular_momentum_magnitude'),
 ('gas', 'specific_angular_momentum_x'),
 ('gas', 'specific_angular_momentum_y'),
 ('gas', 'specific_angular_momentum_z'),
 ('gas', 'sz_kinetic'),
 ('gas', 'szy'),
 ('gas', 'tangential_over_velocity_magnitude'),
 ('gas', 'tangential_velocity'),
 ('gas', 'temperature'),
 ('gas', 'thermal_energy'),
 ('gas', 'total_energy'),
 ('gas', 'total_gravitational_potential'),
 ('gas', 'velocity_cylindrical_radius'),
 ('gas', 'velocity_cylindrical_theta'),
 ('gas', 'velocity_cylindrical_z'),
 ('gas', 'velocity_divergence'),
 ('gas', 'velocity_divergence_absolute'),
 ('gas', 'velocity_magnitude'),
 ('gas', 'velocity_spherical_phi'),
 ('gas', 'velocity_spherical_radius'),
 ('gas', 'velocity_spherical_theta'),
 ('gas', 'velocity_x'),
 ('gas', 'velocity_y'),
 ('gas', 'velocity_z'),
 ('gas', 'vertex_x'),
 ('gas', 'vertex_y'),
 ('gas', 'vertex_z'),
 ('gas', 'vorticity_growth_magnitude'),
 ('gas', 'vorticity_growth_magnitude_absolute'),
 ('gas', 'vorticity_growth_timescale'),
 ('gas', 'vorticity_growth_x'),
 ('gas', 'vorticity_growth_y'),
 ('gas', 'vorticity_growth_z'),
 ('gas', 'vorticity_magnitude'),
 ('gas', 'vorticity_radiation_pressure_growth_magnitude'),
 ('gas', 'vorticity_radiation_pressure_growth_magnitude_absolute'),
 ('gas', 'vorticity_radiation_pressure_growth_timescale'),
 ('gas', 'vorticity_radiation_pressure_growth_x'),
 ('gas', 'vorticity_radiation_pressure_growth_y'),
 ('gas', 'vorticity_radiation_pressure_growth_z'),
 ('gas', 'vorticity_radiation_pressure_magnitude'),
 ('gas', 'vorticity_radiation_pressure_x'),
 ('gas', 'vorticity_radiation_pressure_y'),
 ('gas', 'vorticity_radiation_pressure_z'),
 ('gas', 'vorticity_squared'),
 ('gas', 'vorticity_stretching_magnitude'),
 ('gas', 'vorticity_stretching_x'),
 ('gas', 'vorticity_stretching_y'),
 ('gas', 'vorticity_stretching_z'),
 ('gas', 'vorticity_x'),
 ('gas', 'vorticity_y'),
 ('gas', 'vorticity_z'),
 ('gas', 'x'),
 ('gas', 'xray_emissivity'),
 ('gas', 'y'),
 ('gas', 'z'),
 ('index', 'DM_background_density'),
 ('index', 'DM_background_potential'),
 ('index', 'cell_volume'),
 ('index', 'cylindrical_r'),
 ('index', 'cylindrical_radius'),
 ('index', 'cylindrical_theta'),
 ('index', 'cylindrical_z'),
 ('index', 'disk_angle'),
 ('index', 'dx'),
 ('index', 'dy'),
 ('index', 'dz'),
 ('index', 'grid_indices'),
 ('index', 'grid_level'),
 ('index', 'height'),
 ('index', 'magnitude_cylindrical_radius'),
 ('index', 'magnitude_cylindrical_z'),
 ('index', 'morton_index'),
 ('index', 'ones'),
 ('index', 'ones_over_dx'),
 ('index', 'path_element_x'),
 ('index', 'path_element_y'),
 ('index', 'path_element_z'),
 ('index', 'radius'),
 ('index', 'spherical_phi'),
 ('index', 'spherical_r'),
 ('index', 'spherical_radius'),
 ('index', 'spherical_theta'),
 ('index', 'vertex_x'),
 ('index', 'vertex_y'),
 ('index', 'vertex_z'),
 ('index', 'virial_radius_fraction'),
 ('index', 'x'),
 ('index', 'y'),
 ('index', 'z'),
 ('index', 'zeros'),
 ('io', 'age'),
 ('io', 'birth_mass'),
 ('io', 'creation_time'),
 ('io', 'dynamical_time'),
 ('io', 'mesh_id'),
 ('io', 'metallicity_fraction'),
 ('io', 'particle_As_abundance'),
 ('io', 'particle_As_fraction'),
 ('io', 'particle_As_over_Ba'),
 ('io', 'particle_As_over_C'),
 ('io', 'particle_As_over_Fe'),
 ('io', 'particle_As_over_H'),
 ('io', 'particle_As_over_Mg'),
 ('io', 'particle_As_over_N'),
 ('io', 'particle_As_over_Na'),
 ('io', 'particle_As_over_Ni'),
 ('io', 'particle_As_over_O'),
 ('io', 'particle_Ba_abundance'),
 ('io', 'particle_Ba_fraction'),
 ('io', 'particle_Ba_over_As'),
 ('io', 'particle_Ba_over_C'),
 ('io', 'particle_Ba_over_Ca'),
 ('io', 'particle_Ba_over_Fe'),
 ('io', 'particle_Ba_over_H'),
 ('io', 'particle_Ba_over_Mg'),
 ('io', 'particle_Ba_over_Mn'),
 ('io', 'particle_Ba_over_N'),
 ('io', 'particle_Ba_over_Na'),
 ('io', 'particle_Ba_over_Ni'),
 ('io', 'particle_Ba_over_O'),
 ('io', 'particle_Ba_over_S'),
 ('io', 'particle_Ba_over_Si'),
 ('io', 'particle_Ba_over_Sr'),
 ('io', 'particle_Ba_over_Y'),
 ('io', 'particle_C_abundance'),
 ('io', 'particle_C_fraction'),
 ('io', 'particle_C_over_Ba'),
 ('io', 'particle_C_over_Fe'),
 ('io', 'particle_C_over_H'),
 ('io', 'particle_C_over_Mg'),
 ('io', 'particle_C_over_N'),
 ('io', 'particle_C_over_Na'),
 ('io', 'particle_C_over_Ni'),
 ('io', 'particle_C_over_O'),
 ('io', 'particle_Ca_abundance'),
 ('io', 'particle_Ca_fraction'),
 ('io', 'particle_Ca_over_Ba'),
 ('io', 'particle_Ca_over_C'),
 ('io', 'particle_Ca_over_Fe'),
 ('io', 'particle_Ca_over_H'),
 ('io', 'particle_Ca_over_Mg'),
 ('io', 'particle_Ca_over_N'),
 ('io', 'particle_Ca_over_Na'),
 ('io', 'particle_Ca_over_Ni'),
 ('io', 'particle_Ca_over_O'),
 ('io', 'particle_Fe_abundance'),
 ('io', 'particle_Fe_fraction'),
 ('io', 'particle_Fe_over_As'),
 ('io', 'particle_Fe_over_Ba'),
 ('io', 'particle_Fe_over_C'),
 ('io', 'particle_Fe_over_Ca'),
 ('io', 'particle_Fe_over_H'),
 ('io', 'particle_Fe_over_Mg'),
 ('io', 'particle_Fe_over_Mn'),
 ('io', 'particle_Fe_over_N'),
 ('io', 'particle_Fe_over_Na'),
 ('io', 'particle_Fe_over_Ni'),
 ('io', 'particle_Fe_over_O'),
 ('io', 'particle_Fe_over_S'),
 ('io', 'particle_Fe_over_Si'),
 ('io', 'particle_Fe_over_Sr'),
 ('io', 'particle_Fe_over_Y'),
 ('io', 'particle_H_abundance'),
 ('io', 'particle_H_fraction'),
 ('io', 'particle_He_abundance'),
 ('io', 'particle_He_fraction'),
 ('io', 'particle_Mg_abundance'),
 ('io', 'particle_Mg_fraction'),
 ('io', 'particle_Mg_over_Ba'),
 ('io', 'particle_Mg_over_C'),
 ('io', 'particle_Mg_over_Fe'),
 ('io', 'particle_Mg_over_H'),
 ('io', 'particle_Mg_over_N'),
 ('io', 'particle_Mg_over_Na'),
 ('io', 'particle_Mg_over_Ni'),
 ('io', 'particle_Mg_over_O'),
 ('io', 'particle_Mn_abundance'),
 ('io', 'particle_Mn_fraction'),
 ('io', 'particle_Mn_over_Ba'),
 ('io', 'particle_Mn_over_C'),
 ('io', 'particle_Mn_over_Fe'),
 ('io', 'particle_Mn_over_H'),
 ('io', 'particle_Mn_over_Mg'),
 ('io', 'particle_Mn_over_N'),
 ('io', 'particle_Mn_over_Na'),
 ('io', 'particle_Mn_over_Ni'),
 ('io', 'particle_Mn_over_O'),
 ('io', 'particle_N_abundance'),
 ('io', 'particle_N_fraction'),
 ('io', 'particle_N_over_Ba'),
 ('io', 'particle_N_over_C'),
 ('io', 'particle_N_over_Fe'),
 ('io', 'particle_N_over_H'),
 ('io', 'particle_N_over_Mg'),
 ('io', 'particle_N_over_Na'),
 ('io', 'particle_N_over_Ni'),
 ('io', 'particle_N_over_O'),
 ('io', 'particle_Na_abundance'),
 ('io', 'particle_Na_fraction'),
 ('io', 'particle_Na_over_Ba'),
 ('io', 'particle_Na_over_C'),
 ('io', 'particle_Na_over_Fe'),
 ('io', 'particle_Na_over_H'),
 ('io', 'particle_Na_over_Mg'),
 ('io', 'particle_Na_over_N'),
 ('io', 'particle_Na_over_Ni'),
 ('io', 'particle_Na_over_O'),
 ('io', 'particle_Ni_abundance'),
 ('io', 'particle_Ni_fraction'),
 ('io', 'particle_Ni_over_Ba'),
 ('io', 'particle_Ni_over_C'),
 ('io', 'particle_Ni_over_Fe'),
 ('io', 'particle_Ni_over_H'),
 ('io', 'particle_Ni_over_Mg'),
 ('io', 'particle_Ni_over_N'),
 ('io', 'particle_Ni_over_Na'),
 ('io', 'particle_Ni_over_O'),
 ('io', 'particle_O_abundance'),
 ('io', 'particle_O_fraction'),
 ('io', 'particle_O_over_Ba'),
 ('io', 'particle_O_over_C'),
 ('io', 'particle_O_over_Fe'),
 ('io', 'particle_O_over_H'),
 ('io', 'particle_O_over_Mg'),
 ('io', 'particle_O_over_N'),
 ('io', 'particle_O_over_Na'),
 ('io', 'particle_O_over_Ni'),
 ('io', 'particle_S_abundance'),
 ('io', 'particle_S_fraction'),
 ('io', 'particle_S_over_Ba'),
 ('io', 'particle_S_over_C'),
 ('io', 'particle_S_over_Fe'),
 ('io', 'particle_S_over_H'),
 ('io', 'particle_S_over_Mg'),
 ('io', 'particle_S_over_N'),
 ('io', 'particle_S_over_Na'),
 ('io', 'particle_S_over_Ni'),
 ('io', 'particle_S_over_O'),
 ('io', 'particle_Si_abundance'),
 ('io', 'particle_Si_fraction'),
 ('io', 'particle_Si_over_Ba'),
 ('io', 'particle_Si_over_C'),
 ('io', 'particle_Si_over_Fe'),
 ('io', 'particle_Si_over_H'),
 ('io', 'particle_Si_over_Mg'),
 ('io', 'particle_Si_over_N'),
 ('io', 'particle_Si_over_Na'),
 ('io', 'particle_Si_over_Ni'),
 ('io', 'particle_Si_over_O'),
 ('io', 'particle_Sr_abundance'),
 ('io', 'particle_Sr_fraction'),
 ('io', 'particle_Sr_over_Ba'),
 ('io', 'particle_Sr_over_C'),
 ('io', 'particle_Sr_over_Fe'),
 ('io', 'particle_Sr_over_H'),
 ('io', 'particle_Sr_over_Mg'),
 ('io', 'particle_Sr_over_N'),
 ('io', 'particle_Sr_over_Na'),
 ('io', 'particle_Sr_over_Ni'),
 ('io', 'particle_Sr_over_O'),
 ('io', 'particle_Y_abundance'),
 ('io', 'particle_Y_fraction'),
 ('io', 'particle_Y_over_Ba'),
 ('io', 'particle_Y_over_C'),
 ('io', 'particle_Y_over_Fe'),
 ('io', 'particle_Y_over_H'),
 ('io', 'particle_Y_over_Mg'),
 ('io', 'particle_Y_over_N'),
 ('io', 'particle_Y_over_Na'),
 ('io', 'particle_Y_over_Ni'),
 ('io', 'particle_Y_over_O'),
 ('io', 'particle_age'),
 ('io', 'particle_alpha_5_abundance'),
 ('io', 'particle_alpha_5_over_As'),
 ('io', 'particle_alpha_5_over_Ba'),
 ('io', 'particle_alpha_5_over_C'),
 ('io', 'particle_alpha_5_over_Ca'),
 ('io', 'particle_alpha_5_over_Fe'),
 ('io', 'particle_alpha_5_over_H'),
 ('io', 'particle_alpha_5_over_Mg'),
 ('io', 'particle_alpha_5_over_Mn'),
 ('io', 'particle_alpha_5_over_N'),
 ('io', 'particle_alpha_5_over_Na'),
 ('io', 'particle_alpha_5_over_Ni'),
 ('io', 'particle_alpha_5_over_O'),
 ('io', 'particle_alpha_5_over_S'),
 ('io', 'particle_alpha_5_over_Si'),
 ('io', 'particle_alpha_5_over_Sr'),
 ('io', 'particle_alpha_5_over_Y'),
 ('io', 'particle_alpha_abundance'),
 ('io', 'particle_alpha_over_As'),
 ('io', 'particle_alpha_over_Ba'),
 ('io', 'particle_alpha_over_C'),
 ('io', 'particle_alpha_over_Ca'),
 ('io', 'particle_alpha_over_Fe'),
 ('io', 'particle_alpha_over_H'),
 ('io', 'particle_alpha_over_Mg'),
 ('io', 'particle_alpha_over_Mn'),
 ('io', 'particle_alpha_over_N'),
 ('io', 'particle_alpha_over_Na'),
 ('io', 'particle_alpha_over_Ni'),
 ('io', 'particle_alpha_over_O'),
 ('io', 'particle_alpha_over_S'),
 ('io', 'particle_alpha_over_Si'),
 ('io', 'particle_alpha_over_Sr'),
 ('io', 'particle_alpha_over_Y'),
 ('io', 'particle_angular_momentum'),
 ('io', 'particle_angular_momentum_magnitude'),
 ('io', 'particle_angular_momentum_x'),
 ('io', 'particle_angular_momentum_y'),
 ('io', 'particle_angular_momentum_z'),
 ('io', 'particle_cylindrical_velocity_theta'),
 ('io', 'particle_cylindrical_velocity_z'),
 ('io', 'particle_index'),
 ('io', 'particle_mass'),
 ('io', 'particle_model_E0'),
 ('io', 'particle_model_E1'),
 ('io', 'particle_model_L0'),
 ('io', 'particle_model_L1'),
 ('io', 'particle_model_L_FUV'),
 ('io', 'particle_model_L_LW'),
 ('io', 'particle_model_Q0'),
 ('io', 'particle_model_Q1'),
 ('io', 'particle_model_lifetime'),
 ('io', 'particle_model_luminosity'),
 ('io', 'particle_ones'),
 ('io', 'particle_position'),
 ('io', 'particle_position_cylindrical_radius'),
 ('io', 'particle_position_cylindrical_theta'),
 ('io', 'particle_position_cylindrical_z'),
 ('io', 'particle_position_relative'),
 ('io', 'particle_position_relative_x'),
 ('io', 'particle_position_relative_y'),
 ('io', 'particle_position_relative_z'),
 ('io', 'particle_position_spherical_phi'),
 ('io', 'particle_position_spherical_radius'),
 ('io', 'particle_position_spherical_theta'),
 ('io', 'particle_position_x'),
 ('io', 'particle_position_y'),
 ('io', 'particle_position_z'),
 ('io', 'particle_radial_velocity'),
 ('io', 'particle_radius'),
 ('io', 'particle_specific_angular_momentum'),
 ('io', 'particle_specific_angular_momentum_x'),
 ('io', 'particle_specific_angular_momentum_y'),
 ('io', 'particle_specific_angular_momentum_z'),
 ('io', 'particle_spherical_position_phi'),
 ('io', 'particle_spherical_position_radius'),
 ('io', 'particle_spherical_position_theta'),
 ('io', 'particle_spherical_velocity_phi'),
 ('io', 'particle_spherical_velocity_radius'),
 ('io', 'particle_spherical_velocity_theta'),
 ('io', 'particle_type'),
 ('io', 'particle_velocity'),
 ('io', 'particle_velocity_cylindrical_radius'),
 ('io', 'particle_velocity_cylindrical_theta'),
 ('io', 'particle_velocity_cylindrical_z'),
 ('io', 'particle_velocity_magnitude'),
 ('io', 'particle_velocity_relative'),
 ('io', 'particle_velocity_relative_x'),
 ('io', 'particle_velocity_relative_y'),
 ('io', 'particle_velocity_relative_z'),
 ('io', 'particle_velocity_spherical_phi'),
 ('io', 'particle_velocity_spherical_radius'),
 ('io', 'particle_velocity_spherical_theta'),
 ('io', 'particle_velocity_x'),
 ('io', 'particle_velocity_y'),
 ('io', 'particle_velocity_z'),
 ('io', 'rad_table_T_pos'),
 ('io', 'rad_table_Z_pos'),
 ('io', 'rad_table_g_pos'),
 ('io', 'se_table_M_pos'),
 ('io', 'se_table_Z_pos'),
 ('io', 'sn_mass_ejected'),
 ('io', 'wind_mass_ejected'),
 ('io', 'yield_table_M_pos'),
 ('io', 'yield_table_Z_pos'),
 ('main_sequence_stars', 'age'),
 ('main_sequence_stars', 'birth_mass'),
 ('main_sequence_stars', 'creation_time'),
 ('main_sequence_stars', 'dynamical_time'),
 ('main_sequence_stars', 'mesh_id'),
 ('main_sequence_stars', 'metallicity_fraction'),
 ('main_sequence_stars', 'particle_As_fraction'),
 ('main_sequence_stars', 'particle_Ba_fraction'),
 ('main_sequence_stars', 'particle_C_fraction'),
 ('main_sequence_stars', 'particle_Ca_fraction'),
 ('main_sequence_stars', 'particle_Fe_fraction'),
 ('main_sequence_stars', 'particle_H_fraction'),
 ('main_sequence_stars', 'particle_He_fraction'),
 ('main_sequence_stars', 'particle_Mg_fraction'),
 ('main_sequence_stars', 'particle_Mn_fraction'),
 ('main_sequence_stars', 'particle_N_fraction'),
 ('main_sequence_stars', 'particle_Na_fraction'),
 ('main_sequence_stars', 'particle_Ni_fraction'),
 ('main_sequence_stars', 'particle_O_fraction'),
 ('main_sequence_stars', 'particle_S_fraction'),
 ('main_sequence_stars', 'particle_Si_fraction'),
 ('main_sequence_stars', 'particle_Sr_fraction'),
 ('main_sequence_stars', 'particle_Y_fraction'),
 ('main_sequence_stars', 'particle_angular_momentum'),
 ('main_sequence_stars', 'particle_angular_momentum_magnitude'),
 ('main_sequence_stars', 'particle_angular_momentum_x'),
 ('main_sequence_stars', 'particle_angular_momentum_y'),
 ('main_sequence_stars', 'particle_angular_momentum_z'),
 ('main_sequence_stars', 'particle_cylindrical_velocity_theta'),
 ('main_sequence_stars', 'particle_cylindrical_velocity_z'),
 ('main_sequence_stars', 'particle_index'),
 ('main_sequence_stars', 'particle_mass'),
 ('main_sequence_stars', 'particle_ones'),
 ('main_sequence_stars', 'particle_position'),
 ('main_sequence_stars', 'particle_position_cylindrical_radius'),
 ('main_sequence_stars', 'particle_position_cylindrical_theta'),
 ('main_sequence_stars', 'particle_position_cylindrical_z'),
 ('main_sequence_stars', 'particle_position_relative'),
 ('main_sequence_stars', 'particle_position_relative_x'),
 ('main_sequence_stars', 'particle_position_relative_y'),
 ('main_sequence_stars', 'particle_position_relative_z'),
 ('main_sequence_stars', 'particle_position_spherical_phi'),
 ('main_sequence_stars', 'particle_position_spherical_radius'),
 ('main_sequence_stars', 'particle_position_spherical_theta'),
 ('main_sequence_stars', 'particle_position_x'),
 ('main_sequence_stars', 'particle_position_y'),
 ('main_sequence_stars', 'particle_position_z'),
 ('main_sequence_stars', 'particle_radial_velocity'),
 ('main_sequence_stars', 'particle_radius'),
 ('main_sequence_stars', 'particle_specific_angular_momentum'),
 ('main_sequence_stars', 'particle_specific_angular_momentum_x'),
 ('main_sequence_stars', 'particle_specific_angular_momentum_y'),
 ('main_sequence_stars', 'particle_specific_angular_momentum_z'),
 ('main_sequence_stars', 'particle_spherical_position_phi'),
 ('main_sequence_stars', 'particle_spherical_position_radius'),
 ('main_sequence_stars', 'particle_spherical_position_theta'),
 ('main_sequence_stars', 'particle_spherical_velocity_phi'),
 ('main_sequence_stars', 'particle_spherical_velocity_radius'),
 ('main_sequence_stars', 'particle_spherical_velocity_theta'),
 ('main_sequence_stars', 'particle_type'),
 ('main_sequence_stars', 'particle_velocity'),
 ('main_sequence_stars', 'particle_velocity_cylindrical_radius'),
 ('main_sequence_stars', 'particle_velocity_cylindrical_theta'),
 ('main_sequence_stars', 'particle_velocity_cylindrical_z'),
 ('main_sequence_stars', 'particle_velocity_magnitude'),
 ...]

















In [64]:



    


M            = data['birth_mass']
t_o          = data['creation_time'].convert_to_units('Myr')
MS_lifetime  = data[('io','particle_model_lifetime')].to('Myr')
MS_death     = t_o + MS_lifetime
px           = (data['particle_position_x'] - gal.ds.domain_center[0]).to('pc')
py           = (data['particle_position_y'] - gal.ds.domain_center[1]).to('pc')
pz           = (data['particle_position_z'] - gal.ds.domain_center[2]).to('pc')

recent_death = (MS_death > ds.current_time - dt) * (MS_death <= ds.current_time + 0.001*yt.units.Myr)
alive        = MS_death > ds.current_time + 0.001*yt.units.Myr

AGB           = M < 8.0
massive_star  = (M > 8.0) * (M < 25.0)



    











In [77]:



    


width = 500.0

boxdim = np.array([width*1.25,width*1.25,30.0])*yt.units.pc
region = gal.ds.box(ds.domain_center - boxdim*0.5, ds.domain_center + boxdim*0.5)

proj = yt.ProjectionPlot(gal.ds, 'z', ['number_density','N_over_O','O_over_H'], 
                        weight_field = 'number_density', data_source = region, width = (width,'pc'))
proj.set_unit('number_density','cm**(-3)')
proj.set_cmap('number_density','viridis')

proj.set_cmap('N_over_O','magma')
proj.set_zlim('N_over_O',-2,2)
proj.set_log('N_over_O',False)

proj.set_cmap('O_over_H','cubehelix')
proj.set_zlim('O_over_H',-5,0)
proj.set_log('O_over_H',False)

dt = 5.0 * yt.units.Myr

in_image     = (np.abs(pz) <= boxdim[2]*0.5) * (np.abs(px) <= width*0.5) * (np.abs(py) <= width*0.5)

proj.show()



    


















    


























In [78]:



    


pp = {}
pp['massive_star_winds'] = in_image * alive * massive_star
pp['AGB_winds']          = in_image * recent_death * AGB
pp['SN']                 = in_image * recent_death * massive_star
#pp['other_stars']        = in_image * alive * (np.logical_not(pp['massive_star_winds']))

for k in proj.plots.keys():
    image = proj.plots[k]
    
    #
    # Now select and annotate the points we want
    #
    for s in pp.keys():
        if np.size(px[pp[s]].value) > 0:
            print np.size(px[pp[s]]), 'Particles in ', s, px[pp[s]], py[pp[s]]
            image.axes.scatter(px[pp[s]].value,py[pp[s]].value, s = ps[s], marker = markers[s], color = colors[s])
        else:
            print 'No particles in ', s
proj.show()



    


















    






1 Particles in  AGB_winds [ 218.68010222] pc [-197.64468402] pc
No particles in  SN
4 Particles in  massive_star_winds [  8.63663388  11.28849799  10.37337861  10.24681363] pc [ 10.70303672  12.47341914  14.28900076  14.07335101] pc
1 Particles in  AGB_winds [ 218.68010222] pc [-197.64468402] pc
No particles in  SN
4 Particles in  massive_star_winds [  8.63663388  11.28849799  10.37337861  10.24681363] pc [ 10.70303672  12.47341914  14.28900076  14.07335101] pc
1 Particles in  AGB_winds [ 218.68010222] pc [-197.64468402] pc
No particles in  SN
4 Particles in  massive_star_winds [  8.63663388  11.28849799  10.37337861  10.24681363] pc [ 10.70303672  12.47341914  14.28900076  14.07335101] pc










    


























In [70]:



    


np.min(data['particle_position_x'])



    


















    
Out[70]:








0.491358217208 code_length

















In [ ]:

Content source: aemerick/galaxy_analysis

Similar notebooks:

notebook.community | gallery | about