Abstract
The main purpose of this study is to develop a three dimensional numerical analysis code for the prediction of tsunami drift behavior. First, we developed a numerical method to calculate a shape and external forces for each drift segment divided along the orthogonal grid accurately and stably. The method is based on the vector analyses. Then, we developed the new code by building the method into the three dimensional fluid analysis code which used the FAVOR method and the VOF method. We applied the new code to two types of hydraulic experiments which have different initial conditions. As a result, the code simulated each drift behavior qualitatively and quantitatively.