 Therefore, a transformation of the vertical coordinate to a terrain-following one is performed.Hence, the originally irregularly bounded physical domain is mapped onto one consisting of unit cubes.The discretized equations are solved numerically on a staggered grid, i.e.