These parameters set the size and mesh spacing of the computational domain. They also define boundary conditions.
- Ndim(integer): The number of spatial dimensions in the problem of interest.
- Size_x[idim](real vector): The size of the computational domain. There is one entry for each dimension.
- Esize_x[idim](real vector): The mesh spacing. There is one entry for each dimension (they do not have to be the same).
- Type_bc[0,1][idim](integer array): There are three lines of input associated with Type_bc. Each line should contain two integers.
The different lines designate the boundary conditions in the
different dimensions (some may not be used). The two entries on a given line specify boundary conditions for the two sides
of the domain in a given dimension. The options are:
- -1: IN_WALL: semi-infinite surface : set density=0 when integrating beyond boundary.
- 0: IN_BULK: constant bulk fluid: set density equal to the bulk density when integrating beyond boundary.
- 1: PERIODIC: periodic boundary: set density at node N+k equal to the density at node k (Nodes run from 0 to N).
- 2: REFLECT: reflective boundary: set density at node N+k equal to the density at node N-k.
- 3: LAST_NODE: continuation boundary: set density at node N+k equal to the density at node N.