The numerical computions of tidal differential equations for Jiaozhou Bay are carried out using Flather and Heaps's scheme. The results are satisfactory compared with the observations. The computations show that the shoaling region has much more important influence on the solution than the non-linear advective terms have in Jiaozhou Bay, where the intertidal zone accounts about a quarter of the total area.