Mailman 3 October 2013 - yt-users

Porting .ply files to Blender with colors?
by Jill P. Naiman Aug. 15, 2022

Aug. 15, 2022

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 Dec. 11, 2021

Dec. 11, 2021

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 April 11, 2020

April 11, 2020

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

Light Cone Projection
by Pengfei Chen March 2, 2014

March 2, 2014

Hi all, I downloaded the example script of light cone projectoin( http://yt-project.org/doc/cookbook/cosmological_analysis.html#light-cone-pr…) to my laptop and changed the Enzo simulation parameter file and redshifts. But when I ran it the following error occurred: File "light_cone_projection.py", line 19, in <module> find_outputs=False) File "/Users/madcpf/Documents/dev_yt/src/yt-hg/yt/analysis_modules/cosmological_observation/light_cone/light_cone.py", line 158, in __init__ find_outputs=find_outputs) File "/Users/madcpf/Documents/dev_yt/src/yt-hg/yt/analysis_modules/cosmological_observation/cosmology_splice.py", line 44, in __init__ find_outputs=find_outputs) File "/Users/madcpf/Documents/dev_yt/src/yt-hg/yt/convenience.py", line 131, in simulation find_outputs=find_outputs) File "/Users/madcpf/Documents/dev_yt/src/yt-hg/yt/frontends/enzo/simulation_handling.py", line 81, in __init__ find_outputs=find_outputs) File "/Users/madcpf/Documents/dev_yt/src/yt-hg/yt/data_objects/time_series.py", line 343, in __init__ self._parse_parameter_file() File "/Users/madcpf/Documents/dev_yt/src/yt-hg/yt/frontends/enzo/simulation_handling.py", line 273, in _parse_parameter_file param, vals = (i.strip() for i in line.split("=")) ValueError: need more than 1 value to unpack Does the parameters "near_redshift" and "far_redshift" in LightCone need to be the same with "CosmologyFinalRedshift" and "CosmologyInitialRedshift" in enzo parameter file? Does anyone know how to solve this? Thank you all! Pengfei

2 10

Running Rockstar Halo Finder to create a merger tree
by Hilary Egan Dec. 18, 2013

Dec. 18, 2013

Hi all, I'm quite confused on a number of points related to running the rockstar halo finder, so I hope its alright that I put all these questions into this one email! 1. I can't seem to run the rockstar halo finder at all without getting this error followed by a segmentation fault and crash. [Warning] Network IO Failure (PID XXXXXX): Connection reset by peer [Network] Packet receive retry count at: 1 It sort of seems like this issue ( http://lists.spacepope.org/htdig.cgi/yt-dev-spacepope.org/2012-November/002…) but I couldn't really figure out what the resolution was from the thread. Im attempting to run this on kraken and it doesn't matter if I use a single compute node or multiple, I get the same error. (I hope this isn't the infiniband issue the docs warned about, I couldn't figure out if that is how kraken is connected and I got an error that the suggested flag doesn't exist so I didn't press the issue.) 2. Whenever I finally do get the halo finder to work, I need the results to be in a form that the merger tree can use. It seems as though the MergerTree needs the results in the same form as the other halo finders give, so would getting the halo list and then dumping it as usual be the appropriate strategy? Ie: rh.run() halo_list = rh.halo_list() halo_list.dump('MergerHalos') 2.5. The docs sort of give mixed messages on whether or not I could just be calling MergerTree with the argument halo_finder_function = RockstarHaloFinder. At this point I've pretty thoroughly convinced myself that I can't, but it would be nice if that was clarified. (Just a thoroughly overwhelmed new user's perspective!) 3. I'm a little confused as to whether or not I have to use a TimeSeriesData object rather than the usual single time output when instantiating the halo finder. Under "Rockstar Halo Finding" it uses TimeSeriesData, unlike the rest of the examples, but under the subheading "Output Analysis" it just uses pf. The "Output Analysis" example also doesn't call the run() method, which leads me to believe something else entirely is going on, but its not quite clear. Thanks! -Hilary

7 8

Athena units
by Stephanie Tonnesen Nov. 22, 2013

Nov. 22, 2013

Hi Again, Everyone! So I would also like my figures to be in physical units, instead of in my code units. I *think* the Athena frontend doesn't read in what the unit conversions are (and I am not sure how it would do that). Can I just do something simple like: DensityUnits = 3. density_physunits = grid['density'] * DensityUnits Then, of course, I don't know how I would do a SlicePlot, ProjectionPlot or phaseplot in physical units... Any help/advice/direction is appreciated! Thanks, Stephanie

3 11

chombo/orion2 projection issue
by William Gray Nov. 1, 2013

Nov. 1, 2013

Hello yt-users, I am having an issue in creating off axis projection plots with short depths in yt. I am looking for disks around particles which necessitates using the off axis projection functionality. What I want to do is create a projection where the depth is the same as or similar to the width of the plot. However, I am finding that unless the depth is orders of magnitude larger than the width I will always create a blank plot (for example, the width is ~100 au, the depth to get a plot has to be ~1000xwidth, I want a plot where the width and depth are ~100 au). I have tried using both off_axis_projection and the OffAxisProjection methods with similar (bad) results. As a comparison, I ran a low resolution FLASH simulation that did not seem to have this problem. I also know that such a plot was possible with ORION2 data in the recent past (spring of this year) since a colleague was able produce such a plot. That same script now also fails. Does anyone know of a way around this problem? Has there been any recent changes that broke volume rendering for chombo/orion2 data? Thanks, William Gray

3 4

Loading RAMSES Data into yt.
by k.grisdale＠surrey.ac.uk Oct. 31, 2013

Oct. 31, 2013

Hi, I have been trying to load RAMSES data into yt 2.5. Orginaly it was loading in but the graphs I was able to produce didn’t make sense. So I tried to change to yt3 using the commands cd $YT_DEST/src/yt-hg/ hg pull http://bitbucket.org/yt_analysis/yt-3.0 hg up -C yt-3.0 python2.7 setup.py build_ext -i -f The commands appear to have worked. But when I try to use yt to load in RAMSES data I now get this error: Traceback (most recent call last): File "/Applications/Code/yt/yt-x86_64/src/yt-hg/yt/gui/reason/basic_repl.py", line 81, in execute exec code in self.locals File "", line 1, in File "/Applications/Code/yt/yt-x86_64/src/yt-hg/yt/convenience.py", line 76, in load if c._is_valid(*args, **kwargs): candidates.append(n) File "/Applications/Code/yt/yt-x86_64/src/yt-hg/yt/frontends/art/data_structures.py", line 427, in _is_valid amr_header_vals = read_attrs(fh, amr_header_struct, '>') File "/Applications/Code/yt/yt-x86_64/src/yt-hg/yt/utilities/fortran_utils.py", line 62, in read_attrs vals = list(struct.unpack(net_format, f.read(size))) error: unpack requires a string argument of length 892 Any idea what I need to do? Thanks Kearn

3 13

Re: [yt-users] Scipy installation failure
by Yong Zheng Oct. 30, 2013

Oct. 30, 2013

Dear yt-developers: I am trying to install scipy by using the install script, but it keeps failing all the time. Could you help me to check what's wrong with it? Thank you very much for your help! Here are the last few lines in the log file. Installed /home/yzheng/yt-x86_64/lib/python2.7/site-packages/distribute-0.6.32-py2.7.egg Processing dependencies for distribute==0.6.32 Finished processing dependencies for distribute==0.6.32 After install bootstrap. /home/yzheng/yt-x86_64/lib/python2.7/site-packages/setuptools-0.6c11-py2.7.egg-info already exists Searching for pip Best match: pip 1.4.1 Processing pip-1.4.1-py2.7.egg pip 1.4.1 is already the active version in easy-install.pth Installing pip script to /home/yzheng/yt-x86_64/bin Installing pip-2.7 script to /home/yzheng/yt-x86_64/bin Using /home/yzheng/yt-x86_64/lib/python2.7/site-packages/pip-1.4.1-py2.7.egg Processing dependencies for pip Finished processing dependencies for pip Running from scipy source directory. blas_opt_info: blas_mkl_info: libraries mkl,vml,guide not found in ['/home/yzheng/yt-x86_64/lib', '/usr/local/lib', '/usr/lib', '/usr/lib/x86_64-linux-gnu'] NOT AVAILABLE atlas_blas_threads_info: Setting PTATLAS=ATLAS libraries ptf77blas,ptcblas,atlas not found in ['/home/yzheng/yt-x86_64/lib', '/usr/local/lib', '/usr/lib', '/usr/lib/x86_64-linux-gnu'] NOT AVAILABLE atlas_blas_info: libraries f77blas,cblas,atlas not found in ['/home/yzheng/yt-x86_64/lib', '/usr/local/lib', '/usr/lib', '/usr/lib/x86_64-linux-gnu'] NOT AVAILABLE /home/yzheng/yt-x86_64/lib/python2.7/site-packages/numpy/distutils/system_info.py:1494: UserWarning: Atlas (http://math-atlas.sourceforge.net/) libraries not found. Directories to search for the libraries can be specified in the numpy/distutils/site.cfg file (section [atlas]) or by setting the ATLAS environment variable. warnings.warn(AtlasNotFoundError.__doc__) Best, - Yong

1 0

merger tree crash
by Latif Oct. 29, 2013

Oct. 29, 2013

Hi all, I am running merger tree on Enzo simulations following the recipe given in YT documentation. YT crashes by producing the following error. I am using yt version "b118390aa42c". Any ideas what is going on here? FOF CPU TIME: 2.408634 secs ========== Cycle 00009 (z=9.999997) ========== Parent halo = 0 --> Most massive progenitor == Halo -1 ========== Cycle 00008 (z=14.999989) ========== ========== Cycle 00007 (z=19.999992) ========== ========== Cycle 00006 (z=29.999961) ========== Traceback (most recent call last): File "enzo.py", line 33, in <module> mt.save_halo_evolution('halos.h5') File "/home/uni05/mlatif/ytnew/yt-x86_64/src/yt-hg/yt/analysis_modules/halo_merger_tree/enzofof_merger_tree.py", line 507, in save_halo_evolution halo_id = self.levels[cycle_fin][0].halo_id IndexError: list index out of range Thanks in advance, Cheers Latif

2 2