Mailman 3 June 2018 - yt-users

Porting .ply files to Blender with colors?
by Jill P. Naiman 15 Aug '22

15 Aug '22

Hi All, I was curious if anyone has had any luck porting .ply files generated with yt into Blender with colors? It seems to port in the surfaces just fine, but they are a uniform color. From what I can gather, Blender imports .ply files with face color instead of vertex so I tried the surf.export_ply routine with "face" instead of "vertex" for sample_type, but no luck ("vertex" doesn't work either). Let me know if I can provide any more useful info. Thanks for any help, -J

3 2

Stars on volume rendering
by Elizabeth Harper-Clark 11 Dec '21

11 Dec '21

Hi all, Thanks for all your help over the last couple of days. One more question: - Can I plot particles on a volume rendered image? I have stars and I want to show where they are! Thanks, Libby -- Elizabeth Harper-Clark MA MSci PhD Candidate, Canadian Institute for Theoretical Astrophysics, UofT Sciences and Engineering Coordinator, Teaching Assistants' Training Program, UofT www.astro.utoronto.ca/~h-clark <http://www.astro.utoronto.ca/%7Eh-clark> h-clark(a)cita.utoronto.ca Astronomy office phone: +1-416-978-5759

8 14

Variance in phase/profiles
by david collins 11 Apr '20

11 Apr '20

Hi, Everybody! Does anyone out there have a technique for getting the variance out of a profile object? A profile object is good at getting <X> vs. B, I'd then like to get < (X - <X>)^2 > vs B. Matt and I had spittballed the possibility some time ago, but I was wondering if anyone out there had successfully done it. Thanks, d. -- Sent from my computer.

6 8

3-D Mass power spectra
by 梁家豪 03 Jun '19

03 Jun '19

Dear yt Can current yt calculate 3-D Mass power spectra? I checked the website but I didn't find any information. I think calculating 3-D Mass power spectra is a very useful for cosmological simulations. So I guess maybe yt supports this function now....? Thanks in advance

4 7

extracting regions
by ｔyuta 08 Jul '18

08 Jul '18

Dear yt-users, Hi, I want to know about the function "extract_connected_sets". I want to extract SuperNova enriched bubbles (which I define as Z > -4 or something) and analyze its size or mass or some other values. I'm working with cosmological simulation by enzo, and there are many SN enriched bubbles. What I need is the set of "bubble regions" of (probably) YTRegions. For this purpose, it seemed to me that "extract_connected_sets" would be the perfect function, but I don't understand what this object is, or how I can analyze the object, because when I tried projection plot of the objects, nothing appeared. Could you tell me how I can use "extract_connected_sets"? Or, is there any other way to extract SN enriched region? Best, YT

2 2

Question about star particle
by ｔyuta 28 Jun '18

28 Jun '18

Dear yt users, Hi, I want to ask about star particles in enzo simulation. I want to know whether "star particle filter" only give me "alive" star particles or also counts "dead" particles (which have exploded in the past). I defined the "star particle" with particle filter like this: import yt from yt.data_objects.particle_filters import add_particle_filter def stars(pfilter, data): filter = data[(pfilter.filtered_type, "particle_type")] == 2 return filter def DMparticles(pfilter, data): filter = data[(pfilter.filtered_type, "particle_type")] == 1 return filter yt.add_particle_filter("Stars", function=stars, filtered_type='io', requires=["particle_type"]) yt.add_particle_filter("DMparticles", function=DMparticles, filtered_type='io', requires=["particle_type"]) ds = yt.load("~/IsolatedGalaxy/galaxy0030/galaxy0030") ds.add_particle_filter('Stars') ds.add_particle_filter('DMparticles') ad = ds.all_data() stellar_mass = ad[("Stars", "particle_mass")].in_units('Msun') print(stellar_mass) Does this "stellar_mass" contain both alive and dead stars ? Best, Y.T.

2 2

Annotating BH particles on SlicePlots
by Sam Patrick 19 Jun '18

19 Jun '18

Hi all, I'm trying to annotate BH particles on top of a regular yt SlicePlot. I'm loading from an Enzo dataset and using yt version 3.3.1, the plot field and width in parsecs are given as command line arguments (I tried using density and 3.0). Here is the script I'm working with: import yt import sys from yt.units import pc, Msun, Myr import numpy as np yt.enable_parallelism() #get command line inputs, syntax is [field] [size] field = sys.argv[1] try : size = float(sys.argv[2]) except (TypeError, IndexError): print("Error, argument syntax is 'plot.py [field] [size]'") quit(0) ds = yt.load("DD0121/output_0121") time = '{:.3f}'.format(round(ds.current_time.in_units('Myr') - 194.608*Myr, 4)) val, loc = ds.find_max('density') plt = yt.SlicePlot(ds, 'z', field, center= loc, width= (size, 'pc')) plt.set_zlim(field, 1e-21, 1e-12) plt.annotate_particles(ds, size*pc, ptype= 8) #produces 1st error #plt.annotate_particles(ds, size*pc, p_size= 2.0, col= 'k', marker= 'o', stride= 1, ptype= 8, minimum_mass= None, alpha= 1.0) #produces 2nd error plt.annotate_text((0.05, 0.05), 't=' + time + 'Myr', coord_system= 'axis') plt.save("test.png") The script seems to run ok until attempting to save the plot, when I get the following error: Traceback (most recent call last): File "plot_yt.py", line 35, in <module> plt.save("test.png") File "/opt/apps/yt/3.3.1/yt-conda/lib/python2.7/site-packages/yt/visualization/plot_container.py", line 78, in newfunc args[0]._setup_plots() File "/opt/apps/yt/3.3.1/yt-conda/lib/python2.7/site-packages/yt/visualization/plot_window.py", line 949, in _setup_plots self.run_callbacks() File "/opt/apps/yt/3.3.1/yt-conda/lib/python2.7/site-packages/yt/visualization/plot_window.py", line 1005, in run_callbacks sys.exc_info()[2]) File "/opt/apps/yt/3.3.1/yt-conda/lib/python2.7/site-packages/yt/visualization/plot_window.py", line 999, in run_callbacks callback(cbw) File "/opt/apps/yt/3.3.1/yt-conda/lib/python2.7/site-packages/yt/visualization/plot_modifications.py", line 53, in _check_geometry return func(self, plot) File "/opt/apps/yt/3.3.1/yt-conda/lib/python2.7/site-packages/yt/visualization/plot_modifications.py", line 1575, in __call__ reg = self._get_region((x0,x1), (y0,y1), plot.data.axis, data) File "/opt/apps/yt/3.3.1/yt-conda/lib/python2.7/site-packages/yt/visualization/plot_modifications.py", line 1639, in _get_region LE[zax] = data.center[zax].ndarray_view() - self.width*0.5 yt.utilities.exceptions.YTPlotCallbackError: annotate_particles callback failed with the following error: unsupported operand type(s) for *: 'EnzoDataset' and 'float' I've tried adding all the keyword arguments and their parameters to the annotate_particles call (switching the commented line on my script), but this produces: Traceback (most recent call last): File "plot_yt.py", line 35, in <module> plt.save("test.png") File "/opt/apps/yt/3.3.1/yt-conda/lib/python2.7/site-packages/yt/visualization/plot_container.py", line 78, in newfunc args[0]._setup_plots() File "/opt/apps/yt/3.3.1/yt-conda/lib/python2.7/site-packages/yt/visualization/plot_window.py", line 949, in _setup_plots self.run_callbacks() File "/opt/apps/yt/3.3.1/yt-conda/lib/python2.7/site-packages/yt/visualization/plot_window.py", line 997, in run_callbacks callback = CallbackMaker(*args[1:], **kwargs) TypeError: __init__() got multiple values for keyword argument 'p_size' Any help on getting this to work would be vastly appreciated. I was also wondering whether it was possible to use the annotate_particles callback on OffAxisSlicePlots and OffAxisProjectionPlots, I've tried both and neither seems to support this functionality. Many thanks, Sam Patrick

2 2

Question about radial_velocity
by Yuan Li 15 Jun '18

15 Jun '18

Dear yt-users, I was making a plot with "radial_velocity" and noticed that it is always positive, which is odd (there should be both inflow and outflow). The source code is: >>> print(ds.field_info["gas","radial_velocity"].get_source()) def _radial(field, data): return data[ftype, "%s_spherical_radius" % basename] I find that quite confusing. It also looks the same as (‘gas’, ‘radial_magnetic_field’) according to http://yt-project.org/doc/ reference/field_list.html Thank you!

5 11

add_field r/rvir
by Stephanie Tonnesen 14 Jun '18

14 Jun '18

Hi yt-users! I am trying to add a field that is radius/rvir. This is an idealized galaxy sim with static DM potential (no live halo), so I was planning on putting the virial radius in by hand. I am not totally sure what is causing yt to choke--does it not like that I am putting in a number? See below for details. Thanks in advance for any help! Best, Stephanie The lines of code: def rrvir(field,data): return data['radius'].in_units('kpc')/(218.,'kpc') i =0 while i < len(loop): ds = yt.load("blahblah/DD"+loop[i]+"/sb_"+loop[i]) ds.add_field(('gas','r_rvir'),function=rrvir) ....many unimportant lines.... I get this error message: yt : [INFO ] 2018-06-14 09:02:48,783 Gathering a field list (this may take a moment.) yt_slices_allouts_mli.py:19: UserWarning: Because 'sampling_type' not specified, yt will assume a cell 'sampl ing_type' ds.add_field(('gas','r_rvir'),function=rrvir) Traceback (most recent call last): File "yt_slices_allouts_mli.py", line 19, in <module> ds.add_field(('gas','r_rvir'),function=rrvir) File "/home/stonnesen/yt-conda/src/yt-git/yt/data_objects/static_output.py", line 1221, in add_field deps, _ = self.field_info.check_derived_fields([name]) File "/home/stonnesen/yt-conda/src/yt-git/yt/fields/field_info_container.py", line 366, in check_derived_fi elds fd = fi.get_dependencies(ds = self.ds) File "/home/stonnesen/yt-conda/src/yt-git/yt/fields/derived_field.py", line 210, in get_dependencies e[self.name] File "/home/stonnesen/yt-conda/src/yt-git/yt/fields/field_detector.py", line 108, in __missing__ vv = finfo(self) File "/home/stonnesen/yt-conda/src/yt-git/yt/fields/derived_field.py", line 250, in __call__ dd = self._function(self, data) File "yt_slices_allouts_mli.py", line 13, in rrvir return data['radius'].in_units('kpc')/(300.,'kpc') File "/home/stonnesen/yt-conda/src/yt-git/yt/units/yt_array.py", line 1372, in __array_ufunc__ out=out, **kwargs) TypeError: ufunc 'true_divide' not supported for the input types, and the inputs could not be safely coerced to any supported types according to the casting rule ''safe'' -- Dr. Stephanie Tonnesen Associate Research Scientist CCA, Flatiron Institute New York, NY stonnes(a)gmail.com

3 2

downloading yt using all-in-one script
by Stephanie Tonnesen 11 Jun '18

11 Jun '18

Hi yt-users, I am putting yt on a new machine, and it looks like I will have to use the all-in-one script. I would like to install yt-dev. To do that, what do I need to rename the BRANCH= to? Thanks! Stephanie -- Dr. Stephanie Tonnesen Associate Research Scientist CCA, Flatiron Institute New York, NY stonnes(a)gmail.com

2 4