Editorial Type:
Article
Article Type:
Research Article

A Finite-Volume Icosahedral Shallow-Water Model on a Local Coordinate

Jin-Luen Lee NOAA/Earth System Research Laboratory, Boulder, Colorado

Search for other papers by Jin-Luen Lee in
Current site
Google Scholar
PubMed Close
and
Alexander E. MacDonald NOAA/Earth System Research Laboratory, Boulder, Colorado

Search for other papers by Alexander E. MacDonald in
Current site
Google Scholar
PubMed Close
Print Publication:
01 Apr 2009
DOI:
https://doi.org/10.1175/2008MWR2639.1
Page(s):
1422–1437
Restricted access

Abstract

An icosahedral-hexagonal shallow-water model (SWM) on the sphere is formulated on a local Cartesian coordinate based on the general stereographic projection plane. It is discretized with the third-order Adam–Bashforth time-differencing scheme and the second-order finite-volume operators for spatial derivative terms. The finite-volume operators are applied to the model variables defined on the nonstaggered grid with the edge variables interpolated using polynomial interpolation. The projected local coordinate reduces the solution space from the three-dimensional, curved, spherical surface to the two-dimensional plane and thus reduces the number of complete sets of basis functions in the Vandermonde matrix, which is the essential component of the interpolation. The use of a local Cartesian coordinate also greatly simplifies the mathematic formulation of the finite-volume operators and leads to the finite-volume integration along straight lines on the plane, rather than along curved lines on the spherical surface. The SWM is evaluated with the standard test cases of Williamson et al. Numerical results show that the icosahedral SWM is free from Pole problems. The SWM is a second-order finite-volume model as shown by the truncation error convergence test. The lee-wave numerical solutions are compared and found to be very similar to the solutions shown in other SWMs. The SWM is stably integrated for several weeks without numerical dissipation using the wavenumber 4 Rossby–Haurwitz solution as an initial condition. It is also shown that the icosahedral SWM achieves mass conservation within round-off errors as one would expect from a finite-volume model.

Corresponding author address: Dr. Jin Lee, NOAA/Earth System Research Laboratory, Global Systems Division, 325 Broadway, Boulder, CO 80303. Email: jin.lee@noaa.gov

Abstract

An icosahedral-hexagonal shallow-water model (SWM) on the sphere is formulated on a local Cartesian coordinate based on the general stereographic projection plane. It is discretized with the third-order Adam–Bashforth time-differencing scheme and the second-order finite-volume operators for spatial derivative terms. The finite-volume operators are applied to the model variables defined on the nonstaggered grid with the edge variables interpolated using polynomial interpolation. The projected local coordinate reduces the solution space from the three-dimensional, curved, spherical surface to the two-dimensional plane and thus reduces the number of complete sets of basis functions in the Vandermonde matrix, which is the essential component of the interpolation. The use of a local Cartesian coordinate also greatly simplifies the mathematic formulation of the finite-volume operators and leads to the finite-volume integration along straight lines on the plane, rather than along curved lines on the spherical surface. The SWM is evaluated with the standard test cases of Williamson et al. Numerical results show that the icosahedral SWM is free from Pole problems. The SWM is a second-order finite-volume model as shown by the truncation error convergence test. The lee-wave numerical solutions are compared and found to be very similar to the solutions shown in other SWMs. The SWM is stably integrated for several weeks without numerical dissipation using the wavenumber 4 Rossby–Haurwitz solution as an initial condition. It is also shown that the icosahedral SWM achieves mass conservation within round-off errors as one would expect from a finite-volume model.

Corresponding author address: Dr. Jin Lee, NOAA/Earth System Research Laboratory, Global Systems Division, 325 Broadway, Boulder, CO 80303. Email: jin.lee@noaa.gov

