Statistics of the relative orientation of the magnetic field and density structures for the 3D simulation. Top: HRO diagram computed at 5 Myr for four different number density intervals. Bottom: shape parameter (see Eq. (4.3)) versus number density, where the error bars were computed using Eq. (5) of Planck Collaboration Int. XXXV (2016). In this panel, shape parameter values in black squares are not shown in the HRO at the top panel. It can be seen that the magnetic field is increasingly parallel over the first three intervals, but then the parallel relative orientation decreases towards higher densities.