diff --git a/docs/stretched_grids.md b/docs/stretched_grids.md index 0cf11bf4..90c32988 100644 --- a/docs/stretched_grids.md +++ b/docs/stretched_grids.md @@ -9,13 +9,17 @@ We also store the temperature field `temp` (and if used the humidity field `qhum The velocity components are then discretised on their corresponding cell edges. The boundaries of the domain are specified to lie on the grid of the cell edges, so for example (if the x-domain size is set as `alx3=1`) + $$ x_c[1] = 0, \qquad x_c[n_x] = 1 . $$ + The cell mid-points are always defined as + $$ x_m[k] = \frac{x_c[k] + x_c[k+1]}{2} , $$ + whether the grid is stretched or uniform. The grids in $y$ and $z$ must always be uniformly spaced so that Fourier transforms can be performed in the pressure solver, but the grids in $x$ can be non-uniform. Their stretching is defined by the input parameters `str3` and `istr3`.