yt-users August 2020

yt-users@python.org

6 participants
8 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

comoving scale annotation?
by Salome Mtchedlidze 05 Dec '20

05 Dec '20

Hello again, I am using yt slices and projections with scale annotation. I have my plots in comoving units ("Mpccm/h") but whenever I want to use the same units for annotation: slc.annotate_scale(corner='lower_right', coeff=10, units='Mpccm/h'), yt does not recognize these units. Am I missing something or only the physical scale annotation is possible ? Thanks, Salome

3 10

spatial positions in FLASH code
by Thibault Goudal 26 Aug '20

26 Aug '20

Hi! Working on the FLASH code (dynamic mesh refinement), I am looking for how could I get from the HDF5 Outputs the spatial positions (in physical units). Idealy, if I run a 2D simulation, I would like to obtain x and y arrays at a given time. It is not necessarally a YT issue but it is obvious that some modules could help a lot. Tell me if you have any idea! Thanks -- Thibault Goudal PhD Student CEA CELIA CNRS UMR5107 laboratory +33664160842 thibault.goudal(a)u-bordeaux.fr

2 1

annotate_quiver
by Prateek Gupta 10 Aug '20

10 Aug '20

Dear all, I am trying to annotate quiver over the projection plot of density field using the following code line: * p.annotate_quiver("velocity_x", "velocity_y", factor = 35, normalize = True, plot_args = {"color":'Red'})* but it shows an *error*: File "/home/user1/yt/yt/visualization/plot_window.py", line 997, in run_callbacks callback = CallbackMaker(*args[1:], **kwargs) *TypeError: __init__() got an unexpected keyword argument 'plot_args'* Without *'plot_args'* keyword argument annotation works fine but I need to adjust the color and the width of the vectors. I am using yt Version = 3.5.dev0. It may be a silly question: while annotating vectors (velocity_x and velocity_y, say) over the projection plot, is the annotation done after the weighted line integration of velocity_x and velocity_y over the width of the projection? and for this also the weight field is the same as we specified in the field projection method! If so, Is there any way to change the weight field for the annotation? Many thanks and regards, -Prateek

1 0

Adding divergence field
by Michael Jennings 07 Aug '20

07 Aug '20

Hi all, I am attempting to add a divergence field and am running into the issue apparent in the plot below: The grid artifacts appear to be the meshblocks for mpi parallelization when running the simulation. This doesn’t occur for any other field. I am working with an Athena++ athdf dataset and using yt version 3.7.dev0. Here is the script I used to add the desired field: import yt from yt.fields.derived_field import ValidateSpatial import yt.units as u kappa0 = 9.68e4 ds = yt.load('../pleiades/grf/grf_pr_2_3_nm5_3/TI.uov.00200.athdf',\ units_override=units_override,unit_system="cgs") grad_fields = ds.add_gradient_fields(("gas","temperature")) def _therm_cond(field,data): kappa = kappa0*u.erg/(u.cm*u.K*u.s) sl_left = slice(None, -2, None) sl_right = slice(2, None, None) div_fac = 2.0 ds = div_fac * data['dx'].flat[0] f = data["temperature_gradient_x"][sl_right,1:-1,1:-1]/ds f -= data["temperature_gradient_x"][sl_left ,1:-1,1:-1]/ds if data.ds.dimensionality > 1: ds = div_fac * data['dy'].flat[0] f += data["temperature_gradient_y"][1:-1,sl_right,1:-1]/ds f -= data["temperature_gradient_y"][1:-1,sl_left ,1:-1]/ds if data.ds.dimensionality > 2: ds = div_fac * data['dz'].flat[0] f += data["temperature_gradient_z"][1:-1,1:-1,sl_right]/ds f -= data["temperature_gradient_z"][1:-1,1:-1,sl_left ]/ds new_field = np.zeros(data["temperature_gradient_x"].shape, dtype='float64') new_field = data.ds.arr(new_field,'K/cm**2') new_field[1:-1,1:-1,1:-1] = f return kappa*new_field/data['density'] ds.add_field(('gas','therm_cond'), function=_therm_cond, units="erg/(g*s)", sampling_type="cell", take_log=False, display_name=r'$\nabla \cdot (\kappa\nabla T)/\rho$', force_override=True, validators = [ValidateSpatial(1, ['temperature_gradient_x','temperature_gradient_y','temperature_gradient_z'])]) Has anyone come across this issue before? Is this a known bug or am I just implementing the divergence field incorrectly? Any help is greatly appreciated. Thanks, Michael

2 4

Using 'density' method for smooth() function
by clp3ef＠virginia.edu 06 Aug '20

06 Aug '20

Hi all, I'm trying to use the 'density' method for the smooth() function (which points to the class SmoothedDensityEstimate), but I can't find any documentation about how to use the 'density' method of smooth (i.e. what fields, index fields, etc. I should feed it.) Could someone give me an example of how this method should be used? I'm using yt 3.6 with FLASH particle data. The link to the SmoothedDensityEstimate class code is https://github.com/yt-project/yt/blob/6befef2bc0427250fd62395962599be41b193…. Thanks, Camryn

2 3

FLASH Nearest Neighbor - yt 3.6.0
by clp3ef＠virginia.edu 05 Aug '20

05 Aug '20

Hi! I'm working with a FLASH dataset and am trying to use the add_nearest_neighbor_field() function as briefly described here https://yt-project.org/doc/analyzing/fields.html#computing-the-nth-nearest-…. However, when I run the short script (with the input data replaced with my own FLASH dataset), I get the errors yt.fields.field_exceptions.NeedsGridType: (0, None) and NotImplementedError. Is this function not available for FLASH datasets, and is there a workaround? My final goal is to create a smoothed particle field, so if there is a better way than using the nearest_neighbor_distance_64 field in place of a smoothing_length field, please let me know! Input data: plt file - http://use.yt/upload/a94cdffd part file - http://use.yt/upload/dec62ce5 Test Script: http://paste.yt-project.org/show/259/ Camryn

2 1

Ellipticity/axes of the 3D
by Robin Eappen 05 Aug '20

05 Aug '20

Hi, Just wondering if it is possible to find the ellipticity of the model. I assume one should first start with all the (x,y,z) coordinates and fit them to an ellipsoid. Is it already implemented in yt? Thank you Robin

2 1

2022

2021

2020

2019

2018

2017

2016

2015

2014

2013

2012

2011

2010

2009

2008

2007

yt-users August 2020