added supplement03.f90 to output HDF5 format