yt-users August 2016

yt-users@python.org

23 participants
31 discussions

by david collins

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.

5 months, 1 week

Registration open for the Learn yt Workshop: Oct 10-12, 2016

by Nathan Goldbaum

Hi all, Registration is now open for the Learn yt Workshop at the University of Illinois Urbana-Champaign, October 10-12, 2016. This workshop will introduce attendees to the yt project (yt-project.org), a python toolkit for the analysis and visualization of volumetric data. We have create a website for the workshop with more information: http://yt-project.org/workshop2016/ In addition, if you wold like to register for the workshop, please fill out this online form: https://goo.gl/forms/c6oIzNQywU1YWgOe2 The workshop will cover basic usage of yt, including the yt data model, yt's field system, low-level data inspection, basic visualization workflows, and data analysis and reduction tasks. We will also cover more advanced usages like generating simulated observations, halo finding and halo analysis, volume rendering, advanced 3D visualizations, and advanced data analysis and reduction tasks. Finally, we will cover how to modify and extend yt, as well as the development and contribution process. In addition, there will be time set aside for exploring data you bring to the workshop, along with opportunities to work directly with yt developers to load and explore your data. The workshop will take place at the National Center for Supercomputing Applications building on the north end of the UIUC campus. The NCSA building is about a block away from the conference hotel and is next door to a parking structure that offers metered all-day parking. There are a number of food trucks nearby, a university-run cafeteria about 2 blocks away, and a university business district with many good lunch restaurants about a half a mile away. We are planning to offer funding for hotel and travel for those requesting support. If you request funding, you will be notified of available funds by September 15. Travel awards will be made in the form of arranged lodging and airfare, with reservations being made by the conference organizers. We hope to see you there. On behalf of the organizing committee, Matt Turk Nathan Goldbaum Jill Naiman John ZuHone Kandace Turner

3 years, 12 months

profile object / data container question

by Molly Peeples

Hi folks, I’m trying to construct and save a 2d profile container to make a phase plot. Here’s what I have: nhT_prof_new = yt.create_profile(sph, [('gas', 'H_number_density'),('gas','temperature')], fields='cell_mass', extrema=extrema,fractional=True, n_bins=6) nhT_prof_new.save_as_dataset() prof_ds = yt.load("snap_047_z000p000_Profile2D.h5") (just using a small number of bins for now to play with…) In [115]: prof_ds.data["temperature"] ...: prof_ds.data["H_number_density"] ...: prof_ds.data["cell_mass"] ...: Out[115]: YTArray([[ 0.01685665, 0.0405085 , 0.03575519, 0.03091732, 0.02552936, 0.0230379 ], [ 0.03770747, 0.05441198, 0.04859202, 0.03285325, 0.02385346, 0.04674817], [ 0.03022001, 0.05150457, 0.03591368, 0.01974481, 0.01052408, 0.00331175], [ 0.04278126, 0.01889035, 0.0177686 , 0.01644028, 0.01598129, 0.00578262], [ 0.05010083, 0.06573628, 0.04247398, 0.03914208, 0. , 0. ], [ 0. , 0.11691224, 0. , 0. , 0. , 0. ]]) dimensionless but when I try to create a PhasePlot: In [116]: p = yt.PhasePlot(prof_ds.data, "H_number_density", "temperature","cell_mass",weight_field=None) ...: --------------------------------------------------------------------------- KeyError Traceback (most recent call last) <ipython-input-116-37b81d75b563> in <module>() ----> 1 p = yt.PhasePlot(prof_ds.data, "H_number_density", "temperature","cell_mass",weight_field=None) /Users/molly/anaconda2/envs/astroconda/lib/python2.7/site-packages/yt/visualization/profile_plotter.pyc in __init__(self, data_source, x_field, y_field, z_fields, weight_field, x_bins, y_bins, accumulation, fractional, fontsize, figure_size) 756 757 type(self)._initialize_instance(self, data_source, profile, fontsize, --> 758 figure_size) 759 760 @classmethod /Users/molly/anaconda2/envs/astroconda/lib/python2.7/site-packages/yt/visualization/profile_plotter.pyc in _initialize_instance(cls, obj, data_source, profile, fontsize, figure_size) 773 obj.profile = profile 774 super(PhasePlot, obj).__init__(data_source, figure_size, fontsize) --> 775 obj._setup_plots() 776 obj._initfinished = True 777 return obj /Users/molly/anaconda2/envs/astroconda/lib/python2.7/site-packages/yt/visualization/profile_plotter.pyc in _setup_plots(self) 854 cax = self.plots[f].cax 855 --> 856 x_scale, y_scale, z_scale = self._get_field_log(f, self.profile) 857 x_title, y_title, z_title = self._get_field_title(f, self.profile) 858 /Users/molly/anaconda2/envs/astroconda/lib/python2.7/site-packages/yt/visualization/profile_plotter.pyc in _get_field_log(self, field_z, profile) 813 814 def _get_field_log(self, field_z, profile): --> 815 zfi = profile.field_info[field_z] 816 if self.x_log is None: 817 x_log = profile.x_log KeyError: ('data', 'cell_mass') Suggestions? —Molly

