Аннотация:In this manuscript we discuss a new modification of the grid-characteristic method for calculating the seismic stability of buildings in the Far North, taking into account the presence of soil piles and soil subsidence between them with a topography of an arbitrary shape. We introduce a new method of covering the integration domain with various types of computational grids. This hybrid approach of the use of computational grids to solve this class of problems is discussed in detail in the manuscript. Attention is paid to the issue of the grids generation, the conformity of contacts between them and the possibility of using interpolation, including the introduction of hierarchical coordinate step refinement. The types of grids under consideration are either Cartesian (with a constant step in Cartesian coordinates), or structured curvilinear (body-fitted and boundary-conforming) to take into account soil subsidence of arbitrary shape. A number of test calculations have been performed.