Merge pull request #33 from ThibHlln/fix-multi-div-state-dumps

add netCDF dimension for component state divisions in dump file

changelog.rst

@@ -3,6 +3,11 @@ v0.1.0
 
 Released on 2021, ???.
 
+.. rubric:: Bug fixes
+
+* fix dump file update bug due to missing 'divisions' dimension
+  (`#32 <https://github.com/hydro-jules/cm4twc/issues/32>`_)
+
 .. rubric:: Documentation
 
 * move API reference page to doc tree root

cm4twc/components/_utils/states.py

@@ -119,8 +119,16 @@ def create_states_dump(filepath, states_info, solver_history,
 
         # state variables
         for var in states_info:
-            s = f.createVariable(var, dtype_float(),
-                                 ('time', 'history', *axes))
+            d = states_info[var].get('divisions', 1)
+            if d > 1:
+                f.createDimension('{}_divisions'.format(var), d)
+                s = f.createVariable(var, dtype_float(),
+                                     ('time', 'history', *axes,
+                                      '{}_divisions'.format(var)))
+            else:
+                s = f.createVariable(var, dtype_float(),
+                                     ('time', 'history', *axes))
+
             s.standard_name = var
             s.units = states_info[var]['units']