The discretized equations are solved numerically on a staggered grid, i.e. the scalar quantities \rho, p and \theta are defined at the cell centre while the velocity components u, v and w are defined at the centre of the appropriate interfaces.
Temporal discretization of the prognostic equations is based on the explicit second order Adams-Bashforth scheme.