Predicting the performance of large turbine arrays requires the understanding of many physical factors, such as array geometry, turbine operation, inflow conditions and turbulent wake mixing. Due to the large parameter space that an array may be optimised over, low-order models with low computational cost are often employed. This paper extends one of these models, the inviscid–viscous coupled model, for multi-row turbine modelling. Firstly, an extension to the inviscid actuator disc theory is presented by removing the limit on the number of discrete streamtubes computed. The extended model allows for the quantification of the impact of shear in the bypass and core flows separately. In particular, it is shown that averaging a sheared bypass flow profile can result in a substantial over-prediction of the power of a turbine in a laterally bounded flow as the effective blockage of the flow increases. The model is also used to confirm that an approximation using a limited number of streamtubes in some previous applications of the inviscid–viscous approach has a negligible impact on the results. Secondly, we explore the performance of a multi-row array with either uniform or varying turbine resistance across different rows. Results suggest that by varying resistance across rows, the array may outperform the uniform resistance case. The performance gain is dependent, however, on the arrangement and inter-turbine spacing both in the spanwise and streamwise directions.