A new set of nonlocal boundary conditions are proposed for the higher modes of the 3D inviscid primitive equations. Numerical schemes using the splitting-up method are proposed for these modes. Numerical simulations of the full nonlinear primitive equations are performed on a nested set of domains, and the results are discussed.