[yt-users] Re: Velocity ProjectionPlot with yt-4.0

Hi Daniel, What is your version of NumPy? There is an issue with NumPy 1.18 which causes this: https://github.com/yt-project/yt/issues/2469 https://github.com/yt-project/yt/issues/2469 This is fixed in yt 3.6 but not yet in 4.0. The temporary solution is to downgrade to NumPy 1.17. Also, the problem with TNG halos is that yt won’t recognize them as Arepo datasets. In order to make them work, you have to hack them like this: # Now we have to open up the file we downloaded and fix the header # in the file to include info yt needs to read this as an Arepo file f = h5py.File(saved_filename, "r+") f["Header"].attrs["NumPart_Total"] = np.array(f["Header"].attrs["NumPart_ThisFile"]) f.create_group("Config") f["/Config"].attrs["VORONOI"] = 1 f.flush() f.close() We should fix this in the future. Best, John

On Apr 24, 2020, at 10:11 AM, Daniel DeFelippis wrote:

Hello,

I'm trying to create a projection plot of velocity (e.g. velocity_x) for a halo from IllustrisTNG ('cutout_478037.hdf5' below). I'm using the demeshened version of yt (yt-4.0). When I try to do this it fails with the following error message. This only seems to occur when I try to make a projection of a quantity with positive and negative values. If I instead use 'density' or 'temperature' for example, it works fine. Any suggestions?

Thank you, Daniel

ds = yt.load('cutout_478037.hdf5') yt.ProjectionPlot(ds,'x',('gas','velocity_x'),width=(400,'kpc'))

--------------------------------------------------------------------------- TypeError Traceback (most recent call last) ~/anaconda3/lib/python3.6/site-packages/numpy/core/function_base.py in linspace(start, stop, num, endpoint, retstep, dtype, axis) 116 try: --> 117 num = operator.index(num) 118 except TypeError:

TypeError: 'numpy.float64' object cannot be interpreted as an integer

During handling of the above exception, another exception occurred:

TypeError Traceback (most recent call last) <ipython-input-21-7770240b54f8> in <module> 1 ds = yt.load('cutout_478037.hdf5') ----> 2 yt.ProjectionPlot(ds,'x',('gas','velocity_x'),width=(400,'kpc')) 3 #slc.save('testplot.pdf') 4 #slc.set_zlim(('gas','radial_velocity'),-200,200) 5 #slc.annotate_title(title)

~/yt/yt/visualization/plot_window.py in __init__(self, ds, axis, fields, center, width, axes_unit, weight_field, max_level, origin, right_handed, fontsize, field_parameters, data_source, method, proj_style, window_size, buff_size, aspect) 1609 PWViewerMPL.__init__(self, proj, bounds, fields=fields, origin=origin, 1610 right_handed=right_handed, fontsize=fontsize, window_size=window_size, -> 1611 aspect=aspect, buff_size=buff_size) 1612 if axes_unit is None: 1613 axes_unit = get_axes_unit(width, ds)

~/yt/yt/visualization/plot_window.py in __init__(self, *args, **kwargs) 780 self._plot_type = kwargs.pop("plot_type") 781 self._splat_color = kwargs.pop("splat_color", None) --> 782 PlotWindow.__init__(self, *args, **kwargs) 783 784 def _setup_origin(self):

~/yt/yt/visualization/plot_window.py in __init__(self, data_source, bounds, buff_size, antialias, periodic, origin, oblique, right_handed, window_size, fields, fontsize, aspect, setup) 215 self._field_transform[field] = linear_transform 216 self.setup_callbacks() --> 217 self._setup_plots() 218 219 def __iter__(self):

~/yt/yt/visualization/plot_window.py in _setup_plots(self) 1066 else: # symlog_transform 1067 flinthresh = 10**np.floor( np.log10( self.plots[f].cb.norm.linthresh ) ) -> 1068 mticks = self.plots[f].image.norm( get_symlog_minorticks(flinthresh, vmin, vmax) ) 1069 self.plots[f].cax.yaxis.set_ticks(mticks, minor=True) 1070 else:

~/yt/yt/visualization/plot_container.py in get_symlog_minorticks(linthresh, vmin, vmax) 129 return np.hstack((-get_log_minorticks(linthresh,-vmin)[::-1], 0) ) 130 else: --> 131 return np.hstack((-get_log_minorticks(linthresh,-vmin)[::-1], 0, 132 get_log_minorticks(linthresh, vmax))) 133

~/yt/yt/visualization/plot_container.py in get_log_minorticks(vmin, vmax) 105 expA += 1 106 else: --> 107 lmticks = np.hstack( (lmticks, np.linspace(cofA, cofB, cofB-cofA+1)*10**expA) ) 108 expA += 1 109 return np.array(lmticks)

<__array_function__ internals> in linspace(*args, **kwargs)

~/anaconda3/lib/python3.6/site-packages/numpy/core/function_base.py in linspace(start, stop, num, endpoint, retstep, dtype, axis) 119 raise TypeError( 120 "object of type {} cannot be safely interpreted as an integer." --> 121 .format(type(num))) 122 123 if num < 0:

TypeError: object of type cannot be safely interpreted as an integer. _______________________________________________ yt-users mailing list -- yt-users@python.org To unsubscribe send an email to yt-users-leave@python.org