[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