Changes between Version 145 and Version 146 of doc/app/initialization_parameters
- Timestamp:
- Jan 6, 2011 1:22:31 PM (14 years ago)
Legend:
- Unmodified
- Added
- Removed
- Modified
-
doc/app/initialization_parameters
v145 v146 933 933 }}} 934 934 {{{#!td style="vertical-align:top" 935 1.0E- 6935 1.0E-4 936 936 }}} 937 937 {{{#!td … … 954 954 The user can choose between the following schemes:\\\\ 955 955 '' 'ws-scheme' ''\\\\ 956 956 957 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'. 957 958 … … 959 960 The computation of turbulent fluxes takes place inside the advection routines to get a statistical evaluation consistent to the numerical solution. 960 961 961 '''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 u-component of velocity: CALL exchange_horiz( pt , nbgp ). \\\\ 962 '''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 temerature: CALL exchange_horiz( pt , nbgp ). \\\\ 963 962 964 '' 'pw-scheme' ''\\\\ 963 965 The scheme of Piascek and Williams (1970, J. Comp. Phys., 6, 392-405) with central differences in the form C3 is used.