一、数据收集
选定研究区范围
搜集区域内断层线数据
数据一般为经纬度，建模一般用平面坐标。GMT可以方便的完成转换，对于中央经线为105的研究区域，其UTM带号为48。GMT转换坐标时，代码如下：gmt mapproject fault.txt -Ju+48/1:1 -C -F > fault-utm.txt搜集断层产状数据，并根据断层产状数据将其整理成切向量
详情见此文，断层产状换算成切向量搜集研究区DEM数据
推荐使用GMT官方数据，很方便使用，详情见earth_relief: 全球地形起伏数据。
- 可使用
grd2xyz模块将研究区dem输出成点坐标格式，方便使用 - 使用
mapproject模块将上述数据投影成UTM平面坐标
- 在
crust1.0中搜集研究区的上、中、下地壳层面数据 - 大致了解研究区的莫霍面深度，推荐使用crust1.0，以对整个模型各个层面的厚度有个概念。
二、将数据导入到GOCAD中去
使用GOCAD进行初步的模型构建
-
将上、中、下地壳分层及DEM数据导入GOCAD中，如下所示：
层面信息导入 断层数据导入
-
将断层分类好，如下：
断层数据 -
导入GOCAD，记得导入为断层线
将断层点导入为线 - 如果一些断层成密度不够，适当加密。
3. 将断层线投影到DEM上
投影断层线
查看投影效果，断层已经到DEM上去了
