Monopole Vortex

This directory contains the files needed to run a medium sized problem: a reduced gravity model of a monopole vortex ring on a beta-plane interacting with a western boundary. The model parameters are chosen such that the gravity wave speed is 9m/s, the Rossby deformation radius at mid-basin is 100km, the viscosity is 50m^2/s and the reduced gravity is 0.081m/s^2. The vortex propagates westward, reflects off the western boundary, and generates Kelvin waves and short scale Rossby waves. The Kelvin waves travel around the basin and shed westward propagating Rossby waves as they propagate in the north-south direction on the western and eastern boundaries. More information on this problem can be obtained from Milliff and McWilliams or Iskandarani et al. The list of file is as follow:

Sflags.h The appropriate CPP flags for this problem.
bas11x10.13 A spectral element grid that contains an 11x10 elements of degree 12.
bas11x10.11 A spectral element grid that contains an 11x10 elements of degree 10.
bas11x10.09 A spectral element grid that contains an 11x10 elements of degree 8.
bas11x10.13.dbc
monop.txt The data file for the SEOM code
monop.in The input file to SEOM telling it where to find the grid files and where to find the data file.
pcont.m An matlab script to contour the pressure output of SEOM
pcont.com An NCAR graphics script to contour the pressure output of SEOM using contours
grid.com An NCAR graphics script to plot any of the SE grids using drawcoast.
rectelem.in to construct rectangular SE grids

Two velocity and pressure grid combinations can be used: bas11x10.13 with bas11x10.11; and bas11x10.11 with bas11x10.09. The one with 13 Gauss-Lobatto nodes for the velocity ran stably up to day 160 with a viscosity of 50 m^2/s, but the one with only 11 nodes crashed somewhere between day 140 and 120. The likely culprit are the short waves generated by the reflection and which are not damped sufficiently by the prescribed viscosity.

To run this sample problem:

Build the executable in directory src.
1. copy Sflags.h to src11 directory
2. modify file dimns.F90 in src11 to adjust npts
3. copy initialization file:
```
cp tmonop.F90 tuser.F90
```
4. Build the dependency list:
```
make clean
make depend
```
5. Build the executable:
```
make shallow
```
Move the executable to the monop directory
Run the example problem:
```
 shallow < monop.in
```
Visualize the pressure field by contouring the depth anomaly at different times. This can be done using either NCAR graphics libraries, or using the matlab script "pcont.m"
```
 contours < monop.com 
```