A “Vertically Lagrangian” Finite-Volume Dynamical Core for Global Models

A finite-volume dynamical core with a terrain-following Lagrangian control-volume discretization is described. The vertically Lagrangian discretization reduces the dimensionality of the physical problem from three to two with the resulting dynamical system closely resembling that of the shallow water system. The 2D horizontalto-Lagrangian-surface transport and dynamical processes are then discretized using the genuinely conservative flux-form semi-Lagrangian algorithm. Time marching is split-explicit, with large time steps for scalar transport, and small fractional steps for the Lagrangian dynamics, which permits the accurate propagation of fast waves. A mass, momentum, and total energy conserving algorithm is developed for remapping the state variables periodically from the floating Lagrangian control-volume to an Eulerian terrain-following coordinate for dealing with ‘‘physical parameterizations’’ and to prevent severe distortion of the Lagrangian surfaces. Deterministic baroclinic wave-growth tests and long-term integrations using the Held‐Suarez forcing are presented. Impact of the monotonicity constraint is discussed.

[1]  V. Starr,et al.  A QUASI-LAGRANGIAN SYSTEM OF HYDRODYNAMICAL EQUATIONS , 1945 .

[2]  R. Sadourny.  Conservative Finite-Difference Approximations of the Primitive Equations on Quasi-Uniform Spherical Grids , 1972 .

[3]  William Bourke,et al.  A multi-level spectral model. I. Formulation and hemispheric integrations , 1974 .

[4]  A. Kasahara.  Various Vertical Coordinate Systems Used for Numerical Weather Prediction , 1974 .

[5]  B. Vanleer,et al.  Towards the ultimate conservative difference scheme. IV. A new approach to numerical convection , 1977 .

[6]  B. V. Leer,et al.  Towards the ultimate conservative difference scheme. IV. A new approach to numerical convection , 1977 .

[7]  A. Arakawa,et al.  A Potential Enstrophy and Energy Conserving Scheme for the Shallow Water Equations , 1981 .

[8]  P. Woodward,et al.  The Piecewise Parabolic Method (PPM) for Gas Dynamical Simulations , 1984 .

[9]  P. Woodward,et al.  The numerical simulation of two-dimensional fluid flow with strong shocks , 1984 .

[10]  J. Hack,et al.  Description of the NCAR Community Climate Model (CCM1) , 1987 .

[11]  Richard B. Rood,et al.  Numerical advection algorithms and their role in atmospheric transport and chemistry models , 1987 .

[12]  Akio Arakawa,et al.  Numerical modeling of the atmosphere with an isentropic vertical coordinate , 1990 .

[13]  A. Harten.  ENO schemes with subcell resolution , 1989 .

[14]  S. Osher,et al.  Efficient implementation of essentially non-oscillatory shock-capturing schemes,II , 1989 .

[15]  P. Woodward,et al.  Application of the Piecewise Parabolic Method (PPM) to meteorological modeling , 1990 .

[16]  Richard B. Rood,et al.  Application of a Monotonic Upstream-biased Transport Scheme to Three-Dimensional Constituent Transport Calculations , 1991 .

[17]  Is the Midlatitude Zonal Flow Absolutely Unstable , 1993 .

[18]  M. Suárez,et al.  A proposal for the intercomparison of the dynamical cores of atmospheric general circulation models , 1994 .

[19]  David A. Randall,et al.  Geostrophic Adjustment and the Finite-Difference Shallow-Water Equations , 1994 .

[20]  Shian-Jiann Lin,et al.  A Class of the van Leer-type Transport Schemes and Its Application to the Moisture Transport in a General Circulation Model , 1994 .

[21]  P. Paolucci,et al.  The “Cubed Sphere” , 1996 .

[22]  Shian-Jiann Lin,et al.  Transport-induced interannual variability of carbon monoxide determined using a chemistry and transport model , 1996 .

[23]  Shian‐Jiann Lin,et al.  Multidimensional Flux-Form Semi-Lagrangian Transport Schemes , 1996 .

[24]  James J. Hack,et al.  Description of the NCAR Community Climate Model (CCM3). Technical note , 1996 .

[25]  J. Thuburn.  Multidimensional Flux-Limited Advection Schemes , 1996 .

[26]  H. T. Huynh.  Schemes and constraints for advection , 1997 .

[27]  Shian-Jiann Lin,et al.  An explicit flux‐form semi‐lagrangian shallow‐water model on the sphere , 1997 .

[28]  Shian-Jiann Lin,et al.  A finite‐volume integration method for computing pressure gradient force in general vertical coordinates , 1997 .

[29]  P. Roe.  Approximate Riemann Solvers, Parameter Vectors, and Difference Schemes , 1997 .

[30]  Reply to comments by T. Janjić on ‘A finite‐volume integration method for computing pressure gradient force in general vertical coordinates’ (July B, 1997, 123, 1749–1762) , 1998 .

[31]  S.-J. Lin,et al.  Development of the Joint NASA/NCAR General Circulation Model , 1999 .

[32]  Todd D. Ringler,et al.  Modeling the Atmospheric General Circulation Using a Spherical Geodesic Grid: A New Class of Dynamical Cores , 2000 .

[33]  Shian-Jiann Lin,et al.  The Global Modeling Initiative Assessment Model: Model Description, Integration and Testing of the Transport Shell , 2000 .

[34]  R. LeVeque.  Finite Volume Methods for Hyperbolic Problems: Characteristics and Riemann Problems for Linear Hyperbolic Equations , 2002 .