#!/usr/bin/env python
#
# Implement LAMMPS units style "si" in H5MD.
#
# http://lammps.sandia.gov/doc/units.html
# http://www.bipm.org/en/si/si_brochure/
#
import h5py
import numpy
f = h5py.File("lammps_units_si.h5", "w")
quantities = f.create_group("h5md/quantities")
quantities["mass"] = numpy.dtype(dtype=numpy.float64)
quantities["mass"].attrs["unit"] = "kilogram"
quantities["distance"] = numpy.dtype(dtype=numpy.float64)
quantities["distance"].attrs["unit"] = "metre"
quantities["time"] = numpy.dtype(dtype=numpy.float64)
quantities["time"].attrs["unit"] = "second"
quantities["energy"] = numpy.dtype(dtype=numpy.float64)
quantities["energy"].attrs["unit"] = "joule"
quantities["velocity"] = numpy.dtype(dtype=numpy.float64)
quantities["velocity"].attrs["unit"] = "metre per second"
quantities["force"] = numpy.dtype(dtype=numpy.float64)
quantities["force"].attrs["unit"] = "newton"
quantities["torque"] = numpy.dtype(dtype=numpy.float64)
quantities["torque"].attrs["unit"] = "newton metre"
quantities["temperature"] = numpy.dtype(dtype=numpy.float64)
quantities["temperature"].attrs["unit"] = "kelvin"
quantities["pressure"] = numpy.dtype(dtype=numpy.float64)
quantities["pressure"].attrs["unit"] = "pascal"
quantities["dynamic_viscosity"] = numpy.dtype(dtype=numpy.float64)
quantities["dynamic_viscosity"].attrs["unit"] = "pascal second"
quantities["charge"] = numpy.dtype(dtype=numpy.float64)
quantities["charge"].attrs["unit"] = "coulomb"
quantities["dipole"] = numpy.dtype(dtype=numpy.float64)
quantities["dipole"].attrs["unit"] = "coulomb metre"
quantities["electric_field"] = numpy.dtype(dtype=numpy.float64)
quantities["electric_field"].attrs["unit"] = "volt per metre"
quantities["density"] = numpy.dtype(dtype=numpy.float64)
quantities["density"].attrs["unit"] = "kilogram per cubic metre"