A semiclassical method is used to treat proton tunneling in a planar model of malonaldehyde. Classical trajectories were calculated on a realistic potential-energy surface and WKB tunneling probabilities were calculated at turning points on the barrier separating the two equivalent potential wells. The calculated ground-state splitting, 24.5 cm-1, is in good agreement with the experimental value of 21.6 cm-1. Vibrational mode selectivity was studied by calculating the splitting for all 15 modes for 2 kcal/mol excitation energy (35 kcal/mol total energy). The results show significant mode specific effects for all except some C-H stretching modes. This study demonstrates that multidimensional semiclassical tunneling can be readily and accurately treated.