  In the present paper, the primitive equations, which can be used to simulate
the large scale motion of ocean and atmosphere, are considered in the
three-dimensional domain bounded below by a fixed solid boundary and above by a
free moving boundary. The global existence and uniqueness of strong solutions
are established and the long time convergence to the equilibrium state is
showed either at exponential rate for horizontal periodic domain or at
algebraic rate for horizontal whole space.
