We describe here the Korea ocean prediction system that closely resembles operational numerical weather prediction systems. This prediction system will be served for real-time forecasts. The core of the system is a three-dimensional primitive equation numerical circulation model, based on ${sigma}$-coordinate. Remotely sensed multi-channel sea surface temperature (MCSST) is imposed at the surface. Residual subsurface temperature is assimilated through the relationship between vertical temperature structure function and residual of sea surface height (RSSH) using an optimal interpolation scheme. A unified grid system, named as [K-E-Y], that covers the entire seas around Korea is used. We present and compare hindcasting results during 1990-1999 from a model forced by MCSST without incorporating RSSH data assimilation and the one with both MCSST and RSSH assimilated. The data assimilation is applied only in the East Sea, hence the comparison focuses principally on the mesoscale features prevalent in the East Sea. It is shown that the model with the data assimilation exhibits considerable skill in simulating both the permanent and transient mesoscale features in the East Sea.