ABSTRACT To take into account the strong nonlinearities of vertical fluxes due to small-scale heterogeneities of surface properties, more and more coupled general circulation models compute part of their atmospheric physical parameterizations, either the surface fluxes or the whole package, on the finer grid of their ocean or land model. A modification of a traditional interpolation scheme is presented to calculate the values of atmospheric variables over surface model grid points. In addition to the desirable properties of flux conservation and preservation of a constant field, the new scheme allows discontinuities in the interpolated fields at the surface model's boundaries and orographic jumps, while remaining continuous elsewhere. It can also be tuned separately for each variable. The modified scheme is then evaluated using the circulation model of the Laboratoire de Météorologie Dynamique coupled to the Laboratoire d'Océanographie Dynamique et de Climatologie tropical Pacific Ocean model using the delocalized physics method. The results show a large improvement of heat and humidity fluxes near the focus region of the South American coast in the southeastern equatorial Pacific, and a subsequent westward propagation of significant cold SST anomalies
