Mailman 3 - yt-users

Precise Data
by Daniel Ruvalcaba April 4, 2024

April 4, 2024

Hello I am a new user to yt and am also new to python as a coding language. I do have past coding experiences but that is not important. I have been making graphs based on the files given on yt_data and am now trying to see if I can grab specific data from points within the graph. I saw that yt can do that with a sphere and its radius but I would like to see if there was a way to do the same with projection plots (other graphs as well). Specifically what I want to ask is if I can grab data from … [View More]

3 3

Trouble adding a field to TNG data
by Jacob Morgan April 4, 2024

April 4, 2024

Hey there! I have a relatively extensive plugins file that adds lots of fields. I changed how one of the fields is calculated, but now I am getting an error about an undefined variable. I have the following function: def _get_scalar_func(_ptype, field_name, unit_dict): def _scalar(_ptype, field, data): if field_name in only_gas_fields: _ptype='gas' return data[_ptype, field_name] elif 'velocity_los' in field_name and data.ds.head('RAY')==0: pass else: … [View More]old_data=data[_ptype, field_name] x_old_units=unyt.unyt_array(old_data.value*unit_dict[field_name][0], units=unit_dict[field_name][1], registry=data.ds.unit_registry) return x_old_units.to(unit_dict[field_name][2]) return _scalar it converts some scalar fields for me based on "unit_dict", which holds what units I want for each property. It also adds the field ('PartType0', 'entropy'), which is just a copy of ('gas', 'entropy') with the first "if" statement. However, it says that the variable "_ptype" is referenced before assignment, and none of my plugins stuff is added. When I try to add a field manually outside the plugins file like so: ds.add_field(('PartType0', 'Density_msun_pc3'), function=_get_scalar_func('PartType0', 'Density', gas_scalar_conversions), sampling_type='particle', units='Msun/pc**3') I get the following long error: Traceback (most recent call last): File "file_loader_tester.py", line 35, in <module> new_yt.add_field(('PartType0', 'Density_msun_pc3'), function=_get_scalar_func('PartType0', 'Density', gas_scalar_conversions), sampling_type='particle', units='Msun/pc**3') File "/home/jacob/.local/lib/python3.8/site-packages/yt/data_objects/static_output.py", line 1593, in add_field self.index File "/home/jacob/.local/lib/python3.8/site-packages/yt/data_objects/static_output.py", line 539, in index self.create_field_info() File "/home/jacob/.local/lib/python3.8/site-packages/yt/data_objects/static_output.py", line 613, in create_field_info self.field_info.load_all_plugins(self.default_fluid_type) File "/home/jacob/.local/lib/python3.8/site-packages/yt/fields/field_info_container.py", line 402, in load_all_plugins self.find_dependencies(loaded) File "/home/jacob/.local/lib/python3.8/site-packages/yt/fields/field_info_container.py", line 415, in find_dependencies deps, unavailable = self.check_derived_fields(loaded) File "/home/jacob/.local/lib/python3.8/site-packages/yt/fields/field_info_container.py", line 612, in check_derived_fields fd = fi.get_dependencies(ds=self.ds) File "/home/jacob/.local/lib/python3.8/site-packages/yt/fields/derived_field.py", line 253, in get_dependencies e[self.name] File "/home/jacob/.local/lib/python3.8/site-packages/yt/fields/field_detector.py", line 125, in __missing__ vv = finfo(self) File "/home/jacob/.local/lib/python3.8/site-packages/yt/fields/derived_field.py", line 294, in __call__ dd = self._function(self, data) File "/home/jacob/.config/yt/my_plugins.py", line 69, in _scalar old_data=data[_ptype, field_name] UnboundLocalError: local variable '_ptype' referenced before assignment This has been working for literally years, and I only made minor changes to the inner _scalar function recently. But this almost seems like an issue with the add_field() function, which I don't think I have ever touched. Any ideas? [View Less]

1 1

how to load and plot Particle data output by AMReX
by hyf＠stu.pku.edu.cn April 1, 2024

April 1, 2024

