• Du, Q., V. Faber, and M. Gunzburger, 1999: Centroidal Voronoi tesselations: Applications and algorithms. SIAM Rev., 41 , 637676.

  • Dukowicz, J. K., and J. R. Baumgardner, 2000: Incremental remapping as a transport/advection algorithm. J. Comput. Phys., 160 , 318335.

    • 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 , 702720.

    • Search Google Scholar
    • Export Citation
  • Giraldo, F. X., 1999: Trajectory calculations for spherical geodesic grids in Cartesian space. Mon. Wea. Rev., 127 , 16511662.

  • Heikes, R., 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 , 18621880.

    • Search Google Scholar
    • Export Citation
  • Hunke, E. C., and J. K. Dukowicz, 1997: An elastic–viscous–plastic model for sea ice dynamics. J. Phys. Oceanogr., 27 , 18491867.

  • LeVeque, R. J., 1992: Numerical Methods for Conservation Laws. Birkhauser Verlag, 214 pp.

  • Lipscomb, W. H., and E. C. Hunke, 2004: Modeling sea ice transport using incremental remapping. Mon. Wea. Rev., 132 , 13411354.

  • Masuda, Y., and H. Ohnishi, 1986: An integration scheme of the primitive equations model with an icosahedral-hexagonal grid system and its application to the shallow water equations. Short- and Medium-Range Numerical Weather Prediction, T. Matsuno, Ed., Japan Meteorological Society, 317–326.

    • Search Google Scholar
    • Export Citation
  • McGregor, J. L., 1993: Economical determination of departure points for semi-Lagrangian models. Mon. Wea. Rev., 121 , 221230.

  • Nair, R. D., and B. Machenhauer, 2002: The mass-conservative cell-integrated semi-Lagrangian advection scheme on the sphere. Mon. Wea. Rev., 130 , 649667.

    • Search Google Scholar
    • Export Citation
  • Nair, R. D., J. S. Scroggs, and F. H. M. Semazzi, 2003: A forward-trajectory global semi-Lagrangian transport scheme. J. Comput. Phys., 190 , 275294.

    • Search Google Scholar
    • Export Citation
  • Randall, D. A., T. D. Ringler, R. Heikes, P. Jones, and J. Baumgardner, 2002: Climate modeling with spherical geodesic grids. Comput. Sci. Eng., 4 , 3241.

    • Search Google Scholar
    • Export Citation
  • Ringler, T. D., and D. A. Randall, 2002a: A potential enstrophy and energy conserving numerical scheme for solution of the shallow-water equations on a geodesic grid. Mon. Wea. Rev., 130 , 13971410.

    • Search Google Scholar
    • Export Citation
  • Ringler, T. D., and D. A. Randall, 2002b: The ZM grid: An alternative to the Z grid. Mon. Wea. Rev., 130 , 14111422.

  • Ringler, T. D., R. P. Heikes, and D. A. Randall, 2000: Modeling the atmospheric general circulation using a spherical geodesic grid: A new class of dynamical cores. Mon. Wea. Rev., 128 , 24712490.

    • Search Google Scholar
    • Export Citation
  • Sadourny, R., A. Arakawa, and Y. Mintz, 1968: Integration of the nondivergent barotropic vorticity equation with an icosahedral-hexagonal grid for the sphere. Mon. Wea. Rev., 96 , 351356.

    • Search Google Scholar
    • Export Citation
  • Schär, C., and P. K. Smolarkiewicz, 1996: A synchronous and iterative flux-correction formalism for coupled transport. J. Comput. Phys., 128 , 101120.

    • Search Google Scholar
    • Export Citation
  • Smith, R. D., S. Kortas, and B. Meltz, 1995: Curvilinear coordinates for global ocean models. Tech. Rep. LA-UR-95-1146, Los Alamos National Laboratory, Los Alamos, NM, 53 pp.

  • Smolarkiewicz, P. K., 1991: Nonoscillatory advection schemes. Proc. ECMWF Seminar on Numerical Methods in Atmospheric Models, Reading, United Kingdom, European Centre for Medium-Range Weather Forecasts, 235–256.

  • Stroud, A. H., 1971: Approximate Calculation of Multiple Integrals. Prentice-Hall, 431 pp.

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

  • van Leer, B., 1979: Towards the ultimate conservative difference scheme. V. A second-order sequel to Godunov’s method. J. Comput. Phys., 32 , 101136.

    • Search Google Scholar
    • Export Citation
  • Williamson, D. L., 1968: Integration of the barotropic vorticity equations on a spherical geodesic grid. Tellus, 20 , 642653.

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

    • Search Google Scholar
    • Export Citation
  • Zalesak, S. T., 1979: Fully multidimensional flux-corrected transport algorithms for fluids. J. Comput. Phys., 31 , 335362.

All Time Past Year Past 30 Days
Abstract Views 0 0 0
Full Text Views 2 2 2
PDF Downloads 0 0 0

An Incremental Remapping Transport Scheme on a Spherical Geodesic Grid

View More View Less
  • 1 Group T-3, Los Alamos National Laboratory, Los Alamos, New Mexico
  • | 2 Department of Atmospheric Science, Colorado State University, Fort Collins, Colorado
Restricted access

Abstract

Weather and climate models contain equations for transporting conserved quantities such as the mass of air, water, ice, and associated tracers. Ideally, the numerical schemes used to solve these equations should be conservative, spatially accurate, and monotonicity-preserving. One such scheme is incremental remapping, previously developed for transport on quadrilateral grids. Here the incremental remapping scheme is reformulated for a spherical geodesic grid whose cells are hexagons and pentagons. The scheme is tested in a shallow-water model with both uniform and varying velocity fields. Solutions for standard shallow-water test cases 1, 2, and 5 are obtained with a centered scheme, a flux-corrected transport (FCT) scheme, and the remapping scheme. The three schemes are about equally accurate for transport of the height field. For tracer transport, remapping is far superior to the centered scheme, which produces large overshoots, and is generally smoother and more accurate than FCT. Remapping has a high startup cost associated with geometry calculations but is nearly twice as fast as FCT for each added tracer. As a result, remapping is cheaper than FCT for transport of more than about seven tracers.

Corresponding author address: Dr. William H. Lipscomb, Group T-3, MS B216, Los Alamos National Laboratory, Los Alamos, NM 87545. Email: Lipscomb@lanl.gov

Abstract

Weather and climate models contain equations for transporting conserved quantities such as the mass of air, water, ice, and associated tracers. Ideally, the numerical schemes used to solve these equations should be conservative, spatially accurate, and monotonicity-preserving. One such scheme is incremental remapping, previously developed for transport on quadrilateral grids. Here the incremental remapping scheme is reformulated for a spherical geodesic grid whose cells are hexagons and pentagons. The scheme is tested in a shallow-water model with both uniform and varying velocity fields. Solutions for standard shallow-water test cases 1, 2, and 5 are obtained with a centered scheme, a flux-corrected transport (FCT) scheme, and the remapping scheme. The three schemes are about equally accurate for transport of the height field. For tracer transport, remapping is far superior to the centered scheme, which produces large overshoots, and is generally smoother and more accurate than FCT. Remapping has a high startup cost associated with geometry calculations but is nearly twice as fast as FCT for each added tracer. As a result, remapping is cheaper than FCT for transport of more than about seven tracers.

Corresponding author address: Dr. William H. Lipscomb, Group T-3, MS B216, Los Alamos National Laboratory, Los Alamos, NM 87545. Email: Lipscomb@lanl.gov

Save