Guest User

Untitled

a guest
May 28th, 2013
73
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
  1. In [1]: from yt.mods import *
  2.  
  3. In [2]: pf = load('DD0047/output_0047')
  4. yt : [INFO     ] 2013-05-28 19:38:49,553 Parameters: current_time              = 20.046508365
  5. yt : [INFO     ] 2013-05-28 19:38:49,553 Parameters: domain_dimensions         = [256 256 256]
  6. yt : [INFO     ] 2013-05-28 19:38:49,554 Parameters: domain_left_edge          = [ 0.  0.  0.]
  7. yt : [INFO     ] 2013-05-28 19:38:49,554 Parameters: domain_right_edge         = [ 1.  1.  1.]
  8. yt : [INFO     ] 2013-05-28 19:38:49,554 Parameters: cosmological_simulation   = 1
  9. yt : [INFO     ] 2013-05-28 19:38:49,554 Parameters: current_redshift          = 16.8746585087
  10. yt : [INFO     ] 2013-05-28 19:38:49,554 Parameters: omega_lambda              = 0.724
  11. yt : [INFO     ] 2013-05-28 19:38:49,554 Parameters: omega_matter              = 0.276
  12. yt : [INFO     ] 2013-05-28 19:38:49,554 Parameters: hubble_constant           = 0.703
  13.  
  14. In [3]: pf.h.print_stats()
  15. yt : [INFO     ] 2013-05-28 19:38:52,657 Getting the binary hierarchy
  16. yt : [INFO     ] 2013-05-28 19:38:52,741 Finished with binary hierarchy reading
  17. yt : [INFO     ] 2013-05-28 19:38:52,828 Adding unknown field PotentialField to list of fields
  18. yt : [INFO     ] 2013-05-28 19:38:52,828 Adding unknown field GravitatingMassField to list of fields
  19. level   # grids        # cells       # cells^3
  20. ----------------------------------------------
  21.   0    768        16777216             255
  22.   1    396         4096000             159
  23.   2   1104        22906304             283
  24.   3   1012       147197952             527
  25.   4      7           34976              32
  26.   5      3           13264              23
  27.   6      1            5832              17
  28.   7      1            4608              16
  29.   8      1            3584              15
  30.   9      1            3584              15
  31.  10      1            3584              15
  32.  11      1            3136              14
  33.  12      1            2744              13
  34.  13      1            2744              13
  35.  14      1            2744              13
  36.  15      1            2744              13
  37.  16      1            2744              13
  38.  17      1            2744              13
  39.  18      1            2744              13
  40.  19      1            2744              13
  41.  20      1            2744              13
  42. ----------------------------------------------
  43.       3305       191074736
  44.  
  45.  
  46. t = 2.00465084e+01 = 7.37163789e+15 s = 2.33593109e+08 years
  47.  
  48. Smallest Cell:
  49.     Width: 4.168e-10 mpch
  50.     Width: 5.929e-10 mpc
  51.     Width: 3.725e-09 1
  52.     Width: 3.725e-09 unitary
  53.     Width: 7.451e-09 mpchcm
  54.     Width: 1.060e-08 mpccm
  55.     Width: 3.147e-08 aye
  56.     Width: 4.168e-07 kpch
  57.     Width: 5.929e-07 kpc
  58.     Width: 7.451e-06 kpchcm
  59.     Width: 1.060e-05 kpccm
  60.     Width: 4.168e-04 pch
  61.     Width: 5.929e-04 pc
  62.     Width: 7.451e-03 pchcm
  63.     Width: 1.060e-02 pccm
  64.     Width: 8.598e+01 auh
  65.     Width: 1.223e+02 au
  66.     Width: 1.537e+03 auhcm
  67.     Width: 2.186e+03 aucm
  68.     Width: 1.849e+04 rsunh
  69.     Width: 2.631e+04 rsun
  70.     Width: 3.306e+05 rsunhcm
  71.     Width: 4.702e+05 rsuncm
  72.     Width: 7.992e+09 milesh
  73.     Width: 1.137e+10 miles
  74.     Width: 1.429e+11 mileshcm
  75.     Width: 2.032e+11 milescm
  76.     Width: 1.286e+15 cmh
  77.     Width: 1.830e+15 cm
  78.     Width: 2.299e+16 cmhcm
  79.     Width: 3.270e+16 cmcm
  80.  
  81. In [4]: val, loc = pf.h.find_max('Density')
  82. ---------------------------------------------------------------------------
  83. IndexError                                Traceback (most recent call last)
  84. <ipython-input-4-efca65fda48d> in <module>()
  85. ----> 1 val, loc = pf.h.find_max('Density')
  86.  
  87. /usr/projects/magnetic/jsmidt/yt-x86_64/lib/python2.7/site-packages/yt-2.4-py2.7-linux-x86_64.egg/yt/data_objects/object_finding_mixin.pyc in find_max(self, field, finest_levels)
  88.      59         if (field, finest_levels) in self._max_locations:
  89.      60             return self._max_locations[(field, finest_levels)]
  90. ---> 61         mg, mc, mv, pos = self.find_max_cell_location(field, finest_levels)
  91.      62         self._max_locations[(field, finest_levels)] = (mv, pos)
  92.      63         return mv, pos
  93.  
  94. /usr/projects/magnetic/jsmidt/yt-x86_64/lib/python2.7/site-packages/yt-2.4-py2.7-linux-x86_64.egg/yt/data_objects/object_finding_mixin.pyc in find_max_cell_location(self, field, finest_levels)
  95.      72                     len(source._grids), field)
  96.      73         max_val, maxi, mx, my, mz, mg = \
  97. ---> 74             source.quantities["MaxLocation"]( field, lazy_reader=True)
  98.      75         max_grid = self.grids[mg]
  99.      76         mc = np.unravel_index(maxi, max_grid.ActiveDimensions)
  100.  
  101. /usr/projects/magnetic/jsmidt/yt-x86_64/lib/python2.7/site-packages/yt-2.4-py2.7-linux-x86_64.egg/yt/data_objects/derived_quantities.pyc in __call__(self, *args, **kwargs)
  102.      90                           self._data_source.pf.h.io)
  103.      91         if lazy_reader and not self.force_unlazy:
  104. ---> 92             return self._call_func_lazy(args, kwargs)
  105.      93         else:
  106.      94             return self._call_func_unlazy(args, kwargs)
  107.  
  108. /usr/projects/magnetic/jsmidt/yt-x86_64/lib/python2.7/site-packages/yt-2.4-py2.7-linux-x86_64.egg/yt/data_objects/derived_quantities.pyc in _call_func_lazy(self, args, kwargs)
  109.      97         self.retvals = [ [] for i in range(self.n_ret)]
  110.      98         for gi,g in enumerate(self._get_grids()):
  111. ---> 99             rv = self.func(GridChildMaskWrapper(g, self._data_source), *args, **kwargs)
  112.     100             if not iterable(rv): rv = (rv,)
  113.     101             for i in range(self.n_ret): self.retvals[i].append(rv[i])
  114.  
  115. /usr/projects/magnetic/jsmidt/yt-x86_64/lib/python2.7/site-packages/yt-2.4-py2.7-linux-x86_64.egg/yt/data_objects/derived_quantities.pyc in _MaxLocation(data, field)
  116.     670         maxi = np.argmax(data[field])
  117.     671         ma = data[field][maxi]
  118. --> 672         mx, my, mz = [data[ax][maxi] for ax in 'xyz']
  119.     673         mg = data["GridIndices"][maxi]
  120.     674     return (ma, maxi, mx, my, mz, mg)
  121.  
  122. IndexError: index out of bounds
RAW Paste Data