Changes between Version 254 and Version 255 of doc/app/initialization_parameters
- Timestamp:
- Mar 6, 2015 12:56:32 PM (10 years ago)
Legend:
- Unmodified
- Added
- Removed
- Modified
-
doc/app/initialization_parameters
v254 v255 1082 1082 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 scheme implies a small numerical dissipation that stabilizes the solution. To assure a stable numerical solution the time integration has to be carried out with [#timestep_scheme timestep_scheme] = '' 'runge-kutta-3' '' . This scheme is based on a formulation 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] = '' .T. '' has to be used. 1083 1083 '''Note''': Due to the larger stencil of this scheme vertical grid stretching should be handled with care. 1084 The computation of turbulent fluxes takes place inside the advection routines to get a statistical evaluation consistent to the numerical solution. 1085 1086 '''Important''': The number of ghost layers for 2d and 3d arrays changed. This affects also the user interface. Please adapt the allocation of 2d and 3d arrays in your user interface like [../userint/output#Allocate here]. Furthermore the exchange of ghost layers for 3d variables changed, so calls of exchange_horiz in the user interface have to be modified. Here an example for the potential temperature: CALL exchange_horiz( pt , nbgp ). \\\\ 1084 The computation of turbulent fluxes takes place inside the advection routines to get a statistical evaluation consistent to the numerical solution. \\\\ 1087 1085 '' 'ws-scheme-mono' ''\\\\ 1088 Monotone scheme combining the 5th order upwind scheme of Wicker and Skamarock with 1st order upwind scheme using a monotonic limiter. Please note that this monotone scheme can are very diffusive. Even though this scheme prevents non-physical numerical oscillations at large gradients, the numerical dispersion might be very large, particularly at large gradients. Moreover, it should be noted that in contrast to 'ws-scheme', the computation of turbulent fluxes '''does not''' take place inside the advection routines, resulting in possible kinks in the near-surface profiles of vertical fluxes.\\\\1086 Monotone scheme combining the 5th order upwind scheme of Wicker and Skamarock with 1st order upwind scheme using a monotonic limiter. Please note that this monotone scheme can be very diffusive. Even though the scheme prevents non-physical numerical oscillations in case of strong gradients, the numerical dispersion might be very large, particularly in areas where strong gradients occur. Moreover, it should be noted that in contrast to 'ws-scheme', the computation of turbulent fluxes '''does not''' take place inside the advection routines, resulting in possible kinks in the near-surface profiles of vertical fluxes.\\\\ 1089 1087 '' 'pw-scheme' ''\\\\ 1090 1088 The scheme of Piacsek and Williams (1970, J. Comp. Phys., 6, 392-405) with central differences in the form C3 is used.\\\\