The objective of this paper is twofold. One is to explore extensions of the generalized mean-field empirical Galerkin model, previously developed for wake instabilities 1 to singularly actuated 2D airfoils, including a high lift configuration and a single airfoil at a high angle of attack. We present a minimum order mean field model, explore the role of the mean field as a mediator between actuation and dominant instabilities, and illustrate the need for richer mode-sets, as operating conditions change and during sharp transients. The second objective is to develop computational tools for effective extraction of empirical modes, for such models. We present two such tools. The first is the utilization of temporal harmonic analysis to separate the empirical reference into sub-reference, each representing a single (time varying) coherent flow structures' category. The second is a two-step approximate POD procedure, intended to drastically reduce the computational cost of POD approximations of very long / rich reference ensembles. It is based on a first step reference partitioning and a second step global analysis, with an a priori guaranteed resolution level.