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.
Sent from my computer.
Can current yt calculate 3-D Mass power spectra? I checked the website but
I didn't find any information. I think calculating 3-D Mass power
spectra is a very useful for cosmological simulations. So I guess maybe yt
supports this function now....?
Thanks in advance
I am using arbitrary_grid to access my data without flattening:
> dd=ds.arbitrary_grid(left_edge=[(com_x_1.value-30), (com_y_1.value-30),
> (com_z_1.value-30)]*kpc, right_edge=[(com_x_1.value+30), (com_y_1.value+30),
> (com_z_1.value+30)]*kpc, dims=[128, 128, 128])
and as I attempt to get the gas fields for this object,
> print(dd['gas', 'velocity_x'])
but end up getting the error: (which isn't the case for deposit fields)
yt.units.unit_registry.UnitParseError: Code units have not been defined.
> Try creating the array or quantity using ds.arr or ds.quan instead.
If I'm thinking right, is it because the arbitrary_grid is only applied to
the particles. In that case, is there any equivalent for gas fields where I
can access my 3D data for the gas as it is , without flattening? Thanks
This is a new feature in yt 3.5.0, you’ll need to update to be able to use
it. Let us know if you have questions about updating.
On Tue, Oct 30, 2018 at 2:16 PM Samantha M Benincasa <
> Hi everyone,
> I'm trying to get a halo into a face-on and side-on configuration using
> Is there a way to do this calculation with only the young stars in the
> sphere? The docs suggest that angular_momentum_vector will take a
> "particle_type" keyword, but my version does not recognize this.
> yt-users mailing list -- yt-users(a)python.org
> To unsubscribe send an email to yt-users-leave(a)python.org
I'm trying to get a halo into a face-on and side-on configuration using sp.quantities.angular_momentum_vector.
Is there a way to do this calculation with only the young stars in the sphere? The docs suggest that angular_momentum_vector will take a "particle_type" keyword, but my version does not recognize this.
I would like to access data from Ramses as it is because I would want to
divide my selected region into subgrids and compute some quantity locally
in those subgrids. However, using all_data or ds.r would rearrange the
information as it flattens the data (please correct me if I got this
wrong). Is there a way in yt, where I can simply access the 3D data as it
is and divide into subcubes? Thank you.
We have an open postdoctoral research associate position in computational astrophysics/supernova theory at Michigan State. More info and application portal here: http://careers.msu.edu/cw/en-us/job/500160/research-associatefixed-term. Please pass along to anyone you think might be interested.
Apologies if you receive this announcement multiple times.
Sean M. Couch, Ph.D.
Department of Physics and Astronomy
Department of Computational Mathematics, Science, and Engineering
National Superconducting Cyclotron Laboratory/Facility for Rare Isotope Beams
Michigan State University
567 Wilson Rd, 3250 BPS
East Lansing, MI 48824
(517) 884-5035 —— couch(a)pa.msu.edu<mailto:email@example.com> —— www.pa.msu.edu/~couch<http://www.pa.msu.edu/~couch>
I am still a bit new to yt and am trying to compute local velocity
dispersions of the gas in my simulation by dividing the entire volume into
‘n’ number of sub grids. Does yt offer any functionality that could perform
this in a simple fashion? (I looked up the documentation and couldn’t find
any such function). Thank you.
I am creating a derived field for new_vorticity_magnitude. But when I take
the SlicePlot of this derived field I am getting the following error:
slc = yt.SlicePlot(ds, 'z', "Vorticity_mag", center= center_old , width =
File "/home/dumbo/yt/yt/visualization/plot_window.py", line 1953, in
return AxisAlignedSlicePlot(ds, normal, fields, *args, **kwargs)
File "/home/dumbo/yt/yt/visualization/plot_window.py", line 1266, in
File "/home/dumbo/yt/yt/data_objects/data_containers.py", line 1336, in
File "/home/dumbo/yt/yt/data_objects/data_containers.py", line 1356, in
fd = self._generate_field(field)
File "/home/dumbo/yt/yt/data_objects/data_containers.py", line 318, in
tr = self._generate_fluid_field(field)
File "/home/dumbo/yt/yt/data_objects/data_containers.py", line 338, in
rv = finfo(gen_obj)
File "/home/dumbo/yt/yt/fields/derived_field.py", line 237, in __call__
dd = self._function(self, data)
File "vorticity_test.py", line 63, in _newvorticity_x_
f = (data[ftype, "zVel"][sl_center,sl_right,sl_center] - data[ftype,
"zVel"][sl_center,sl_left,sl_center]) / (div_fac*just_one(data["index",
File "/home/dumbo/yt/yt/units/yt_array.py", line 1036, in __getitem__
ret = super(YTArray, self).__getitem__(item)
IndexError: too many indices for array
The current version of yt is:
Version = 3.5.dev0
Changeset = 0cb7ad5befb2
Please let me know, why I am getting this error?
I am trying to get multiple profile plots of density pdf for ramses data
and I encounter the following error:
> FutureWarning: Conversion of the second argument of issubdtype from `float`
> to `np.floating` is deprecated. In future, it will be treated as
> `np.float64 == np.dtype(float).type`.
> from ._conv import register_converters as _register_converters
> Traceback (most recent call last):
> File "pdf.py", line 8, in <module>
> sim = yt.simulation(filename, "Ramses")
> "/home/samvad/yt-conda/lib/python3.6/site-packages/yt/convenience.py", line
> 112, in simulation
> raise YTSimulationNotIdentified(simulation_type)
> yt.utilities.exceptions.YTSimulationNotIdentified: Simulation time-series
> type ramses not defined.
and the following is the code I am using to do so:
> import yt
> from yt.units import kpc
> import matplotlib.pyplot as plt
> import numpy as np
> import yt.visualization.eps_writer as eps
> sim = yt.simulation(filename, "Ramses")
> for ds in sim:
> ds.define_unit("hpc", (1.674*10**(-24), "g/cm**3"))
> dens=ad.cut_region(["obj['density'].in_units('hpc') > 0"])
> profiles.append(yt.create_profile(dens, "density", "cell_mass",
> weight_field=None, fractional=True))
> labels.append("time=%.2f" % ds.current_time.in_units('Myr'))
> plot_specs.append(dict(linewidth=2, alpha=0.7))
> plot = yt.ProfilePlot.from_profiles(profiles, labels=labels,
> plot.set_log("density", True)
> plot.set_unit("density", "hpc")
> plot.set_xlim(1e-7, 1e7)
could someone please let me know how I can get the simulation_type for