Hi all,
I feel like this should be relatively straightforward, but I have a lot of
rather small clouds that I have identified already. I can get info about
them by using their position and then creating a sphere around it.
I already have the maximum density and minimum temperature values for these
clouds. All I want to do is see if these are the same cell, and if not,
what the difference in temp (or density) is. So I know how to find a
maximum quantity, but not how to then find values of other variables in
that cell. I was thinking that I basically want an array with max density,
min temperature, density at min temperature, temperature at max density.
So can anyone point me to how to find the temperature of the cell with the
highest density in a sphere object?
Thanks!
Stephanie
--
--
Dr. Stephanie Tonnesen
Associate Research Scientist
CCA, Flatiron Institute
New York, NY
stonnes(a)gmail.com
Dear all,
Registration is open for the upcoming "yAC: yt at CCA" workshop, to take
place at the Flatiron Institute from Monday March 23rd through Wednesday
March 25th:
https://forms.gle/56bmuYKcX5NCDF5aA
This workshop will introduce the yt visualization and analysis tool (
https://yt-project.org/), provide tutorials and example workflows, and time
to interact with the developers. Additionally, there will be ample free
time for practice and hacking in the afternoons.
Participation is limited to ~40, so please register early if you are
interested. More details and a schedule will come in a few weeks.
local organizers:
Stephanie Tonnesen
Michael Zingale
--
Michael Zingale
Associate Professor
Dept. of Physics & Astronomy • Stony Brook University • Stony Brook, NY
11794-3800
*phone*: 631-632-8225
*e-mail*: Michael.Zingale(a)stonybrook.edu
*web*: http://www.astro.sunysb.edu/mzingale
github: http://github.com/zingale
Dear Yt users,
I am trying to install YT with the attached installing script. I got the
following error message.
*Error::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::*
Installing yt
Testing that yt can be imported
install_script.sh: line 414: 45562 Illegal instruction (
${DEST_DIR}/bin/${PYTHON_EXEC} -c "import yt" 2>&1 ) >> ${LOG_FILE}
********************************************
FAILURE REPORT:
********************************************
yt conda-forge/linux-64::yt-3.5.1-py37hc1659b7_0
Downloading and Extracting Packages
yt-3.5.1 | 9.4 MB | ########## | 100%
Preparing transaction: ...working... done
Verifying transaction: ...working... done
Executing transaction: ...working... done
********************************************
********************************************
Failure. Check /export/home/minato/yt2020/yt-conda/yt_install.log. The
last 10 lines are above.
*:::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::*
I also attached the log file for your reference. Please help me.
Many thanks,
-Prateek
Greetings yt community,
I am a graduate student in Cornell's Department of Computer Science seeking suggestions for semester long software development projects. I wonder, is there anyone in the yt community with a yt development project in mind that might benefit from a dedicated team of (aspiring) scientific software developers? The guidelines of the project require that we have a "client" to interface with on a weekly basis. To learn more about the course and the projects, see https://www.cs.cornell.edu/courses/cs5150/2020sp/ and follow the Projects links.
Note especially that the projects are production systems, not prototypes or (purely) academic research. A project can be a small, complete system, or part of a larger system; students undertake projects in teams of 6 to 8. These teams develop real software applications for clients who intend to place them into production. At the end of the course, students are expected to hand over operational software, which has been thoroughly tested, with good documentation for future maintenance.
Thanks for considering.
Sincerely,
Ariel Kellison
Hi all,
I'm trying to make slice plots for 3D chombo data sets created with pluto version 4.3, but I'm consistently (not just this dataset) running into the issue where some boxes are empty. I've also made a slice plot using VisIt where these empty boxes do not show up, so it's probably an issue in yt.
Is pluto v4.3 supported? Is there maybe an easy work around?
yt-version: 4.0.dev0
The issue was also present in the latest stable version.
The code to produce it is simple:
fn = "data.0125.hdf5"
var = 'density'
fig = plt.figure()
ds = yt.load(fn)
p = yt.SlicePlot(ds, 'z', var)
p.zoom(2)
p.pan_rel((0.,-0.1))
p.set_zlim(var, 1, 2)
p.set_cmap(field=var, cmap='rainbow')
p.set_width((1., 1.5))
p.set_log(var, True)
p.show()
An image of the problem (red squares):
https://i.imgur.com/llVsvn0.png
data file:
http://use.yt/upload/4a0db3c5
pluto.ini file:
http://use.yt/upload/781478c2
Cheers,
Richard
I would like to annotate Mass weighted Magnetic field quivers on density projection plot so that i can display the magnetic field features.
Do annotate_magnetic_field() , annotate_quiver() work correctly on projection plot or do they work only on slice plot correctly.
If they work on projection plot then against which quantity magnetic field quivers are weighed or not weighted?
Hello All,
I have a question about the Projection plot function in yt. I have a script that plots the smallest mass darkmatter particles of a simulation in all 3 of the 2-d coordinate projections, and overplots the halos found by a halo finder. This script works, and currently plots the colorbar based on the mass in a given 2-d bin as done in some of the basic yt plots included in the website. I am wondering if it is possible to have the color scheme of the particles be based on the squared mass or squared density of the particles in any given 2-d bin, so that the higher density regions are highlighted compared to the background lower density regions. I have included the plotting part of my script below, in case it is helpful.
I hope everyone is having a great New Year!
ds = yt.load(VELA_snaps[position[0]])
domain_width = float(ds.domain_width.in_units('Mpc/h')[0])
ad = ds.all_data()
masses = yt.np.unique(ad[('darkmatter', 'particle_mass')])
#filter out the darkmatter0 particles
def mass_filter(pfilter, data):
filter = data[(pfilter.filtered_type, 'particle_mass')] == masses[0]
return filter
yt.add_particle_filter('darkmatter0', function=mass_filter, filtered_type='darkmatter', requires=['particle_mass'])
ds.add_particle_filter('darkmatter0')
scale = float(ds.scale_factor)
fig = plt.figure()
grid = AxesGrid(fig, (0.075,0.075,10,5),
nrows_ncols = (3, 1),
axes_pad = 1.0,
label_mode = "L",
share_all = False,
cbar_location="right",
cbar_mode="each",
cbar_size="3%",
cbar_pad="0%")
zoom = 10
x0 = float(consistent.halo_data_sorted[index][0][17]) * scale / domain_width
y0 = float(consistent.halo_data_sorted[index][0][18]) * scale / domain_width
z0 = float(consistent.halo_data_sorted[index][0][19]) * scale / domain_width
center = [x0, y0, z0]
a = yt.ParticlePlot(ds, ('darkmatter0', 'particle_position_x'), ('darkmatter0', 'particle_position_y'),\
('darkmatter0', 'particle_mass'), center=center)
a.set_unit(('darkmatter0','particle_mass'), 'Msun')
a.zoom(zoom)
b = yt.ParticlePlot(ds, ('darkmatter0', 'particle_position_y'), ('darkmatter0', 'particle_position_z'),\
('darkmatter0', 'particle_mass'), center=center)
b.set_unit(('darkmatter0','particle_mass'), 'Msun')
b.zoom(zoom)
c = yt.ParticlePlot(ds, ('darkmatter0', 'particle_position_z'), ('darkmatter0', 'particle_position_x'),\
('darkmatter0', 'particle_mass'), center=center)
c.set_unit(('darkmatter0','particle_mass'), 'Msun')
c.zoom(zoom)
for halos in consistent.halo_data_sorted[index]:
x = float(halos[17]) * scale / domain_width
y = float(halos[18]) * scale / domain_width
z = float(halos[19]) * scale / domain_width
r = float(halos[11]) * scale / .7
center = [x, y, z]
a.annotate_sphere(center, radius=(r, 'kpc'), circle_args={'color':'red'})
b.annotate_sphere(center, radius=(r, 'kpc'), circle_args={'color':'red'})
c.annotate_sphere(center, radius=(r, 'kpc'), circle_args={'color':'red'})
index = 0
for letter in [a,b,c]:
plot = letter.plots[('darkmatter0', 'particle_mass')]
plot.figure = fig
plot.axes = grid[index].axes
plot.cax = grid.cbar_axes[index]
letter._setup_plots()
index = index + 1
plt.savefig('{}/CatalogProjectionYT_VELA{}_Scale{}.png'.format(out_dir, VELA_number, str(scale)[2:5]), bbox_inches='tight')
plt.close()
Dear all,
We are pleased to announce the upcoming "yAC: yt at CCA" workshop, to take
place at the Flatiron Institute in NY from Monday March 23rd through
Wednesday March 25th, focusing on the yt visualization and analysis tool (
https://yt-project.org/). The format will be tutorials in the morning and
early afternoon + practice / hacking in the afternoons. A website with
registration links will be announced early in the new year, but we wanted
to share the preliminary announcement and date now. We hope you can join
us.
local organizers:
Stephanie Tonnesen
Michael Zingale
feel free to forward!
--
Michael Zingale
Associate Professor
Dept. of Physics & Astronomy • Stony Brook University • Stony Brook, NY
11794-3800
*phone*: 631-632-8225
*e-mail*: Michael.Zingale(a)stonybrook.edu
*web*: http://www.astro.sunysb.edu/mzingale
github: http://github.com/zingale
Hi guys,
I am very new to Ytini's mailing list, mailing lists as well as Houdini.
So please bear with me while I pick up the etiquette around this medium.
With that out of the way here's my first question:
After having succeeded installing Ytini for Houdini on a Mac OS X 10.14 I
would now like to achieve the same with my Houdini install on a Windows 10
OS.
I have tried to somewhat reproduce all steps described here
http://www.ytini.com/getstarted.html#collapseTwoOut (Houdini 17.5 Tested on
Mac OS X) on my Windows machine but wasn't surprised to learn that of
course Houdini could not find the ytini module when I tried to import it in
the python shell.
Is there somebody, that could point me . in the right direction of how to
install Ytini for Houdini on a Windows 10 OS environment?
Not afraid of reading so even hints a re more than welcome.
Many Thanks in advance
Ur4n
