In [1]:

    

%matplotlib inline

import numpy as np
import matplotlib.pyplot as plt
import galaxy_analysis as ga


    

    




yt : [INFO     ] 2017-02-21 16:25:06,046 Loading plugins from /home/emerick/.yt/my_plugins.py
/home/emerick/anaconda/lib/python2.7/site-packages/matplotlib/__init__.py:1350: 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.

  warnings.warn(_use_error_msg)

In [7]:

    

home = "/home/emerick/work/enzo_runs/leo_p_stampede/fullres/run13/3parsec/"

galaxy = ga.Galaxy("DD0040", wdir = home)


    

    




yt : [INFO     ] 2017-02-21 16:33:18,461 Parameters: current_time              = 156.999907313
yt : [INFO     ] 2017-02-21 16:33:18,462 Parameters: domain_dimensions         = [128 128 128]
yt : [INFO     ] 2017-02-21 16:33:18,463 Parameters: domain_left_edge          = [ 0.  0.  0.]
yt : [INFO     ] 2017-02-21 16:33:18,463 Parameters: domain_right_edge         = [ 1.  1.  1.]
yt : [INFO     ] 2017-02-21 16:33:18,464 Parameters: cosmological_simulation   = 0.0
yt : [INFO     ] 2017-02-21 16:33:18,478 Parameters: current_time              = 156.999907313
yt : [INFO     ] 2017-02-21 16:33:18,478 Parameters: domain_dimensions         = [128 128 128]
yt : [INFO     ] 2017-02-21 16:33:18,479 Parameters: domain_left_edge          = [ 0.  0.  0.]
yt : [INFO     ] 2017-02-21 16:33:18,479 Parameters: domain_right_edge         = [ 1.  1.  1.]
yt : [INFO     ] 2017-02-21 16:33:18,480 Parameters: cosmological_simulation   = 0.0
Parsing Hierarchy : 100%|██████████| 302/302 [00:00<00:00, 45088.81it/s]
yt : [INFO     ] 2017-02-21 16:33:18,500 Gathering a field list (this may take a moment.)

In [8]:

    

galaxy.compute_everything()
galaxy.save()


    

    




yt : [INFO     ] 2017-02-21 16:34:00,253 Projection completed
yt : [INFO     ] 2017-02-21 16:34:00,386 Projection completed

In [4]:

    

galaxy.time_data['time']


    

    
Out[4]:





YTArray([  1.27787509e+08,   1.29787509e+08,   1.31787509e+08,
         1.33787509e+08,   1.35787509e+08,   1.37787509e+08,
         1.39787509e+08,   1.41787509e+08,   1.43787509e+08,
         1.45787509e+08,   1.47787509e+08,   1.49787509e+08,
         1.51787509e+08,   1.53787509e+08,   1.55787509e+08]) yr

In [5]:

    

galaxy.time_data['SFR']


    

    
Out[5]:





array([  0.00000000e+00,   5.05809196e-05,   4.29914982e-04,
         3.73939524e-04,   3.25628158e-04,   1.24878338e-03,
         2.07376278e-03,   2.81310595e-03,   5.93158125e-03,
         4.16779286e-03,   3.28089148e-03,   3.50902450e-03,
         2.70662064e-03,   2.44202003e-03,   3.64011393e-03])

In [ ]:

    




    

In [ ]:

    




    

In [ ]:

    




    

In [ ]: