Changes between Version 254 and Version 255 of doc/app/initialization_parameters


Ignore:
Timestamp:
Mar 6, 2015 12:56:32 PM (10 years ago)
Author:
raasch
Comment:

--

Legend:

Unmodified
Added
Removed
Modified
  • doc/app/initialization_parameters

    v254 v255  
    10821082      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.
    10831083'''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 ). \\\\
     1084The computation of turbulent fluxes takes place inside the advection routines to get a statistical evaluation consistent to the numerical solution. \\\\
    10871085'' '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.\\\\
    10891087'' 'pw-scheme' ''\\\\
    10901088      The scheme of Piacsek and Williams (1970, J. Comp. Phys., 6, 392-405) with central differences in the form C3 is used.\\\\