We implement an adaptive isostable reduction strategy to obtain a data-driven reduced order model that captures the dynamics of observables in a computational model for fluid flow over an airfoil at moderate Reynolds numbers. The resulting model characterizes the response to both time-varying inflow conditions and Dirichlet boundary conditions on the surface of the airfoil meant to represent suction or blowing through the action of surface jets. The resulting reduced order model behaviors agree well with the dynamics of the full order model simulations in response to both open loop and closed loop inputs. This study provides a proof of concept that reduced order modeling techniques adaptive isostable coordinates can be successfully used in realistic fluid flow models using geometries with practical relevance.