The wave-induced velocity and pressure fields beneath a large amplitude internal solitary wave of depression propagating over a smooth, flat, horizontal, and rigid boundary in a shallow two-layer fluid are computed numerically. A numerical ocean model is utilised, the set-up of which is designed and tuned to replicate the previously published experimental results of Carr and Davies (Phys Fluids 18(1):016,601-1-016,601-10, 2006). Excellent agreement is found between the two data sets and, in particular, the numerical simulation replicates the finding of a reverse flow along the bed aft of the wave. The numerically computed velocity and pressure gradients confirm that the occurrence of the reverse flow is a consequence of boundary layer separation in the adverse pressure