Recent work has shown the significance of a nonlocal dielectric response and surface correction for nanometer-length-scale plasmonic structures. In this paper, we propose a new surface hydrodynamic model which incorporates such nonlocality and surface correction. Our approach, which is based on the hydrodynamic model, uses a numerical approach based on full-wave numerical computation and thus enables calculation of surface plasmon polariton (SPP) fields. The model also introduces a discontinuity in the normal component of the electric displacement at the interface by taking into account the change in the electron distribution at the interface. The method makes use of the Feibelman d parameter approach for the surface correction by equating the corrected Fresnel reflection of a p-polarized incident field across a planar interface. To demonstrate our method, we first examine numerical calculation of SPP propagation at a single interface structure; this work is then followed by a demonstration for a set of more complex thin-film structures. The latter demonstrates that our method is suitable for use in numerical modeling of complex nanophotonic structures.