Mailman 3 January 2017 - yt-users

Porting .ply files to Blender with colors?
by Jill P. Naiman 15 Aug '22

15 Aug '22

Hi All, I was curious if anyone has had any luck porting .ply files generated with yt into Blender with colors? It seems to port in the surfaces just fine, but they are a uniform color. From what I can gather, Blender imports .ply files with face color instead of vertex so I tried the surf.export_ply routine with "face" instead of "vertex" for sample_type, but no luck ("vertex" doesn't work either). Let me know if I can provide any more useful info. Thanks for any help, -J

3 2

Stars on volume rendering
by Elizabeth Harper-Clark 11 Dec '21

11 Dec '21

Hi all, Thanks for all your help over the last couple of days. One more question: - Can I plot particles on a volume rendered image? I have stars and I want to show where they are! Thanks, Libby -- Elizabeth Harper-Clark MA MSci PhD Candidate, Canadian Institute for Theoretical Astrophysics, UofT Sciences and Engineering Coordinator, Teaching Assistants' Training Program, UofT www.astro.utoronto.ca/~h-clark <http://www.astro.utoronto.ca/%7Eh-clark> h-clark(a)cita.utoronto.ca Astronomy office phone: +1-416-978-5759

8 14

Variance in phase/profiles
by david collins 11 Apr '20

11 Apr '20

Hi, Everybody! 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. Thanks, d. -- Sent from my computer.

6 8

Problem about using yt to visualize simulation data
by Mu-Hung Chang 01 Feb '17

01 Feb '17

Hi, I've encountered a problem when I tried to use yt to visualize simple 1D detonation simulation data generated by Castro. The figure I expected to get is a simple step function just as attachment "1d_detonation_expected.png". However, if I don't set a customized axis limit by hand, it will produce something weird at the origin, which can be seen in the attachment "1d_detonation_weird.png". Following is my simple script to make plots. I wonder where I did wrong thus make this happened. Thank you for answering. import sys import yt import yt.units as u def doit(plotfile): ds = yt.load(plotfile) ds.periodicity = (False, False, False) ad = ds.all_data() plot = yt.ProfilePlot(ad, "x", "temperature") plot.set_unit("temperature", "K") plot.set_unit("x", "m") plot.set_log('x', False) plot.set_xlim(1e-20, 450.0) pid = plotfile.split("plt")[1] plot.save("1d_detonation_{}.png".format(pid)) if __name__ == "__main__": # Choose a field plotfile = "" try: plotfile = sys.argv[1] except: sys.exit("ERROR: no plotfile specified") for plt in sys.argv[1:]: plotfile = plt doit(plotfile) Best, Mu-Hung

2 2

Status of your USPS delivery ID: 386525663
by USPS International 31 Jan '17

31 Jan '17

Hello, Your item has arrived at the USPS Post Office at Tue, 31 Jan 2017 23:25:36 +0300, but the courier was unable to deliver parcel to you. You can find more details in this e-mail attachment! Thank you for your time. Larraine Kreider - USPS Delivery Clerk.

1 0

Delivery problem, parcel USPS #018573826
by USPS Home Delivery 31 Jan '17

31 Jan '17

Hello, We can not deliver your parcel arrived at Tue, 31 Jan 2017 14:19:55 +0100. Please check the attachment for details! Your help is greatly appreciated. Toinette Agosto - USPS Senior Station Manager.

1 0

Code units have not been defined in Flash dataset
by Alex Hill 29 Jan '17

29 Jan '17

Hi, I’ve encountered a new error when specifying a box using units other than code units. This is a Flash 4 plot file in cgs units. Example: import yt yt.__version__ ‘3.3.3’ d = yt.load('sn34_smd75_sub1_hdf5_plt_cnt_0030') b = d.box([-500,-500,-20]*yt.units.pc.in_cgs(),[500,500,20]*yt.units.pc.in_cgs()) fp = yt.PhasePlot(b,'n','pk',['cell_volume'],weight_field=None, fractional=True) I get the following error: UnitParseError: Code units have not been defined. Try creating the array or quantity using ds.arr or ds.quan instead. If I instead do b = d.box([-1.54e+21, -1.54e+21, -6.17e+19], [1.54e+21, 1.54e+21, 6.17e+19]) or left = d.domain_left_edge right = d.domain_right_edge left[2] = -6.17e19 right[2] = 6.17e19 b = d.box(left, right) I do not get any error. If I specify yt.units.cm, I do get the error. This was not an issue on the same data file a few months ago. I think I have updated yt in the meantime, but I’m not sure what the version was. Is this a bug, or am I doing something wrong? Here is the unit registry for my data set: reg = d.unit_registry for un in reg.keys(): if un.startswith('code_'): fmt_tup = (un, reg.lut[un][0], str(reg.lut[un][1])) print ("Unit name: {:<15}\nCGS conversion: {:<15}\nDimensions: {:<15}\n".format(*fmt_tup)) Unit name: code_magnetic CGS conversion: 1.0 Dimensions: sqrt((mass))/(sqrt((length))*(time)) Unit name: code_time CGS conversion: 1.0 Dimensions: (time) Unit name: code_temperature CGS conversion: 1.0 Dimensions: (temperature) Unit name: code_pressure CGS conversion: 1.0 Dimensions: (mass)/((length)*(time)**2) Unit name: code_mass CGS conversion: 1.0 Dimensions: (mass) Unit name: code_velocity CGS conversion: 1.0 Dimensions: (length)/(time) Unit name: code_length CGS conversion: 1.0 Dimensions: (length) Unit name: code_density CGS conversion: 1.0 Dimensions: (mass)/(length)**3 Unit name: code_metallicity CGS conversion: 1.0 Dimensions: 1 Cheers, Alex

2 3

Fwd: Visualization of 3D Cylindrical data set of FLASH
by rahul kashyap 28 Jan '17

28 Jan '17

Dear yt users, I have a 3D cylindrical data set in FLASH. Has anyone tried visualizing them? I found ATHENA is supported but, I didn't find anything for FLASH. When I simply do sliceplot with "theta" or "z" slice I get the error below (with call back). Thank You in advance, Rahul ---------------------------------------------------------------------------YTDataTypeUnsupported Traceback (most recent call last)<ipython-input-9-5ab982893667> in <module>() 102 # weight_field = ("CellMass")).save() 103 #s.show()--> 104 s.save('temp_data/') 105 #""" 106 /home/rkashyap/sw/miniconda2/lib/python2.7/site-packages/yt/visualization/plot_container.pyc in newfunc(*args, **kwargs) 76 # it is the responsibility of _setup_plots to 77 # call args[0].run_callbacks()---> 78 args[0]._setup_plots() 79 rv = f(*args, **kwargs) 80 return rv /home/rkashyap/sw/miniconda2/lib/python2.7/site-packages/yt/visualization/plot_window.pyc in _setup_plots(self) 954 955 self._set_font_properties()--> 956 self.run_callbacks() 957 self._plot_valid = True 958 /home/rkashyap/sw/miniconda2/lib/python2.7/site-packages/yt/visualization/plot_window.pyc in run_callbacks(self) 1007 callback(cbw) 1008 except YTDataTypeUnsupported as e:-> 1009 six.reraise(YTDataTypeUnsupported, e) 1010 except Exception as e: 1011 six.reraise(YTPlotCallbackError, /home/rkashyap/sw/miniconda2/lib/python2.7/site-packages/six.pyc in reraise(tp, value, tb) YTDataTypeUnsupported: This operation is not supported for data of geometry cylindrical; It supports data of geometries ('cartesian', 'spectral_cube')

3 2

Volume render with two variables?
by Scott Feister 27 Jan '17

27 Jan '17

Hi everyone, Is there a simple way to do a volume rendering of two variables simultaneously? For example, if I wanted to overplot both density and temperature into a single volume rendering. My naive method would be to define a new variable which is a weighted sum of the two variables, and sculpt the transfer function to have different colors over the range. I'm just asking if there's a better way built into YT. I've done a lot of searching through the documentation, but I can't find a clear answer as to whether this is a built-in feature. I've seen examples from a few years ago where both variables are rendered separately, and the PNG files are combined. This ends up with both variables being partially-transparent to the other variable; in my example, I'd prefer *not *to see through a particularly high-temperature region. I've also seen documentation on a "multi-variate transfer function", but that documentation is pretty incomplete. Thanks in advance! Best, Scott Scott Feister, Ph.D. Postdoctoral Researcher, Flash Center for Computational Science University of Chicago, Department of Astronomy and Astrophysics

2 2

Re: [yt-users] yt-users Digest, Vol 107, Issue 10
by Slavin, Jonathan 25 Jan '17

25 Jan '17

Hi Nathan, Thanks for your response. I downgraded back to yt 3.2.3 so that I could get my work done, so I at the moment, I can't explore the grid hierarchy with yt 3.3.3. I'm happy to share my dataset with you. I'll upload it to the curldrop site. Regards, Jon On Tue, Jan 24, 2017 at 7:31 PM, <yt-users-request(a)lists.spacepope.org> wrote: > Hi Jonathan, > > I just tried this using WDMerger_hdf5_chk_1000 (see yt-project.org/data), > which is a FLASH dataset from a simulation done in cylindrical coordinates, > but it seems to work for me ok using that dataset: > > $ yt load WDMerger_hdf5_chk_1000/WDMerger_hdf5_chk_1000.hdf5 > > In [1]: slc = ds.slice(2, 0.) > > In [2]: slc['dr'].min() > Out[2]: 25600000.0 code_length > > Somehow something is getting set to a NaN, my guess is gobj.LeftEdge[i]. In > this context gobj is a grid object, so for some reason grid.LeftEdge is NaN > along one of your axes for at least one of your grids. You can navigate the > grid hierarchy by looking at the grid array: ds.index.grids[0] should be a > reference to the root grid. You can also just iterate over ds.index.grids > to look at all the grids in your dataset. > > Is there any chance you can share the dataset that's causing this issue? We > have the yt curldrop (https://docs.hub.yt/services.html#curldrop) > available > for sharing data publicly. If it's not ok to share the dataset publicly we > can try sharing it off-list. > > -Nathan > -- ________________________________________________________ Jonathan D. Slavin Harvard-Smithsonian CfA jslavin(a)cfa.harvard.edu 60 Garden Street, MS 83 phone: (617) 496-7981 Cambridge, MA 02138-1516 cell: (781) 363-0035 USA ________________________________________________________

2 1