[yt-users] finding resolution of a hotspot

Nathan Goldbaum nathan12343 at gmail.com
Thu Sep 15 21:53:54 PDT 2016


On Thu, Sep 15, 2016 at 11:33 PM, tazkera haque <h.tazkera at 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 at lists.spacepope.org
> http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org
>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.spacepope.org/pipermail/yt-users-spacepope.org/attachments/20160915/d9dae7c4/attachment.htm>


More information about the yt-users mailing list