The Explicit Planetary Isentropic-Coordinate (EPIC) Atmospheric Model

Abstract We describe a new general circulation model (GCM) designed for planetary atmospheric studies called the EPIC model. This is a finite-difference model based on the isentropic-coordinate scheme of Hsu and Arakawa (1990. Mon. Wea. Rev. 118, 1933–1959). We report on previously undocumented modifications, additions, and key practical issues that experience running the model has revealed to be important. The model integrates the hydrostatic primitive equations, which are valid for large-scale atmospheric dynamics and include gravity waves (buoyancy waves), planetary waves (Rossby waves), and horizontally propagating sound waves (Lamb waves), but not vertically propagating sound waves because of the hydrostatic approximation. The vertical coordinate is entropy in the form of potential temperature, which coincides with material surfaces for adiabatic motion. This means that there is no vertical velocity except where there is heating, which improves accuracy and helps the model maintain conservation properties over long integrations. An isentropic vertical coordinate is natural for the atmospheres of Jupiter, Saturn, Uranus, and Neptune, which are believed to have essentially adiabatic interiors that match up with the bottom of the model and is also excellent for middle-atmosphere studies on any planet. Radiative processes are parameterized by Newtonian cooling, and the latent heat of ortho–para hydrogen conversion is included when appropriate, with a suitably defined mean potential temperature. The model is written with general map factors that make it easy to configure in oblate spherical, cylindrical, or Cartesian coordinates. The code includes optional Message Passing Interface (MPI) library calls and hence runs on any Unix-based parallel computer or network cluster. An optional graphical user interface to commercial visualization software facilitates control of the model and analysis of output. Memory is allocated dynamically such that the user does not recompile to change horizontal or vertical resolution or range. Applications to date include comet impact forecasts and hindcasts for Jupiter, meridional circulation studies of Uranus and Neptune, and the accompanying paper on three-dimensional simulations of Neptune's Great Dark Spot (1998. Icarus 132, 239–265).

[1]  Stephen S. Leroy,et al.  Temperature and circulation in the stratosphere of the outer planets , 1990 .

[2]  A. E. Gill.  Atmosphere-Ocean Dynamics , 1982 .

[3]  C. Fletcher.  Computational techniques for fluid dynamics , 1992 .

[4]  Joseph B. Klemp,et al.  Numerical Simulation of Hydrostatic Mountain Waves , 1978 .

[5]  R. LeBeau,et al.  EPIC Simulations of Time-Dependent, Three-Dimensional Vortices with Application to Neptune's Great Dark Spot , 1997 .

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

[7]  J. Holton.  An introduction to dynamic meteorology , 2004 .

[8]  Dale R. Durran,et al.  The Third-Order Adams-Bashforth Method: An Attractive Alternative to Leapfrog Time Differencing , 1991 .

[9]  M. Carter.  Computer graphics: Principles and practice , 1997 .

[10]  Akio Arakawa,et al.  A locally one-dimensional semi-implicit scheme for global gridpoint shallow-water models , 1994 .

[11]  Theodore G. Shepherd,et al.  Symmetries, Conservation Laws, and Hamiltonian Structure in Geophysical Fluid Dynamics , 1990 .

[12]  Edward N. Lorenz,et al.  Available Potential Energy and the Maintenance of the General Circulation , 1955 .

[13]  T. Dowling.  Dynamics of Jovian atmospheres , 1995 .

[14]  Anthony Skjellum,et al.  Using MPI - portable parallel programming with the message-parsing interface , 1994 .

[15]  P. Gierasch,et al.  Global variation of the para hydrogen fraction in Jupiter's atmosphere and implications for dynamics on the outer planets , 1984 .

[16]  G. J. Haltiner.  Numerical Prediction and Dynamic Meteorology , 1980 .

[17]  Leonard Kollender Nash,et al.  Elements of classical and statistical thermodynamics , 1970 .

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

[19]  T. Dowling.  Estimate of Jupiter's Deep Zonal-Wind Profile from Shoemaker-Levy 9 Data and Arnol'd's Second Stability Criterion , 1995 .

[20]  Jack Dongarra,et al.  MPI: The Complete Reference , 1996 .

[21]  R. LeBeau,et al.  Dynamic response of Jupiter's atmosphere to the impact of comet Shoemaker–Levy 9 , 1994, Nature.

[22]  Akio Arakawa,et al.  Computational Design of the Basic Dynamical Processes of the UCLA General Circulation Model , 1977 .