4 years

phase plot background color

by Molly Peeples

Hi all, Another phase plot question: is there a way to get the background to be white instead of this weird purpley color? FWIW, the colormap goes down to a very pale grey, not white. [cid:3A7656EE-7769-4312-843B-2CD6051619ED] Thanks, —Molly

4 years

And Google Summer of Code has ended this year!!!

by DVD PS

This year has been an amazing one! This year was the first one in which we've participated as an independent organisation covering many of the astronomy open source efforts under the Open Astronomy umbrella organisation. http://openastronomy.org/members/ Under the Open Astronomy organisation we have been accepted by Google directly, while in previous years the different organisations had participated under the Python Software Foundation or NumFocus. This is great because it gives us a higher visibility for students and in the open source communities, and in addition other important advantages: we get some compensation from Google and two of us will represent the organisation in the Google Summer of Code mentor summit happening at the end of October. We got 8 slots from Google, of which 5 went to Astropy and 3 to SunPy. Sadly the other sub-organisations (yt, chiantipy, juliaAstro, ...) didn't get any qualified student applicants this time but surely it will change in the future. From these 8 slots, 7 of them passed the mid-term and all the 7 passed the final evaluation. The students, their projects and the mentors are the following: For the Astropy project: Micky Costa: Astropy: Bridge Sherpa and astropy fitting <http://myopensauceadventure.blogspot.co.uk/2016/08/dessert.html> (Tom Aldcroft, Omar Laurino, and Moritz Guenther) Olga Vorokh: Image processing and source detection in Gammapy <https://github.com/search?p=1&q=author%3AOlgaVorokh+created%3A%222016-02-10…> (Johannes King and Christoph Deil) Zé Vinícius: Implement PSF photometry for fitting several overlapping objects at once <http://mirca.github.io/gsoc-astropy-final/> (Moritz Guenter, Brigitta Sipocz, and Erik Tollerud) Karl Vyhmeister: Scheduling capabilities for Astroplan <http://kvyhastroplan.blogspot.com/2016/08/cleaning-up-and-documenting.html?> (Brett Morris and Erik Tollerud) For SunPy: Tessa Wilkinson: Implementing AIA response function in Sunpy <http://tdwilkinson.blogspot.co.uk/2016/08/aia-response-summary-this-is-my-f…> (Drew Leonard and Will Barnes) Punyaslok Pattnaik: Improvements to the SunPy Database <https://punyaslokpattnaik.wordpress.com/2016/08/22/its-already-time/> (Stuart Mumford and Simon Liedtke) Sudarshan Konge: Real Time Data Access and Visualisation tools <https://sudonymousblog.wordpress.com/2016/08/19/gsoc-summary-post/> (David Perez-Suarez and Jack Ireland) You can still see all their trip that brought them till the end reading their blogs, available at: http://openastronomy.org/Universe_OA/ If you would like to give us suggestions on how to improve the program (to us as an umbrella, or to Google itself) or you want to get more involved either as a mentor or as a sub-organisation, do not hesitate to contact us. There's a lot that we can do better in the coming years!! Thanks for all the summer! David Perez-Suarez, Tom Aldcroft, Erik Tollerud, Stuart Mumford

4 years

enzo embedded parallel yt problem with gradient field

by Bodo Schwabe

