he3_configs_hdf5.py

#!/usr/bin/env python
# -*- coding: utf-8 -*-

from __future__ import division, print_function

import sys

import numpy as np
import h5py


def main():
    try:
        outfile = sys.argv[1]
    except IndexError:
        print('usage: {} output_hdf5_file'.format(sys.argv[0]))
        exit(1)

    with h5py.File(outfile, 'w') as f:
        for n, l in enumerate(sys.stdin.buffer):
            positions = np.array(l.split()[:9], dtype=float).reshape(3, -1)
            f.create_dataset('config_{}'.format(n), data=positions)


if __name__ == "__main__":
    main()