Julien Le Sommer

Slicing/isosurfaces methods for 3D grids.

5

Eventually, 3D grids should come with methods for projecting 3D fields onto 2D isosurfaces (ex : temperature on an isopycnal). This issue gathers informations / ideas as to how to...

Problem with boundary values in .core.grids._horizontal_gradient

4

numpy.gradient compute the spatial derivatives with a second order scheme except near boundaries where the derivatives are evaluated with a first order scheme. oocgcm.core.grids._horizontal_gradient wraps numpy.gradient for xarray.DataArray with dask....

use geojson for defining spatial diagnostics.

1

following a discussion with @toddringler : it would be nice to have a internal representation of the geometric features used for defining spatial integrals (horizontal averages, flux across section). this...

Subclassing DataArray to create a GriddedArray structure that refers to the grid ?

3

In the long term, it might be more pythonic for methods related to e.g. a scalar field (eg : its gradient, smoothing etc...) to be attributes of the dataArray itself....

Porting equation of state related quantities from PyDom.

10

see [PyDom.eos](http://servforge.legi.grenoble-inp.fr/projects/PyDom/browser/trunk/lib/eos.py) - PyDom.eos to oocgcm to oocgcm.oceanfuncs. - possibly in a oocgcm.oceanmodels.eos first for model specific versions of eos. - possibly with numba for non-trivial equations.

What tools for facilitating the creation of xarrays from NEMO output files?

3

we could think of using something similar to : [mpas xarray wrapper](https://github.com/pwolfram/mpas_xarray_wrapper)

Problem when slicing a NEMO 2d grid objects defined with dask chunks.

4

a slice of a grid object defined with chunks appears to be disfunctionnal. eg. > from oocgcm.oceanmodels.nemo import grids > _grd = grids.nemo_2d_grid(...,chunks=xr_chunks) > grd = _grd[500:1500,500:1500] this works :...

Should we use a external package to deal with units ?

7

We have decided to set automatically the unit of the quantities computed with modelgrids (ex : derivatives). Should we use a external package for computing the new units ? which...