Mailman 3 April 2012 - yt-users

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

Re: [yt-users] Luminosity within r
by Massimo Gaspari 18 Jul '12

18 Jul '12

Matt, Sam, Eric, Britton, Thank you for all the previous infos (projections and 2D profiles). Now I have a clearer vision of these features. Sam, you are right, the derived field is created before the radial binning.I will let you know if I have further doubts (for sure). As I said to Matt previously, I want to extend my congratulations to all the YT developers.You are doing a great job. Cheers, Max

2 1

problem with fields that use ghost zones
by Britton Smith 30 May '12

30 May '12

Hi everyone, I'm having a problem using fields that use ghost zones. The following simple script: http://paste.yt-project.org/show/2010/ gives this error: http://paste.yt-project.org/show/bOikDPScBBtDiUGvH11X/ I am working from the tip, but I get the same behavior from yt/2.3. In yt/2.2, everything is working. I am working now to narrow that range down a bit, but does anyone have an idea? Britton

8 19

yt Google+ Hangout Tomorrow
by Sam Skillman 01 May '12

01 May '12

Hi all, Tomorrow we are going to try out a Google Hangout hosted by +yt. This is a chance for you to ask questions or iterate on a particular piece of analysis you'd like to perform. Also feel free to give feedback or just stop by and say hi to us. We'll be starting around 2pm EST tomorrow and there should be several of the yt developers hanging out. Depending on how it goes, we'll consider running one ever so often as a form of office hours or discussion of new capabilities. Anyways, if you have Google+ enabled and you head over and add +yt to your circles at https://plus.google.com/107728486871834552760/posts, you should see the hangout pop up in your feed. We'll also probably send out the link on the irc channel when we get started. I hope to see you there! Sam

3 3

Light Cone + All Sky
by Richard P Wagner 29 Apr '12

29 Apr '12

Hi, I'm looking at the light cone and all sky examples, and I'm wondering how much effort (in human units, not computer time) is required to produce an all sky map based on the "shift and stack" method used for light cone generation. Naively, the answer seems to be based on whether or not the light cone generator projection function could use the HEALpix camera. Is there a simple way to create an all sky map combining several data sets? Thanks, Rick

3 3

Making custom plots consistent with yt's styles
by Christopher Moody 29 Apr '12

29 Apr '12

Hi all, Aesthetically, it'd be nice to make matplotlib plots with a consistent style as yt-made plots so that fonts, sizes, color palettes and so on agree. I imagine that this means I have to find yt's default matplotlib rc, or something analogous, somewhere. Many thanks for the help! chris

4 3

Radial profile of 2D map
by Massimo Gaspari 23 Apr '12

23 Apr '12

Hi guys, suppose I performed a 2D data projection of this form: frb = proj.to_frb(width, res, center=c) how can I obtain the radial profile of this 2D "map" (say frb["Temperature"])?(I have found only how to perform radial profiles of 3D data, setting a certain sphere) Do you have some template? Thanks in advance! Max

5 4

Luminosity within r
by Massimo Gaspari 23 Apr '12

23 Apr '12

Hi guys, I am trying to calculate the (X-ray) luminosity within a sphere of radius r:Lx = sum emissivity *dvol (within r) Suppose I have defined a new derived field XRayEmissivity, dependent on data['Density'] and data['Temperature']. How would you carry out that? I am trying to use BinnedProfile1D. However, it seems that the derived fields are computed only *after* binning, e.g.XRayEmissivity = <Density>**2 * <Temperature>^0.5 I don't want to use the averaged quantities (in each shell). I would like to compute the emissivity point by point, for each leaf block(very important for the luminosity). Is it possible? Thanks, as always! Max

2 1

volume rendering job dies
by Geoffrey So 17 Apr '12

17 Apr '12

Hi all, I'm starting to do some volume rendering of my 800 cube simulation, but the job in the super computer just dies without any kind of error in the logs. The last few lines are: yt : [INFO ] 2012-04-15 16:15:22,264 Warning: no_ghost is currently True (default). This may lead to artifacts at grid boundaries. yt : [INFO ] 2012-04-15 16:15:33,307 Max Value is 2.46548e-24 at 0.9743750000000000 0.0193750000000000 0.1781250000000000 in grid EnzoGrid_0442 at level 0 (79, 15, 42) yt : [INFO ] 2012-04-15 16:15:33,309 Warning: no_ghost is currently True (default). This may lead to artifacts at grid boundaries. yt : [INFO ] 2012-04-15 16:15:47,133 Max Value is 2.46548e-24 at 0.9743750000000000 0.0193750000000000 0.1781250000000000 in grid EnzoGrid_0442 at level 0 (79, 15, 42) yt : [INFO ] 2012-04-15 16:15:47,135 Warning: no_ghost is currently True (default). This may lead to artifacts at grid boundaries. yt : [INFO ] 2012-04-15 16:15:49,839 Max Value is 2.46548e-24 at 0.9743750000000000 0.0193750000000000 0.1781250000000000 in grid EnzoGrid_0442 at level 0 (79, 15, 42) yt : [INFO ] 2012-04-15 16:15:49,841 Warning: no_ghost is currently True (default). This may lead to artifacts at grid boundaries. I'm wondering if it is again an memory issue. I ran this job with 16 nodes with 16 cores per node and 4GB per core. This setup is enough memory to do the halo finding, I'm wondering what the memory requirements is for volume rendering, and what's the biggest volume people have done up to date? The script works on my tiny test problem of 64 cube: http://paste.yt-project.org/show/2289/ If this isn't a problem people have with rendering larger volumes than mine, then it might be a machine issue. I haven't moved the data to another super computer, but will try that next if others have done larger volume rendering with no problems at all. From G.S.

4 5