GprMax2D ——一般命令

1. #title:
2. #domain:
3. #dx_dy:
4. #time_step_stability_factor:
5. #time_window:
6. #messages:
7. #number_of_media:
8. #nips_number:
9. #media_file:
10. #geometry_file:

#title:

使用#title命令，您可以包含模型的标题。此标题保存在输出文件中。该命令的语法是：

#title：str

这里的参数str可以包含空格字符以分隔单个单词，标题必须包含在一行中。

#domain:

命令#domain：用于指定模型的大小（以米为单位），它的语法是：

#domain：f1 f2

参数f1和f2分别是模型在x和y方向上的大小（以米为单位），例如：

#domain：1.0 1.5

将模型的大小设置为 1. 0 × 1 .5 米

#dx_dy:

该命令#dx_dy用于指定在空间离散化X和Y方向（即Δ X和Δ Y），该命令的语法是：

#dx_dy：f1 f2

其中F1是在空间步骤x方向（Δ X），F2是在空间步骤y方向（Δ Y），因此命令为

#dx_dy：0.1 0.1

将设置的离散化步骤Δ x和Δ y在每个方向上10厘米

离散化的选择步骤Δ x和Δ y对构建模型是非常重要的，此命令与#domain命令结合使用可确定将在模型中使用的单元格数，从而确定模型计算机内存的最大部分要求。例如，使用命令：

#domain：1.0 1.5

#dx_dy：0.1 0.1

模型中的单元格数量为 10 × 15 ，将上面改为

#domain：1.0 1.5

#dx_dy：0.01 0.01

模型中的单元格数量将增加到100 150，此外，空间离散控制的最大允许时间步长Δ t，利用该最大允许时间步长t达到所需的模拟时间窗口。Δx、Δy和Δt之间的关系：

其中c是光速。 在GprMax2D中，上式中使用相等来确定Δt，中可以明显看出，Δx和Δy的值较小将导致t的值较小，这也
意味着达到给定的模拟时间要进行更多的迭代，但要注意Δx，Δy和Δt的值越小，你的模型就越准确。

#time_step_stability_factor:

用命令time_step_stability_factor：您可以更改GprMax2D计算的时间步长t的值。但是，新值应在上述公式确定的允许范围内。如上所述，GprMax2D使用公式中的相等，那是最大允许时间步长。如果需要较小的时间步长则使用下述命令

#time_step_stability_factor：f1

参数f1可以取值0 <f1<=1 ，GprMax2D将使用时间步长f1·Δt（Δ t是在使用相等计算出的）

#time_window:

命令#time_window：用于指定所需的总模拟时间，它的语法是

#time_window：f1 或者 #time_window：i1

在第一种情况下，f1参数确定所需的模拟时间（以秒为单位）。例如，如果要模拟20纳秒的GPR轨迹使用

#time_window：20e-9

GprMax2D将执行必要的迭代次数，以达到所需的模拟时间。

或者，使用i1指定命令，则GprMax2D会将此值解释为指定的总迭代次数。因此命令

#time_window：100

意味着将执行100次迭代，迭代次数和总模拟时间窗口相关

 tw = Δt × Nit

其中TW是时间窗口（秒），Δ t 时间步长和 Nit 迭代次数

#number_of_media:

当您的模型需要使用超过10 种不同的介质时，必须使用命令#number_of_media。最初，GprMax2D分配空间来存储10个媒体的组成参数（free_space，pec和8个用户定义介质）。这个空间通常适用于大多数模型，但如果需要更多介质，则应发出命令

#number_of_media：i1

其中参数i1应大于10且至少等于您需要的不同介质的数量。即使是1000种不同的媒体分配空间也没有问题。但是，有价值的计算机存储器可能被浪费。如果GprMax2D需要更多空间来存储媒体参数，则会终止执行并发出相应的错误消息。

#media_file:

使用命令#media_file：您可以指定描述包含常用介质的构成参数的文件的位置和文件名。通过将这些参数存储在这样的文件中，您不必在要使用它们的每个模型的输入文件中指定它们。该命令的语法是

#media_file：file1

其中参数file1是媒体文件的文件名（必要时包括路径）。

#geometry_file:

使用命令#geometry_file：您可以指定一个文件，其中有关模型几何的信息以二进制格式存储。 此信息可用于创建模型的图像并检查其是否正确构造。该命令的语法是

#geometry_file：file1

参数file1是几何文件的文件名。例如命令

#geometry_file：model.geo

将指示GprMax2D存储有关文件模型的几何信息model.geo。

#messages:

在输入文件中使用命令#messages：可以部分控制运行时屏幕上显示的信息量。该命令的语法是：

#messages：c1

参数C1可以是 y或 n，它可以打开或关闭屏幕上的消息。默认值为n， 当消息打开时，GprMax2D将在屏幕上显示空间和时间值到单元坐标和迭代数（整数）值的转换。此信息可用于错误检查。

#nips_number:

只有当GprMax2D请求使用时，命令#nips_number：才应包含在输入文件中。该命令的语法是

#nips_number：i1

其中参数i1是GprMax2D建议您使用的整数。 此命令控制用于存储有关模型的重要信息的数组的大小。尽管可以在内部计算这些阵列的大小，但是这种计算考虑到所有可能所需空间，结果通常比大多数情况下实际需要的大得多。因此，使用较小的数字以节省计算机存储器。当GprMax2D检测到需要更多空间时，它会发出错误消息并提示您在输入文件中使用此命令，这将为分配留出更多空间。

终于把一般命令学习完了，通过这次学习，我对这些命令的理解更加深刻了，我想这也为我以后的建模打下了坚实的基础。