Articles | Volume 11, issue 5
https://doi.org/10.5194/gmd-11-1785-2018
https://doi.org/10.5194/gmd-11-1785-2018
Methods for assessment of models
 | 
08 May 2018
Methods for assessment of models |  | 08 May 2018

Impacts of the horizontal and vertical grids on the numerical solutions of the dynamical equations – Part 2: Quasi-geostrophic Rossby modes

Celal S. Konor and David A. Randall
Abstract

We use a normal-mode analysis to investigate the impacts of the horizontal and vertical discretizations on the numerical solutions of the quasi-geostrophic anelastic baroclinic and barotropic Rossby modes on a midlatitude β plane. The dispersion equations are derived for the linearized anelastic system, discretized on the Z, C, D, CD, (DC), A, E and B horizontal grids, and on the L and CP vertical grids. The effects of various horizontal grid spacings and vertical wavenumbers are discussed. A companion paper, Part 1, discusses the impacts of the discretization on the inertia–gravity modes on a midlatitude f plane.

The results of our normal-mode analyses for the Rossby waves overall support the conclusions of the previous studies obtained with the shallow-water equations. We identify an area of disagreement with the E-grid solution.

Dates
1 Introduction

In a companion paper (Konor and Randall, 2018; hereafter Part 1), we discuss the horizontal discretization of the linearized anelastic equations on the Z, C, D, CD, (DC), A, E and B grids, and vertical discretization on the L and CP grids. We introduced the DC grid in Part 1 to test the hypothesis that the CD-grid (and DC-grid) solutions are dominated by the corrector step and the grid used with it. Part 1 focuses on the dispersion of nonhydrostatic inertia–gravity modes on an f plane. The present paper gives a corresponding analysis of the dispersion of three-dimensional Rossby modes on a midlatitude β plane. Previous studies (e.g., Neta and Williams, 1989; Dukowicz, 1995) have mostly used the discrete shallow-water equations on a midlatitude β plane. Thuburn (2008) analyzed the inaccuracies of the Rossby modes on the hexagonal C grid and proposed a discretization that minimizes these inaccuracies.

We use the quasi-geostrophic and quasi-static equations in our analysis because Rossby waves are not significantly influenced by ageostrophic or nonhydrostatic effects. Furthermore, the quasi-hydrostatic equations produce an exact solution on the β plane while the basic dynamical equations, including fully compressible and anelastic equations, produce exact solutions only for particular cases. A more detailed discussion is given in the Supplement.

In Sect. 2, we present the continuous linearized anelastic equations with the quasi-geostrophic and quasi-static approximations on the midlatitude β plane and discuss the dispersion of the Rossby modes. Section 3 discusses the discretization of these equations on the seven horizontal grids listed above and the discrete dispersion of the modes. At the end of Sect. 3, we present a comparison of the performance of the grids in simulating the Rossby modes. The vertical discretization using the L and CP grids is discussed in Sect. 4. Finally, a summary and conclusions are provided in Sect. 6. Additional details are given in the Supplement.

2 Linearized anelastic equations with an isothermal basic state

In this section, we derive the basic linearized equations with the quasi-geostrophic (and quasi-static) approximations, referring to the equations of Part 1 when possible, for brevity.

Basic equations

Following Arakawa and Konor (2009), we assume quasi-geostrophic (and quasi-static) balance with the midlatitude β-plane approximation to obtain the dispersion relationship for the baroclinic and barotropic Rossby waves.

Baroclinic Rossby modes. Baroclinic modes involve vertical motions and are influenced by the static stability (B≠0 and w≠0). The equations for this case can be obtained by assuming D/t=0 and w/t=0 in Eqs. (2)–(7) of Part 1, adding a β term in the form of -β/f0P/x to the vorticity equation and replacing f with f0 in the divergence equation. The results are

(1)ωzt=-f0D-βf0Px,(2)0=f0ωz-H2P,

and

(3) N 2 D = 2 z 2 - 1 2 ρ 0 ρ 0 z 2 P t .

