TDEFNODE 编译安装这事说起来真是一把辛酸泪。做地壳形变研究的朋友应该都有体会手里的 GPS 速度场数据攒了好几年测站覆盖也够密就是想找一个能把块体运动、断层闭锁和应变积累综合分析的工具。商业软件要么授权费劝退要么功能不符合需求。开源方案里 TDEFNODE 几乎是绕不开的选择它是用 Fortran 写的经典地球物理反演工具专门做 GPS、地震滑动数据的块体模型建模。但这玩意儿的安装过程对习惯 Python 生态的朋友绝对是个不小的冲击没有 setup.py没有 pip install连个 conda 包都没有拿到手就是一个源码压缩包全靠 gfortran 和 make 手动编译。我第一次编译的时候卡了整整一天从 LAPACK 链接问题到隐式变量类型报错踩坑踩到怀疑人生。这篇文章就是把我从零开始编译、配置、运行 TDEFNODE 的完整过程记录下来包括环境准备、源码编译、输入文件格式说明、跑通一个最小案例以及我在实操中遇到的各种报错和解法。目标是让任何一个有 Linux 基础的人能顺着这篇文章把 TDEFNODE 从源码变成能用的可执行文件并且成功跑出一个像样的结果。不管你是在 Ubuntu、CentOS 还是 macOS 上思路都是一样的。1. TDEFNODE 是干什么的为什么值得花时间装1.1 三步讲清这个工具的定位TDEFNODE 全称是 Tectonic DeforMation and NODE由哈佛大学 Brendan Meade 团队开发和维护本质是一个基于 Fortran 的地壳形变反演程序。它能干什么我用大白话给你讲清楚。你手里有一堆 GPS 站点坐标和速度数据这些站点分布在若干个地质块体上块体之间是断层。TDEFNODE 就是你把这些数据丢进去它会自动调整块体旋转参数、断层滑动参数、应变率参数让模型预测的站点速度尽量拟合你实测的 GPS 速度场。它内部用的是弹性半空间位错理论块体旋转用欧拉矢量描述断层闭锁部分用滑动亏损来模拟最终输出每个块体的旋转参数、断层滑动速率、闭锁深度、应变率分布以及每个站点的观测残差。听起来是不是有点像玩拼图你把每个块体的运动方式想象成一块刚体拼图断层就是拼图之间的接缝。TDEFNODE 在干的事就是告诉你这块拼图怎么转、接缝处怎么滑才能让整个拼图面的形变场和你看到的最吻合。1.2 为什么非要用源码编译说实话我第一次知道 TDEFNODE 只能源码编译时也懵了一下。不过后来了解了原因就理解了这个软件是学术课题组开发的主要面向专业研究用户加上 Fortran 科学计算代码本身可移植性就没那么好直接提供源码让用户在自己机器上编译是最省事也最稳妥的发布方式。另一个原因TDEFNODE 需要链接 LAPACK 数值计算库不同系统的 LAPACK 安装路径和版本差异很大只有本地编译才能检测出正确的头文件和库路径。你可以把 LAPACK 理解成数学公式库TDEFNODE 里的矩阵求逆、最小二乘求解都依赖它。如果你的 gfortran 编译器版本和 LAPACK 版本不匹配链接阶段就会出现各种 undefined reference 报错这部分我在第 3 节会详细展开。1.3 谁适合读这篇文章如果你属于下面这几类人这篇文章就是为你准备的刚拿到 TDEFNODE 源码在 Linux 或 macOS 上编译遇到困难代码下载了但不知道从哪下手。编译成功了但看不懂输入文件格式不知道怎么准备数据或者跑了半天程序报错退出。想做形变反演但完全没接触过 Fortran 和命令行工具想知道这套流程该走通。当然如果你已经是个熟练工可以跳着看第 6 节的常见问题表格和第 7 节的进阶技巧那部分是我踩坑最密集的地方。2. 编译前的环境准备依赖装不好后面全白搭2.1 Linux 系统下的基础依赖清单在开始编译之前你需要确保系统里已经装好了这几样东西依赖作用Ubuntu/Debian 安装命令CentOS/RHEL 安装命令gfortranFortran 编译器核心工具apt install gfortranyum install gcc-gfortranmake构建工具apt install makeyum install makeliblapack-devLAPACK 数值库头文件和链接库apt install liblapack-devyum install lapack-devellibblas-devBLAS 基础线性代数库apt install libblas-devyum install blas-develwget/curl下载源码apt install wgetyum install wget检查是否已经装好可以在终端执行gfortran --version make --version ldconfig -p | grep lapack如果 gfortran 命令找不到说明还没装。如果 ldconfig 那一步没有输出说明 LAPACK 还没装。有个小细节值得注意很多人只装了 lapack但忘记了 lapack-dev 包导致编译时找不到头文件。Dev 包和运行库是两回事Dev 包里才带链接时需要的 .so 符号文件。2.2 gfortran 版本的选择逻辑这一步非常关键也是我踩过最大的一个坑。TDEFNODE 是带了多年历史的 Fortran 代码里面有很多早期 Fortran 77 风格的语法。gfortran 新版本默认启用了更严格的语法检查某些老语法在 12.x 以上版本会直接报错。我个人的经验是gfortran 9.x 到 11.x 这个区间编译 TDEFNODE 最顺滑。如果你系统里默认装的是 12 以上版本也不要慌有两个办法安装旧版本Ubuntu 可以用apt install gfortran-10然后指定make FCgfortran-10。关闭严格检查在 Makefile 里给 FFLAGS 加上-stdlegacy -fallow-argument-mismatch这两个选项。-stdlegacy的作用是告诉编译器“我这个代码可能不是那么符合现代标准你宽容一点”-fallow-argument-mismatch是为了允许函数调用时实参和形参类型不一致老代码里这种问题特别多。2.3 macOS 用户的特殊注意点macOS 用户建议直接装 Homebrew 然后执行brew install gcc lapack这里有个容易踩的坑macOS 自带的 clang 编不了 Fortran你装完 brew 后gfortran 命令通常链接到的是 brew 里的 gcc-12 或类似的特定版本需要确认which gfortran gfortran --version如果提示 command not found就检查一下/opt/homebrew/bin是否在 PATH 里。Apple Silicon 芯片的 Mac 用的是这个路径Intel 芯片的老 Mac 是/usr/local/bin。还有一个独有的大坑macOS 的链接器对动态库路径特别敏感如果在编译的最后链接阶段报错找不到 liblapack.dylib建议用export DYLD_LIBRARY_PATH/opt/homebrew/lib:$DYLD_LIBRARY_PATH试试。2.4 Windows 用户怎么办Windows 确实没法直接编译因为整个流程依赖 Unix 风格的 make 和链接器。我的建议是装 WSL2Windows Subsystem for Linux然后在 WSL 的 Ubuntu 环境里按 Linux 流程走。WSL2 的性能损失很小使用体验和一个真实的 Linux 虚拟机几乎一样。别折腾 Cygwin 和 MSYS2前者可能会出一些奇奇怪怪的路径问题后者则容易在 LAPACK 链接时找不到库浪费时间。3. 源码获取与编译全过程从 tar.gz 到可执行文件3.1 源码下载与目录结构认识TDEFNODE 的源码可以从项目官网的“Downloads”页面直接获取通常是一个类似tdefnode_x.x.tar.gz的压缩包。用 wget 下载下来后执行解压wget 官网下载链接 tar -zxvf tdefnode_x.x.tar.gz cd tdefnode_x.x这时候你用ls命令看一下目录会发现一个典型的 Fortran 项目结构。以我解压的版本为例LICENSE Makefile README.txt doc/ # 使用文档和 manual src/ # 源码文件 examples/ # 官方示例数据再进去看一眼 src 目录ls src/你会看到一堆.f90、.f结尾的文件比如tdefnode.f90是主程序和各个子模块还有sacio.f90、lapack_wrapper.f90这类功能模块。3.2 编译前必须修改的 Makefile 参数开编之前一定要打开 Makefile 看一眼不做这一步直接 make 的人十有八九会在链接阶段卡住。用 vim 或者 nano 打开 Makefilevim Makefile重点关注下面几个变量FC gfortran FFLAGS -O3 -fopenmp LAPACK_LIBS -llapack -lblasFC是编译器名如果系统默认 gfortran 版本不对这里可以改成gfortran-10。FFLAGS是编译选项-O3是最高级优化-fopenmp是开启 OpenMP 并行支持TDEFNODE 支持多核并行就是靠这个选项。LAPACK_LIBS是链接 LAPACK 库的指令-llapack -lblas的意思是链接名为 liblapack 和 libblas 的动态库。如果你的 LAPACK 安装路径比较特殊比如装在/usr/local/lib那么就需要在 Makefile 里添加链接路径LDFLAGS -L/usr/local/lib怎么确认 LAPACK 到底装在哪执行find /usr -name liblapack* 2/dev/null找到liblapack.so或liblapack.dylib的所在目录把路径填进LDFLAGS。这里有个更容易踩的细节如果是用 apt 安装的通常库文件在/usr/lib/x86_64-linux-gnu/Debian 系的这个路径经常被忽视导致链接不到。所以如果默认参数编译报错找不到 lapack首选检查这里。3.3 开始编译实际执行与日志观察确认 Makefile 没问题后在根目录执行make你会看到屏幕上刷过一大串编译信息这个过程一般需要 1 到 5 分钟取决于你机器的性能期间 gfortran 会把 .f90 文件逐行编译成目标文件.o最后再把它们链接成可执行文件。编译正常结束时根目录下会生成一个名为tdefnode的可执行文件。验证一下ls -l tdefnode file tdefnodefile命令会告诉你这是一个 ELF 可执行文件。如果这一步成功了恭喜你最难的编译阶段就过去了。如果中间报错了别慌跳到第 6 节我整理了几乎所有的典型报错和对应解法。3.4 编译成功后的标准动作编译成功后有几个推荐的标准动作能帮你少走不少弯路。先把源码里的示例数据拷贝到一个工作目录单独存放mkdir ~/tdefnode_work cp -r examples/* ~/tdefnode_work/ cd ~/tdefnode_work然后把可执行文件的路径加入环境变量或者干脆创建个软链接方便在任意目录下调用sudo ln -s /path/to/tdefnode/tdefnode /usr/local/bin/tdefnode这样下次你想运行 TDEFNODE 时直接tdefnode就能启动不用每次都在源码目录里干活。养成好习惯源码目录保持原样工作目录单独建这样调试起来不混乱。4. 输入文件格式解析官方示例跑通的关键一步4.1 控制文件.in 文件的结构编译通过只是第一关真正让人头大的是准备输入文件。TDEFNODE 的运行需要至少一个控制文件通常以.in结尾它告诉程序数据文件在哪儿、要反演哪些参数、迭代多少次、输出文件前缀等。我建议直接用官方 examples 目录里自带的输入文件来跑先别急着用自己的数据。打开官方示例的控制文件它长这样tdefnode.in文件内容是按固定格式排列的文本核心是几个命令块stations_file stations.dat fault_file faults.dat grid_file grid.dat out_prefix example_out ...stations_file指定站点文件路径。fault_file指定断层文件路径。grid_file指定网格文件路径网格用于计算应变率等衍生量。out_prefix所有输出文件的前缀比如example_out.def、example_out.sum。控制文件还有一些参数是控制反演流程的比如niter 1000 temp_start 100.0 temp_end 1.0niter是迭代次数temp_start和temp_end是模拟退火算法的起止温度。模拟退火是 TDEFNODE 使用的全局优化算法温度从高到低逐渐冷却类似金属退火过程。这个类比很贴切高温时系统接受较差解的概率高用于大范围搜索低温时逐渐收敛到最优区域。第一次跑示例不需要动这些参数用默认值就行。4.2 站点文件格式经纬度、速度与误差站点文件用来提供 GPS 测站数据。每一行通常对应一个测站格式大约是这个样子站名 经度 纬度 东向速度 北向速度 东向误差 北向误差 块体编号 S001 116.32 39.95 23.12 -5.34 1.23 1.45 1 S002 117.05 40.11 24.56 -4.98 1.10 1.30 1这里有几个容易搞混的点经纬度用的是十进制度数不是度分秒。如果你的数据是度分秒格式记得先转换不然位置全偏到海里去。速度单位通常是毫米/年mm/yr注意和你的原始数据统一。误差也必须是毫米/年而且要有东向和北向两个分量程序计算卡方时会用到。块体编号是整数告诉你这个站点属于哪个块体用于块体旋转拟合。有一个细节很容易忽略TDEFNODE 对站点文件的行格式非常敏感列与列之间可以用空格或者制表符分隔但千万不要混用而且像 S001 这种字符串不要带特殊符号否则读取时会直接中断。4.3 断层文件与网格文件断层文件定义了断层几何参数。和站点文件不同断层文件不是直接给一张地表断层线的经纬度列表而是定义断层分段segment的端点位置、滑动力学参数、闭锁深度等。以官方示例为例断层文件的格式大致是fault_name segment_id x1 y1 x2 y2 depth lock_depth slip_rate ...这里x1 y1和x2 y2是断层段的两个端点坐标depth是断层底部深度lock_depth是闭锁深度slip_rate是滑动速率。这个文件头信息很严格少一列或多一列都会导致读取错位然后计算结果全部变成 NaN。网格文件则是你在研究区域内设置的网格点TDEFNODE 会用这些点来计算应变率或绘图结果的插值。网格文件格式通常是每行一个点经度 纬度 116.0 39.0 116.5 39.0 ...网格密度取决于你的研究区域大小和计算力示例里默认的密度已经足够验证整个流程了。如果网格文件缺失有些功能会跑不了所以第一次别删它。4.4 修改示例文件改成自己的项目名理解了文件格式后我先做了一个最简单的修改把控制文件中的out_prefix改成我自己项目的名字比如test_run。同时把站点文件、断层文件的路径确认了一遍。查看一下文件是否存在路径是否对ls -l stations.dat faults.dat grid.dat确认没问题后就可以运行了。5. 正式跑通第一个案例5.1 运行命令与执行过程在示例目录下执行tdefnode tdefnode.in如果用的是绝对路径就写/path/to/tdefnode/tdefnode tdefnode.in这时终端会打出一串运行日志包括当前迭代步数、目标函数值、每个块体的旋转参数估计等。如果配置了 OpenMP还会显示正在使用多少个线程。这一步比较耗时示例数据规模不大一般几分钟内就能跑完。如果数据量很大去泡杯咖啡接着等。跑完后你会发现工作目录下多了一批文件文件名前缀都是你设置的那个out_prefix比如test_run.lonlat test_run.def test_run.out test_run.sum test_run.vel5.2 关键输出文件解读这一步是让我看懂 TDEFNODE 在干嘛的关键也是我觉得全网教程里讲得不够细的地方。test_run.out这是最主要的输出文件里面记录了拟合的最终参数包括所有块体的欧拉极、断层滑动速率、闭锁深度、应变率等。块体欧拉极一般用经度、纬度和旋转速率deg/Myr表示你可以理解成一个想象中的旋转轴和角速度用来描述整个块体的刚性运动。test_run.vel这是站点速度的预测值和观测值对比文件。每一行有测站名、观测速度、预测速度、残差。这个文件是检验拟合质量的第一手资料。理想情况下残差应该很小而且没有系统性偏差。test_run.sum这是拟合优度统计文件包括加权均方根残差WRMS、卡方值Chi-square等。WRMS 越小说明模型和数据越接近。这个文件是判断反演是否收敛的根本依据。test_run.def断层滑动分布的定义文件记录了每个断层段的深度、滑动矢量、闭锁程度等信息。如果你打开这些文件看到的数值都比较合理没有出现满屏的 NaN 或者无穷大说明你的安装和输入文件格式已经基本没问题了。5.3 怎么判断结果“对不对劲”判断拟合质量的几个经验法则WRMS 值如果残差在 1~3 mm/yr 量级说明拟合得不错。如果残差大得离谱比如几十 mm/yr要么是数据本身有问题要么是模型参数初值给得不好。残差方向如果所有站的残差都偏向同一个方向说明模型存在系统性偏差比如某个块体旋转速率被高估了。断层的滑动速率检查一下得到的滑动速率是否符合该地区的地质背景常识。比如一条倾滑断层如果反演出纯走滑那你就要回头看看断层几何参数是不是写错了。第一次跑通不代表你理解了结果但至少证明从安装到数据准备这条路是通的。接下来你才能放心地把自己的数据填进去去解决实际科学问题。6. 常见问题与避坑实录这些坑我替你踩过了6.1 编译阶段报错速查表报错信息原因解决方法gfortran: command not found编译器没装apt install gfortranundefined reference todgesv_LAPACK 链接失败检查 Makefile 的 LAPACK_LIBS 和 LDFLAGSError: Type mismatch老代码踩到新版语法检查在 FFLAGS 里加-stdlegacy -fallow-argument-mismatchmake: *** No rule to make targetMakefile 路径配置问题检查在源码根目录执行的 make不是 src 目录Error: Cannot open module file ‘xxx.mod’编译依赖顺序问题先make clean再重新执行make我在编译时遇到最多的就是undefined reference to dgesv_。dgesv 是 LAPACK 里解线性方程组的函数很多 Fortran 程序都依赖它。这个报错说明你的编译器找到了代码但链接器找不到 LAPACK 库。那时候我突然想起当时安装 Dev 包时偷懒用了最小安装后来补上liblapack-dev再执行make clean make问题就解决了。6.2 运行阶段常见问题编译过了运行时报错同样让人抓狂。这里把运行时报错也整理成一个速查表报错信息原因解决方法Error opening file输入文件路径不对检查控制文件里的文件名和当前目录NaN铺满所有输出站点文件或断层文件格式错位用head -5检查每个输入文件确认列数和顺序Segmentation fault数组越界或文件行数不足检查输入文件行数、块体编号是否越界程序卡住不退出反演参数设置不合理模拟退火温度过高调低temp_start或增加niter所有站点的速度预测一模一样块体编号没对上检查站点文件的块体编号是否超出块体数量范围特别注意这个 NaN 的问题。我首次拿到自己数据时输出文件里整个 vel 文件全是 NaN排查了一下午才发现是断层文件里多打了一个制表符。TDEFNODE 的读取程序是非常机械的格式上一个字符不对后面的数据全都会错位然后计算得到 NaN。这个问题排查起来非常反直觉因为你不知道是哪一行哪一列错位了只能一行行对齐。6.3 实战排查流程演示假如你运行tdefnode tdefnode.in后输出文件全是 NaN按这个顺序排查第一步看输入文件字符编码。如果文件是 Windows 下编辑的很可能带了回车符\rTDEFNODE 读取时会把它当成文件名的一部分直接报错打不开文件。用sed -i s/\r$// stations.dat faults.dat grid.dat清理一下。第二步用cat -n看文件每行末尾是否有多余空格或制表符。命令行检查比编辑器直观。第三步检查控制文件里各输入文件路径是否写错大小写是否敏感。第四步将断层文件的行数与站点文件的行数简化到最小比如只保留一条断层和 10 个站点看是否能跑通。如果简化版能跑通说明问题出在某个特定数据行上二分法逐段排查。实测下来90% 以上的运行报错都能通过这四步定位。6.4 其他实战技巧工作文件尽量不要放在中文路径或带空格的路径下Fortran 程序内部处理文件名可能有问题我吃过这个亏路径里有空格直接读不出来。修改完输入文件后先跑一次小规模快速测试。建议把网格文件里的网格点删掉一部分只留少量点这样跑一轮只要几秒钟能快速检验格式是否正确然后再跑完整数据。备份/lib里的动态库依赖。tdefnode 可执行文件如果要在其他机器上复制使用需要目标机器也安装匹配版本的 LAPACK否则会报error while loading shared libraries这是我在换服务器时踩到的额外坑。尽可能用静态链接编译。Makefile 里LDFLAGS加上-static把必要的库都打进可执行文件里这样换机器运行就不用再纠结动态库问题了。不过要注意-static和-fopenmp在某些系统上可能会有兼容性问题需要实际测试。7. 进阶使用与二次开发建议7.1 多核并行与效率调优TDEFNODE 编译时如果带了-fopenmp运行时会自动使用所有能用的 CPU 核但你也可以用环境变量来控制线程数。export OMP_NUM_THREADS4我建议第一次跑的时候不要把所有核都用满先限制到 2 到 4 个核试试因为反演程序不是简单的并行任务线程越多部分变量同步的开销也越大。实测中我发现某些环节比如模拟退火里的随机扰动并不适合大规模并行线程数超过 8 后加速比下降明显。如果反演参数多、数据量大建议先跑短迭代观察目标函数收敛趋势再放长跑不要一上来就十万次迭代否则你没法判断参数初值是否合理。7.2 如何把自己的数据接入 TDEFNODE这是把示例跑通后必然会遇到的事情。核心思路先用官方示例把整个流程跑通然后逐文件替换。替换站点文件把你的 GPS 站点数据整理成示例站点文件一样的格式块体编号按你的研究区划分填好。替换断层文件断层几何参数可以从地震地质资料里提取但要注意断层分段的端点坐标必须和你的站点文件坐标系统一致别一个用 WGS84 一个用别的椭球。调整控制文件把块体数量、初始值、迭代次数相应改掉。替换后先跑一个低迭代次数的快速测试确认程序不报错再慢慢加次数跑正式反演。很多新手一上来就指望一次跑成功结果出了问题也不知道是自己数据格式的问题还是参数初值的问题反而多花时间。7.3 和 GMT 绘图工具的联动TDEFNODE 的输出文件格式是文本表格很方便用 GMTGeneric Mapping Tools做可视化。用 GMT 把站点速度观测值和预测值画出来可以直观检验拟合效果。推荐的做法gmt begin velocity_map png gmt coast -R115/119/38/42 -JM15c -Ba -W0.5p gmt plot test_run.vel -Sc0.1c -W0.5p gmt end这里test_run.vel包含观测和预测速度你可以按列提取分别画箭头用不同颜色表示观测和预测叠加在断层图上。可视化这件事对结果的判断帮助极大肉眼比任何统计指标都容易发现问题。7.4 源码二次修改的切入点如果你要做更深度的定制比如修改反演算法或者增加新的约束条件可以从src/tdefnode.f90的主程序作为起点找到对应的模块修改。建议先通读doc/目录下的 manual 文档理解每个子函数的作用再动手改。二次开发时配套建议用gfortran -g -O0编译一个调试版本这个版本没有优化适合用 GDB 逐行打断点排查问题。8. 写在最后的个人经验TDEFNODE 这套工具说难装也难装说简单也简单核心就三件事装对编译器、配好 LAPACK、老老实实按格式准备输入文件。我头一回编译时一天之内从对 Fortran 一无所知到在终端里把几十条报错一条条啃下来整个过程非常磨人但回头看其实思路很清晰先查 Makefile 里的变量再查依赖库路径最后查输入文件格式一步步排查总能解决。最后分享一个小技巧TDEFNODE 官方文档doc 目录里的 manual写得非常详实但很多人习惯跳过它直接编译。我建议你在等编译的过程中把 manual 里关于输入文件格式的章节从头到尾读一遍大概花掉一个多小时但能帮你少走一整天的弯路。我就是没读文档直接上结果在输入文件格式上卡了最久绕了很大一圈才明白。源码编译类工具就是这样一个门槛跨过去之后你会觉得它的数据组织和反演思路都极其干净直接真正把数据丢进去得到合理结果的那一刻前面所有的折腾都值了。
TDEFNODE编译安装实战：从源码配置到地壳形变反演全流程
TDEFNODE 编译安装这事说起来真是一把辛酸泪。做地壳形变研究的朋友应该都有体会手里的 GPS 速度场数据攒了好几年测站覆盖也够密就是想找一个能把块体运动、断层闭锁和应变积累综合分析的工具。商业软件要么授权费劝退要么功能不符合需求。开源方案里 TDEFNODE 几乎是绕不开的选择它是用 Fortran 写的经典地球物理反演工具专门做 GPS、地震滑动数据的块体模型建模。但这玩意儿的安装过程对习惯 Python 生态的朋友绝对是个不小的冲击没有 setup.py没有 pip install连个 conda 包都没有拿到手就是一个源码压缩包全靠 gfortran 和 make 手动编译。我第一次编译的时候卡了整整一天从 LAPACK 链接问题到隐式变量类型报错踩坑踩到怀疑人生。这篇文章就是把我从零开始编译、配置、运行 TDEFNODE 的完整过程记录下来包括环境准备、源码编译、输入文件格式说明、跑通一个最小案例以及我在实操中遇到的各种报错和解法。目标是让任何一个有 Linux 基础的人能顺着这篇文章把 TDEFNODE 从源码变成能用的可执行文件并且成功跑出一个像样的结果。不管你是在 Ubuntu、CentOS 还是 macOS 上思路都是一样的。1. TDEFNODE 是干什么的为什么值得花时间装1.1 三步讲清这个工具的定位TDEFNODE 全称是 Tectonic DeforMation and NODE由哈佛大学 Brendan Meade 团队开发和维护本质是一个基于 Fortran 的地壳形变反演程序。它能干什么我用大白话给你讲清楚。你手里有一堆 GPS 站点坐标和速度数据这些站点分布在若干个地质块体上块体之间是断层。TDEFNODE 就是你把这些数据丢进去它会自动调整块体旋转参数、断层滑动参数、应变率参数让模型预测的站点速度尽量拟合你实测的 GPS 速度场。它内部用的是弹性半空间位错理论块体旋转用欧拉矢量描述断层闭锁部分用滑动亏损来模拟最终输出每个块体的旋转参数、断层滑动速率、闭锁深度、应变率分布以及每个站点的观测残差。听起来是不是有点像玩拼图你把每个块体的运动方式想象成一块刚体拼图断层就是拼图之间的接缝。TDEFNODE 在干的事就是告诉你这块拼图怎么转、接缝处怎么滑才能让整个拼图面的形变场和你看到的最吻合。1.2 为什么非要用源码编译说实话我第一次知道 TDEFNODE 只能源码编译时也懵了一下。不过后来了解了原因就理解了这个软件是学术课题组开发的主要面向专业研究用户加上 Fortran 科学计算代码本身可移植性就没那么好直接提供源码让用户在自己机器上编译是最省事也最稳妥的发布方式。另一个原因TDEFNODE 需要链接 LAPACK 数值计算库不同系统的 LAPACK 安装路径和版本差异很大只有本地编译才能检测出正确的头文件和库路径。你可以把 LAPACK 理解成数学公式库TDEFNODE 里的矩阵求逆、最小二乘求解都依赖它。如果你的 gfortran 编译器版本和 LAPACK 版本不匹配链接阶段就会出现各种 undefined reference 报错这部分我在第 3 节会详细展开。1.3 谁适合读这篇文章如果你属于下面这几类人这篇文章就是为你准备的刚拿到 TDEFNODE 源码在 Linux 或 macOS 上编译遇到困难代码下载了但不知道从哪下手。编译成功了但看不懂输入文件格式不知道怎么准备数据或者跑了半天程序报错退出。想做形变反演但完全没接触过 Fortran 和命令行工具想知道这套流程该走通。当然如果你已经是个熟练工可以跳着看第 6 节的常见问题表格和第 7 节的进阶技巧那部分是我踩坑最密集的地方。2. 编译前的环境准备依赖装不好后面全白搭2.1 Linux 系统下的基础依赖清单在开始编译之前你需要确保系统里已经装好了这几样东西依赖作用Ubuntu/Debian 安装命令CentOS/RHEL 安装命令gfortranFortran 编译器核心工具apt install gfortranyum install gcc-gfortranmake构建工具apt install makeyum install makeliblapack-devLAPACK 数值库头文件和链接库apt install liblapack-devyum install lapack-devellibblas-devBLAS 基础线性代数库apt install libblas-devyum install blas-develwget/curl下载源码apt install wgetyum install wget检查是否已经装好可以在终端执行gfortran --version make --version ldconfig -p | grep lapack如果 gfortran 命令找不到说明还没装。如果 ldconfig 那一步没有输出说明 LAPACK 还没装。有个小细节值得注意很多人只装了 lapack但忘记了 lapack-dev 包导致编译时找不到头文件。Dev 包和运行库是两回事Dev 包里才带链接时需要的 .so 符号文件。2.2 gfortran 版本的选择逻辑这一步非常关键也是我踩过最大的一个坑。TDEFNODE 是带了多年历史的 Fortran 代码里面有很多早期 Fortran 77 风格的语法。gfortran 新版本默认启用了更严格的语法检查某些老语法在 12.x 以上版本会直接报错。我个人的经验是gfortran 9.x 到 11.x 这个区间编译 TDEFNODE 最顺滑。如果你系统里默认装的是 12 以上版本也不要慌有两个办法安装旧版本Ubuntu 可以用apt install gfortran-10然后指定make FCgfortran-10。关闭严格检查在 Makefile 里给 FFLAGS 加上-stdlegacy -fallow-argument-mismatch这两个选项。-stdlegacy的作用是告诉编译器“我这个代码可能不是那么符合现代标准你宽容一点”-fallow-argument-mismatch是为了允许函数调用时实参和形参类型不一致老代码里这种问题特别多。2.3 macOS 用户的特殊注意点macOS 用户建议直接装 Homebrew 然后执行brew install gcc lapack这里有个容易踩的坑macOS 自带的 clang 编不了 Fortran你装完 brew 后gfortran 命令通常链接到的是 brew 里的 gcc-12 或类似的特定版本需要确认which gfortran gfortran --version如果提示 command not found就检查一下/opt/homebrew/bin是否在 PATH 里。Apple Silicon 芯片的 Mac 用的是这个路径Intel 芯片的老 Mac 是/usr/local/bin。还有一个独有的大坑macOS 的链接器对动态库路径特别敏感如果在编译的最后链接阶段报错找不到 liblapack.dylib建议用export DYLD_LIBRARY_PATH/opt/homebrew/lib:$DYLD_LIBRARY_PATH试试。2.4 Windows 用户怎么办Windows 确实没法直接编译因为整个流程依赖 Unix 风格的 make 和链接器。我的建议是装 WSL2Windows Subsystem for Linux然后在 WSL 的 Ubuntu 环境里按 Linux 流程走。WSL2 的性能损失很小使用体验和一个真实的 Linux 虚拟机几乎一样。别折腾 Cygwin 和 MSYS2前者可能会出一些奇奇怪怪的路径问题后者则容易在 LAPACK 链接时找不到库浪费时间。3. 源码获取与编译全过程从 tar.gz 到可执行文件3.1 源码下载与目录结构认识TDEFNODE 的源码可以从项目官网的“Downloads”页面直接获取通常是一个类似tdefnode_x.x.tar.gz的压缩包。用 wget 下载下来后执行解压wget 官网下载链接 tar -zxvf tdefnode_x.x.tar.gz cd tdefnode_x.x这时候你用ls命令看一下目录会发现一个典型的 Fortran 项目结构。以我解压的版本为例LICENSE Makefile README.txt doc/ # 使用文档和 manual src/ # 源码文件 examples/ # 官方示例数据再进去看一眼 src 目录ls src/你会看到一堆.f90、.f结尾的文件比如tdefnode.f90是主程序和各个子模块还有sacio.f90、lapack_wrapper.f90这类功能模块。3.2 编译前必须修改的 Makefile 参数开编之前一定要打开 Makefile 看一眼不做这一步直接 make 的人十有八九会在链接阶段卡住。用 vim 或者 nano 打开 Makefilevim Makefile重点关注下面几个变量FC gfortran FFLAGS -O3 -fopenmp LAPACK_LIBS -llapack -lblasFC是编译器名如果系统默认 gfortran 版本不对这里可以改成gfortran-10。FFLAGS是编译选项-O3是最高级优化-fopenmp是开启 OpenMP 并行支持TDEFNODE 支持多核并行就是靠这个选项。LAPACK_LIBS是链接 LAPACK 库的指令-llapack -lblas的意思是链接名为 liblapack 和 libblas 的动态库。如果你的 LAPACK 安装路径比较特殊比如装在/usr/local/lib那么就需要在 Makefile 里添加链接路径LDFLAGS -L/usr/local/lib怎么确认 LAPACK 到底装在哪执行find /usr -name liblapack* 2/dev/null找到liblapack.so或liblapack.dylib的所在目录把路径填进LDFLAGS。这里有个更容易踩的细节如果是用 apt 安装的通常库文件在/usr/lib/x86_64-linux-gnu/Debian 系的这个路径经常被忽视导致链接不到。所以如果默认参数编译报错找不到 lapack首选检查这里。3.3 开始编译实际执行与日志观察确认 Makefile 没问题后在根目录执行make你会看到屏幕上刷过一大串编译信息这个过程一般需要 1 到 5 分钟取决于你机器的性能期间 gfortran 会把 .f90 文件逐行编译成目标文件.o最后再把它们链接成可执行文件。编译正常结束时根目录下会生成一个名为tdefnode的可执行文件。验证一下ls -l tdefnode file tdefnodefile命令会告诉你这是一个 ELF 可执行文件。如果这一步成功了恭喜你最难的编译阶段就过去了。如果中间报错了别慌跳到第 6 节我整理了几乎所有的典型报错和对应解法。3.4 编译成功后的标准动作编译成功后有几个推荐的标准动作能帮你少走不少弯路。先把源码里的示例数据拷贝到一个工作目录单独存放mkdir ~/tdefnode_work cp -r examples/* ~/tdefnode_work/ cd ~/tdefnode_work然后把可执行文件的路径加入环境变量或者干脆创建个软链接方便在任意目录下调用sudo ln -s /path/to/tdefnode/tdefnode /usr/local/bin/tdefnode这样下次你想运行 TDEFNODE 时直接tdefnode就能启动不用每次都在源码目录里干活。养成好习惯源码目录保持原样工作目录单独建这样调试起来不混乱。4. 输入文件格式解析官方示例跑通的关键一步4.1 控制文件.in 文件的结构编译通过只是第一关真正让人头大的是准备输入文件。TDEFNODE 的运行需要至少一个控制文件通常以.in结尾它告诉程序数据文件在哪儿、要反演哪些参数、迭代多少次、输出文件前缀等。我建议直接用官方 examples 目录里自带的输入文件来跑先别急着用自己的数据。打开官方示例的控制文件它长这样tdefnode.in文件内容是按固定格式排列的文本核心是几个命令块stations_file stations.dat fault_file faults.dat grid_file grid.dat out_prefix example_out ...stations_file指定站点文件路径。fault_file指定断层文件路径。grid_file指定网格文件路径网格用于计算应变率等衍生量。out_prefix所有输出文件的前缀比如example_out.def、example_out.sum。控制文件还有一些参数是控制反演流程的比如niter 1000 temp_start 100.0 temp_end 1.0niter是迭代次数temp_start和temp_end是模拟退火算法的起止温度。模拟退火是 TDEFNODE 使用的全局优化算法温度从高到低逐渐冷却类似金属退火过程。这个类比很贴切高温时系统接受较差解的概率高用于大范围搜索低温时逐渐收敛到最优区域。第一次跑示例不需要动这些参数用默认值就行。4.2 站点文件格式经纬度、速度与误差站点文件用来提供 GPS 测站数据。每一行通常对应一个测站格式大约是这个样子站名 经度 纬度 东向速度 北向速度 东向误差 北向误差 块体编号 S001 116.32 39.95 23.12 -5.34 1.23 1.45 1 S002 117.05 40.11 24.56 -4.98 1.10 1.30 1这里有几个容易搞混的点经纬度用的是十进制度数不是度分秒。如果你的数据是度分秒格式记得先转换不然位置全偏到海里去。速度单位通常是毫米/年mm/yr注意和你的原始数据统一。误差也必须是毫米/年而且要有东向和北向两个分量程序计算卡方时会用到。块体编号是整数告诉你这个站点属于哪个块体用于块体旋转拟合。有一个细节很容易忽略TDEFNODE 对站点文件的行格式非常敏感列与列之间可以用空格或者制表符分隔但千万不要混用而且像 S001 这种字符串不要带特殊符号否则读取时会直接中断。4.3 断层文件与网格文件断层文件定义了断层几何参数。和站点文件不同断层文件不是直接给一张地表断层线的经纬度列表而是定义断层分段segment的端点位置、滑动力学参数、闭锁深度等。以官方示例为例断层文件的格式大致是fault_name segment_id x1 y1 x2 y2 depth lock_depth slip_rate ...这里x1 y1和x2 y2是断层段的两个端点坐标depth是断层底部深度lock_depth是闭锁深度slip_rate是滑动速率。这个文件头信息很严格少一列或多一列都会导致读取错位然后计算结果全部变成 NaN。网格文件则是你在研究区域内设置的网格点TDEFNODE 会用这些点来计算应变率或绘图结果的插值。网格文件格式通常是每行一个点经度 纬度 116.0 39.0 116.5 39.0 ...网格密度取决于你的研究区域大小和计算力示例里默认的密度已经足够验证整个流程了。如果网格文件缺失有些功能会跑不了所以第一次别删它。4.4 修改示例文件改成自己的项目名理解了文件格式后我先做了一个最简单的修改把控制文件中的out_prefix改成我自己项目的名字比如test_run。同时把站点文件、断层文件的路径确认了一遍。查看一下文件是否存在路径是否对ls -l stations.dat faults.dat grid.dat确认没问题后就可以运行了。5. 正式跑通第一个案例5.1 运行命令与执行过程在示例目录下执行tdefnode tdefnode.in如果用的是绝对路径就写/path/to/tdefnode/tdefnode tdefnode.in这时终端会打出一串运行日志包括当前迭代步数、目标函数值、每个块体的旋转参数估计等。如果配置了 OpenMP还会显示正在使用多少个线程。这一步比较耗时示例数据规模不大一般几分钟内就能跑完。如果数据量很大去泡杯咖啡接着等。跑完后你会发现工作目录下多了一批文件文件名前缀都是你设置的那个out_prefix比如test_run.lonlat test_run.def test_run.out test_run.sum test_run.vel5.2 关键输出文件解读这一步是让我看懂 TDEFNODE 在干嘛的关键也是我觉得全网教程里讲得不够细的地方。test_run.out这是最主要的输出文件里面记录了拟合的最终参数包括所有块体的欧拉极、断层滑动速率、闭锁深度、应变率等。块体欧拉极一般用经度、纬度和旋转速率deg/Myr表示你可以理解成一个想象中的旋转轴和角速度用来描述整个块体的刚性运动。test_run.vel这是站点速度的预测值和观测值对比文件。每一行有测站名、观测速度、预测速度、残差。这个文件是检验拟合质量的第一手资料。理想情况下残差应该很小而且没有系统性偏差。test_run.sum这是拟合优度统计文件包括加权均方根残差WRMS、卡方值Chi-square等。WRMS 越小说明模型和数据越接近。这个文件是判断反演是否收敛的根本依据。test_run.def断层滑动分布的定义文件记录了每个断层段的深度、滑动矢量、闭锁程度等信息。如果你打开这些文件看到的数值都比较合理没有出现满屏的 NaN 或者无穷大说明你的安装和输入文件格式已经基本没问题了。5.3 怎么判断结果“对不对劲”判断拟合质量的几个经验法则WRMS 值如果残差在 1~3 mm/yr 量级说明拟合得不错。如果残差大得离谱比如几十 mm/yr要么是数据本身有问题要么是模型参数初值给得不好。残差方向如果所有站的残差都偏向同一个方向说明模型存在系统性偏差比如某个块体旋转速率被高估了。断层的滑动速率检查一下得到的滑动速率是否符合该地区的地质背景常识。比如一条倾滑断层如果反演出纯走滑那你就要回头看看断层几何参数是不是写错了。第一次跑通不代表你理解了结果但至少证明从安装到数据准备这条路是通的。接下来你才能放心地把自己的数据填进去去解决实际科学问题。6. 常见问题与避坑实录这些坑我替你踩过了6.1 编译阶段报错速查表报错信息原因解决方法gfortran: command not found编译器没装apt install gfortranundefined reference todgesv_LAPACK 链接失败检查 Makefile 的 LAPACK_LIBS 和 LDFLAGSError: Type mismatch老代码踩到新版语法检查在 FFLAGS 里加-stdlegacy -fallow-argument-mismatchmake: *** No rule to make targetMakefile 路径配置问题检查在源码根目录执行的 make不是 src 目录Error: Cannot open module file ‘xxx.mod’编译依赖顺序问题先make clean再重新执行make我在编译时遇到最多的就是undefined reference to dgesv_。dgesv 是 LAPACK 里解线性方程组的函数很多 Fortran 程序都依赖它。这个报错说明你的编译器找到了代码但链接器找不到 LAPACK 库。那时候我突然想起当时安装 Dev 包时偷懒用了最小安装后来补上liblapack-dev再执行make clean make问题就解决了。6.2 运行阶段常见问题编译过了运行时报错同样让人抓狂。这里把运行时报错也整理成一个速查表报错信息原因解决方法Error opening file输入文件路径不对检查控制文件里的文件名和当前目录NaN铺满所有输出站点文件或断层文件格式错位用head -5检查每个输入文件确认列数和顺序Segmentation fault数组越界或文件行数不足检查输入文件行数、块体编号是否越界程序卡住不退出反演参数设置不合理模拟退火温度过高调低temp_start或增加niter所有站点的速度预测一模一样块体编号没对上检查站点文件的块体编号是否超出块体数量范围特别注意这个 NaN 的问题。我首次拿到自己数据时输出文件里整个 vel 文件全是 NaN排查了一下午才发现是断层文件里多打了一个制表符。TDEFNODE 的读取程序是非常机械的格式上一个字符不对后面的数据全都会错位然后计算得到 NaN。这个问题排查起来非常反直觉因为你不知道是哪一行哪一列错位了只能一行行对齐。6.3 实战排查流程演示假如你运行tdefnode tdefnode.in后输出文件全是 NaN按这个顺序排查第一步看输入文件字符编码。如果文件是 Windows 下编辑的很可能带了回车符\rTDEFNODE 读取时会把它当成文件名的一部分直接报错打不开文件。用sed -i s/\r$// stations.dat faults.dat grid.dat清理一下。第二步用cat -n看文件每行末尾是否有多余空格或制表符。命令行检查比编辑器直观。第三步检查控制文件里各输入文件路径是否写错大小写是否敏感。第四步将断层文件的行数与站点文件的行数简化到最小比如只保留一条断层和 10 个站点看是否能跑通。如果简化版能跑通说明问题出在某个特定数据行上二分法逐段排查。实测下来90% 以上的运行报错都能通过这四步定位。6.4 其他实战技巧工作文件尽量不要放在中文路径或带空格的路径下Fortran 程序内部处理文件名可能有问题我吃过这个亏路径里有空格直接读不出来。修改完输入文件后先跑一次小规模快速测试。建议把网格文件里的网格点删掉一部分只留少量点这样跑一轮只要几秒钟能快速检验格式是否正确然后再跑完整数据。备份/lib里的动态库依赖。tdefnode 可执行文件如果要在其他机器上复制使用需要目标机器也安装匹配版本的 LAPACK否则会报error while loading shared libraries这是我在换服务器时踩到的额外坑。尽可能用静态链接编译。Makefile 里LDFLAGS加上-static把必要的库都打进可执行文件里这样换机器运行就不用再纠结动态库问题了。不过要注意-static和-fopenmp在某些系统上可能会有兼容性问题需要实际测试。7. 进阶使用与二次开发建议7.1 多核并行与效率调优TDEFNODE 编译时如果带了-fopenmp运行时会自动使用所有能用的 CPU 核但你也可以用环境变量来控制线程数。export OMP_NUM_THREADS4我建议第一次跑的时候不要把所有核都用满先限制到 2 到 4 个核试试因为反演程序不是简单的并行任务线程越多部分变量同步的开销也越大。实测中我发现某些环节比如模拟退火里的随机扰动并不适合大规模并行线程数超过 8 后加速比下降明显。如果反演参数多、数据量大建议先跑短迭代观察目标函数收敛趋势再放长跑不要一上来就十万次迭代否则你没法判断参数初值是否合理。7.2 如何把自己的数据接入 TDEFNODE这是把示例跑通后必然会遇到的事情。核心思路先用官方示例把整个流程跑通然后逐文件替换。替换站点文件把你的 GPS 站点数据整理成示例站点文件一样的格式块体编号按你的研究区划分填好。替换断层文件断层几何参数可以从地震地质资料里提取但要注意断层分段的端点坐标必须和你的站点文件坐标系统一致别一个用 WGS84 一个用别的椭球。调整控制文件把块体数量、初始值、迭代次数相应改掉。替换后先跑一个低迭代次数的快速测试确认程序不报错再慢慢加次数跑正式反演。很多新手一上来就指望一次跑成功结果出了问题也不知道是自己数据格式的问题还是参数初值的问题反而多花时间。7.3 和 GMT 绘图工具的联动TDEFNODE 的输出文件格式是文本表格很方便用 GMTGeneric Mapping Tools做可视化。用 GMT 把站点速度观测值和预测值画出来可以直观检验拟合效果。推荐的做法gmt begin velocity_map png gmt coast -R115/119/38/42 -JM15c -Ba -W0.5p gmt plot test_run.vel -Sc0.1c -W0.5p gmt end这里test_run.vel包含观测和预测速度你可以按列提取分别画箭头用不同颜色表示观测和预测叠加在断层图上。可视化这件事对结果的判断帮助极大肉眼比任何统计指标都容易发现问题。7.4 源码二次修改的切入点如果你要做更深度的定制比如修改反演算法或者增加新的约束条件可以从src/tdefnode.f90的主程序作为起点找到对应的模块修改。建议先通读doc/目录下的 manual 文档理解每个子函数的作用再动手改。二次开发时配套建议用gfortran -g -O0编译一个调试版本这个版本没有优化适合用 GDB 逐行打断点排查问题。8. 写在最后的个人经验TDEFNODE 这套工具说难装也难装说简单也简单核心就三件事装对编译器、配好 LAPACK、老老实实按格式准备输入文件。我头一回编译时一天之内从对 Fortran 一无所知到在终端里把几十条报错一条条啃下来整个过程非常磨人但回头看其实思路很清晰先查 Makefile 里的变量再查依赖库路径最后查输入文件格式一步步排查总能解决。最后分享一个小技巧TDEFNODE 官方文档doc 目录里的 manual写得非常详实但很多人习惯跳过它直接编译。我建议你在等编译的过程中把 manual 里关于输入文件格式的章节从头到尾读一遍大概花掉一个多小时但能帮你少走一整天的弯路。我就是没读文档直接上结果在输入文件格式上卡了最久绕了很大一圈才明白。源码编译类工具就是这样一个门槛跨过去之后你会觉得它的数据组织和反演思路都极其干净直接真正把数据丢进去得到合理结果的那一刻前面所有的折腾都值了。