*Hi all,** ** **I currently try to run enzo with embedded parallelized python/yt with the following user_script.py :* /import yt// // //from yt.frontends.enzo.api import EnzoDatasetInMemory// // //from yt.fields.field_plugin_registry import \// // register_field_plugin// //from yt.fields.fluid_fields import \// // setup_gradient_fields// // //@register_field_plugin// //def setup_my_fields(registry, ftype="enzo", slice_info=None):// // setup_gradient_fields(registry, ('enzo','AxImaginary'), '', slice_info)// // //yt.enable_parallelism()// // //def main():// // // ds = EnzoDatasetInMemory()// // ds.index// // dd = ds.all_data()// // print dd.quantities.total_quantity('AxImaginary_gradient_z')/ *It works fine on a single processor (mpirun -n 1 ...) but gives my the following error message for multiple processors:* /Global Dir set to .// //ENZO_layout 1 x 1 x 2// //Successfully read in parameter file Bosonstar.enzo.// //INITIALIZATION TIME = 3.75330448e-02// //yt : [INFO ] 2016-08-30 15:02:40,868 Global parallel computation enabled: 0 / 2// //yt : [INFO ] 2016-08-30 15:02:40,868 Global parallel computation enabled: 1 / 2// //Continuation Flag = 1// //TopGrid dt = 2.000000e-05 time = 0 cycle = 0// //P000 yt : [INFO ] 2016-08-30 15:02:41,001 Parameters: current_time = 2e-05// //P000 yt : [INFO ] 2016-08-30 15:02:41,001 Parameters: domain_dimensions = [32 32 32]// //P000 yt : [INFO ] 2016-08-30 15:02:41,002 Parameters: domain_left_edge = [-0.256 -0.256 -0.256]// //P000 yt : [INFO ] 2016-08-30 15:02:41,002 Parameters: domain_right_edge = [ 0.256 0.256 0.256]// //P000 yt : [INFO ] 2016-08-30 15:02:41,003 Parameters: cosmological_simulation = 0.0// //P000 yt : [INFO ] 2016-08-30 15:02:41,005 Gathering a field list (this may take a moment.)// // File "<string>", line 1, in <module>// // File "./user_script.py", line 21, in main// // print dd.quantities.total_quantity('AxImaginary_gradient_z')// // File "/gfs2/work/nipbschw/yt-shared-smp/src/yt-hg/yt/data_objects/derived_quantities.py", line 176, in __call__// // rv = super(TotalQuantity, self).__call__(fields)// // File "/gfs2/work/nipbschw/yt-shared-smp/src/yt-hg/yt/data_objects/derived_quantities.py", line 67, in __call__// // sto.result = self.process_chunk(ds, *args, **kwargs)// // File "/gfs2/work/nipbschw/yt-shared-smp/src/yt-hg/yt/data_objects/derived_quantities.py", line 182, in process_chunk// // for field in fields]// // File "/gfs2/work/nipbschw/yt-shared-smp/src/yt-hg/yt/data_objects/data_containers.py", line 279, in __getitem__// // File "<string>", line 1, in <module>// // File "./user_script.py", line 21, in main// // print dd.quantities.total_quantity('AxImaginary_gradient_z')// // File "/gfs2/work/nipbschw/yt-shared-smp/src/yt-hg/yt/data_objects/derived_quantities.py", line 176, in __call__// // self.get_data(f)// // File "/gfs2/work/nipbschw/yt-shared-smp/src/yt-hg/yt/data_objects/data_containers.py", line 1291, in get_data// // self._generate_fields(fields_to_generate)// // File "/gfs2/work/nipbschw/yt-shared-smp/src/yt-hg/yt/data_objects/data_containers.py", line 1311, in _generate_fields// // fd = self._generate_field(field)// // File "/gfs2/work/nipbschw/yt-shared-smp/src/yt-hg/yt/data_objects/data_containers.py", line 316, in _generate_field// // tr = self._generate_fluid_field(field)// // File "/gfs2/work/nipbschw/yt-shared-smp/src/yt-hg/yt/data_objects/data_containers.py", line 334, in _generate_fluid_field// // rv = self._generate_spatial_fluid(field, ngt_exception.ghost_zones)// // File "/gfs2/work/nipbschw/yt-shared-smp/src/yt-hg/yt/data_objects/data_containers.py", line 364, in _generate_spatial_fluid// // gz[field][ngz:-ngz, ngz:-ngz, ngz:-ngz],// // rv = super(TotalQuantity, self).__call__(fields)// // File "/gfs2/work/nipbschw/yt-shared-smp/src/yt-hg/yt/data_objects/derived_quantities.py", line 67, in __call__// // File "/gfs2/work/nipbschw/yt-shared-smp/src/yt-hg/yt/data_objects/data_containers.py", line 279, in __getitem__// // sto.result = self.process_chunk(ds, *args, **kwargs)// // File "/gfs2/work/nipbschw/yt-shared-smp/src/yt-hg/yt/data_objects/derived_quantities.py", line 182, in process_chunk// // for field in fields]// // File "/gfs2/work/nipbschw/yt-shared-smp/src/yt-hg/yt/data_objects/data_containers.py", line 279, in __getitem__// // self.get_data(f)// // File "/gfs2/work/nipbschw/yt-shared-smp/src/yt-hg/yt/data_objects/construction_data_containers.py", line 628, in get_data// // self.get_data(f)// // File "/gfs2/work/nipbschw/yt-shared-smp/src/yt-hg/yt/data_objects/data_containers.py", line 1291, in get_data// // self._generate_fields(fields_to_generate)// // File "/gfs2/work/nipbschw/yt-shared-smp/src/yt-hg/yt/data_objects/data_containers.py", line 1311, in _generate_fields// // fd = self._generate_field(field)// // File "/gfs2/work/nipbschw/yt-shared-smp/src/yt-hg/yt/data_objects/data_containers.py", line 316, in _generate_field// // tr = self._generate_fluid_field(field)// // File "/gfs2/work/nipbschw/yt-shared-smp/src/yt-hg/yt/data_objects/data_containers.py", line 334, in _generate_fluid_field// // rv = self._generate_spatial_fluid(field, ngt_exception.ghost_zones)// // File "/gfs2/work/nipbschw/yt-shared-smp/src/yt-hg/yt/data_objects/data_containers.py", line 364, in _generate_spatial_fluid// // if len(fill) > 0: self._fill_fields(fill)// // File "/gfs2/work/nipbschw/yt-shared-smp/src/yt-hg/yt/data_objects/construction_data_containers.py", line 947, in _fill_fields// // gz[field][ngz:-ngz, ngz:-ngz, ngz:-ngz],// // File "/gfs2/work/nipbschw/yt-shared-smp/src/yt-hg/yt/data_objects/data_containers.py", line 279, in __getitem__// // for chunk in ls.data_source.chunks(fields, "io"):// // File "/gfs2/work/nipbschw/yt-shared-smp/src/yt-hg/yt/data_objects/data_containers.py", line 1190, in chunks// // self.get_data(f)// // File "/gfs2/work/nipbschw/yt-shared-smp/src/yt-hg/yt/data_objects/construction_data_containers.py", line 628, in get_data// // self.get_data(fields)// // File "/gfs2/work/nipbschw/yt-shared-smp/src/yt-hg/yt/data_objects/data_containers.py", line 1279, in get_data// // fluids, self, self._current_chunk)// // File "/gfs2/work/nipbschw/yt-shared-smp/src/yt-hg/yt/geometry/geometry_handler.py", line 245, in _read_fluid_fields// // chunk_size)// // File "/gfs2/work/nipbschw/yt-shared-smp/src/yt-hg/yt/frontends/enzo/io.py", line 380, in _read_fluid_selection// // if len(fill) > 0: self._fill_fields(fill)// // File "/gfs2/work/nipbschw/yt-shared-smp/src/yt-hg/yt/data_objects/construction_data_containers.py", line 947, in _fill_fields// // for chunk in ls.data_source.chunks(fields, "io"):// // File "/gfs2/work/nipbschw/yt-shared-smp/src/yt-hg/yt/data_objects/data_containers.py", line 1190, in chunks// // self.get_data(fields)// // File "/gfs2/work/nipbschw/yt-shared-smp/src/yt-hg/yt/data_objects/data_containers.py", line 1279, in get_data// // fluids, self, self._current_chunk)// // File "/gfs2/work/nipbschw/yt-shared-smp/src/yt-hg/yt/geometry/geometry_handler.py", line 245, in _read_fluid_fields// // data_view = self.grids_in_memory[g.id][fname][self.my_slice].swapaxes(0,2)// //P001 yt : [ERROR ] 2016-08-30 15:02:41,783 KeyError: 1// // chunk_size)// // File "/gfs2/work/nipbschw/yt-shared-smp/src/yt-hg/yt/frontends/enzo/io.py", line 380, in _read_fluid_selection// // data_view = self.grids_in_memory[g.id][fname][self.my_slice].swapaxes(0,2)// //P000 yt : [ERROR ] 2016-08-30 15:02:41,784 KeyError: 2// //[hsmp16:59147] 1 more process has sent help message help-mpi-btl-openib.txt / default subnet prefix// //[hsmp16:59147] Set MCA parameter "orte_base_help_aggregate" to 0 to see all help / error messages// //[hsmp16:59147] 1 more process has sent help message help-mpi-runtime.txt / mpi_init:warn-fork// //[hsmp16:59147] 1 more process has sent help message help-mpi-api.txt / mpi-abort /*Thank you very much for any helpful suggestions,** **Bodo*/ /

