We describe numerical experiments with a limited domain (15¿--67¿N, 65¿ west to east) coarse-resolution two-and-a-half-layer upper ocean model. The model consists of two active variable density layers: a Niiler and Kraus (1977) type mixed layer and a pycnocline layer, which overlays a semipassive deep ocean. The mixed layer is forced with a cosine wind stress and Haney type heat and precipitation-evaporation fluxes, which were derived from zonally averaged climatological (Levitus, 1982) surface temperatures and salinities for the North Atlantic. The second layer is forced from below with (1) Newtonian cooling to climatological temperatures and salinities at the lower boundary, (2) convective adjustment, which occurs whenever the density of the second layer is unstable with respect to climatology, and (3) mass entrainment in areas of strong upwelling, when the deep ocean ventilates through the bottom surface. The sensitivity of this model to changes in its internal (mixed layer) and external (e.g., a Newtonian coupling coefficient) parameters is investigated and compared to the results from a control experiment. We find that the model is not overly sensitive to changes in most of the parameters that were tested, albeit these results may depend to some extent on the choice of the control experiment. ¿ American Geophysical Union 1990