Note that for an isothermal atmosphere 1/ρ0ρ0/z=-1/H and N2=gκ/H. By using Eq. (10) of Part 1 in Eqs. (1)–(3), we obtain the continuous dispersion relation for the baroclinic Rossby waves as

(4) ν = - β k k 2 + 2 + f 0 2 N 2 m 2 + 1 4 H 2 .



Barotropic Rossby modes. Barotropic modes involve purely horizontal motion and are not affected by the static stability (B=0). They also satisfy w=0, D=0 and /z-1/2ρ0ρ0/zP=0. The equations that govern the barotropic motion can be obtained by using these assumptions in Eqs. (1)–(3) as

(5) ω z t = - β f 0 P x ,

and

(6) 0 = f 0 ω z - H 2 P .

Similarly, the continuous dispersion relation for the barotropic modes is given by

(7) ν = - β k k 2 + 2 .
3 Horizontal discretization on different grids and discrete dispersion equation

In this section, we discuss the discretization of the basic equations and derive the discrete dispersion relation on each horizontal grid. At the end of this section, we present an illustrative discussion of the dispersion equations showing frequency plots that are similar to the ones presented in Part 1.

3.1 Solutions for the Z grid

Baroclinic Rossby modes. We horizontally discretize Eqs. (1)–(3) on the Z grid shown in Fig. 1a of Part 1 as

(8)tωzi,j=-f0Di,j-βf0Pi+1,j-Pi-1,j2d,(9)0=f0ωzi,j-1d2Pi+1,j+Pi-1,j+Pi,j+1+Pi,j-1-4Pi,j,

and

(10) N 2 D i , j = 2 z 2 - 1 2 ρ 0 ρ 0 z 2 P i , j t ,

respectively. By using Eq. (16) of Part 1 in Eqs. (1)–(3), we obtain the discrete dispersion relation as

(11) ν = - β ξ ̃ k ξ 2 k 2 + η 2 2 + f 0 2 N 2 m 2 + 1 4 H 2 ,

where

(12)ξ̃sinkdkd,ξsin12kd12kdandηsin12d12d.



Barotropic Rossby modes. We horizontally discretize Eqs. (5) and (6) on the Z grid as

(13) t ω z i , j = - β f 0 P i + 1 , j - P i - 1 , j 2 d ,

and

0=f0ωzi,j(14)-1d2Pi+1,j+Pi-1,j+Pi,j+1+Pi,j-1-4Pi,j,

respectively. Equations (12) and (14) can also be obtained by assuming D=0 in Eqs. (8) and (9), respectively. The discrete dispersion equation for the barotropic modes is

(15) ν = - β ξ ̃ k ξ 2 k 2 + η 2 2 ,

where ξ̃, ξ and η are given by Eq. (12), which have the same definitions in Part 1. For d→0, both Eqs. (11) and (15) become identical to their continuous counterparts given by Eqs. (4) and (7), respectively. This confirms that the discrete solutions are consistent and that they correspond to the solutions of the continuous equations. On the other hand, as the zonal scale approaches the shortest resolvable zonal scale (hereafter SRZS), i.e., kdπ and ξ̃0, the discrete modes lose their ability to recognize the β effect, and the frequency of the modes becomes zero at the SRZS. This result has been derived using the β-plane approximation. It is not immediately clear whether or not it holds in true spherical geometry. This could be studied through a discrete normal-mode analysis on the sphere and/or numerical integrations of the linearized equations on the sphere.

As in Part 1, we present plots of the discrete dispersion of the Rossby modes generated by using the Z, C, D, CD, A, E and B grids. The basic state and plot design are the same as Part 1. We use β=1.62×10-11 m−1 s−1, which is typical for a midlatitude plane.

