957 | 956 | The 5th order upwind scheme of Wicker and Skamarock (2002, Mon. Wea. Rev, 130, 2088-2097) is used. The dispersion error is much smaller than the dispersion error of 'pw-scheme'. The 5th order schemes implies a small numerical dissipation that stabilized the solution. To assure a stable numerical solution the time integration has to be carry out with [#timestep_scheme timestep_scheme] = 'runge-kutta-3' . This scheme is based on a formultion of the advection term in flux form, which requires a vanishing divergence of the flow field, else a stable numerical solution is not given. So [#call_psolver_at_all_substeps call_psolver_at_all_substeps] = 'True' has to be used. When [#psolver psolver] = 'multigrid' the [#residual_limit residual_limit] should be smaller than 1.0E-6 for a sufficient reduction of the divergences. For the moment the scheme can only be used in conjunction with [#topography topography] = 'flat'. So far [#bc_lr bc_lr] /= 'cyclic' or [#bc_ns bc_ns] /= 'cyclic' are not allowed with [#loop_optimization loop_optimization] = 'vector' in combination with 'ws-scheme'. |