Issue #933: linear bins in create_profile crash (yt_analysis/yt)
New issue 933: linear bins in create_profile crash https://bitbucket.org/yt_analysis/yt/issue/933/linear-bins-in-create_profile... Stephanie Tonnesen: I am trying to make a profile of column density vs radius of a galaxy disk using create_profile. I would like linear bins, but cannot make it work. After some back-and-forth, Nathan Goldbaum thinks that it is a bug. Here is the code, using the RamPressureStripping available data. import yt units_override = {"length_unit":(26,"kpc"),"time_unit":(3.086e14,"s"),"mass_unit":(5.165e39,"g")} ds = yt.load("Downloads/RamPressureStripping/id0/rps.0062.vtk", units_override=units_override) alld = ds.all_data() gal = alld.cut_region("obj['density'] > 1.01e-25") logs = {('index', 'cylindrical_r'): False} profile = yt.create_profile(gal, [('index', 'cylindrical_r')], # the bin field [('gas', 'cell_mass')], # profile field weight_field=None, n_bins=17,logs=logs) coldens = profile["gas","cell_mass"]/1.6733e-24 coldens[0] = coldens[0]/(3.14159*profile.x[0]*profile.x[0]) i = 1 while i < len(coldens): coldens[i] = coldens[i]/((3.14159*profile.x[i]*profile.x[i]) - (3.14159*profile.x[i-1]*profile.x[i-1])) i = i+1 mp.plot(profile.x/3.086e21,coldens) mp.show() print profile.x
participants (1)
-
Stephanie Tonnesen