Re: [yt-users] Luminosity within r
by Massimo Gaspari
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
6 years, 8 months
yt Google+ Hangout Tomorrow
by Sam Skillman
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
6 years, 10 months
Light Cone + All Sky
by Richard P Wagner
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
6 years, 10 months
Making custom plots consistent with yt's styles
by Christopher Moody
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
6 years, 10 months
Radial profile of 2D map
by Massimo Gaspari
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
6 years, 11 months
Luminosity within r
by Massimo Gaspari
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
6 years, 11 months
volume rendering job dies
by Geoffrey So
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.
6 years, 11 months
(no subject)
by Nathan Goldbaum
Hi all,
I'm currently analyzing a simulation that contains a disk which is inclined with respect to the simulation box. I've figured out how to make slices through the disk midplane using the cutting plane object but when I try to slice vertically through the disk, I get back images with the disk inclined with respect to the viewing angle.
I could fix this if I could control the position angle of the cutting plane but it doesn't look like there is any way to do that currently. Looking at AMRCuttingPlaneBase the docstring says that the cutting plane tries to guess an 'up' direction and that the user cannot set it themselves. Is there a reason for that?
Thanks very much for your help,
-Nathan Goldbaum
6 years, 11 months
Enzo file visualization problem
by Abul Latif
Dear YT users,
I am running enzo test problem with MPI. I am unable to visualize the data
with YT and it produces following error. If i do not run enzo with MPI
(MPICH2) then it works fine. Any ideas, what is going wrong here.
Thanks in advance
Cheers
Latif
t : [INFO ] 2012-04-05 17:42:42,770 Parameters:
current_time = 7
yt : [INFO ] 2012-04-05 17:42:42,771 Parameters:
domain_dimensions = [16 16 16]
yt : [INFO ] 2012-04-05 17:42:42,771 Parameters:
domain_left_edge = [ 0. 0. 0.]
yt : [INFO ] 2012-04-05 17:42:42,771 Parameters:
domain_right_edge = [ 1. 1. 1.]
yt : [INFO ] 2012-04-05 17:42:42,771 Parameters:
cosmological_simulation = 0.0
yt : [INFO ] 2012-04-05 17:42:42,775 Getting the binary hierarchy
Traceback (most recent call last):
File "./enzo.py", line 7, in <module>
pc=PlotCollection(pf)
File
"/home/astro73/mlatif/YT/yt-x86_64/src/yt-hg/yt/visualization/plot_collection.py",
line 108, in __init__
v,self.c = pf.h.find_max("Density") # @todo: ensure no caching
File
"/home/astro73/mlatif/YT/yt-x86_64/src/yt-hg/yt/data_objects/static_output.py",
line 172, in hierarchy
self, data_style=self.data_style)
File
"/home/astro73/mlatif/YT/yt-x86_64/src/yt-hg/yt/frontends/enzo/data_structures.py",
line 217, in __init__
AMRHierarchy.__init__(self, pf, data_style)
File
"/home/astro73/mlatif/YT/yt-x86_64/src/yt-hg/yt/data_objects/hierarchy.py",
line 71, in __init__
self._parse_hierarchy()
File
"/home/astro73/mlatif/YT/yt-x86_64/src/yt-hg/yt/frontends/enzo/data_structures.py",
line 286, in _parse_hierarchy
if self._parse_binary_hierarchy(): return
File
"/home/astro73/mlatif/YT/yt-x86_64/src/yt-hg/yt/frontends/enzo/data_structures.py",
line 357, in _parse_binary_hierarchy
self.grid_dimensions[:] = f["/ActiveDimensions"][:]
ValueError: operands could not be broadcast together with shapes (32,3)
(8,3)
6 years, 11 months
