Designing dynamical cores that meet the challenges imposed by simulating the continuous equations that govern geophysical flows has a long history (Williamson, 2007). Various numerical methods are employed to achieve accuracy, efficiency, and stability. However, careful compromises are required because these goals are often in conflict: significant dissipation helps with stability at the cost of accuracy, and high-order schemes deliver accuracy at the expense of computing cost. This work explores the discontinuous Galerkin (DG) method for simulating atmospheric motions.