Keaton J. Burns

While loading tasks to DataArrays seems to work great, I think loading full save files to DataSets in general might be bad -- if the variables have different grids for...

In the last commit I've also added the xarray interface to the post-processing tutorial notebook alongside the current h5py instructions. Maybe keeping them both makes sense until we have more...

Rebased onto latest master, and I think good to merge. This does not include lazy loading, but we could add that after if it's really a problem for big datasets,...

**Original comment by Keaton Burns (Bitbucket: [kburns](https://bitbucket.org/%7B3d3e64f1-bf12-45df-b655-4543d8fb34c4%7D/)).** ---------------------------------------- Seems like adding a method to the IVP solver would be good – but have to be careful about invalidating the L.X...

Yes, good catch, I think you're exactly right -- `solver.domain.grid_dtype` should be passed when the `GlobalArrayReducer` is built. Would you like to put in a PR for that fix?

Sorry I missed this update! Hmm, this is tricky -- in a complex domain, Dedalus will treat everything (including fields that are the output of the `abs()` operator) as complex....

Hi Gökhan, sorry for the long delay addressing this. I think it's a good question and solution approach. But we're generally moving away from the meta keywords (they're gone completely...

It looks like the zero_max function in AMD is implementing a hard max rather than a soft max. This may be part of the problem -- computing the hard max...

Well I'm not sure which closure @wenegrat was referring to, but it doesn't look like we currently have any regularization on the constant smagorinsky eddy viscosity, right?

I think it's only guaranteed to be positive on the dealiased (scales=3/2) grid points, where it is calculated, but the spectral interpolant may end up being negative at other points,...