[yt-dev] Issue #1300: yt.create_profile fails to find the innermost cells (yt_analysis/yt)

Hsi-Yu Schive issues-reply at bitbucket.org
Tue Nov 29 08:18:40 PST 2016


New issue 1300: yt.create_profile fails to find the innermost cells
https://bitbucket.org/yt_analysis/yt/issues/1300/ytcreate_profile-fails-to-find-the

Hsi-Yu Schive:

The following script and dataset fail to find the innermost cells and return zeros in all profile bins. Data are spherical symmetric with a center at (7.13 Mpc, 7.13 Mpc, 7.13 Mpc). The innermost 8 cells have the same peak density about 1.23651e-25 g/cm^3.

[Dataset URL](http://use.yt/upload/ca211222)

Script:
```
#!python

ds = yt.load('./Data_000000')
sp = ds.sphere( 'c', (15, "kpc") )

print 'peak density = ', sp['density'].max()

profile = yt.create_profile( sp, 'radius', 'density', n_bins=5 )
print 'radial bins = ', profile.x.value
print 'density profile = ', profile['density'].value
```





More information about the yt-dev mailing list