注意,以下的说明都是根据SEM2D自带的xmeshfem2D实现的网格化

此页面对应程序中interfaces文件的修改

设置介质的界面

在二维的SEM中,一般会选取X-Z平面作为研究区域,所以对研究区域作为分层,决定分层的这一界面实际上是一条线。

1. 水平界面

由于两点在二维平面内可以确定的生成一条直线,所以确定水平界面只需要给出界面的起点和终点的X,Z坐标

对应一个两层的水平分层的介质来说,相应的界面示例文件:interfaces_flat.dat为(中文仅是说明,不需要包括在此文件内):

# number of interfaces
 3
# 划分出两层,需要有3个界面
# for each interface below, we give the number of points and then x,z for each point
# interface number 1 (bottom of the mesh)
 2 
#决定这一界面的点的个数
 0 0 
#起点的X,Z坐标,单位是米
 5000 0 
# 终点的X,Z坐标,单位是米
# interface number 2
 2
 0 1000
 5000 1000
# interface number 3 (top of the mesh)
 2
 0 3000
 5000 3000
# for each layer, we give the number of spectral elements in the vertical direction
# 对于每一层,给出垂向上计算单元的个数,这里对于垂向上的1000m会出现20个计算单元,每个单元Z方向的长度为50m
# layer number 1 (bottom layer)
 20
# layer number 2 (top layer)
 40

2. 弯曲界面

在二维平面内确定一条曲线需要多个点,点数越多就可以模拟越复杂的界面,与确定水平界面同理,需要给出每个点的X,Z坐标。相应的界面示例文件为EXAMPLES/simple_topography_and_also_a_simple_fluid_layer/DATA中的interfaces_simple_topo_curved.dat