Second harmonic generation (SHG) is a nonlinear optical phenomenon where two photons at the frequency ω combine to form a single photon at the second-harmonic frequency 2ω. Since that second-order process is very weak in bulk isotropic media, optical SHG responses of interfaces provide a powerful and versatile technique to probe the molecular structure and dynamics of liquid interfaces. Both local dipole contributions and non-local quadrupole contributions can be interesting to investigate different properties of the interface, such as the molecular orientation or the charge density. However, a major difficulty is to comprehend the link between the S-SHG intensity and molecular details. This article reports a numerical approach to model the polarization-resolved SHG intensities of a model vapor/liquid interface of pure water. The influence of the interfacial local environment on the hyperpolarizability is taken into account using quantum mechanical/molecular mechanics calculations. The numerical predictions are in very good agreement with experiments. We detail the hypotheses made during the modeling steps and discuss the impact of various factors on the modeled SHG intensities, including the description of the exciting field in the interfacial layer, the effect of neighboring molecules on the second-harmonic polarization, and the presence of an additional static electric field at the interface.