Creates a Finite Difference Grid over a Two-Dimensional Rectangular Domain

Share:

Description

Creates a finite difference grid over a rectangular two-dimensional model domain starting from two separate one-dimensional grids (as created by setup.grid.1D).

Usage

1
setup.grid.2D(x.grid = NULL, y.grid = NULL)

Arguments

x.grid

list containing the one-dimensional grid in the vertical direction - see setup.grid.1D for the structure of the list

y.grid

list containing the one-dimensional grid in the horizontal direction - see setup.grid.1D for the structure of the list

Value

a list of type grid.2D containing:

x.up

position of the upstream interface in x-direction (i.e. if x is vertical, the upper boundary); one value

x.down

position of the downstream interface in x-direction (i.e. if x is vertical, the lower boundary); one value

x.mid

position of the middle of the grid cells in x-direction; vector of length x.N

x.int

position of the interfaces of the grid cells in x-direction; vector of length x.N+1

dx

distance between adjacent cell interfaces in x-direction (thickness of grid cells); vector of length x.N

dx.aux

auxiliary vector containing the distance between adjacent cell centers; at the upstream and downstream boundary calculated as (x[1]-x.up) and (x.down-x[x.N]) respectively; vector of length x.N+1

x.N

total number of grid cells in the x direction; one value

y.left

position of the upstream interface in y-direction (i.e. if y us the horizontal, the left boundary); one value

y.right

position of the downstream interface in y-direction (i.e. if y us the horizontal, the right boundary); one value

y.mid

position of the middle of the grid cells in y-direction; vector of length y.N

y.int

position of the interfaces of the grid cells in y-direction; vector of length y.N+1

dy

distance between adjacent cell interfaces in y-direction (thickness of grid cells); vector of length y.N

dy.aux

auxiliary vector containing the distance between adjacent cell centers; at the upstream and downstream boundary calculated as (y[1]-y.up) and (y.down-y[y.N]) respectively; vector of length y.N+1

y.N

total number of grid cells in the y direction; one value

Author(s)

Filip Meysman <filip.meysman@nioz.nl>, Karline Soetaert <karline.soetaert@nioz.nl>

See Also

tran.2D, for a discretisation of the general transport equation in 2-D

setup.grid.1D, for the creation of grids in 1-D

setup.prop.2D for defining properties on the 2-D grid.

Examples

1
2
3
4
# test of the setup.grid.2D functionality
x.grid  <- setup.grid.1D(x.up = 0, L = 10, N = 5)
y.grid  <- setup.grid.1D(x.up = 0, L = 20, N = 10)
(grid2D <- setup.grid.2D(x.grid, y.grid))