I noticed your library uses output to numpy arrays in a few places
In the Fortran 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: fortran-lang/stdlib#581