The dispersion plots for baroclinic and barotropic Rossby modes with the Z grid are presented in Fig. 1. The most striking feature is that the frequencies of all modes, for all vertical scales and horizontal grid spacings, approach zero at the SRZS. We use k=ℓ to plot these results. This is a consequence of the use of the centered finite difference to approximate the zonal pressure gradient at cell centers. As a result, the β effect cannot be recognized by any of the modes at the SRZS. Consequently, a dynamically inert mode is generated. Again, it should be checked whether or not this conclusion carries over to the linearized equations on the sphere.

https://www.geosci-model-dev.net/11/1785/2018/gmd-11-1785-2018-f01

Figure 1Plots of the absolute value of frequency of the baroclinic (red lines) and barotropic (dashed red lines) Rossby modes obtained on the Z grid for the grid spacings (a) 2 km, (b) 10 km, (c) 25 km and (d) 100 km, and for the various vertical wavenumbers. The thin blue and thick green dashed lines are the corresponding true baroclinic and barotropic frequencies, respectively.

Download

https://www.geosci-model-dev.net/11/1785/2018/gmd-11-1785-2018-f02

Figure 2Same as Fig. 1 but on the C grid.

Download

3.2 Solutions for the C grid

Baroclinic Rossby modes. We horizontally discretize Eqs. (1) and (2) on the C grid shown in Fig. 1b of Part 1 as

ωzi+1/2,j+1/2t=-f014Di,j+Di+1,j+Di,j+1+Di+1,j+1(16)-βf0Pi+3/2,j+1/2-Pi-1/2,j+1/22d,

where

(17)Pi+1/2,j+1/214Pi,j+Pi+1,j+Pi,j+1+Pi+1,j+1,and0=f014ωzi+1/2,j+1/2+ωzi+1/2,j-1/2+ωzi-1/2,j+1/2+ωzi-1/2,j-1/2-1d2Pi+1,j+Pi-1,j+Pi,j+1+Pi,j-1-4Pi,j.

Equations (16), (17) and (10) complete the set of discrete equations for the C grid. By using Eqs. (16) and (22) of Part 1, we obtain the discrete dispersion relation as

(18) ν = - μ 2 ξ ̃ β k ξ 2 k 2 + η 2 2 + μ 2 f 0 2 N 2 m 2 + 1 4 H 2 ,

where ξ̃, ξ and η are given by Eq. (12), and

(19) μ cos 1 2 k d cos 1 2 d .

The definition of μ is identical to that used in Part 1.

Barotropic Rossby modes. By using D=0 in Eqs. (16) and (17) and then using Eqs. (16) and (22) of Part 1, we obtain the discrete dispersion relation for barotropic Rossby modes on the C grid as

(20) ν = - μ 2 ξ ̃ β k ξ 2 k 2 + η 2 2 .

The discrete baroclinic and barotropic dispersion relations, Eqs. (18) and (20), for the C grid include an averaging factor μ2. This is a difference from their Z-grid counterparts, Eqs. (11) and (15). Averaging of pressure term P from the cell centers to the corners in Eq. (17) leads to the factor of μ2 in the numerators of Eqs. (18) and (20). A factor of μ2 also appears in the inertia term at the denominator of Eq. (18), due to the averaging of divergence and vorticity to each other's grid points. Since μ and ξ̃ are both equal to zero at the SRZS, dynamically inert modes exist for both the baroclinic and barotropic Rossby modes on the C grid, similar to those that exist in the Z-grid solutions.

The C-grid solutions shown in Fig. 2 are qualitatively similar to the Z-grid solutions, but the C-grid solution deviates slightly because the dispersion relation for the C grid given by Eq. (18) contains an averaging factor μ2 in the numerator. Since μ also approaches zero at the SRZS, and ξ̃ approaches zero, the small-scale modes on the C grid move or oscillate more slowly than on the Z grid. As mentioned above, at the SRZS, a dynamically inert mode is generated with the C grid.

3.3 Solutions for the D grid

Baroclinic Rossby modes. We horizontally discretize Eqs. (1) and (2) on the D grid shown in Fig. 1c of Part 1 as

