An ensemble adjustment Kalman filter system is developed to assimilate Argo profiles into the Northwest Pacific MASNUM wave-circulation coupled model, which is based on the Princeton Ocean Model (POM). This model was recoded in FORTRAN-90 style, and some new data types were defined to improve the efficiency of system design and execution. This system is arranged for parallel computing by using UNIX shell scripts: it is easier with single models running separately with the required information exchanged through input/output files. Tests are carried out to check the performance of the system: one for checking the ensemble spread and another for the performance of assimilation of the Argo data in 2005. The first experiment shows that the assimilation system performs well. The comparison with the Satellite derived sea surface temperature (SST) shows that modeled SST errors are reduced after assimilation; at the same time, the spatial correlation between the simulated SST anomalies and the satellite data is improved because of Argo assimilation. Furthermore, the temporal evolution/trend of SST becomes much better than those results without data assimilation. The comparison against GTSPP profiles shows that the improvement is not only in the upper layers of ocean, but also in the deeper layers. All these results suggest that this system is potentially capable of reconstructing oceanic data sets that are of high quality and are temporally and spatially continuous.