A numerical isopycnal ocean model has been designed and applied to model the Indian Ocean north of 25¿S. Vertical normal modes are used in the open boundary conditions and for selections of initial layer depths. A 21-year integration with a reduced Hellerman--Rosenstein monthly averaged wind stress has been made with 3.5-layer and 1.5-layer versions of the model. Both solutions reproduce the main features of the observed wind-driven seasonal circulation in the Indian Ocean above the main thermocline. The transient semiannual equatorial surface jets are more intense, more coherent, and in better phase agreement with observations when three layers are active. The associated undercurrents below the main thermocline are also included in the 3.5-layer-model solution. Second baroclinic-mode, reflecting, equatorial Kelvin and Rossby waves combined to give a semiannual, resonant basin mode. Experiments with an equatorial band of semiannual zonal winds suggest a very strong response of the Indian Ocean to wind forcing with this period. Further, the amplitudes of the 28--30 day oscillations in the western equatorial model region are found to be strongly damped with depth; they have upward phase propagation and downward energy propagation. ¿ American Geophysical Union 1993