[yt-users] plot with source region and field parameters

[Nathan Goldbaum]

Hey Dave,

Have you tried passing the field_parameters in as a kwargs dict?  Something like this:

field_parameters = {'avg_vx',7}
proj = oober.pf.h.proj(2,'linewidth',source=region, **field_parameters)

Although as Matt notes, ideally it should be possible to avoid the kwargs construction.

-Nathan

On May 9, 2013, at 3:20 PM, david collins <dcollins4096 at gmail.com> wrote:

> Hi, Everybody--
> 
> I have a field that I have defined that uses a field_parameter, but I don't seem to be able to get it to work with a projection.
> 
> >>> region
> AMRRegion (time0050): center=[ 0.5         0.03369141  0.68457031], left_edge=[ 0.          0.00195312  0.63378906], right_edge=[ 1.          0.06542969  0.73535156]
> 
> def linewidth(field,data):
>     avg_vx = data.get_field_parameter('avg_vx')
>     return (data['x-velocity']-avg_vx)**2
> add_field('linewidth',function=linewidth, validators=[ValidateParameter('avg_vx')])
> avg_vx = 7
> region.set_field_parameter('avg_vx',avg_vx)
> print region['linewidth'].min()
> print region.get_field_parameter('avg_vx')
> 
> which yields expected output,
> 1.15751883639e-07
> 7
> 
> But trying to make a projection of that source gives the error below.  Is it obvious what I'm doing wrong? I can plot 'Density' just fine.
> 
> Thanks!
> 
> >>> proj = oober.pf.h.proj(2,'linewidth',source=region)
> Initializing tree  0 /  4100% |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||| Time: 00:00:00 
> Traceback (most recent call last):                                                                                                                     | ETA:  --:--:-- 
>   File "<stdin>", line 1, in <module>
>   File "/nics/b/home/collins/local_2013_03_12/src/yt-hg/yt/data_objects/data_containers.py", line 1758, in __init__
>     self._refresh_data()
>   File "/nics/b/home/collins/local_2013_03_12/src/yt-hg/yt/data_objects/data_containers.py", line 318, in _refresh_data
>     self.get_data()
>   File "/nics/b/home/collins/local_2013_03_12/src/yt-hg/yt/data_objects/data_containers.py", line 1859, in get_data
>     self._add_level_to_tree(tree, level, fields)
>   File "/nics/b/home/collins/local_2013_03_12/src/yt-hg/yt/data_objects/data_containers.py", line 1981, in _add_level_to_tree
>     self._add_grid_to_tree(tree, grid, fields, zero_out, dls)
>   File "/nics/b/home/collins/local_2013_03_12/src/yt-hg/yt/data_objects/data_containers.py", line 1925, in _add_grid_to_tree
>     for fd in self._get_data_from_grid(grid, fields)]
>   File "/nics/b/home/collins/local_2013_03_12/src/yt-hg/yt/data_objects/data_containers.py", line 93, in save_state
>     tr = func(self, grid, field, *args, **kwargs)
>   File "/nics/b/home/collins/local_2013_03_12/src/yt-hg/yt/data_objects/data_containers.py", line 2003, in _get_data_from_grid
>     return [grid[field] * bad_points for field in fields]
>   File "/nics/b/home/collins/local_2013_03_12/src/yt-hg/yt/data_objects/grid_patch.py", line 157, in __getitem__
>     self.get_data(key)
>   File "/nics/b/home/collins/local_2013_03_12/src/yt-hg/yt/data_objects/grid_patch.py", line 200, in get_data
>     self._generate_field(field)
>   File "/nics/b/home/collins/local_2013_03_12/src/yt-hg/yt/data_objects/grid_patch.py", line 132, in _generate_field
>     self.pf.field_info[field].check_available(self)
>   File "/nics/b/home/collins/local_2013_03_12/src/yt-hg/yt/data_objects/field_info_container.py", line 364, in check_available
>     validator(data)
>   File "/nics/b/home/collins/local_2013_03_12/src/yt-hg/yt/data_objects/field_info_container.py", line 440, in __call__
>     raise NeedsParameter(doesnt_have)
> yt.data_objects.field_info_container.NeedsParameter: (['avg_vx'])
> 
> 
> _______________________________________________
> yt-users mailing list
> yt-users at lists.spacepope.org
> http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org