--- you can reply above this line ---
New issue 531: "RadiusCode" field is ill-defined for periodic domains https://bitbucket.org/yt_analysis/yt/issue/531/radiuscode-field-is-ill-defin...
Kacper Kowalik:
As per ml [thread](http://lists.spacepope.org/pipermail/yt-dev-spacepope.org/2013-March/002979....) following script:
``` #!python
from yt.mods import * from yt.testing import * from yt.data_objects.universal_fields import get_radius from yt.data_objects.field_info_container import ValidateParameter
@derived_field(name="invrad", validators = [ValidateParameter("center")]) def _invrad(field, data): data.set_field_parameter("center", [1.0, 1.0, 1.0]) return 1.0/data["RadiusCode"]
pf = fake_random_pf(64) sp = pf.h.sphere([1.0, 1.0, 1.0], 0.25) r = get_radius(sp, "")
if hasattr(sp, "size"): v = 3 p1 = pf.h.proj("invrad", 0, data_source = sp, style="mip") p2 = pf.h.proj("Density", 0, data_source = sp) else: v = 2 p1 = pf.h.proj(0, "invrad", source = sp, proj_style="mip") p2 = pf.h.proj(0, "Density", source = sp)
pw1 = p1.to_pw() pw1.save("yt_v%s" % v)
pw2 = p2.to_pw() pw2.set_zlim("Density", 0, 0.32) pw2.save("yt_v%s" % v) ``` produces   while it should yield  
--
This is an issue notification from bitbucket.org. You are receiving this either because you are the owner of the issue, or you are following the issue.
-
Kacper Kowalik