The constructed 3D numerical model of thermal convection in the upper mantle of the Earth is based on applying the vector potential-vorticity variables and the fractional step method. Detailed testing of the model is carried out. Results of 3D modeling of convection beneath cratons of central Asia are presented. The calculation results illustrate the structure of convection flows.
