Abstract.A new set of nonlocal boundary conditions is 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.Mathematics Subject Classification. 35L50, 65M06, 76B99, 86A05.