Active gust load alleviation is an important technology for designing future passenger airplanes to be lighter and thus more environmentally friendly. Unsteady Reynolds-averaged Navier–Stokes (URANS) simulations are typically used to accurately calculate gust loads, but because of their high computational cost, they can only be performed at a few selected operating points. In simpler potential theory models, stall is neglected, resulting in loss of accuracy. In this paper, a low-order unsteady aerodynamics wing model is presented, which is able to represent well compressible flow with stall. Furthermore, the model offers the possibility to modularly incorporate actuators, which allows the design and evaluation of active load alleviation systems. The model is based on a conventional unsteady 2D airfoil model including a dynamic stall model. The dynamic stall model requires viscous steady coefficients, e.g. from 2D steady RANS computations. This 2D airfoil model is coupled with a 3D steady-state lifting line model. The model is applied to the LEISA research airplane and extensively validated with URANS results. It performs well in calculating gust loads with and without simultaneous flap deflections, and provides significantly more accurate results in the case of stall than when stall is neglected.