4 years

volume rendering annotation error for a derived field

by tazkera haque

Hi yt people, I just updated to yt 3.3.1 to use all the new features of volume rendering. I am trying to render ejected density which is a derived field and draw the coordinate triads, domain boundary along with labels showing the log range of ejected density. my code is attached. While my code worked fine with default field 'density', the following figure was produced for 'density', but I am getting error when the source field is set to 'ejected_density'. Also the 'density' figure still does not show either the coordinate triads or the domain boundary. I would appreciate your help to fix the error and thanks in advance. *Traceback (most recent call last):* * File "volume_render_test.py", line 117, in <module>* * text_annotate=[[(.1, 1.05), text_string]])* * File "/work/03858/thaque56/sw/yt-new-3.3/yt-conda/lib/python2.7/site-packages/yt/visualization/volume_rendering/scene.py", line 395, in save_annotated* * label = rs.data_source.ds._get_field_info(rs.field).get_label()* * File "/work/03858/thaque56/sw/yt-new-3.3/yt-conda/lib/python2.7/site-packages/yt/fields/derived_field.py", line 225, in get_label* * units = Unit(self.units, registry=self.ds.unit_registry)* *AttributeError: 'NoneType' object has no attribute 'unit_registry'* *Thanks* *Tazkera*

4 years

