Keaton J. Burns

Just to add a reference, [Kevlahan et al. 2015](http://doi.org/10.5194/gmd-8-3891-2015) doesn't discuss multiple layers, but seems to do a good job incorporating bathymetry and immersed-boundary continents that have proper wave reflection...

Yes seeing the same here, also pinning sphinx-gallery < 0.11 for now.

**Original comment by Tomás Chor (Bitbucket: [tomchor](https://bitbucket.org/%7B9844151e-c1d7-470a-a61b-af043d145a6a%7D/)).** ---------------------------------------- I’m new to Dedalus but I can say that this feature would be very much appreciated. HDF5 isn’t the most intuitive format....

**Original comment by J. S. Oishi (Bitbucket: [jsoishi](https://bitbucket.org/%7Bab7c2db2-7642-47a8-9264-0939f601d904%7D/)).** ---------------------------------------- Hi Tomas, Thanks for your comment; it’s really helpful to hear from users! The primary disadvantage of outputting NETCDF4 is that...

@jsoishi it looks like we may be able to implement an xarray backend that lazily reads Dedalus HDF5 files: https://xarray.pydata.org/en/stable/internals/how-to-add-new-backend.html#rst-lazy-loading It might also be worth revisiting the idea of making...

Looks like netcdf does not support complex numbers, so that's out.

A basic xarray interface has been added with #213. I'll close this for now, but would definitely be open to further discussions about enhancing this interface with lazy loading, etc.,...

Thanks for the report. Here's a minimum working example showing the failure: [implicit_identity_mwe.py.txt](https://github.com/DedalusProject/dedalus/files/9261569/implicit_identity_mwe.py.txt)

Looks like the issue is that the intertwiners are not applied to the components of constant fields (a choice we actively made so that spherical BC fields were in spin...

Maybe the best we can do right now is raise an error and force you to use NCCs with bases in curvilinear problems.