Manage chunksizes for 2d or more vector in netcdf

Author: AntSimi

src/py_eddy_tracker/observations/observation.py

@@ -839,6 +839,17 @@ def create_variable(
         scale_factor=None,
         add_offset=None,
     ):
+        dims = kwargs_variable.get('dimensions', None)
+        # Manage chunk in 2d case
+        if dims is not None and len(dims) > 1:
+            chunk = [1]
+            cum = 1
+            for dim in dims[1:]:
+                nb = len(handler_nc.dimensions[dim])
+                chunk.append(nb)
+                cum *= nb
+            chunk[0] = int(400000 / cum)
+            kwargs_variable['chunksizes'] = chunk
         var = handler_nc.createVariable(zlib=True, complevel=1, **kwargs_variable)
         attrs = list(attr_variable.keys())
         attrs.sort()