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?
Hi all, I'm trying to use "gas", "radial velocity" to compare inflows and outflows in a simulation, but the results don't really look right.
using radial_velocity as a binary field, I defined new fields as follows
def inflow_density(field,data):
tr = data['gas','density']*((data['gas','radial_velocity'])<0).astype(float)
return tr
def outflow_density(field,data):
tr = data['gas','density']*((data['gas','radial_velocity'])>0).astype(float)
return tr
A ProjectionPlot gives one side which is all inflow, one side which is all outflow. This resembles what you get if you have an incorrect bulk_velocity field, but when I correct for it like this
sp = ds.sphere(center, (Rvir, "kpc"))
bulk_velocity = sp.quantities.bulk_velocity()
plot = yt.ProjectionPlot(ds, 0, center=center, fields=('gas','inflow_density'), width = (w,"kpc"), field_parameters={'bulk_velocity':bulk_velocity})
it doesn't change at all. Am I using this right?
Thanks,
Clayton
Hi yt users!
I'm not sure if this is the place to ask this but I was following
The example to export to RADMC3D
https://yt-astro-analysis.readthedocs.io/en/latest/radmc3d_export.html
Number of processors: 4
Number of threads in use: 1
Reading global frequencies/wavelengths...
Reading grid file and prepare grid tree...
At line 3196 of file ./amr_module.f90 (unit = 1, file = 'amr_grid.inp')
Fortran runtime error: Bad integer for item 5 in list input
I've run examples from radmc3dpy webpage and those worked. I'm not sure if
the problem is the format in which yt writes the amr input file or radmc3d
Saludos!
José Utreras
Greetings!
I'm very pleased to announce the release of Trident version 1.2!
Trident is an extension of yt for creating absorption spectra from
simulations.
This is the first major Trident release in nearly two years. Trident
1.2 contains a number of new features, bug fixes, and minor
enhancements. All are listed below.
*This is the last version of Trident that will support Python 2.*
All versions to follow will require the use of Python 3.5 or later.
*New Features*
- Add support for creating spectra in velocity space
- Add ability to set wavelength limits to 'auto'
- Add ability to make spectrum from data container
*Bug Fixes*- off by one error in subgrid index calculation
- fixing error in _ion_mass
- H_p0_number_density to default, H_number_density to alias
- Implements a better way of calculating the redshift on the light ray
- Assures that ion_fraction field reflects on-disk fields
- Fix atomic data for Si II 1260 according to Morton 2003
- A check to avoid failure when no continuum absorbers were found in the ray
- Auto-addition of H_nuclei_density to a LightRay when present in the base
dataset
- Adding max_box_fraction as a kwarg to make_compound_ray
- Updated trident_path() to be OS Independent
- simplify setting up ion fields using the "local" field type
- split and join filenames using os.sep instead of assuming unix
- updated oscillator strengths and gamma's for Si II 1206 and Si III 1260
*Minor Enhancements*- Calculating LOS velocity with relative velocities to
account for bulk motion
- Enabling use of `output_absorbers_file` kwarg in SpectrumGenerator
- Switching imports from yt.analysis_modules to yt_astro_analysis
- Enable passing in-memory LineDatabase to SpectrumGenerator
- Added equivalent width calculation to line_observables_dict
- Numerous documentation updates
- Updates and fixes to testing
The latest version of Trident can be installed with pip:
pip install trident
Documentation: https://trident.readthedocs.io/en/stable/
Repository: https://github.com/trident-project/trident
Trident has received contributions from 23 people. 9 people
contributed changes to the latest release (underline) with 3 of those
people doing so for the first time (bold). Thanks to everyone for
their hard work!
- Trident development team
*Brendan Boyd*
Lauren Corlies
*Bili Dong*Hilary Egan
*Nathan GoldbaumCameron HummelsBharat Kambalur*BW Keller
Kacper Kowalik
Erwin Lau
Andrew Myers
*Molly S. Peeples*Douglas Rudd
Anthony Scopatz
*Devin Silvia*Stephen Skory
*Britton SmithClayton Strawn*Stephanie Tonnesen
Matthew Turk
Miguel de Val-Borro
John Wise
John ZuHone
Hi, I would like to use yt to make selections of data from Enzo files (e.g. gas with certain properties in a particular volume), and then output the data to an external file like an ascii or fits table, to be analysed by other codes. Can you please tell me if such a tool is already available, or how I might be able to do it otherwise? Thanks.
Hello all,
I want to use a user-defined function (say, Maxwell distribution function)
with a variable (say, p--momentum, because I want to integrate the function
on that variable), and want to use density and temperature field
information of each grid. This function applies individually on each grid.
Please suggest me the way to do so. This could be very naive to ask but I
stuck in it. That function can't be a yt-field (I think) because it has a
variable, and if I want to use density and temperature field even after
removing units inside the function, it will show an error, because the
field is in vector form (3Dim) and the variable is scaler (basically shape
miss-match error).
I have to use many such functions as post-processing, please suggest me a
way to achieve this goal.
Many thanks and regards,
-Prateek Gupta
Hello!
I'm having issues running a yt 3.5.1 script to process an
OffAxisProjectionPlot of some Enzo 2.5 data. Running mpi 4.0.1 I get this
error:
File "plot_time.py", line 28, in <module>
val, loc = ds.find_max('density')
File "/opt/apps/pkgs/yt/3.5.1/src/yt/data_objects/static_output.py", line
826, in find_max
source.quantities.max_location(field)
File "/opt/apps/pkgs/yt/3.5.1/src/yt/data_objects/derived_quantities.py",
line 635, in __call__
rv = super(MaxLocation, self).__call__(field, sample_fields)
File "/opt/apps/pkgs/yt/3.5.1/src/yt/data_objects/derived_quantities.py",
line 591, in __call__
rv = super(SampleAtMaxFieldValues, self).__call__(field, sample_fields)
File "/opt/apps/pkgs/yt/3.5.1/src/yt/data_objects/derived_quantities.py",
line 68, in __call__
for sto, ds in parallel_objects(chunks, -1, storage = storage):
File
"/opt/apps/pkgs/yt/3.5.1/src/yt/utilities/parallel_tools/parallel_analysis_interface.py",
line 520, in parallel_objects
to_share, datatype = 'dict', op = 'join')
File
"/opt/apps/pkgs/yt/3.5.1/src/yt/utilities/parallel_tools/parallel_analysis_interface.py",
line 284, in passage
return func(self, *args, **kwargs)
File
"/opt/apps/pkgs/yt/3.5.1/src/yt/utilities/parallel_tools/parallel_analysis_interface.py",
line 723, in par_combine_object
data.update(self.comm.recv(source=i, tag=0))
File "mpi4py/MPI/Comm.pyx", line 1173, in mpi4py.MPI.Comm.recv
File "mpi4py/MPI/msgpickle.pxi", line 302, in mpi4py.MPI.PyMPI_recv
File "mpi4py/MPI/msgpickle.pxi", line 261, in mpi4py.MPI.PyMPI_recv_match
P000 yt : [ERROR ] 2019-09-14 12:26:53,382 Exception: MPI_ERR_INTERN:
internal error
P000 yt : [ERROR ] 2019-09-14 12:26:53,383 Error occured on rank 0.
I know that someone posted a similar error recently but I don't think there
was a solid solution to the problem. If anyone can help please let me know.
Also i'm sorry if I didn't explain this well enough, please let me know if
anything is not clear.
Many thanks
Andy
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
Hello,
I want to compute the magnetic flux over a disk which length is the minimum possible i.e. the resolution. I wanted to approximate this quantity using :
sp=ds.disk(center,normal,(1,'pc'),(1*res,'pc'))
flux=sp.calculate_flux("sqrt_magx","sqrt_magy","sqrt_magz","sqrt_magx")
Where I defined the quantities sqrt_magx (root square of the magnetic field along the x direction) with units of Gauss^0.5, so the output of calculate_flux has Gauss*cm^2 units. However, it is not possible to use calculate_flux on YTDisks and only calculate_isocontour_flux is available for this data container.
Do you know how I can obtain exactly the same result for the magnetic flux using the available calculate_isocontour_flux instead of the calculate_flux function?
Cheers,
Hi all,
I've been encountering an error recently whenever I try to save a slice
plot, such as:
data_sets = yt.load("DD*/output_[0-9][0-9][0-9][0-9]")
for ds in data_sets:
filename =
"density_t={:.3f}_ho.png".format(ds.current_time.in_units('Myr').value)
val, loc = ds.find_max('density')
print("{:.3f}".format(val))
plt = yt.SlicePlot(ds, L, field, center= loc, width= (size, 'pc'))
plt.save(filename)
YT seems to think issue is with the call to 'find_max' but a little
experimentation shows this is not the case. I'm unable to work out any
of what comes back to me asides that:
File "plot.py", line 67, in <module>
val, loc = ds.find_max('density')
File "/opt/apps/pkgs/yt/3.5.1/src/yt/data_objects/static_output.py",
line 826, in find_max
source.quantities.max_location(field)
File
"/opt/apps/pkgs/yt/3.5.1/src/yt/data_objects/derived_quantities.py",
line 635, in __call__
rv = super(MaxLocation, self).__call__(field, sample_fields)
File
"/opt/apps/pkgs/yt/3.5.1/src/yt/data_objects/derived_quantities.py",
line 591, in __call__
rv = super(SampleAtMaxFieldValues, self).__call__(field, sample_fields)
File
"/opt/apps/pkgs/yt/3.5.1/src/yt/data_objects/derived_quantities.py",
line 68, in __call__
for sto, ds in parallel_objects(chunks, -1, storage = storage):
File
"/opt/apps/pkgs/yt/3.5.1/src/yt/utilities/parallel_tools/parallel_analysis_interface.py",
line 520, in parallel_objects
to_share, datatype = 'dict', op = 'join')
File
"/opt/apps/pkgs/yt/3.5.1/src/yt/utilities/parallel_tools/parallel_analysis_interface.py",
line 284, in passage
return func(self, *args, **kwargs)
File
"/opt/apps/pkgs/yt/3.5.1/src/yt/utilities/parallel_tools/parallel_analysis_interface.py",
line 723, in par_combine_object
data.update(self.comm.recv(source=i, tag=0))
File "mpi4py/MPI/Comm.pyx", line 1173, in mpi4py.MPI.Comm.recv
File "mpi4py/MPI/msgpickle.pxi", line 302, in mpi4py.MPI.PyMPI_recv
File "mpi4py/MPI/msgpickle.pxi", line 261, in mpi4py.MPI.PyMPI_recv_match
Many thanks,
Sam Patrick
