Britton, Devin, and I are writing up the Trident method paper, and we
wanted to be explicit about how exactly yt deposits SPH field quantities on
to an octree for particle-based fields. Are there any formal equations for
this written down anywhere? I realize SPH support was added since the
first yt method paper, so it isn't in there, and I cannot seem to find
anything in the docs. I wonder if anyone has anything written down for
this or if I'm missing it somewhere?
NSF Postdoctoral Fellow
Department of Astronomy
California Institute of Technology
New issue 1275: Doc notebooks are not properly evaluated
There's no difference between ["raw"](http://yt-project.org/docs/dev/analyzing/analysis_modules/Particle_Trajectories.ipynb) and ["evaluated"](http://yt-project.org/docs/dev/analyzing/analysis_modules/Particle_Trajectories_evaluated.ipynb)
New issue 1274: RAMSES temperature units wrong
I tried newest pip yt today after a long break and plotted RAMSES isolated galaxy. The temperatures I got were ~1e10 K, which is incorrect. I think in the conversion of temperature you forgot to check the boxsize (in my case it was 60. in code units), so I'd see the problem in this.
New issue 1271: Derived fields that return data with incorrect shape should raise errors
Here's a test script that illustrates the problem:
import numpy as np
def wrong_field(field, data):
return np.array([1, 2, 3])
ds = yt.load('IsolatedGalaxy/galaxy0030/galaxy0030')
ad = ds.all_data()
On my machine, this prints:
YTArray([1, 2, 3]) (dimensionless)
The fact that this works without erroring out right now leads to broknness later and user confusion. See Issue #1270.
New issue 1270: unable to plot created field
I have created some fields. It is able to plot some but others (particularly my Enclosed mass field and every other field that uses that field) cannot be plotted. When I try to use Sliceplot the warning "Plot image for field ('gas', 'Enclosed_Mass') has zero dynamic range. Min = Max = 0.000000.". However, when the field is retrieved using dd['Enclosed_Mass'], it is not all zero, in fact min=6.01427540665e+32 g and max=2.11538515805e+33 g.
This field is not created in the conventional way, with one to one matrix math, and it is a bit slow to calculate because I used radial bins to approximate the value. But I think the main issue is that, the field is non-zero everywhere, but Sliceplot is unable to plot the field or other fiends that use it.