import numpy as np, h5py

h5 = h5py.File('skymap.h5', mode='w')
h5.create_dataset('ne2001/ddm', data=np.load('datacube_ddm_ne2001.npy'))
h5.create_dataset('ymw16/ddm', data=np.load('datacube_ddm_ymw16.npy'))
h5.create_dataset('ne2001/ddm_tau', data=np.load('datacube_ddm_tau_ne2001.npy'))
h5.create_dataset('ymw16/ddm_tau', data=np.load('datacube_ddm_tau_ymw16.npy'))


h5.create_dataset('ne2001/dmd', data=np.load('datacube_dmd_ne2001.npy'))
h5.create_dataset('ymw16/dmd', data=np.load('datacube_dmd_ymw16.npy'))
h5.create_dataset('ne2001/dmd_tau', data=np.load('datacube_dmd_tau_ne2001.npy'))
h5.create_dataset('ymw16/dmd_tau', data=np.load('datacube_dmd_tau_ymw16.npy'))

h5.create_dataset('ne2001/xyz', data=np.load('ne2001_xyz.npy'))
h5.create_dataset('ymw16/xyz', data=np.load('ymw16_xyz.npy'))

h5.create_dataset('gl', data=np.load('datadim_gl.npy'))
h5.create_dataset('gb', data=np.load('datadim_gb.npy'))
h5.create_dataset('dist', data=np.load('datadim_dist.npy'))
h5.create_dataset('dm', data=np.load('datadim_dm.npy'))
h5.create_dataset('x', data=np.load('xyz_x.npy'))
h5.create_dataset('y', data=np.load('xyz_y.npy'))
h5.create_dataset('z', data=np.load('xyz_z.npy'))

h5.close()