eigenv
- 简介:
计算面波频散曲线
语法
grt eigenv -Mmodel -Cpath -SR|L -Ff1[/f2/df][+p] [ -N[max_order] ] [ -Xfreq[+ccmin/cmax][+iiref] ] [ -T[+ttol][+ccgap][+rthrd][+vvgap][+udc] ] [ -Pnthreads ] [ -s ] [ -h ]
描述
通过搜索由广义反射透射系数矩阵法定义的久期函数的零点来计算层状介质的面波相速度频散曲线，
结果保存为 NetCDF 网格格式。
如果使用 ncdump -h 命令查看该文件的基本信息，结果类似于：
netcdf phase_R {
dimensions:
freq = 500 ;
mode = 54 ;
variables:
double freq(freq) ;
int mode(mode) ;
int cnum(freq) ;
double c(freq, mode) ;
c:_FillValue = 0. ;
int ciref(freq, mode) ;
ciref:_FillValue = -1 ;
// global attributes:
:command = "grt eigenv -Mmod2 -SR -F0/5/0.01 -N -Cphase_R.nc" ;
:model = "mod2" ;
:isRayl = 1 ;
}
维度（dimensions）
freq- 频率点数mode- 阶数点数
变量（variables）
freq- 频率数组 (Hz)mode- 阶数数组cnum- 不同频率下的相速度点数c- 不同频率不同阶的相速度 (km/s)ciref- 搜索对应相速度使用的久期函数层位
全局属性（global attributes）
command- 生成该文件的命令model- 模型路径isRayl- 是否为Rayleigh波结果
必选选项
- -Mmodel
半空间层状模型文件的路径。模型格式如下:
Thickness(km) Vp(km/s) Vs(km/s) Rho(g/cm^3) [Qp] [Qs]
后两列的 Qp, Qs 可省略。
如果首列首行的值为正数，则首列表示 每层厚度 (km) ，最后一层为半空间，其厚度设置为 0 即可（不会被使用）。
如果首列首行的值为 0，则首列表示 每层层顶深度 (km) 。
如果某层设置 Vs = 0，则该层为液体层。
- -Cpath
eigenv 模块计算的相速度频散结果保存文件路径。
- -SR|L
计算 Rayleigh 波频散还是 Love 波频散。
- -Ff1[/f2/df][+p]
设置要计算的频率，有以下几种用法：
-Ff1 - 单个频率 f1 ，单位 Hz
-Ff1/f2/df - 从频率 f1 到 f2 ，间距 df 的一系列等距频率，单位 Hz
在以上用法后加上 +p ，则 f1,f2,df 均变为周期，单位 secs。 例如 -F0/10/2+p 表示计算周期 2s, 4s, 6s, 8s, 10s 处的频散。
注： 在使用 f1/f2/df 定义等距频率点时，为了方便，用户可以设置 f1=0， 此时程序会自动跳到下一个频率点 (df)。
可选选项
- -N[max_order]
设置计算的最大阶数，0 对应基阶。
如果不使用 -N ，默认仅计算基阶，相当于设置 -N0
如果省略 max_order ，即仅给定 -N ，则尽可能计算所有阶。
- -Xfreq[+ccmin/cmax][+iiref]
该选项常用于 debug，可输出频率为 freq Hz，相速度范围在 [cmin, cmax] 之间的第 iref 层的久期函数。 此时 -F 和 -C 不必再设置，久期函数结果将输出到标准输出，包含三列，分别为相速度、久期函数实部、久期函数虚部。 采样间隔受搜根方法控制，详见 -T 。
- -T[+ttol][+ccgap][+rthrd][+vvgap][+udc]
设置搜根过程中的一些参数，绝大多数情况下默认参数即可。 默认情况下， eigenv 模块使用自适应采样久期函数的方式进行搜根 (Zhu et al., 2026)， 其依赖参数：
+ttol - 自适应策略的收敛容差 \(10^{-\it \text{tol}}\) ， tol 默认为 3。
作为测试，也可使用等相速度间隔的方式进行搜根，其依赖参数：
+udc - 等间隔搜索的相速度间隔，单位 km/s 。
搜根过程也依赖以下参数：
+ccgap - 同频率下相邻零点的最小相速度间隔 \(10^{-\it \text{cgap}}\) ， cgap 默认为 7。
+rthrd - 识别为零点的久期函数振幅阈值 \(10^{-\it \text{thrd}}\) ， thrd 默认为 4。
+vvgap - 零点位置与模型速度的最小间隔 \(10^{-\it \text{vgap}}\) ， vgap 默认为 7。
- -Pnthreads
多线程数，默认使用全部线程。
eigenv 模块会自动并行搜索不同频率的频散结果，提高计算效率。
- -s
静默输出，不在终端打印参数和进度条。
- -h
打印帮助文档。