Save
Cover Monthly Weather Review
Monthly Weather Review

  • Arakawa, A., and V. Lamb, 1981: A potential enstrophy and energy conserving scheme for the shallow-water equations. Mon. Wea. Rev., 109 , 18–36.

    • Search Google Scholar
    • Export Citation

  • Bonaventura, L., and T. Ringler, 2005: Analysis of discrete shallow-water model on geodesic Delaunay grids with C-type staggering. Mon. Wea. Rev., 133 , 2351–2373.

    • Search Google Scholar
    • Export Citation

  • Durran, D. R., 1991: The third-order Adams–Bashforth method: An attractive alternative to leapfrog time differencing. Mon. Wea. Rev., 119 , 702–720.

    • Search Google Scholar
    • Export Citation

  • Haltiner, G., and R. Williams, 1980: Numerical Prediction and Dynamic Meteorology. John Wiley and Sons, 477 pp.

  • Heikes, R. H., and D. A. Randall, 1995: Numerical integration of the shallow-water equations on a twisted icosahedral grid. Part I: Basic design and results of tests. Mon. Wea. Rev., 123 , 1862–1880.

    • Search Google Scholar
    • Export Citation

  • Lee, J-L., G. L. Browning, and Y-F. Xie, 1995: Estimating divergence and vorticity from the wind profiler network hourly wind measurements. Tellus, 47A , 892–910.

    • Search Google Scholar
    • Export Citation

  • Lin, S. J., and R. Rood, 1996: Multidimensional flux-form semi-Lagrangian transport schemes. Mon. Wea. Rev., 124 , 2046–2070.

  • Lin, S. J., and R. Rood, 1997: An explicit flux-form semi-Langrangian shallow-water model on the sphere. Quart. J. Roy. Meteor. Soc., 123 , 2477–2498.

    • Search Google Scholar
    • Export Citation

  • MacDonald, A. E., J. L. Lee, and S. Sun, 2000: QNH: Design and test of a quasi-nonhydrostatic model for mesoscale weather prediction. Mon. Wea. Rev., 128 , 1016–1036.

    • Search Google Scholar
    • Export Citation

  • Majewski, D., D. LierMann, P. Prohl, B. Ritter, M. Buchhold, T. Hanisch, G. Paul, and W. Wergen, 2002: The operational global icosahedral-hexagonal gridpoint model GME: Description and high-resolution tests. Mon. Wea. Rev., 130 , 319–338.

    • Search Google Scholar
    • Export Citation

  • Orszag, S. A., 1970: Transform method for calculation of vector-coupled sums: Application to the spectral form of the vorticity equation. J. Atmos. Sci., 27 , 890–895.

    • Search Google Scholar
    • Export Citation

  • Philips, N., 1973: Principle of large-scale numerical weather prediction. Dynamic Meteorology, P. Morel, Ed., Riedel, 3–96.

  • Ringler, T. D., and D. A. Randall, 2002: A potential enstrophy and energy conserving numerical scheme for solution of the shallow-water equations on a geodesic grid. Mon. Wea. Rev., 130 , 1397–1410.

    • Search Google Scholar
    • Export Citation

  • Sadourny, R., A. Arakawa, and Y. Mintz, 1968: Integration of the non-divergent barotropic vorticity equation with an icosahedral-hexagonal grid for the sphere. Mon. Wea. Rev., 96 , 351–356.

    • Search Google Scholar
    • Export Citation

  • Stuhne, G., and W. Peltier, 1996: Vortex erosion and amalgamation in a new model of large scale flow on the sphere. J. Comput. Phys., 128 , 58–81.

    • Search Google Scholar
    • Export Citation

  • Stuhne, G., and W. Peltier, 1999: New icosahedral grid-point discretizations of the shallow water equations on the sphere. J. Comput. Phys., 148 , 23–58.

    • Search Google Scholar
    • Export Citation

  • Thuburn, J., 1997: A PV-based shallow-water model on a hexagonal-icosahedral grid. Mon. Wea. Rev., 125 , 2328–2347.

  • Thuburn, J., and Y. Li, 2000: Numerical simulations of Rossby–Haurwitz waves. Tellus, 52A , 181–189.

  • Tomita, H., M. Tsugawa, M. Satoh, and K. Goto, 2001: Shallow water model on a modified icosahedral geodesic grid by using spring dynamics. J. Comput. Phys., 174 , 579–613.

    • Search Google Scholar
    • Export Citation

  • Tomita, H., M. Satoh, and K. Goto, 2002: An optimization of the icosahedral grid modified by spring dynamics. J. Comput. Phys., 183 , 307–331.

    • Search Google Scholar
    • Export Citation

  • Tomita, H., M. Satoh, and K. Goto, 2004: A new dynamical framework of global nonhydrostatic model using the icosahedral grid. Fluid Dyn. Res., 34 , 357–400.

    • Search Google Scholar
    • Export Citation

  • van Leer, B., 1977: Towards the ultimate conservative difference scheme. III. Upstreamcentered finite-difference schemes for ideal compressible flow. J. Comput. Phys., 23 , 263–275.

    • Search Google Scholar
    • Export Citation

  • Williamson, D., 1968: Integration of the barotropic vorticity equation on a spherical geodesic grid. Tellus, 20 , 642–653.

  • Williamson, D., J. B. Drake, J. J. Hack, R. Jakob, and P. N. Swarztrauber, 1992: A standard test set for numerical approximations to the shallow water equations in spherical geometry. J. Comput. Phys., 102 , 221–224.

    • Search Google Scholar
    • Export Citation
All Time Past Year Past 30 Days
Abstract Views 0 0 0
Full Text Views 298 88 7
PDF Downloads 215 59 2