Mailman 3 July 2019 - yt-users

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

14 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

Ytini installation for Houdini on Windows
by uran san 16 Dec '19

16 Dec '19

Hi guys, I am very new to Ytini's mailing list, mailing lists as well as Houdini. So please bear with me while I pick up the etiquette around this medium. With that out of the way here's my first question: After having succeeded installing Ytini for Houdini on a Mac OS X 10.14 I would now like to achieve the same with my Houdini install on a Windows 10 OS. I have tried to somewhat reproduce all steps described here http://www.ytini.com/getstarted.html#collapseTwoOut (Houdini 17.5 Tested on Mac OS X) on my Windows machine but wasn't surprised to learn that of course Houdini could not find the ytini module when I tried to import it in the python shell. Is there somebody, that could point me . in the right direction of how to install Ytini for Houdini on a Windows 10 OS environment? Not afraid of reading so even hints a re more than welcome. Many Thanks in advance Ur4n

6 23

HOP vs Rockstar virial radii
by Nick Gnedin 07 Aug '19

07 Aug '19

Folks, I ran both HOP and Rockstar halo finders on the same ART simulations, and the mass functions match very well (attached hmf_m.pdf file). However, when I compare virial radii (attached hmf_r.pdf file), they don't match. Even if I try rescaling the Rockstar virial radius by 1.5, the match is much worse than for mass functions. How are the 'virial_radius' fields defined for both halo finders? This is a z=7 simulations, so all mass definitions are the same, Mvir = M200m = M200c. I also attach the script to plot virial radii if you want to check it, although it is trivial. Many thanks (as usual), Nick

2 1

Question about YT Units Conversion
by Sean Larkin 07 Aug '19

07 Aug '19

Hello All, I am currently attempting to compare a catalog with Mvir and Rvir units of Mpc/h and Kpc/h comoving to a catalog with Mpc and Kpc proper. I am currently using yt version 3.5.1, which I just checked for updates. When doing this conversion of units within yt, there appears to be a problem where an amount of h*h is not handeled properly. I have included a screenshot of a jupyter session that outlines the problem as well as a .py and jupyter script so it can be tested. The problem arises when I try to calculate delta_vir, which is the density of a galaxy/ critical density. The correct value for both catalogs I am comparing is around 170. When I use yt to calculate the critical density at z = 2, it get the correct value, which I convert to Mpc*h*h//kpc/kpc/kpc to match with my /h catalog. Then, when I divide the densities, even thought their units are identical, their division is different than if just the float values are divided, by an exact factor of h*h. I wonder if there is something wrong with how I am initializing my values with the ytquantity, or if there is something wrong going on in the code. [image: Screen Shot 2019-07-31 at 2.36.01 PM.png]

2 1

Lacking `particle_type` in yt 3.4.1 under Python 3.7
by Yi's Toy Repositories 28 Jul '19

28 Jul '19

I recently installed Anaconda3 with Python 3.7, as well as the latest yt 3.4.1. But when I ran my code, calculating the center of mass and angular momentum, an error saying: count_values() got an unexpected keyword argument 'particle_type' came up. Is there anything I can do to fix it except installing an Anaconda with Python 3.6?

4 7

Increasing number of CPU cores can not speed up the works
by Tseng, Po-Hsun 18 Jul '19

18 Jul '19

Dear yt developers, I found that increasing number of CPU cores in a single node that have 16 cores and 16 additional hyper-threadings can dramatically deteriorate speed-up. figure [3] shows this result. You may reproduce figure [3] by the following steps: Running shell script [1] can automatically invoke yt script [2] and record elapsed time on files, called core_* [4]. [1] shell script: http://paste.yt-project.org/show/150/ [2] yt script: http://paste.yt-project.org/show/151/ [3] # of cores vs. elapsed time http://i.imgur.com/mFk6lKu.png [4] all scripts, data and log files: http://use.yt/upload/d8efa692 Thanks for your help. Rocky Tseng

3 4

Python 2 support removed from yt-4 branch
by Matthew Turk 13 Jul '19

13 Jul '19

Hi all, With the merger today of https://github.com/yt-project/yt/pull/2163 , Python 2 support has been removed from the yt-4 branch. Those of you using the master branch or any release will not be affected, but if you're on yt-4, keep this in mind. Best, Matt

1 0

cut_region() not working as intended
by Andrew Wu 11 Jul '19

11 Jul '19

Hi yt, When making projection plots of total gas density and cold gas density, averaged over the projection length (i.e. with a uniform weighting), and using cut_region() to filter out hot gas to make the cold gas plots, it appears that some areas in the cold gas plot have a higher density than the corresponding areas in the total gas plot, which shouldn't happen. If you look at the top-middle of the two plots (link below) and compare the densities, the cold gas appears to have a higher density there. The code is below: (link to plots: https://drive.google.com/drive/folders/1UJ6fmKNFQzGr8cfjov5-GwOzfS8nt2e8?us…) import yt ds = yt.load("/tigress/cen/CZ15/redshift0100") z = 0.62 # the box size is 120/0.7/(1+z) Mpc, where z is redshift boxSize = 120e3/0.7/(1+z) # virial radius in code units, r_v is virial radius in kpc codeWidth = r_v / boxSize center = [xCoord, yCoord, zCoord] left_edge = [xCoord-2.1*codeWidth, yCoord-2.1*codeWidth, zCoord-2.1*codeWidth] right_edge = [xCoord+2.1*codeWidth, yCoord+2.1*codeWidth, zCoord+2.1*codeWidth] # want a box with full width about 4 times the virial radius reg = ds.region(center, left_edge,right_edge) # total gas plot prj = yt.ProjectionPlot(ds, "x", "density", center=center, width=(4*r_v, 'kpc'), data_source=reg, method='integrate', weight_field='ones') u = 100 / 0.048 l = 0.07 / 0.048 prj.set_zlim('density', l,u) prj.set_unit('density','mbd') prj.set_cmap(field='density', cmap='CMRmap') prj.show() # cold gas plot cold_reg = reg.cut_region(["obj['temperature'] < 1e5"]) prj = yt.ProjectionPlot(ds, "x", "density", center=center, width=(4*r_v, 'kpc'), data_source=cold_reg, method='integrate', weight_field='ones') prj.set_zlim('density', l,u) prj.set_unit('density','mbd') prj.set_cmap(field='density', cmap='CMRmap') prj.show() Thanks, Andrew

2 4