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

[David Collins]

If you do the projection of the region, it correctly gets the region right?
>
>
pf.h.proj(2,'Density',source=region)

does the correct thing for that region, if that's what you mean.

d.


> If so, file a bug, assign it to me, and I will correct it when I can
> get back to the computer.  I think I know what's going on -- I suspect
> during the _get_data_from_grid call (or whatever it is in the
> quadtree) it's not passing field parameters through.  As a backup for
> the moment you might be able to use overlap_proj, but this is 100% a
> bug and needs to be fixed, so file a bug assigned to me and I'll get
> at it asap.
>
> -Matt
>
> On Thu, May 9, 2013 at 6: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
> >
> _______________________________________________
> yt-users mailing list
> yt-users at lists.spacepope.org
> http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.spacepope.org/pipermail/yt-users-spacepope.org/attachments/20130509/3c85a607/attachment.html>