yt 3.3.1 not working with yt 3.1 python scripts

by tazkera haque

Hi Nathan, I have just updated yt 3.3 and now I have this weird issue with the scripts , I am uploading my script for sliceplots, which should just analyse the super3d_hdf5_plt_cnt files but I have got issues which also pulls up the super3d_hdf5_part files which are also located in the same path and produces the error. *yt : [INFO ] 2016-08-26 13:23:36,817 Particle file found: super3d_hdf5_part_0000* *yt : [INFO ] 2016-08-26 13:23:36,840 integer runtime parameter checkpointfilenumber overwrites a simulation scalar of the same name* *yt : [INFO ] 2016-08-26 13:23:36,862 Parameters: current_time = 0.0* *yt : [INFO ] 2016-08-26 13:23:36,862 Parameters: domain_dimensions = [128 128 128]* *yt : [INFO ] 2016-08-26 13:23:36,863 Parameters: domain_left_edge = [ -2.80000000e+10 -2.80000000e+10 -2.80000000e+10]* *yt : [INFO ] 2016-08-26 13:23:36,863 Parameters: domain_right_edge = [ 2.80000000e+10 2.80000000e+10 2.80000000e+10]* *yt : [INFO ] 2016-08-26 13:23:36,864 Parameters: cosmological_simulation = 0.0* *yt : [INFO ] 2016-08-26 13:23:38,500 xlim = -28000000000.000000 28000000000.000000* *yt : [INFO ] 2016-08-26 13:23:38,500 ylim = -28000000000.000000 28000000000.000000* *yt : [INFO ] 2016-08-26 13:23:38,501 xlim = -28000000000.000000 28000000000.000000* *yt : [INFO ] 2016-08-26 13:23:38,501 ylim = -28000000000.000000 28000000000.000000* *yt : [INFO ] 2016-08-26 13:23:38,502 Making a fixed resolution buffer of (('flash', 'temp')) 800 by 800* *yt : [INFO ] 2016-08-26 13:23:39,016 Making a fixed resolution buffer of (('flash', 'temp')) 800 by 800* *yt : [INFO ] 2016-08-26 13:23:39,252 Saving plot Slice_z_temp/super3d_hdf5_plt_cnt_0000_Slice_z_temp.png* *yt : [INFO ] 2016-08-26 13:23:39,885 Particle file found: super3d_hdf5_part_0001* *Traceback (most recent call last):* * File "allSlices.py", line 115, in <module>* * ds=yt.load(filename)* * File "/work/03858/thaque56/sw/yt-new-3.3/yt-conda/lib/python2.7/site-packages/yt/convenience.py", line 86, in load* * return candidates[0](*args, **kwargs)* * File "/work/03858/thaque56/sw/yt-new-3.3/yt-conda/lib/python2.7/site-packages/yt/frontends/flash/data_structures.py", line 215, in __init__* * raise IOError('%s and %s are not at the same time.' % (self.particle_filename, filename))* *IOError: /work/03858/thaque56/runs/run_176/super3d_hdf5_part_0001 and /work/03858/thaque56/runs/run_176/super3d_hdf5_plt_cnt_0001 are not at the same time.* *[c558-502.stampede.tacc.utexas.edu:mpispawn_0][child_handler] MPI process (rank: 0, pid: 13788) exited with status 1* *TACC: MPI job exited with code: 1* Please take a look at the error and let me know how to avoid the particle files in this circumstance. Best Tazkera