(21)ωzi,jt=-f014Di-1/2,j-1/2+Di+1/2,j-1/2+Di-1/2,j+1/2+Di+1/2,j+1/2-βf0Pi+1,j-Pi-1,j2d,

and

(22)0=f014ωzi+1,j+1+ωzi,j+1+ωzi+1,j+ωzi,j-1d2Pi+3/2,j+1/2+Pi+1/2,j+3/2+Pi+1/2,j-1/2+Pi-1/2,j+1/2-4Pi+1/2,j+1/2,

respectively. In Eq. (22), Pi+1/2,j+1/214Pi,j+Pi+1,j+Pi,j+1+Pi+1,j+1. By adding the discrete version of Eq. (3) given by

(23)N214Di-1/2,j-1/2+Di+1/2,j-1/2+Di-1/2,j+1/2+Di+1/2,j+1/2-2z2-12ρ0ρ0z2tPi,j=0

to Eqs. (21) and (22), we complete the discrete equations for the D grid. The resulting discrete dispersion relation is

(24) ν = - ξ ̃ β k ξ 2 k 2 + η 2 2 + f 0 2 N 2 m 2 + 1 4 H 2 .

Barotropic Rossby modes. By using D=0 in Eqs. (21) and (22), and using Eqs. (16) and (22) of Part 1, we obtain the discrete dispersion relation for the discrete barotropic modes as

(25) ν = - ξ ̃ β k ξ 2 k 2 + η 2 2 .

The dispersion equation for the discrete baroclinic and barotropic Rossby modes on the D grid is identical to that of the Z-grid solution. In the linear system, every averaging introduces a factor μ. For nontrivial solutions of Eqs. (21)–(23), the factors of μ cancel each other. As a result, the dispersion equation is identical to that of the Z grid.

Figure 1 is effectively a plot of the frequencies for the D grid because the dispersion equations for the Z grid given by Eqs. (11) and (15) are identical to those for the D grid, as given by Eqs. (24) and (25), respectively.

3.4 Solutions for the CD grid

Baroclinic Rossby modes. By dropping the finite-difference time derivatives of divergence and vertical velocity in Eqs. (32)–(42) of Part 1 and adding -i12τμβ/f0ξ̃kP^ and -iτβ/f0ξ̃kP^ to Eqs. (32) and (38) of Part 1, respectively, we write the CD-grid equations for a midlatitude β plane as

Predictor step on the C grid:

(26)ω^z(*)=μω^zn-12τf0D^-i12τμβ/f0ξ̃kP^,(27)0=f0w^z(n)+L2P^,(28)0=-im+12HP^+B^(n),(29)B^(*)=B^n-12τN2w^,(30)D^+im-12Hw^=0.

Corrector step on the D grid:

(31)ω^zn+1=ω^zn-τf0μD^-iτβ/f0ξ̃kP^,(32)0=f0w^z(n)+μL2P^,(33)0=-im+12HP^+B^(*),(34)B^n+1=B^n-τN2w^,(35)μD^+im-12Hw^=0.

In these equations, ξ̃ is given by Eq. (12), L2ξ2k2+η22, ξ and η are given by Eq. (12) and μ is given by Eq. (15).

In this system, the divergence is a diagnostic variable, defined on the cell corners. This is why the divergence D^ is multiplied by the averaging factor μ in Eq. (31) but not in Eq. (26). Using Scheme I, as discussed in Part 1, we eliminate ω^z() by using Eq. (26) in Eq. (32) and eliminate B^() by using Eq. (29) in Eq. (33). Then Eq. (43) of Part 1 is used to obtain the real frequency and amplification factor equations as follows:

(36)e2νiτμ2N2L2+f02σm2sin(2νrτ)+12τe2νiτμ2N2βξ̃cos(2νrτ)+τeνiτμ2N2cosνrτβξ̃k+2f02eνiτμ2-1σm2sinνrτ+12τμ2N2βξ̃k=0,

and

(37) e ν i τ = - b + b 2 - 4 a c 2 a ,

where

