The finite-element method with B splines is used for definition of vertical operators in the nonhydrostatic fully compressible dynamical core of the ALADIN system. It represents a generalization of the same method used in the hydrostatic dynamical core shared by the ALADIN system and the global forecast system ARPEGE/IFS. The method is shown to be robust enough in idealized academic tests and real simulations. Its theoretical superiority is shown when compared with the finite-difference method.