As an example, the second order () divergence damping can be expressed as
A balance between minimizing the smoothing and including the damping effects (i.e. prescribing ) is achieved through a combination of using values from the literature and testing. For the example of second order damping (Eq.s 87-88), Shapiro (1991) tested values of in the range 210 to 610 m s derived assuming a grid spacing km. For the current model, the value = 0.028 was used with explicit forward-in time differencing of the prognostic momentum equations, which results in = 9.310 m s for km (assuming a Courant number of 1), which obviously falls about in the middle of values that Shapiro tested. Model results can be quite sensitive to divergence damping when an explicit time integration method for the momentum equations is used, depending on the chosen value for . Slightly larger-than-optimum values damp the vertical velocity fields quite a bit in time, and too low values can result in numerical instabilities. Even reasonable values will tend to dampen the vertical velocity fields such that the amplitudes decrease noticeably in time, but again, use of higher order derivatives can greatly alleviate this problem at fairly low cost.A forward-backward time-splitting scheme for treating gravity waves is used to integrate the horizontal momentum equations which reduces the need for divergence damping, so that a lower value for can be used. For example, Skamarock et al. (2005) use a value of =0.01 for the WRF model using second order damping. Testing using this value together with the forward-backward scheme was found to give good results, however there was still a general damping over time of wavelengths larger than . For this reason, a very high order damping is currently used (=4) in order to more selectively damp small scale noise while leaving the overall divergence field relatively intact. The fourth order divergence damping terms can be expressed as
The use of fourth order damping maintains a very dynamic pattern with good structure in the vertical velocity field during the entire integration period. See the section on time discretization for more details on the forward-backward time-split integration method.