Yup, just about -- but try instead doing all_data_level_max['elec_den'].  I *think* it will be zero in the places not included in the cutout.

On Sun, Jan 31, 2021 at 4:10 PM Salome Mtchedlidze <salomchedlidze@gmail.com> wrote:
Hi Matthew,
Did you mean something like this?

def _create_elec_den(ds):
    def _elec_den(field,data):
        return ()

all_data_level_max = ds.covering_grid(----)
cutout_obj = sph2 - sph1
all_data_level_max._data_source = cutout_obj

If it is so, then
gives me again flattened array..
yt-users mailing list -- yt-users@python.org
To unsubscribe send an email to yt-users-leave@python.org
Member address: matthewturk@gmail.com