Thanks @salhus. I pushed a few updates to this PR:

• updated error handling for Khs
• combined the code that searches the .nc and .dat files for Khs into one segment for more robust error handling
• added the same Khs workflow for the other hydrostatics parameters that are read from Hydrostatics.dat

If a user properly adds all hydrostatics data to their Capytaine dataset (hydrostatic_stiffness, center_of_mass, center_of_buoyancy, volume), users should not need any custom .dat files.

@mancellin It would be much more robust for BEMIO to read center_of_mass, center_of_buoyancy, and volume if they were added to the output dataset by default.

Noting that right now the Khs will not read from the .nc file properly if bodies have <6 DOFs. This is an ongoing problem and I'm working on a broader fix for readCapytaine in this case.

I will merge once the tests are passing

It would be much more robust for BEMIO to read center_of_mass, center_of_buoyancy, and volume if they were added to the output dataset by default.

I was delaying this until I have a proper way to handle also the multiple body cases, but maybe we can implement something already for the single rigid body case.

thanks @mancellin, great to hear that it is in progress. No rush from my perspective. When this is added on the Capytaine side we will update the BEMIO function accordingly.