We discuss the 2+1 flavor Polyakov loop enhanced Nambu--Jona-Lasinio model in
a finite volume. The main objective is to check the volume scaling of
thermodynamic observables for various temperatures and chemical potentials. We
observe the possible violation of the scaling with system size in a
considerable window along the whole transition region in the T−μq​ plane
