I noticed your library uses output to numpy arrays in a few places In the Fortran [stdlib](https://github.com/fortran-lang/stdlib) project, we recently got support for numpy arrays of general rank and kind. See: https://stdlib.fortran-lang.org/page/specs/stdlib_io.html#save_npy The PR where this was implemented: https://github.com/fortran-lang/stdlib/pull/581