In this paper we present a new tool which can be used to simplify and speed up the reconstruction of real Earth surfaces, cake-layered domains and planar fault sources for numerical simulations. The tool makes use of the CUBIT-Python interface in order to directly ’communicate’ and to allow for maximum portability across different operating systems. We will focus the use of the tool to earthquake simulations, although many other types of simulation are able to make use of the same output. Indeed the features created using this software can be applied to numerous other model scenarios, including but not limited to: water flow models, avalanche models, subsurface cavity effects on travelling waves, and earthquake simulations.