This study focuses on the compilation of bathymetric data for the Bay of Bengal domain, specifically for the area between 15°N to 23°N latitudes and 85°E to 95°E longitudes, to ensure compatibility with the Finite Element Method (FEM). A high-resolution color map of the region was generated using ArcGIS software, which served as the basis for discretizing the domain into a triangular mesh. Bathymetric data for the Bay of Bengal, including the Meghna River, was then compiled using a combination of recent data from the Bangladesh Inland Water Transport Authority (BIWTA) and generated data, which were interpolated using a two-way cubic spline interpolation method. The Modified Inverse Distance Weighted Interpolation (MIDWI) method was employed to map this bathymetric data to the nodal points of the triangular mesh, ensuring accurate representation of the seafloor’s varying depths. This work is essential for developing a reliable storm surge prediction model for the Bay of Bengal using FEM, as it provides the necessary geometric and bathymetric details required for precise numerical modeling.