A Fortran subroutine that generates harmonic oscillator matrix elements of the nucleon-nucleon interaction Daejeon16, is provided.

We note that that there was a mistake in the fortran file Daejeon16_public_v1.f90 downloadable before June 21, 2017. We recommend to use the current version Daejeon16_public_v2.f90, where this bug is corrected and which works much faster.