ordering occurs when the glass susceptibility, which is the spatial average of correlation function \[20\] $$G_G(md, {\bf R})=d \sum_j \int d^2r {\overline {|< \, \psi^*_j({\bf
r}) \psi_{j+m}({\bf r}+{\bf R}) \, >|^2}} \eqno(1.1)$$ expressed in terms of the pair-field (superconducting order parameter) $\psi_j({\bf r})$ at $j$-th layer, becomes divergent. In eq.(1.1), the angular bracket and the overbar denote, respectively, the thermal and the random averages, and $d$ is the interlayer spacing. This approach can be formulated \[17-19\] as a natural extension of the nonGaussian superconducting fluctuation theory \[6\] to lower temperature at which the vortex pinning due to structural disorder is not negligible [*even in clean systems*]{}. Previously, this approach was criticized in a review paper \[21\] because it was not easy \[16\] for this approach to justify the phenomenological guess \[20,21\] that, in thermodynamic limit, the vortex solid in the pinning-free case with nonzero vortex flow resistance should be replaced in clean limit of real systems by a vortex glass with zero linear resistance. This obstacle was overcome \[17,22\] at least in high field case by combining properties of the pinning-free Abrikosov solid in 2D limit with the framework of vortex glass fluctuation based on eq.(1.1). Further, since eq.(1.1) is 