In the data set we provide both mantle velocity and maximum principal stress orientation resulting from a geodynamical model. The data are calculated with use of the ProSpher 3D code in a spectral domain by spherical harmonics decomposition. The resolution of the model is of 120 spherical harmonics laterally and 50 km in depth. For velocity data (file set: Petrunin-etal19-Vel_XXX.dat), the 1st column represents longitude, 2nd column – latitude, 3d, 4th , 5th – longitudinal, latitudinal, and radial components of velocity in mm/yr, correspondingly. For maximum principal stress orientation data (file set: Petrunin-etal19-SH_XXX.dat), the 1st column represents longitude, 2nd column – latitude, 3d, 4th – longitudinal and latitudinal components of the unit vector representing maximum principal stress direction.