Data from a coupled ROMS+BEC (Shchepetkin and McWilliams, 2006; Moore et al. 2013; Frischknecht et al. 2018) hindcast simulation for February-June 2018. The coupled model was run on a northern Benguela upwelling system regular grid of horizontal resolution of 3.3 km, using a vertical sigma-layer coordinate. Original model output was saved in the form of 3D daily means. These NetCDF files contain model output that was post-processed for the creation of figures included in a submitted manuscript (preliminary title "Mesopelagic particle layers of the hypoxic northern Benguela are shaped by diel vertical migration, zooplankton communities and near shore currents." by Lovecchio et al.). A description of the model setup and an evaluation of the model are provided in the manuscript and its supplement. This model run was created to complement and upscale in-situ glider data collected in the same period in the northern Benguela (published here https://doi.pangaea.de/10.1594/PANGAEA.938221).