(38a)aμ2N2L2cos2νrτ-12τμ2N2βξ̃ksin2νrτ+f02σm2cos2νrτ,(38b)b2f02μ2-1σm2cosνrτ-τμ2N2βξ̃ksinνrτ,

and

(38c) c f 0 2 1 - 2 μ 2 σ m 2 - μ 2 N 2 L 2 .

Barotropic Rossby modes. By eliminating the divergence, vertical velocity and buoyancy in Eqs. (26)–(35), we obtain the two-part dispersion equation for the barotropic modes as

(39)0=L22-12τβξ̃k2sinνrτ+τβξ̃kL2cosνrτ,

and

(40) e ν i τ = L 2 L 2 cos ν r τ - 1 2 τ β ξ ̃ k sin ν r τ .

At the SRZS, for which ξ̃=0 and μ=0, the real frequency νr becomes 0 in Eqs. (36) and (39), and the amplification factor eνiτ becomes 1 in Eqs. (37) and (40). The Supplement gives a more detailed derivation of the discrete equations.

The CD-grid solution shown by Fig. 3 is virtually identical to that for the Z-grid solutions (and D-grid solutions) shown in Figs. 1 and 2, respectively.

https://www.geosci-model-dev.net/11/1785/2018/gmd-11-1785-2018-f03

Figure 3Same as Fig. 1 but on the CD grid.

Download

DC grid

As stated above, the CD grid behaves similarly to the D grid rather than the C grid in the numerical solution of the Rossby waves on a midlatitude β plane. The normal-mode analysis of the Rossby waves with the DC grid produce a solution that is very close to the C-grid solution. A detailed discussion and frequency plots are presented in the Supplement. This is consistent with the findings of Part 1 that the correction step dominates the solutions with the CD and DC grids.

3.5 Solutions for the A grid

Baroclinic Rossby modes. We horizontally discretize Eqs. (1)–(3) on the A grid shown in Fig. 1e of Part 1 as

(41) t ω z i , j = - f 0 D i , j - β f 0 P i + 1 , j - P i - 1 , j 2 d ,

and

(42)0=fωzi,j-14d2Pi+2,j+Pi,j+2+Pi,j-2+Pi-2,j-4Pi,j,

respectively. Similarly, we obtain the discrete dispersion relation for the baroclinic Rossby modes as

(43) ν = - β ξ ̃ k ξ ̃ 2 k 2 + η ̃ 2 2 + f 0 2 N 2 m 2 + 1 4 H 2 ,

where the definition of ξ̃ is given by Eq. (12) and

(44) η ̃ sin d d .

The frequency becomes zero at the SRZS because ξ̃ is zero in the nominator of Eq. (43). This indicates the existence of a non-moving and non-oscillating computational mode. Moreover, the factor of ξ̃2 in the denominator causes the frequency to behave badly near the smallest resolvable horizontal scale.

Barotropic Rossby modes. By dropping f02/N2m2+1/4H2 in Eq. (43), the discrete dispersion relation for the barotropic Rossby modes can be obtained as

(45) ν = - β ξ ̃ k ξ ̃ 2 k 2 + η ̃ 2 2 .

The frequency of the barotropic modes becomes strongly negative (retrogressing) at the SRZS. This means that small-scale barotropic Rossby modes can behave very badly. We discuss the behavior of these modes in connection with the plots below.

https://www.geosci-model-dev.net/11/1785/2018/gmd-11-1785-2018-f04

Figure 4Same as Fig. 1 but on the A grid.

Download

https://www.geosci-model-dev.net/11/1785/2018/gmd-11-1785-2018-f05

Figure 5Same as Fig. 1 but on the E grid.

Download

https://www.geosci-model-dev.net/11/1785/2018/gmd-11-1785-2018-f06

Figure 6Same as Fig. 1 but on the B grid.

Download

Figure 4 shows the frequency of the Rossby modes obtained on the A grid. The A grid produces very fast retrogression speeds of the barotropic mode at the SRZS. The baroclinic modes with short vertical scales retrograde faster than the true solution near the SRZS, but right at the SRZS, they do not move at all.

