The recently developed mixed quantum-classical propagation method is extended to treat tunneling effects in multidimensional systems. Formulated for systems consisting of a quantum primary part and a classical bath of heavier particles, the method employs a frozen Gaussian description for the bath degrees of freedom, while the dynamics of the quantum subsystem is governed by a corrected propagator. The corrections are defined in terms of matrix elements of zeroth-order propagators. The method is applied to a model system of a double-well potential bilinearly coupled to a harmonic oscillator. The extension of the method, which includes nondiagonal elements of the correction propagator, enables an accurate treatment of tunneling in an antisymmetric double-well potential.