Hydrodynamic model outputs are used in urban flood risk modelling, flood alert systems, and Monte Carlo hazard assessment. This study tackles an under-explored challenge wherein regular updates to the spatial characteristics of the watershed – due to factors such as changing land use – alter the watershed’s response to rainfall forcing, thus rendering existing model outputs obsolete. Because state-of-the-art hydrodynamic models are computationally expensive, frequently re-running simulations can be costly. Modularization addresses this problem by requiring re-computation only for a limited domain affected by the land use changes. This article introduces a novel approach by modularizing the 2D domain into independent sub-domains before (‘discrete’) and after (‘abstract’) the numerical computations. Using the Hydrologic Engineering Center River Analysis System (HEC-RAS) 2D model of a large urban watershed in Houston as an illustrative and generalizable testbed, we show that both the discrete and abstract modularization closely approximates the results from re-running the entire model. The computational cost of modularization scales linearly with model size for memory requirements as storing the solution on the interior boundaries (discrete) or throughout the domain (abstract) are necessary. This trade-off of memory for computation may facilitate advances in surrogate modelling or Monte Carlo flood risk assessment.