4 years

yt update and error

by tazkera haque

Hi people, I would like to upgrade yt to the latest version, I tried typing *yt update *in a folder on cluster /work/03858/thaque56/sw/yt-3.1 and it gives me following error: Traceback (most recent call last): File "/work/03858/thaque56/sw/yt-3.1/yt-x86_64/bin/yt", line 9, in <module> load_entry_point('yt==3.2.3', 'console_scripts', 'yt')() File "build/bdist.linux-x86_64/egg/pkg_resources/__init__.py", line 558, in load_entry_point File "build/bdist.linux-x86_64/egg/pkg_resources/__init__.py", line 2682, in load_entry_point File "build/bdist.linux-x86_64/egg/pkg_resources/__init__.py", line 2355, in load File "build/bdist.linux-x86_64/egg/pkg_resources/__init__.py", line 2361, in resolve File "/work/03858/thaque56/sw/yt-3.1/yt-x86_64/src/yt-hg/yt/__init__.py", line 123, in <module> from yt.data_objects.api import \ File "/work/03858/thaque56/sw/yt-3.1/yt-x86_64/src/yt-hg/yt/data_objects/api.py", line 16, in <module> from .grid_patch import \ File "/work/03858/thaque56/sw/yt-3.1/yt-x86_64/src/yt-hg/yt/data_objects/grid_patch.py", line 19, in <module> from yt.data_objects.data_containers import \ File "/work/03858/thaque56/sw/yt-3.1/yt-x86_64/src/yt-hg/yt/data_objects/data_containers.py", line 55, in <module> from yt.utilities.parallel_tools.parallel_analysis_interface import \ File "/work/03858/thaque56/sw/yt-3.1/yt-x86_64/src/yt-hg/yt/utilities/parallel_tools/parallel_analysis_interface.py", line 33, in <module> from yt.utilities.lib.quad_tree import \ ImportError: No module named quad_tree and now I cannot excecute any of my python scripts there, please help me upgrade and fix the problem. thanks Tazkera

4 years

Loading particles

by Ariadna Murguia

Dear yt users, I'm trying to load a set of particles (a bunch of particles, so I'm running it in the hyades cluster). I have their positions and velocities. What I'm doing is just defining a box as the domain of my simulation, the data as a dictionary of velocities and positions, and just saying yt.load_particles(data=data, bbox=bbox). I tried it for some particle files that I'm using, and it works, but for others, it crashes. I'm getting the following error: /pfs/sw/python/yt/yt-x86_64/src/yt-hg/yt/units/yt_array.py:793: RuntimeWarning: invalid value encountered in divide return YTArray(super(YTArray, self).__div__(ro)) /pfs/sw/python/yt/yt-x86_64/src/yt-hg/yt/utilities/math_utils.py:861: RuntimeWarning: invalid value encountered in arccos return np.arccos( JdotCoords / np.sqrt(np.sum(coords**2,axis=0)) ) I tried it in my computer with less particles,and it works perfectly. Am I loading the particles wrong or is there a way to fix this error? Thank you very much in advance All the best, Ariadna ---------------- Ariadna Murguia Berthier University of California, Santa Cruz Department of Astronomy and Astrophysics

4 years

2020

2019

2018

2017

2016

2015

2014

2013

2012

2011

2010

2009

2008

2007

yt-users August 2016