A set of fundamental equations that govern atmospheric dynamics are referred to as the primitive equations. Solving this system numerically is a very complex and time intensive task. A library of numeric tools called PETSc (Portable, Extensible, Toolkit for Numeric Computation) is utilized to aid in this task. Applying PETSc to help solve this system provides opportunities for atmospheric scientists to acquire and use advanced numerical methods utilized in the PETSc libraries and gives PETSc developers an opportunity to test and optimize their library's ability to solve complex nonlinear systems of equations.
To demonstrate PETSc's ability to help solve challenges in the atmospheric sciences a simplified version of the 2-Dimensional atmospheric primitive equations was developed. A crude surface energy balance model was also developed to obtain the diabatic changes in the thermodynamic energy equation. Several meteorological scenarios are available for users to manipulate in order to provide a full picture of the dynamical system and how PETSc solves this system. Several numerical solvers are also available for use that can be applied to the system. Such solvers are Backward Euler, GL, and an implicit variable time-stepping Sundials solver. Other options such as number of grid points, visualizations, and running in parallel are also available as runtime command line arguments to interested users.
Supplementary URL: http://www.mcs.anl.gov/petsc/petsc-as/