yt-users October 2018

yt-users@python.org

14 participants
26 discussions

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

ACCESSING GAS FILEDS THROUGH ARBITRARY_GRID
by Vadlamani Samhitha 05 Nov '18

05 Nov '18

Hi, I am using arbitrary_grid to access my data without flattening: > dd=ds.arbitrary_grid(left_edge=[(com_x_1.value-30), (com_y_1.value-30), > (com_z_1.value-30)]*kpc, right_edge=[(com_x_1.value+30), (com_y_1.value+30), > (com_z_1.value+30)]*kpc, dims=[128, 128, 128]) > and as I attempt to get the gas fields for this object, > print(dd['gas', 'velocity_x']) > but end up getting the error: (which isn't the case for deposit fields) yt.units.unit_registry.UnitParseError: Code units have not been defined. > Try creating the array or quantity using ds.arr or ds.quan instead. If I'm thinking right, is it because the arbitrary_grid is only applied to the particles. In that case, is there any equivalent for gas fields where I can access my 3D data for the gas as it is , without flattening? Thanks

2 2

Re: using young stars to calculate the angular momentum vector
by Nathan Goldbaum 31 Oct '18

31 Oct '18

Hi Sam, This is a new feature in yt 3.5.0, you’ll need to update to be able to use it. Let us know if you have questions about updating. Nathan On Tue, Oct 30, 2018 at 2:16 PM Samantha M Benincasa < smbenincasa(a)ucdavis.edu> wrote: > Hi everyone, > > I'm trying to get a halo into a face-on and side-on configuration using > sp.quantities.angular_momentum_vector. > > Is there a way to do this calculation with only the young stars in the > sphere? The docs suggest that angular_momentum_vector will take a > "particle_type" keyword, but my version does not recognize this. > > Cheers, > Sam > _______________________________________________ > yt-users mailing list -- yt-users(a)python.org > To unsubscribe send an email to yt-users-leave(a)python.org >

2 1

using young stars to calculate the angular momentum vector
by Samantha M Benincasa 30 Oct '18

30 Oct '18

Hi everyone, I'm trying to get a halo into a face-on and side-on configuration using sp.quantities.angular_momentum_vector. Is there a way to do this calculation with only the young stars in the sphere? The docs suggest that angular_momentum_vector will take a "particle_type" keyword, but my version does not recognize this. Cheers, Sam

1 0

ACCESSING DATA AS IT IS
by Vadlamani Samhitha 24 Oct '18

24 Oct '18

Hi, I would like to access data from Ramses as it is because I would want to divide my selected region into subgrids and compute some quantity locally in those subgrids. However, using all_data or ds.r would rearrange the information as it flattens the data (please correct me if I got this wrong). Is there a way in yt, where I can simply access the 3D data as it is and divide into subcubes? Thank you.

2 3

Postdoc in Comp. Astro. at MSU
by Sean M. Couch 22 Oct '18

22 Oct '18

Hi yt-users, We have an open postdoctoral research associate position in computational astrophysics/supernova theory at Michigan State. More info and application portal here: http://careers.msu.edu/cw/en-us/job/500160/research-associatefixed-term. Please pass along to anyone you think might be interested. Apologies if you receive this announcement multiple times. Thanks! Sean ---------------------------------------------------------------------------------------------------------- Sean M. Couch, Ph.D. Assistant Professor Department of Physics and Astronomy Department of Computational Mathematics, Science, and Engineering National Superconducting Cyclotron Laboratory/Facility for Rare Isotope Beams Michigan State University 567 Wilson Rd, 3250 BPS East Lansing, MI 48824 (517) 884-5035 —— couch(a)pa.msu.edu<mailto:couch@pa.msu.edu> —— www.pa.msu.edu/~couch<http://www.pa.msu.edu/~couch>

1 0

BUILT-IN VELOCITY DISPERSION COMPUTATION
by Vadlamani Samhitha 22 Oct '18

22 Oct '18

Hi, I am still a bit new to yt and am trying to compute local velocity dispersions of the gas in my simulation by dividing the entire volume into ‘n’ number of sub grids. Does yt offer any functionality that could perform this in a simple fashion? (I looked up the documentation and couldn’t find any such function). Thank you.

1 0

Error while deriving vorticity field
by Prateek Gupta 22 Oct '18

22 Oct '18

Dear yt-users I am creating a derived field for new_vorticity_magnitude. But when I take the SlicePlot of this derived field I am getting the following error: ############################################### slc = yt.SlicePlot(ds, 'z', "Vorticity_mag", center= center_old , width = (4,'Mpccm/h')) File "/home/dumbo/yt/yt/visualization/plot_window.py", line 1953, in SlicePlot return AxisAlignedSlicePlot(ds, normal, fields, *args, **kwargs) File "/home/dumbo/yt/yt/visualization/plot_window.py", line 1266, in __init__ slc.get_data(fields) File "/home/dumbo/yt/yt/data_objects/data_containers.py", line 1336, in get_data self._generate_fields(fields_to_generate) File "/home/dumbo/yt/yt/data_objects/data_containers.py", line 1356, in _generate_fields fd = self._generate_field(field) File "/home/dumbo/yt/yt/data_objects/data_containers.py", line 318, in _generate_field tr = self._generate_fluid_field(field) File "/home/dumbo/yt/yt/data_objects/data_containers.py", line 338, in _generate_fluid_field rv = finfo(gen_obj) File "/home/dumbo/yt/yt/fields/derived_field.py", line 237, in __call__ dd = self._function(self, data) File "vorticity_test.py", line 63, in _newvorticity_x_ f = (data[ftype, "zVel"][sl_center,sl_right,sl_center] - data[ftype, "zVel"][sl_center,sl_left,sl_center]) / (div_fac*just_one(data["index", "dy"])) File "/home/dumbo/yt/yt/units/yt_array.py", line 1036, in __getitem__ ret = super(YTArray, self).__getitem__(item) IndexError: too many indices for array ############################################### The current version of yt is: --- Version = 3.5.dev0 Changeset = 0cb7ad5befb2 --- Please let me know, why I am getting this error? Best regards -Prateek

3 4

2022

2021

2020

2019

2018

2017

2016

2015

2014

2013

2012

2011

2010

2009

2008

2007

yt-users October 2018