3.3.1. Mesh file

The mesh defines the model region and has the following structure:

../../_images/mesh3D.png
  • \(NE\): Number of cells in the East direction.

  • \(NN\): Number of cells in the North direction

  • \(NZ\): Number of cells in the vertical direction

  • \(E_o, N_o, Z_o\): Coordinates, in meters, of the southwest top corner, specified in (Easting, Northing, Elevation). The elevation can be relative to a reference elevation other than the sea level, but it needs to be consistent with the elevation used to specify the locations, observations, and topography files.

  • \(\Delta E_n\): \(n^{th}\) cell width in the easting direction (ordered W to E).

  • \(\Delta N_n\): \(n^{th}\) cell width in the northing direction (ordered S to N).

  • \(\Delta Z_n\): \(n^{th}\) cell thickness (ordered top to bottom).

The mesh can be designed in accordance with the area of interest and the spacing of the data available in the area. In general, the mesh consists of a core region which is directly beneath the area of available data, and a padding zone surrounding this core mesh. Within the core mesh, the size of the cells should be comparable with the spacing of the data. There is no restriction on the relative position of data location and nodal points in horizontal direction. The cell width in this area is usually uniform. Beyond the core region, the mesh should be padded with cells that increase (typically no more than 40% of the previous length).

The vertical position of the mesh is specified in elevation. This is to accommodate the inversion of a data set acquired over a topographic surface. When there is strong topographic relief, which the user wishes to incorporate it into the inversion, special care should be taken to design the mesh. A conceptually simple approach is first to design a rectangular mesh whose top (specified by \(Z_o\)) is just below the highest elevation point, and then to strip off cells that are above the topographic surface. This is the approach taken in . The number of cells to be stripped off in each column is determined by the user-supplied topography file. Only the remaining cells will be used in the forward modelling or included in the inversion as model parameters.

3.3.1.1. Example

This example shows a mesh that consists of 26 cells in easting, 27 cells in the northing, and 23 cells in the vertical directions. The top of the mesh is located at 0 m of elevation and the southwest corner is at -350 m easting and -400 m northing. The cells in the core portion of the mesh are all 50 m \(\times\) 50 m \(\times\) 25 m. There are three cells in the padding zone in every direction except the top of the core mesh.

../../_images/mesh3Dex.png