Hi Elizabeth,<div><br></div><div>So the problem is that when you add_field and ask to validate the height_vector, it doesn't exist yet. I believe you can get the desired behavior with:</div><div><br></div><div><div>from yt.mods import *</div>
<div><br></div><div>def _RayTest(field, data):</div><div><br></div><div> center = data.get_field_parameter("center")</div><div> h_vec = na.array(data.get_field_parameter("height_vector"))</div><div>
<br></div><div> print "center :", center, "h_vec: ", h_vec</div><div><br></div><div> return data["Density"]</div><div><br></div><div>add_field("RayTest", function=_RayTest)</div>
<div><br></div><div># Grab data</div><div>fn = "GravPotential/DD0301/GT_BTAccel_256AMR4_PeHeat_sf5_SNe_0301"<br></div><div>pf = load(fn)</div><div>ray1 = pf.h.ray((2.0, 2.0 ,2.0), (12.0, 12.0, 12.0)) </div><div>
ray1.set_field_parameter('height_vector',na.array([0.,0.,1.]))</div><div>ray1.set_field_parameter('center',[0.5]*3)</div><div>ray1["RayTest"]</div><div><br></div><div>When I run this I get: </div>
<div><div>center : [0.5, 0.5, 0.5] h_vec: [ 0. 0. 1.]</div></div><div><br></div><div>I'm not sure what dd was in your dd.center since dd wasn't defined in that script, but you should be able to set center to whatever you want here. I also changed density to Density. </div>
<div><br></div><div>Best,</div><div>Sam</div><div><br></div><div class="gmail_quote">On Wed, Aug 10, 2011 at 9:36 AM, Elizabeth Tasker <span dir="ltr"><<a href="mailto:taskere@mcmaster.ca">taskere@mcmaster.ca</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;">Hi Matt,<br>
<br>
I think this script explains what I mean. When I run it, I get:<br>
<br>
/1/home/taskere/yt/yt/data_<u></u>objects/field_info_container.<u></u>py in __call__(self, data)<br>
357 doesnt_have.append(p)<br>
358 if len(doesnt_have) > 0:<br>
--> 359 raise NeedsParameter(doesnt_have)<br>
360 return True<br>
361<br>
<br>
NeedsParameter: (['height_vector'])<br>
<br>
<br>
As a side note, I think 'ray' sets its own value for the field parameter, "center", making it the start of the ray. If that's true, does it have problems for calling radius? e.g.<br>
<br>
ray["Radius"]<br>
<br>
where radius would need to be the center of the object, not necessarily the start of the ray? (This particular issue wouldn't affect me, but I thought it might bite someone later).<br>
<br>
Elizabeth<br>
<br>
<br>
<br>
Matthew Turk wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div class="im">
<br>
Hi Elizabeth,<br>
<br>
Could you elaborate on what you mean, and provide an example script?<br>
<br>
Thanks,<br>
<br>
-Matt<br>
<br></div><div class="im">
On Aug 10, 2011 7:22 AM, "Elizabeth Tasker" <<a href="mailto:taskere@mcmaster.ca" target="_blank">taskere@mcmaster.ca</a> <mailto:<a href="mailto:taskere@mcmaster.ca" target="_blank">taskere@mcmaster.ca</a>>> wrote:<br>
> Hi,<br>
><br>
> I'm having problems passing parameters to a ray object. I'm trying:<br>
><br>
> ray1 = pf.h.ray( com[c], com[nearestcore], field_parameters =<br>
> {'height_vector':na.array([0,<u></u>0,1]), 'center':dd.center})<br>
> Omega0 = ray1["Omega"]<br>
><br>
> and I also tried out:<br>
><br>
> ray1.set_field_parameter("<u></u>height_vector", na.array([0,0,1]))<br>
> ray1.set_field_parameter("<u></u>center", dd.center)<br>
><br>
> But it can't recognise them.<br>
><br>
> Elizabeth<br>
> ______________________________<u></u>_________________<br>
> yt-users mailing list<br></div>
> <a href="mailto:yt-users@lists.spacepope.org" target="_blank">yt-users@lists.spacepope.org</a> <mailto:<a href="mailto:yt-users@lists.spacepope.org" target="_blank">yt-users@lists.<u></u>spacepope.org</a>><div class="im">
<br>
> <a href="http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org" target="_blank">http://lists.spacepope.org/<u></u>listinfo.cgi/yt-users-<u></u>spacepope.org</a><br></div>
------------------------------<u></u>------------------------------<u></u>------------<div class="im"><br>
<br>
______________________________<u></u>_________________<br>
yt-users mailing list<br>
<a href="mailto:yt-users@lists.spacepope.org" target="_blank">yt-users@lists.spacepope.org</a><br>
<a href="http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org" target="_blank">http://lists.spacepope.org/<u></u>listinfo.cgi/yt-users-<u></u>spacepope.org</a><br>
<br>
</div></blockquote>
<br>
<br>_______________________________________________<br>
yt-users mailing list<br>
<a href="mailto:yt-users@lists.spacepope.org">yt-users@lists.spacepope.org</a><br>
<a href="http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org" target="_blank">http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org</a><br>
<br></blockquote></div><br></div>