Numerical Simulation of Baroclinic Waves with a Parameterized Boundary Layer