The following presents a method for determining topographic elevation error for overland unstructured finite element meshes derived from bare earth LiDAR for use in a shallow water equations model. This thesis investigates the development of an optimal interpolation method to produce minimal error for a given element size. In hydrodynamic studies, it is vital to represent the floodplain as accurately as possible since terrain is a critical factor that influences water flow. An essential step in the development of a coastal inundation model is processing and resampling dense bare earth LiDAR to a DEM and ultimately to the mesh nodes; however, it is crucial that the correct DEM grid size and interpolation method be employed for an accurate representation of the terrain. The following research serves two purposes: 1) to assess the resolution and interpolation scheme of bare earth LiDAR data points in terms of its ability to describe the bare earth topography and its subsequent performance during relevant tide and storm surge simulations