Mesh Parameters
These parameters set the size and mesh spacing of the computational domain. They also define boundary conditions.
 Ndim(integer): The number of spatial dimensions to be solved numerically 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 need not 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: semiinfinite 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 Nk.
 3: LAST_NODE: continuation boundary: set density at node N+k equal to the density at node N.
 4: LAST_NODE_RESTART: continuation boundary: set density at node N+k equal to the density at node N, but keep solution values obtained from a restart file for these nodes.
This can be useful for pinning a free interface away from some kind of surface or inclusion.
