Simulation data of the density and phase modulation shown in Fig use : import h5py h5_file = h5py.File(path, 'r') image = h5_file['image'] # get meta_data experiment data meta_data = h5_file.attrs['key'], keys = {t_max,conversion_to_um,pixel_size,magnification,I} units of the image are atoms / pixel , to transform in atoms / µm use image = image / pixel_size * magnification # get meta_data simulation data meta_data = h5_file.attrs['key'], keys = {x_step_size,t_step_size,x_grid,t_grid,I} units of the image are atoms / pixel , to transform in atoms / µm use image = image / x_step_size see python file beside. The Shapiro curves are stored in csv files as in Fig_1.