A three-dimensional isopycnic-coordinate ocean model for the study of internal tides is presented. In this model, the ocean interior is viewed as a stack of isopycnic layers, each characterized by a constant density. The isopycnic coordinate performs well at tracking the depth variance of the thermocline, and is suitable for simulation of internal tides. This model consists of external and internal modes, and barotropic and baroclinic motions are calculated in the two modes, respectively. The capability of simulating internal tides was verified by comparing model results with an analytical solution. The model was then applied to the simulation of internal tides in the South China Sea (SCS) with the forcing of M2 and K1 tidal constituents. The results show that internal tides in the SCS are mainly generated in the Luzon Strait. The generated M2 internal tides propagate away in three different directions (branches). The branch with the widest tidal beam propagates eastward into the Pacific Ocean, the most energetic branch propagates westward toward Dongsha Island, and the least energetic branch propagates southwestward into the basin of the SCS. The generated K1 internal tides propagate in two different directions (branches). One branch propagates eastward into the Pacific Ocean, and the other branch propagates southwestward into the SCS basin. The steepening process of internal tides due to shoaling effects is described briefly. Meridionally integrated westward energy fluxes into the SCS are comparable to the meridionally integrated eastward energy fluxes into the Pacific Ocean. A three-dimensional isopycnic-coordinate ocean model for the study of internal tides is presented. In this model, the ocean interior is viewed as a stack of isopycnic layers, each characterized by a constant density. The isopycnic coordinate performs well at tracking the depth variance of the thermocline, and is suitable for simulation of internal tides. This model consists of external and internal modes, and barotropic and baroclinic motions are calculated in the two modes, respectively. The capability of simulating internal tides was verified by comparing model results with The analytical model was then applied to the simulation of internal tides in the South China Sea (SCS) with the forcing of M2 and K1 tidal constituents. The results show that internal tides in the SCS are mainly generated in the Luzon Strait. The generated M2 internal tides propagate away in three different directions (branches). The branch with the widest tidal beam propagates eastward into the Pacif ic Ocean, the most energetic branch propagates westward toward Dongsha Island, and the least energetic branch propagates southwestward into the basin of the SCS. The generated K1 internal tides propagate in two different directions (branches). One branch propagates eastward into the Pacific Ocean, and the other branch propagates southwestward into the SCS basin. The steepening process of internal tides due to shoaling effects is briefly briefly.