[yt-users] Over plot particles
REMUDIN RESHID MEKURIA
remudin at yahoo.com
Thu May 24 14:04:49 PDT 2012
Dear Matt, Casey,
In the list of field, I do not have any of the particle positions, I am not sure if I need to. But these is the result of your quick test:
print pf.h.field_list ['density', 'xmom', 'ymom', 'zmom', 'rho_E', 'rho_e', 'rho_K', 'Temp', 'rho_H', 'rho_He', 'phi_grav', 'grav_x', 'grav_y', 'grav_z', 'prod_sgs', 'diss_sgs', 'turb_src', 'pressure', 'MachNumber', 'magvort', 'particle_count', 'particle_mass_density']
print pf.use_particles False
Ok, I haven't seen some thing like "particles.write_in_plotfile" in inputs. I have attached the inputsI used here within.
Thank you again for taking much of your time.
Cheers,
Remudin
On 2012 Cam 24, at 22:40, Matthew Turk wrote:
> Hi Remudin,
>
> Ah, bummer. Revisiting what Casey said, it sounds like you're trying
> to plot and it doesn't know that there are particles. What is the
> output of these commands:
>
> print pf.h.field_list
> print pf.use_particles
>
> And, can you check the inputs file to see what the value of
> "particles.write_in_plotfile" is? That parameter dictates (for Nyx)
> whether particles are written out in the plotfiles, and without it yt
> cannot read them.
>
> -Matt
>
> On Fri, May 25, 2012 at 5:36 AM, REMUDIN RESHID MEKURIA
> <remudin at yahoo.com> wrote:
>> Dear Matt,
>>
>> I have included the changes in the script and the error doesn't
>> change. KeyError: 'particle_position_y'.
>>
>> Right, it will be of great simplification when those Units are clarified in
>> the recipes by you guys soon, thanks.
>>
>>
>> Cheers,
>>
>> Remudin
>>
>> On 2012 Cam 24, at 22:15, Matthew Turk wrote:
>>
>> Hi Remudin,
>>
>> The units between Enzo and Nyx are different -- and to get the recipe
>> to reflect that there are a few places that have to be modified. In
>> response to the community survey results that came in, just yesterday
>> (coincidentally enough) Britton Smith and I have been going through
>> trying to make changes to the recipes to de-enzo-ify them and to make
>> it more clear what are units and what aren't.
>>
>> I'm not sure this will fix the problem, but it might help!
>>
>> On Fri, May 25, 2012 at 5:08 AM, REMUDIN RESHID MEKURIA
>> <remudin at yahoo.com> wrote:
>>
>> Thanks Casey,
>>
>>
>> This is the script I used from yt-cookbook Overplot particles
>>
>>
>> #!/usr/users/rmekuri/data/yt-x86_64/bin/python
>>
>>
>>
>> from yt.mods import * # set up our namespace
>>
>>
>>
>>
>> fn = "plt00736" # parameter file to load
>>
>>
>>
>>
>> pf = load(fn) # load data
>>
>>
>>
>> pc = PlotCollection(pf, [0.5, 0.5, 0.5])
>>
>>
>> Try changing [0.5, 0.5, 0.5] to "c" like this:
>>
>> pc = PlotCollection(pf, "c")
>>
>> this will put you at the domain center.
>>
>> p = pc.add_projection("Density", 0)
>>
>> p.modify["particles"](1.0) # 1.0 is the 'width' we want for our slab of
>>
>>
>> try changing this 1.0 to the full domain:
>>
>> 1.0/pf['unitary']
>>
>> 'unitary' is the unit meaning the full width of the domain.
>>
>>
>>
>> # particles -- this governs the allowable
>>
>> locations
>>
>>
>> # of particles that show up on the image
>>
>>
>>
>> # NOTE: we can also supply a *ptype* to cut
>>
>> based
>>
>>
>> # on a given (integer) particle type
>>
>>
>>
>> pc.set_width(1.0, '1') # change width of our plot to the full domain
>>
>>
>> Let's set this also to the full domain:
>>
>> pc.set_width(1.0, 'unitary')
>>
>>
>> -Matt
>>
>>
>>
>> pc.save(fn) # save all plots
>>
>>
>> Cheers,
>>
>>
>> Remudin
>>
>>
>>
>> On 2012 Cam 24, at 22:03, Casey W. Stark wrote:
>>
>>
>> Hi Remudin.
>>
>>
>> Could you include the script that is generating the error?
>>
>>
>> I'm not sure without the script, but it looks like you are trying to access
>>
>> the "particle_position_y" field and it is not in the plotfile you are
>>
>> working with.
>>
>>
>> Best,
>>
>> Casey
>>
>>
>>
>> On Thu, May 24, 2012 at 11:15 AM, REMUDIN RESHID MEKURIA <remudin at yahoo.com>
>>
>> wrote:
>>
>>
>> Dear all,
>>
>>
>> Thank you for your help in advance.
>>
>> I am currently running cosmological simulations using Nyx code.
>>
>> And while trying to over plotting particles on the density slice, I
>>
>> encountered the following error which I couldn't resolve. I will greatly
>>
>> appreciate your help.
>>
>>
>> Traceback (most recent call last):
>>
>> File "./OverPlotingParticles.py", line 14, in <module>
>>
>> pc.set_width(1.0, '1') # change width of our plot to the full domain
>>
>> File
>>
>> "/home/uni06/cosmo/rmekuri/yt-x86_64/src/yt-hg/yt/visualization/plot_collection.py",
>>
>> line 267, in set_width
>>
>> plot.set_width(width, unit)
>>
>> File
>>
>> "/home/uni06/cosmo/rmekuri/yt-x86_64/src/yt-hg/yt/visualization/plot_types.py",
>>
>> line 441, in set_width
>>
>> self._refresh_display_width()
>>
>> File
>>
>> "/home/uni06/cosmo/rmekuri/yt-x86_64/src/yt-hg/yt/visualization/plot_types.py",
>>
>> line 459, in _refresh_display_width
>>
>> self._redraw_image()
>>
>> File
>>
>> "/home/uni06/cosmo/rmekuri/yt-x86_64/src/yt-hg/yt/visualization/plot_types.py",
>>
>> line 406, in _redraw_image
>>
>> self._run_callbacks()
>>
>> File
>>
>> "/home/uni06/cosmo/rmekuri/yt-x86_64/src/yt-hg/yt/visualization/plot_types.py",
>>
>> line 259, in _run_callbacks
>>
>> cb(self)
>>
>> File
>>
>> "/home/uni06/cosmo/rmekuri/yt-x86_64/src/yt-hg/yt/visualization/plot_modifications.py",
>>
>> line 985, in __call__
>>
>> & ( reg[field_y] >= y0 ) & ( reg[field_y] <= y1 ) )
>>
>> File
>>
>> "/home/uni06/cosmo/rmekuri/yt-x86_64/src/yt-hg/yt/data_objects/data_containers.py",
>>
>> line 319, in __getitem__
>>
>> self.get_data(key)
>>
>> File
>>
>> "/home/uni06/cosmo/rmekuri/yt-x86_64/src/yt-hg/yt/data_objects/data_containers.py",
>>
>> line 2410, in get_data
>>
>> if self._generate_field(field):
>>
>> File
>>
>> "/home/uni06/cosmo/rmekuri/yt-x86_64/src/yt-hg/yt/data_objects/data_containers.py",
>>
>> line 352, in _generate_field
>>
>> raise KeyError(field)
>>
>> KeyError: 'particle_position_y'
>>
>>
>> With best regards,
>>
>>
>>
>> Cheers,
>>
>>
>> Remudin
>>
>>
>> ------
>>
>>
>> Remudin Reshid Mekuria phone: +49 (0)551
>>
>> 39-13805
>>
>>
>> IAG, Friedrich-Hund-Plaz 1 fax: +49 (0)551
>>
>> 39-5043
>>
>>
>> D-37077 Goettingen
>>
>>
>> Germany
>>
>>
>>
>>
>>
>>
>> _______________________________________________
>>
>> 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
>>
>>
>>
>> ------
>>
>>
>> Remudin Reshid Mekuria phone: +49 (0)551
>>
>> 39-13805
>>
>>
>> IAG, Friedrich-Hund-Plaz 1 fax: +49 (0)551 39-5043
>>
>>
>> D-37077 Goettingen
>>
>>
>> Germany
>>
>>
>>
>>
>>
>>
>> _______________________________________________
>>
>> 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
>>
>>
>> ------
>>
>> Remudin Reshid Mekuria phone: +49 (0)551
>> 39-13805
>>
>> IAG, Friedrich-Hund-Plaz 1 fax: +49 (0)551 39-5043
>>
>> D-37077 Goettingen
>>
>> Germany
>>
>>
>>
>>
>>
>> _______________________________________________
>> 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
------
Remudin Reshid Mekuria phone: +49 (0)551 39-13805
IAG, Friedrich-Hund-Plaz 1 fax: +49 (0)551 39-5043
D-37077 Goettingen
Germany
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.spacepope.org/pipermail/yt-users-spacepope.org/attachments/20120524/ec222b37/attachment.html>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: inputs
Type: application/octet-stream
Size: 3082 bytes
Desc: not available
URL: <http://lists.spacepope.org/pipermail/yt-users-spacepope.org/attachments/20120524/ec222b37/attachment-0001.obj>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.spacepope.org/pipermail/yt-users-spacepope.org/attachments/20120524/ec222b37/attachment-0001.html>
More information about the yt-users
mailing list