The numerical method of lines (MOLs) in coordination with the classical fourth-order Runge?Kutta (RK(4, 4)) method is used to solve shallow water equations (SWEs) for foreseeing water levels owing to the nonlinear interaction of tide and surge accompanying with a storm along the coast of Bangladesh. The SWEs are developed by extending the body forces with tide generating forces (TGFs). Spatial variables of the SWEs along with the boundary conditions are approximated by means of finite difference technique on an Arakawa C-grid to attain a system of ordinary differential equations (ODEs) of initial valued in time, which are being solved with the aid of the RK(4, 4) method. Nested grid technique is adopted to solve coastal complexities closely with least computational cost. A stable tidal solution in the region of our choice is produced by applying the tidal forcing with the major tidal constituent M2 (lunar semi-dial) along the south open-sea boundary of the outer scheme. Numericalexperimentations are carried out to simulate water levels generated by the cyclonic storm AILA along the coast of Bangladesh. The model simulated results are found to be in a reasonable agreement with the limited available reported data and observations.