This paper describes a modelling approach suitable for assessments of future scenarios for renewable energy integration in large and interconnected power systems, based on sequential optimal power flow computations that take into account variability in power consumption, in renewable power production, energy storage, and flexible demand. The approach and the implementation as an open source Python package called Power Grid And Market Analysis is described in some detail. Particular emphasis is put on the modelling of energy storage systems, and the use of storage values as a means to define storage utilisation strategies. A case study representing a 2030 scenario for the Western Mediterranean region is then analysed using this approach. The main aim of this study is to assess the benefit for the system of adding flexibility in terms of storage associated with concentrated solar power or flexible demand. But other results are also presented, such as the resulting energy mix, generation costs, price variations, and grid congestion.