Describe the bug
According to the multiscene documentation, the property all_same_area does:

Determine if all contained Scenes have the same ‘area’.

However, I have created a multiscene where all scenes have the same area (they just differ between datasets), yet the property returns Fa

https://earthenginepartners.appspot.com/science-2013-global-forest

There are various examples of how to define an AreaDefinition throughout the documentation and in test files like https://github.com/pytroll/pyresample/blob/master/pyresample/test/test_files/areas.yaml. However, some of these have projection definitions with only a defined (equatorial radius). This is technically incomplete and although it works in PROJ, rasterio/GDAL are a little more strict an

#747 flagged that the basemap package doesn't exist on Sandbox. I thought at first this was just a version discrepancy but actually it turns out that basemap has been deprecated and replaced with cartopy. The climate driver rainfall notebook (and any others that use basemap) should be changed to use cartopy.

I've been using nb_black lately and it's wonderful. https://github.com/dnanhkhoa/nb_black

We can add this with %load_ext lab_black for our notebook examples and it will auto-format every cell to fit black to make it more clean and readable without us having to worry about it.

Looping over the data_vars in fit/predict/transform/... should be done with parallel.

It's really impressive how much pyomeca has progressed! Have you ever considered including some form of simple 3D marker visualization?

This should use the __getitem__ logic and return None if __getitem__ files so that the two are consistent (as expected from dicts and other mappings)