Guest User

Untitled

a guest
May 28th, 2013
83
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

Adblocker detected! Please consider disabling it...

We've detected AdBlock Plus or some other adblocking software preventing Pastebin.com from fully loading.

We don't have any obnoxious sound, or popup ads, we actively block these annoying types of ads!

Please add Pastebin.com to your ad blocker whitelist or disable your adblocking software.

×