The development of realistic ocean circulation models may improve our understanding of the ocean processes and enable the prediction of the ocean response to the external forcing. The reanalysis and satellite data may be selected as control data for the simulated variations of the temperature fields on the surface. The 3D finite element North Pacific Circulation Model developed in the ICMMG is used to study the variability of the North Pacific circulation under the varying in time boundary conditions at the sea surface associated with the period of the El-Nino 1982, 1987 and La-Nina 1988 events. The results of the simulation for the period of 1981-1991 with wind-stress and a heat flux at the sea surface, adopted from the reanalysis data of the European Center of the Medium-Range Weather Forecast were analyzed. The anomalies of the circulation, temperature distributions and heat fluxes are discussed in comparison with the reanalysis and satellite measurements for this period.