A new and alternative numerical approach towards efficient and effective studying of the physics involved in the earthquake phenomena was developed. LSMearth, capable of simulating the fracture and failure of the rocks in microscopic scale was coupled with the macroscopic, elasto-static element-free Galerkin model. A coupler was developed to transfer the physical values between the two to describe the deformation. Since the two are written in different programming languages, C++ and F90, all data are transferred to one another through the coupler by initiating a function call by the job controller and/or EFG software, written in F90. This preliminary report illustrates the methodologies implemented for this coupling process and a simple analysis case, which shows to be a promising approach for the future studies.