On Thu, Sep 15, 2016 at 11:33 PM, tazkera haque <h.tazkera@gmail.com> wrote:
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")

From here, try:

point = ds.point(c)

print(point['x'])
print(point['y'])
print(point['z'])
print(point['dx'])
print(point['dy'])
print(point['dz'])



 

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 

_______________________________________________
yt-users mailing list
yt-users@lists.spacepope.org
http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org