Forced by different wind stress as well as different surface heat flux, three 10 year integrations using L30T63 oceanic GCM are carried out. The annual mean state, seasonal cycle and interannual variability of the upper ocean temperature over the tropical Pacific are examined. In the case of annual mean state of sea surface temperature, relults from D15c and D15e are alike and have good agreement with NCEP data. Defects in D15d include the excessively higher SST over warm pool region and higher SST along the eastern coast. The corresponding westward displacement of cold water is also noticeable. All the experiments represent the annual cycle of SST in the eastern equatorial Pacific while the amplitudes are weaker compared to NCEP data. The intensity of interannual variability of SST in eastern equatorial Pacific is markedly improved D15d and is perfectly comparative with that from NCEP data. The eastward propagations of the interannual variability of upper ocean heat content along the equator in 1982/1983 E1 Nio are well represented in all experiments. On the onset stage of 1982/1983 E1 Nio, the modeled interannual variabilities of SST and heat content are generally the same among all experiments. It is the weakest interannual variability of surface heat flux, which damps SST variability during the mature phase of E1 Nio, that results in a most intensified SST variability in D15d. Forced by different wind stress as well as different surface heat flux, three 10-year integrations using L30T63 oceanic GCM are carried out. The annual mean state, seasonal cycle and interannual variability of the upper ocean temperature over the tropical Pacific are examined. In the case of annual mean state of sea surface temperature, relults from D15c and D15e are alike and have good agreement with NCEP data. Defects in D15d include the excessively higher SST over warm pool region and higher SST along the eastern coast. The corresponding westward displacement of cold All the experiments represent the annual cycle of SST in the eastern equatorial Pacific while the amplitudes are weaker compared to NCEP data. The intensity of interannual variability of SST in eastern equatorial Pacific is markedly improved D15d and is perfectly comparative with that from NCEP data. The eastward propagations of the interannual variability of upper ocean heat content along t he equator in 1982/1983 E1 Nio are well represented in all experiments. On the onset stage of 1982/1983 E1 Nio, the modeled interannual variabilities of SST and heat content are generally the same among among all experiments. It is the weakest interannual variability of surface heat flux, which damps SST variability during the mature phase of E1 Nio, that results in a most intensified SST variability in D15d.