3.6 Baroclinic and barotropic Rossby modes with the E grid

Part 1 discusses in detail the horizontal discretization on the E grid. There it is pointed out that the E grid can be viewed as the superposition of the two C grids, in which the cell centers of one C grid are placed at the corners of a second C grid. It is also shown that, from the vorticity and divergence point of view, the E grid can be viewed as a superposition of two independent and non-interacting Z grids, as shown in Fig. 1f of Part 1. The dispersion relation for the E grid is identical to that for the Z grid, but the smallest resolvable zonal scale extends to kmax=2π/d (and max=2π/d) for the E grid. Therefore, the dispersions of baroclinic and barotropic modes on the E grid are governed by Eqs. (11) and (15) with kmax and max as described above. Recall that we use a grid spacing of 2d with the E grid to maintain the same cell density as with the other grids.

The E grid produces the wildest solutions, as shown in Fig. 5. It is the only grid that generates prograding Rossby modes. The modes with all vertical scales and horizontal grid spacings used in the models generate prograding solutions near the SRZS. The deeper the mode is, the faster the progradation speed is. The prograding modes are generated near the SRZS because the factor ξ̃ yields negative values for k>π/d. A interpretation is that the finite-difference pressure gradient determined over the two-grid distance is subject to aliasing errors for zonal waves with k>π/d, which causes the system to recognize the pressure gradient with the wrong sign.

3.7 Solutions for the B grid

Baroclinic Rossby modes. We can obtain the equations for the B grid by ignoring D/t in Eq. (63) of Part 1, replacing f with f0 and using Eqs. (8) and (9). Similarly, we obtain the discrete dispersion relation for the baroclinic Rossby modes on the B grid as

(46) ν = - ξ ̃ β k ξ 2 k 2 + η 2 2 - 1 2 d 2 ξ 2 k 2 η 2 2 + f 0 N 2 m 2 + 1 4 H 2 ,

where the factors ξ, η and ξ̃ are defined by Eq. (12). The frequency becomes zero for the SRZS because ξ̃ is zero. The Laplacian term ξ2k2+η22-12d2ξ2k2η22 also approaches zero in the numerator of Eq. (46) as the zonal wavenumber approaches the SRZS. This makes the frequency behave similarly to that of the A grid.

Barotropic Rossby modes. By dropping f02/N2m2+1/4H2 in Eq. (46), we obtain the discrete dispersion relation of the barotropic Rossby modes as

(47) ν = - ξ ̃ β k ξ 2 k 2 + η 2 2 - 1 2 d 2 ξ 2 k 2 η 2 2 .

The denominator approaches zero at the SRZS, which yields an infinite retrogression speed for these modes.

Figure 6 shows the frequency of the Rossby modes on the B grid. As with the A-grid solutions, the B grid produces infinitely fast retrogression speeds for the barotropic mode at the SRZS, and the shallow baroclinic modes retrograde faster than the true solution near the SRZS and do not move at all at the SRZS.

https://www.geosci-model-dev.net/11/1785/2018/gmd-11-1785-2018-f07

Figure 7Plots of (a, d) true and discrete frequencies for the baroclinic and barotropic Rossby modes obtained on the (b, e) L and (c, f) CP grids. The thick blue and green dashed lines on the left panels indicate the true baroclinic and barotropic frequencies, respectively. The thick red and red dashed lines on the center and right panels indicate the discrete baroclinic and barotropic frequencies, respectively. The upper and lower panels show the plots for the maximum vertical integer wavenumbers of nmax=320 (δz=250 m) and nmax=80 (δz=1 km), respectively.

Download

As discussed in Sect. 3.8 of Part 1, the A, E and B grids generate multiple (or non-unique) solutions and dynamically inert modes. Here, we see that the impact of the dynamically inert modes on the short Rossby waves is very severe.

