Real-time flood prediction is crucial for coastal urban cities prone to flooding; enabling communities, emergency services, and transportation authorities to prepare and take necessary precautions. Traditional approaches for flood prediction involve high-fidelity physics-based hydrodynamic models, which are computationally expensive, particularly in terms of time. To address these limitations, data-driven machine learning models have been developed for real-time predictions, focusing primarily on specific spatial dimensions such as street segments, floodplain areas, or geographic locations. In this study, we propose a generalized machine learning model based on Fourier Neural Operators (FNO), capable of operating across geographic locations not seen during training. For this research, we analyze eleven storm events in Norfolk, VA, spanning from 2017 to 2022, each lasting approximately four days with a 15-minute time interval. Our FNO model utilizes water depth maps generated by the TUFLOW hydrodynamic model, with a spatial resolution of 2.5m x 2.5m, along with rainfall data from seven observation sites maintained by the Hampton Roads Sanitation District. We employ inverse-distance weighted interpolation to account for geographic variations of rainfall across the study area. For model evaluation and generalization, we implement a k-fold cross-validation approach, randomly dividing the study area into five folds. Our findings show that training for 24 time steps (360 minutes) result in a model capable of accurate predictions in the next 6 time steps (90 minutes). Additionally, sequential experiments using k-fold with eleven storm event folds demonstrate that a look-back period of 360 minutes yields similarly low error rate predictions in the 15 minutes look-ahead interval. These findings underscore the efficiency and generalizability of our FNO-based model, demonstrating its capability to effectively handle predictions over extended periods for unseen geographic locations. Additionally, our machine learning model achieves an order of magnitude speed up compared to traditional physics-based hydrodynamic models.