A model for the time and space variation of the internal boundary-layer height over a land area with an irregular coastline is presented. It is based on the analytical model of the boundary-layer height proposed by Gryning and Batchvarova (1990) and Batchvarova and Gryning (199 1). The model accounts for the temperature jump and the mean vertical air motion at the top of the internal boundarylayer. Four cases from experiments in Nanticoke and Vancouver are used for model validation. The agreement between the calculated and measured internal boundary layer height at the observational sites is fairly good. The input information for the model consist of wind speed and direction, friction velocity and kinematic heat flux in time and space for the area, and the potential temperature gradient and the mean vertical air motion above the internal boundary layer. For the experiments used in the validation the effect of subsidence is relatively important in the afternoon under low wind speed high pressure conditions, lowering the height of the internal boundary layer by up to lo%, and it is negligible in the morning hours. The effect of the mixing height over the sea is found to be negligible.