The results of our normal-mode analysis of the nonhydrostatic anelastic barotropic and baroclinic Rossby waves on a midlatitude β plane the C, D, A, E and B grids overall agree with the results of Dukowicz's (1995) normal-mode analysis with the shallow-water equations. An exception is that we include the prograding modes with the E-grid solutions, whereas Dukowicz (1995) excludes them as “inadmissible”.

Table 1A summary of the continuous and discrete dispersion relations with various horizontal and vertical grids.

Download Print Version | Download XLSX

3.8 Vertical discretization of the linear anelastic equations on the L and CP grids and discrete dispersion equation

Part 1 presents a discussion on the vertical grids, including a historical perspective, used in atmospheric models. Our purpose in this section is to assess and compare the performance of the L and CP grids in simulating Rossby modes on a midlatitude β plane through a normal-mode analysis.

3.9 The L grid

By replacing f with f0 in Eqs. (65) and (66) of Part 1, adding the β term -β/f0Pk/x to the right-hand side of Eq. (65) of Part 1 and dropping Dk/t and wk+1/2/t in Eqs. (66) and (67) of Part 1, respectively, and using Eq. (70) of Part 1, we obtain (after some manipulations) the discrete dispersion relation for the baroclinic Rossby modes as

(48) ν = - μ z 2 β k μ z 2 k 2 + 2 + f 0 2 N 2 ζ 2 m 2 + μ z 2 1 4 H 2 ,

where

(49) ζ 1 1 2 m δ z sin 1 2 m δ z and μ z cos 1 2 m δ z .

By dropping f02/N2ζ2m2+μz2/4H2 in Eq. (48), we obtain the discrete dispersion relation for the barotropic Rossby mode as

(50) ν = - β k k 2 + 2 .

In Eq. (48), the numerator is proportional to μz2, which is zero for the smallest resolvable vertical scale (SRVS), for which mδz=π. This means that, for all horizontal scales, the modes with the SRVS cannot propagate. They are dynamically inert (computational) modes. The pressure in the β term cannot recognize the SRVS buoyancy perturbation in the vertical velocity equation Eq. (67) of Part 1 with the quasi-static assumption (wk+1/2/t0). The frequency of the discrete barotropic mode given by Eq. (50) is identical to the true frequency in Eq. (7), which is expected because the barotropic mode has no vertical structure and therefore is not affected by the vertical discretization.

3.10 The CP grid

We now derive the discrete dispersion relation for the baroclinic and barotropic Rossby modes on the CP grid, following the same strategy used with the L grid. The results are

(51) ν = - β k k 2 + 2 + f 0 2 N 2 ζ 2 m 2 + μ z 2 1 4 H 2 ,

and

(52) ν = - β k k 2 + 2 ,

respectively. The dispersion equation for the baroclinic Rossby modes on the CP grid given by Eq. (51) does not have an averaging factor in the numerator, and therefore it does not allow a dynamically inert mode with zero frequency at the SRZS.

Figure 7 shows the frequencies as functions of composite horizontal wavenumber of barotropic and baroclinic Rossby modes obtained with the L and CP grids. The true frequencies are also shown in separate panels of the figure. The figure shows the results for two vertical wavenumbers (or number of layers), namely nmax=320 and 80. We included additional frequency lines corresponding to more vertical wavenumbers than were used in the plots of Sect. 3 (indicated by thinner solid lines in the plots). In the L-grid solutions shown in Fig. 7b and e, the frequency of the smallest vertical resolvable mode, identified by nmax , deviates greatly from the true frequency, which yields zero values. Similar to the case of the inertia–gravity modes, as the vertical scale approaches the smallest resolvable scale, the modes gradually lose their ability to recognize the effects of buoyancy and therefore baroclinicity. For the mode with the smallest scale, the buoyancy and baroclinicity are completely decoupled from the wind field; for that mode, the buoyancy is dynamically inert. In contrast, the frequency of the CP-grid solutions shown in Fig. 7c and f is generally close to the true frequency but slightly higher.

4 Summary and conclusions

