A model is developed for use in integrated modelling codes to predict the height, width and shape of the H-mode pedestal as well as the frequency and width of edge localized modes (ELMs). The model for the H-mode pedestal in tokamak plasmas is based on flow shear reduction of anomalous transport, while the periodic ELM crashes are triggered by MHD instabilities. The formation of the pedestal and the L-H transition in this model are the direct result of E→r × B flow shear suppression of transport. Suppression of the anomalous transport enhances the role of neoclassical transport in the pedestal region. The ratio of suppression of anomalous thermal transport in electron and ion channels controls the ratio of electron to ion temperature at the top of the pedestal. Two mechanisms for triggering ELMs are considered. ELMs are triggered by ballooning modes if the pressure gradient exceeds the ballooning limit or by peeling modes if the edge current density exceeds the peeling mode criterion. The models for the pedestal and ELMs are used in a predictive integrated modelling code to follow the time evolution of tokamak discharges from L-mode through the transition from L-mode to H-mode, with the formation of the H-mode pedestal, and, subsequently, the triggering of ELMs. The objective is to produce self-consistent predictions of the width, height and shape of the H-mode pedestal and the frequency of ELMs. The dependences of pedestal temperature, pedestal width and ELM frequency as a function of plasma heating power, magnetic field and density are discussed.