The proper orthogonal decomposition (POD) is successfully employed in a variety of projection-based methods for parametric model order reduction (pMOR) of large dynamical systems. It extracts the most energetic modes describing the dynamics of a system from time snapshots of the solution. In practice, these snapshots are computed at user-defined parameters of the system with a high fidelity model. Then either all the snapshots are concatenated to extract a global POD basis, or a different POD basis is extracted for each parameter value. In the latter case, for unseen target parameters, the POD bases need to be adapted with a dedicated technique. An established method addressing this task is the interpolation on the tangent space to the Grassmann manifold (ITSGM). It interpolates the linear subspaces spanned by the known POD bases, discarding by construction the knowledge on the modes ordering by energy levels. In this paper, we formalize an ordered reduced basis interpolation (ORBI) technique which preserves such ordering. This approach improves the adaptation accuracy of the resulting POD basis as more information is taken into account in the construction of the interpolation operator. Numerical investigations are conducted on the parametric reduced order model of a dynamical system representing a small-scale gas-bearings supported rotor. Results show that the proposed method is more accurate than the ITSGM method at a similar computation cost. Trained at only three points of the parameters space, the developed hyper reduced order model (h-ROM) performs 30x to 200x faster simulations with satisfactory accuracy, even far from the training points.