We have discussed the effects on the dispersion of middle-latitude Rossby waves of the horizontal and vertical discretizations of the quasi-geostrophic (quasi-static) linearized equations on the A, B, C, CD, (DC), D, E and Z horizontal grids and the L and CP vertical grids. We present a summary of the discrete dispersions of Rossby modes for the horizontal and vertical grids in Table 1 for an easy comparison.

The Z, C, D and CD (DC) grids generate similar dispersion of the baroclinic and barotropic Rossby modes. All have a dynamically inert mode at the SRZS because these scales cannot recognize the β effect. The dispersion equations for the A and B grids give infinite frequencies at the SHZS. Among all horizontal grids, the E grid produces the wildest solutions. The Rossby modes of all vertical scales near the SHZS prograde, while the true modes retrograde. The A, E and B grids generate multiple (non-unique) solutions, including dynamically inert (computational) modes. The impact of the computational modes on the short Rossby modes appears very severe on these grids.

The results of our normal-mode analysis of the Rossby waves for the C, D, A, E and B grids overall agree with the results of Dukowicz's (1995) normal-mode analysis with the shallow-water equations. Dukowicz (1995) considers the prograding modes with the E-grid solutions “inadmissible”, however, while we include them.

The selection of the vertical grid impacts the dispersion of the Rossby modes as much as the horizontal grid selection. The modes with the smallest resolvable vertical scale on the L grid do not retrograde. The CP-grid solutions are much more accurate than the L-grid solutions.

Code and data availability

Fortran codes that are used to compute and plot the frequencies for the CD grid will be provided by the corresponding author upon request. Related files can also be found in http://doi.org/10.5281/zenodo.1117930.

Supplement

The supplement related to this article is available online at: https://doi.org/10.5194/gmd-11-1785-2018-supplement.

Competing interests

The authors declare that they have no conflict of interest.

Acknowledgements

We are grateful to Bill Skamarock for his comments and suggestions to improve the manuscript. We thank the reviewer Almut Gaßmann and an anonymous reviewer for their constructive and helpful comments. This research was supported by the National Science Foundation (NSF) under AGS-1500187, the US Department of Energy Office of Science DE-SC07050 (SciDAC), DE-SC00016273 (ACME) and DE-SC00016305 (CMDV).

Edited by: Paul Ullrich
Reviewed by: Almut Gaßmann and one anonymous referee

References

Arakawa, A. and Konor, C. S.: Unification of the anelastic and quasi-hydrostatic systems of equations, Mon. Weather Rev., 137, 710–726, 2009. 

Dukowicz, J. K.: Mesh effects for Rossby waves, J. Comput. Phys., 160, 336–368, 1995. 

Konor, C.: Codes and data used for plotting some figures in Konor and Randall (2017), Zenodo, https://doi.org/10.5281/zenodo.1117930, 2017. 

Konor, C. S. and Randall, D. A.: Impacts of the horizontal and vertical grids on the numerical solutions of the dynamical equations – Part 1: Nonhydrostatic inertia–gravity modes, Geosci. Model Dev., 11, 1753–1784, https://doi.org/10.5194/gmd-11-1753-2018, 2018. 

Neta, B. and Williams, R. T.: Rossby wave frequencies and group velocities for finite element and finite difference approximations to the vorticity-divergence and primitive forms of the shallow-water equations, Mon. Weather Rev., 117, 1439–1457, 1989. 

Thuburn, J.: Numerical wave propagation on the hexagonal C grid, J. Comput. Phys., 227, 5826–5858, 2008. 

Short summary
We have discussed the effects of discretization on dispersion of midlatitude Rossby waves on the A, B, C, CD, (DC), D, E and Z horizontal grids, and the L and CP vertical grids. The Z, C, D and CD (DC) grids generate similar dispersion of the baroclinic and barotropic Rossby modes. The A, E and B grids generate multiple (non-unique) solutions. The modes with the smallest resolvable vertical scale on the L grid do not retrograde. The CP-grid solutions are much more accurate than the L grid.