Hi, The stars will be flat arrays unless you deposit them -- you might want to check out the particle deposition routines, like cic and deposit. -Matt On Mon, Apr 29, 2019 at 5:56 AM Vadlamani Samhitha <vadlamani.samhitha@gmail.com> wrote:
Hi,
I am trying to add value of a particle field to the gas field (same units) by accessing a data region of fixed resolution. I expect the fields to have same sizes but it doesn't seem to be the case. Could someone please shed some light on this? Thank you.
DIM=256
dd=ds.arbitrary_grid(left_edge= ds.arr([(com_x_1.value-10), (com_y_1.value-10), (com_z_1.value-10)], 'kpc'), right_edge= ds.arr([(com_x_1.value+10), (com_y_1.value+10), (com_z_1.value+10)], 'kpc'), dims=[DIM,DIM,DIM])
m=dd["gas", "cell_mass"].in_units('Msun')
starmass=dd[("new_stars", "particle_mass")].in_units('Msun')
b=m.value+starmass.value
Error:
Traceback (most recent call last):
File "gasinflow.py", line 95, in <module>
b=m.value+starmass.value
ValueError: operands could not be broadcast together with shapes (16777216,) (8978,)
_______________________________________________ yt-users mailing list -- yt-users@python.org To unsubscribe send an email to yt-users-leave@python.org