[yt-users] finding resolution of a hotspot

tazkera haque h.tazkera at gmail.com
Thu Sep 15 21:33:15 PDT 2016


Hi people,

I have a hotspot in one of my sliceplots [see attached image] and I want to
find the resolution/cell size of the hotspot, in other words, the ['x'],
['y']. ['z'] co-ordinates  and ['dx'] of the cells those contain that
maximum temperature. I have written a code like this:


*import yt*
*ds = yt.load("~/super3d_hdf5_plt_cnt_0018")*

*dd = ds.all_data()*
*v, c = ds.find_max("temperature")*

*print dd['dx']*
*print dd['dy']*
*print dd['dz']*
*print dd['x']*
*print dd['y']*
*print dd['z']*

and I have an output like this:

*yt : [INFO     ] 2016-09-15 22:49:18,265 Parameters: current_time
     = 28.501984704*
*yt : [INFO     ] 2016-09-15 22:49:18,265 Parameters: domain_dimensions
    = [128 128 128]*
*yt : [INFO     ] 2016-09-15 22:49:18,293 Parameters: domain_left_edge
     = [ -2.80000000e+10  -2.80000000e+10  -2.80000000e+10]*
*yt : [INFO     ] 2016-09-15 22:49:18,294 Parameters: domain_right_edge
    = [  2.80000000e+10   2.80000000e+10   2.80000000e+10]*
*yt : [INFO     ] 2016-09-15 22:49:18,294 Parameters:
cosmological_simulation   = 0.0*
*yt : [INFO     ] 2016-09-15 22:49:25,099 Max Value is 4.81587e+09 at
-184570312.5000000000000000 116210937.5000000000000000
6835937.5000000000000000*
*[  4.37500000e+08   4.37500000e+08   4.37500000e+08 ...,   4.37500000e+08*
*   4.37500000e+08   4.37500000e+08] code_length*
*[  4.37500000e+08   4.37500000e+08   4.37500000e+08 ...,   4.37500000e+08*
*   4.37500000e+08   4.37500000e+08] code_length*
*[  4.37500000e+08   4.37500000e+08   4.37500000e+08 ...,   4.37500000e+08*
*   4.37500000e+08   4.37500000e+08] code_length*
*[ -2.77812500e+10  -2.77812500e+10  -2.77812500e+10 ...,   2.77812500e+10*
*   2.77812500e+10   2.77812500e+10] code_length*
*[ -2.77812500e+10  -2.77812500e+10  -2.77812500e+10 ...,   2.77812500e+10*
*   2.77812500e+10   2.77812500e+10] code_length*
*[ -2.77812500e+10  -2.73437500e+10  -2.69062500e+10 ...,   2.69062500e+10*
*   2.73437500e+10   2.77812500e+10] code_length*

but I believe my code is wrong in this context, because my highest
resolution is 136 km and this is an AMR code , so I should see ['dx'] very
close to the order of 136 km. please help me get the code. thanks

Tazkera
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.spacepope.org/pipermail/yt-users-spacepope.org/attachments/20160916/8d89dad4/attachment-0001.htm>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: super3d_hdf5_plt_cnt_0018_3slice_temperature.png
Type: image/png
Size: 78479 bytes
Desc: not available
URL: <http://lists.spacepope.org/pipermail/yt-users-spacepope.org/attachments/20160916/8d89dad4/attachment-0001.png>


More information about the yt-users mailing list