Hi all,

I'm trying to do a volume render of the temperature field inside a
sphere with the following code:

val, loc = ds.find_max('density')
sp = ds.sphere(loc, SP_DIST)
source = VolumeSource(sp, 'temperature')
scene = Scene()
scene.add_source(source)

# v is calculated based on angular momentum
camera = scene.add_camera(sp, lens_type= 'perspective')
cam_position = loc + v * CAM_DIST
camera.set_position(cam_position)
camera.set_focus(loc)

source.set_field('temperature')
source.set_log(True)
bounds = (UPPER, LOWER)
tf = yt.ColorTransferFunction(np.log10(bounds))
tf.add_layers(12, colormap='jet')
source.tfh.tf = tf
source.tfh.bounds = bounds
source.tfh.plot('transfer_function.png', profile_field='temperature')
scene.save("render_t={:.3f}".format(times[i]), sigma_clip= 1)

Upon saving though I encounter the error:

Traceback (most recent call last):
   File "plot_3d.py", line 51, in <module>
     scene.save("render_t={:.3f}".format(times[i]), sigma_clip= 1)      
#last sigma = 2
   File
"/opt/apps/pkgs/cpython/3.7.1/intel64/intel_2019.2/lib/python3.7/site-packages/yt/visualization/volume_rendering/scene.py",
line 310, in save
yt : [INFO     ] 2020-09-15 10:34:02,783 Saving render render_t=0.000.png
     self._last_render.write_png(fname, sigma_clip=sigma_clip)
   File
"/opt/apps/pkgs/cpython/3.7.1/intel64/intel_2019.2/lib/python3.7/site-packages/yt/data_objects/image_array.py",
line 289, in write_png
     scaled = self.rescale(inline=False)
   File
"/opt/apps/pkgs/cpython/3.7.1/intel64/intel_2019.2/lib/python3.7/site-packages/yt/data_objects/image_array.py",
line 242, in rescale
     np.clip(out, 0.0, 1.0, out)
   File "<__array_function__ internals>", line 6, in clip
   File
"/opt/apps/pkgs/cpython/3.7.1/intel64/intel_2019.2/lib/python3.7/site-packages/numpy/core/fromnumeric.py",
line 2084, in clip
     return _wrapfunc(a, 'clip', a_min, a_max, out=out, **kwargs)
   File
"/opt/apps/pkgs/cpython/3.7.1/intel64/intel_2019.2/lib/python3.7/site-packages/numpy/core/fromnumeric.py",
line 61, in _wrapfunc
     return bound(*args, **kwds)
   File
"/opt/apps/pkgs/cpython/3.7.1/intel64/intel_2019.2/lib/python3.7/site-packages/numpy/core/_methods.py",
line 132, in _clip
     um.clip, a, min, max, out=out, casting=casting, **kwargs)
   File
"/opt/apps/pkgs/cpython/3.7.1/intel64/intel_2019.2/lib/python3.7/site-packages/numpy/core/_methods.py",
line 85, in _clip_dep_invoke_with_casting
     return ufunc(*args, out=out, **kwargs)
   File
"/opt/apps/pkgs/cpython/3.7.1/intel64/intel_2019.2/lib/python3.7/site-packages/yt/units/yt_array.py",
line 1401, in __array_ufunc__
     "added to YTArray." % (str(ufunc), len(inputs)))
RuntimeError: Support for the <ufunc 'clip'> ufunc with 3 inputs has not
beenadded to YTArray.

I would also like to eventually add a colorbar to the side of the
created image which represents the transfer function and annotate text
as well. Does the render functionality allow me to do this? Any help on
these topics would be much appreciated.

Many thanks,
Sam Patrick
_______________________________________________
yt-users mailing list -- yt-users@python.org
To unsubscribe send an email to yt-users-leave@python.org
https://mail.python.org/mailman3/lists/yt-users.python.org/
Member address: matthewturk@gmail.com