The filling factor φ against pressure P in [Pa]. This figure is same as Fig. 12, but plotted with a linear scale of φ and reversal of xy axis to compare with previous studies (see Fig. 4 in Seizinger et al. 2012). The dotted line is the result of numerical simulations in the high-density region (φ ≳ 0.1) in Seizinger et al. (2012) and the thin solid line is the fitting formula proposed by Güttler et al. (2009). Our results consistently connect to the previous simulations in the high-density region.