Changes between Version 148 and Version 149 of doc/app/initialization_parameters
- Timestamp:
- Jan 12, 2011 12:22:26 PM (14 years ago)
Legend:
- Unmodified
- Added
- Removed
- Modified
-
doc/app/initialization_parameters
v148 v149 758 758 The user can choose between the following schemes:\\\\ 759 759 '' 'ws-scheme' ''\\\\ 760 The 5 th 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'.760 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. When [#psolver psolver] = '' 'multigrid' '' the [#residual_limit residual_limit] should be smaller than 1.0E-6 for a sufficient reduction of the divergences. So far 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' ''. 761 761 762 762 '''Note''': Due to the larger stencil of this scheme vertical grid stretching should be handled with care. 763 763 The computation of turbulent fluxes takes place inside the advection routines to get a statistical evaluation consistent to the numerical solution. 764 764 765 '''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( u, nbgp ). \\\\765 '''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 ). \\\\ 766 766 '' 'pw-scheme' ''\\\\ 767 767 The scheme of Piascek and Williams (1970, J. Comp. Phys., 6, 392-405) with central differences in the form C3 is used. … … 954 954 The user can choose between the following schemes:\\\\ 955 955 '' 'ws-scheme' ''\\\\ 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'.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 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. When [#psolver psolver] = '' 'multigrid' '' the [#residual_limit residual_limit] should be smaller than 1.0E-6 for a sufficient reduction of the divergences. So far 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 957 958 958 '''Note''': Due to the larger stencil of this scheme vertical grid stretching should be handled with care. 959 959 The computation of turbulent fluxes takes place inside the advection routines to get a statistical evaluation consistent to the numerical solution. 960 960 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 potential tem erature: CALL exchange_horiz( pt , nbgp ). \\\\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 potential temperature: CALL exchange_horiz( pt , nbgp ). \\\\ 962 962 '' 'pw-scheme' ''\\\\ 963 963 The scheme of Piascek and Williams (1970, J. Comp. Phys., 6, 392-405) with central differences in the form C3 is used.
