Numerical code

As a starting point, we identified a serial shallow-water Fortran code written by Paul Swarztrauber of NCAR. His code was based on a paper by Robert Sadourny [1].

Arakawa C staggered grid.

The code is a second-order finite-difference code that is calculated on a staggered grid [2]. The code prognosticates wind velocities U and V, and potential pressure P, and diagnoses potential voticity, free surface height, and the U, V components of momentum. It will be gradually modified to demonstrate how to modernize the code and allow it to be run on parallel architecture machines.

[1] R. Sadourny, The Dynamics of Finite-Difference Models of the Shallow-Water Equations, 1975. Journal of the Atmospheric Sciences, v. 32, p. 680-689.

http://ams.allenpress.com/perlserv/?request=get-pdf&doi=10.1175%2F1520-0469%281975%29032%3C0680%3ATDOFDM%3E2.0.CO%3B2

[2] A. DeCaria, Millersville University, ESCI 445 - Numerical Modeling of the Atmospherie and Oceans, Lesson 9 - Staggered Grids. 2006.

http://www.atmos.millersville.edu/~adecaria/ESCI445/esci445_09_staggered_grids.html

If you have questions about this document, please contact SCD Customer Support. You can also reach us by telephone 24 hours a day, seven days a week at 303-497-1278. Additional contact methods: consult1@ucar.edu and during business hours in NCAR Mesa Lab Suite 39.