Atmospheric general circulation models (AGCMs) are often “coupled” with time varying observations of boundary conditions or some other aspect of the climate system. A typical example is the Atmospheric Model Intercomparison Project (AMIP) experimental protocol, which required the specification of sea surface temperature and sea-ice extent from observed monthly means. AGCMs ordinarily incorporate the prescribed conditions by evaluating an interpolating function at each time step. Typical schemes, such as that used in the second generation GCM (GCM2) of the Canadian Centre for Climate Modelling and Analysis (CCC), do not preserve monthly means and have a smoothing effect on the interpolated time series which tends to reduce the amplitude of annual cycle and interannual variability of sea surface temperature (SST). By solving a large set of linear equations, a simple linear time-interpolation scheme that preserves the observed monthly mean SST and hence its variability can be obtained. The new scheme improves upon that used previously in CCC GCM2 by eliminating the substantial loss of interannual variability (up to 20%) and the small attenuation of the annual cycle (less than 4% on average) incurred with the old scheme. The improved linear interpolation scheme is easily adapted to other quantities.