Hi everyone, I have output particle data by the function 'pc.Checkpoint' of AMREX after 'WriteMultiLevelPlotfile' you can see details here: https://amrex-codes.github.io/amrex/docs_html/Particle.html#particle-io however, I cannot find the documention about how to load and handle with the output particle data. could you give me some references or suggestions, thank a lot! best wishes

3 9

how to load and plot Particle data output by AMReX
by hyf＠stu.pku.edu.cn March 27, 2024

March 27, 2024

Hi everyone, I have output particle data by the function 'pc.Checkpoint' of AMREX after 'WriteMultiLevelPlotfile' you can see details here: https://amrex-codes.github.io/amrex/docs_html/Particle.html#particle-io however, I cannot find the documention about how to load and handle with the output particle data. could you give me some references or suggestions, thank a lot! best wishes

1 0

outputting unigrid Enzo data cubes with yt
by Daniel Whalen March 16, 2024

March 16, 2024

Hi, everyone: Quick question: Muhammad Latif and I would like to output a unigrid data cube of a Pop III supernova remnant in an Enzo AMR simulation in either HDF5 or ASCII format for post processing with a radiative transfer code. We need to include all primordial, metal and molecular mass fractions in addition to densities and energies. We went through the yt docs but didn't immediately see a way to do this, does anyone have any suggestions? Thanks! Dan Whalen

2 2

Off-Axis SlicePlot Error with TNG data in yt
by 缘边 March 15, 2024

March 15, 2024

Hi I've recently started using yt for visualizing the gas density of TNG galaxies with an OffAxisSlicePlot. I want the top of the image to align with the galaxy's velocity direction [0.9196489, 0.15788311, -0.35960928], so I set north_vector = [0.9196489, 0.15788311, -0.35960928]. The projection direction is the normal vector perpendicular to the galaxy's velocity direction. Here's my code snippet: c = ds.arr([Subhalo["SubhaloPos"][0], Subhalo["SubhaloPos"][1], Subhalo["SubhaloPos"][2]], "… [View More]

2 1

Spherical harmonic decomposition of sphere surface
by Samuel Patrick March 15, 2024

March 15, 2024

Hi all, I would like to perform a spherical harmonic decomposition of the density within a spherical shell, essentially a Fourier transform of the density (or any scalar quantity) on the surface of a sphere object . I'm aware of python libraries like Healpy which have functions such as map2alm which take a healpix map and return the corresponding Alm coefficients; are there any utilities within yt that can either do this directly or return a healpix map from a sphere object? Thanks and best, … [View More]

2 1

Change the center position while using plot_2d function.
by vashishthanitin5＠gmail.com March 4, 2024

March 4, 2024

Hi yt Team, I am not able to change the centre position of the selection while plotting the density from a simulation output file using yt.plot_2d function. It is taking the domain centre as the centre value, even if I change it. I am also not able to change the origin of the output plot as well. Following is the code, please look into it. The output file is from amrvac code. It is the output of a 2.5 MHD code in spherical geometry. My domain centre values are [5.5 1.5707963267948966 0.5] … [View More]

2 1

FLIP IN A MULTI PLOT
by FERNANDO ILLACANCHI GUERRA Feb. 27, 2024

Feb. 27, 2024

Hi Everyone, I have a simulation with an axis symmetry so i want to plor full simulation using the flip function. However, both plots are flipped in te resulting imagen. Is there a way to flip just one of them? Thanks in advance. for i in range(0,2): print(i) # Load the data and create a single plot ds = yt.load(fn) # load data p = yt.plot_2d(ds, ("Temp")) p.set_width((4,1)) if i == 1: p.flip_vertical() print("vertical") # This forces the … [View More]

2 1

yt notebook
by Daniel Ruvalcaba Feb. 27, 2024

Feb. 27, 2024

Hello I am currently trying to install and setup yt on my laptop using the miniconda/anaconda/pip/whatever will work and I keep running into the issue of not being able to open my yt notebook. For pip/git it does not accept the command and for miniconda/anaconda it keeps giving me this error: Traceback (most recent call last): File "C:\ProgramData\miniconda3\Scripts\yt-script.py", line 10, in <module> sys.exit(run_main()) ^^^^^^^^^^ File "C:\Users\.\AppData\Roaming\… [View More]

2 1