qq_45174849的博客

gprMax使用的是将一个.in文件导入在gprMax的运行环境下进行仿真的，首先要提供一个.in文。在输入.in文件中，指令都是通过哈希字符(#)作为开头的，即表示传递给gprMax的命令的开头。件作为输入，这个.in文件包括运行GPR模型所需的所有必要信息，.in文件是一个ASCII文本文件，哈希字符(#)，如果一行以哈希字符(#)开头，程序会默认是一个有效的命令，如果没有哈希字符(#)f：是一个实数，可以用小数，例如1.5，可以用科学记数法，例如15e-1或0.15e1.如果命令的名称不正确，

在gprMax中大多数命令都是可选的，但有一些基本命令是构建任何模型所必需的，例如模型的大小，gprMax将终止执行并发出错误消息。

其中f1可以取0 < f1≤1，那么实际使用的时间步长是f1 *∆𝑡，其中∆𝑡是使用CFL条件中的等式计。使用Python代码块中创建命令更容易，其中gprMax中有内置模块，包含一些最常用的函数形式。file1可以是在与输入文件相同的目录中包含命令的文件名称，也可以是包含命令的文件的完整。功能：从指定文件中包含的命令插入到输入文件中放置#include_file命令的位置。其中c1可以是y(是)或n(否)，打开或关闭显示屏幕上的信息量，默认值为y。功能：更改gprMax使用的时间步长∆𝑡的值。

式中，∆𝜖𝑟𝑝=𝜖𝑟𝑠𝑝−𝜖𝑟∞，其中，𝜖𝑟𝑠𝑝为极点的零频率相对介电常数，𝜖𝑟∞为无穷频下的相对介电常数，𝜏𝑝为极点弛豫时间。式中，∆𝜖𝑟𝑝 = 𝜖𝑟𝑠𝑝 −𝜖𝑟∞，𝜖𝑟𝑠𝑝其为极点的零频率相对介电常数，𝜖𝑟∞其为无限频率时的相对介电。f1：为第一个Debye poles极零频率相对介电常数与无限频率相对介电常数之差，即∆𝜖𝑟𝑝=𝜖𝑟𝑠𝑝1−𝜖𝑟∞。f3：为第二个Debye poles极零频率相对介电常数与无限频率相对介电常数之差，即∆𝜖𝑟𝑝=𝜖𝑟𝑠𝑝2−𝜖𝑟∞。与极点频率和松弛时间相关的时间值应始。

对象构造命令在输入文件是按照命令行的顺序进行处理的，在进行对象构建的时候，可以把模型中的空间当成一个画布，其中引入了物体，并且可以将一个物体覆盖在另一个物体上，从而覆盖其属性，以产生所需的几何形状，对象构造命令可创建复杂的形状和配置。f4、f5、f6：为#fractal_box上表面的右上角(x,y,z)坐标。例如，指定一个其轴在y方向，长度为0.7 m，半径为100mm，并且是一个完美的电导体的圆柱体。例如，#sand: 3 0 10 0是材料，文件中的第一个材料，它将与HDF5文件中的任何为零的整数。

例如，一个y方向，内阻为50欧姆，幅度为5，中心频率为1.2 GHz高斯波形的电压源。例如，一个z向传输线源，其电阻为75欧姆，幅度为5，中心频率为1.2 GHz高斯波形。式中𝐽𝑠为电流密度，𝐼为电流，∆𝑙为无穷小电偶极子的长度，∆x，∆y，∆z为网格的空间分辨。功能：指定一个ASCII文件，该文件包含幅度值列，这些幅度值列指定可以与模型中的源一起。file1：为与输入文件在同一目录下包含指定波形的文件的名称，或者file1可以是包含指定波形的文。如果它的电阻为零，即指定电场分量的时间变化是规定的，

gprMax的默认行为是使用一阶CFS PML，该PML在模型域的六面各有10个单元格的厚度。这个厚度可以通过以下命令进行更改。