Topography Response to Horizontal Slab Tearing and Oblique Continental Collision: 3D Thermomechanical Numerical Modelling

This repository submission provides the 3D thermo-mechanical numerical modelling code I3ELVIS developed by the workgroup of Prof. Taras Gerya (ETH Zürich). The code is based on finite-difference and marker-in-cell methods (Gerya & Yuen, 2003, 2007; Gerya, 2019). The code solves the momentum, continuity, and energy equations on the fixed Eulerian grid and transports the physical properties by Lagrangian markers using the velocity field. The code also accounts for the major phase transitions in the Earth’s mantle and internal heat sources arising from adiabatic, radiogenic, and frictional heating. Partial melting and melt extraction processes are neglected for the sake of simplicity. A more detailed description of the code, including the governing equations and the adopted rheological model can be found in various published literature (Andrić-Tomašević et al., 2023; Boonma et al., 2023; Maiti et al., 2024). Interested users are recommended to contact Prof. Taras Gerya ([e-mail]). The repository also provides input and output files to run and reproduce model results and manuscript figures of Maiti et al., 2024 (JGR Solid Earth). Paraview States for processing .vtr result files and visualizing the model output data. Matlab script to visualise the topography from .grd files.本仓库提交的是由苏黎世联邦理工学院（ETH Zürich）塔拉斯·格里亚（Taras Gerya）教授课题组开发的三维热-力学数值模拟代码I3ELVIS。该代码基于有限差分（finite-difference）与质点-网格（marker-in-cell）方法（Gerya & Yuen, 2003, 2007; Gerya, 2019），可在固定欧拉（Eulerian）网格上求解动量、连续性与能量方程，并借助速度场通过拉格朗日（Lagrangian）质点输运物理属性。此外，该代码还考虑了地球地幔中的主要相变，以及由绝热加热、放射性加热和摩擦加热产生的内部热源；为简化计算，未纳入部分熔融与熔体提取过程。有关该代码的更详细说明（含控制方程与采用的流变模型），可参阅已发表的相关文献（Andrić-Tomašević et al., 2023; Boonma et al., 2023; Maiti et al., 2024）。有兴趣的使用者可联系塔拉斯·格里亚教授（邮箱：[e-mail]）。 本仓库同时提供了可运行并复现Maiti等人2024年发表于《JGR Solid Earth》（《地球物理学研究杂志：固体地球》）的模型结果与论文配图所需的输入、输出文件；用于处理.vtr格式结果文件、可视化模型输出数据的帕拉维尤（Paraview）状态文件；以及用于从.grd格式文件中可视化地形的Matlab脚本。