Raytracing: sound waves in the ocean, gravity waves approaching a beach
Geostrophic adjustment-1: geostrophic_adjustment.m (also need this: waves_2modes.m)
Script to plot the Eady solution: eady_wave.m, eady_wave.pdf
Shallow Water Equations: Fortran files
Everything as a single tar file: swe.tar
Precompiled binary for Mac OS X: swe2d
This is a fairly complete implementation of the linear shallow water equations on the beta-plane for a fluid of variable depth. The numerics, while by no means state-of-the-art, are adequate and the code is modular and easy to follow. You can use the code to study everything from geostrophic adjustment to tides. (But watch out for the bit that switches evaluation of the Coriolis term at even and odd time steps. This avoids setting up a nasty computational mode.)
To compile: type "make" (on Linux, this expects pgf77 or pgf90 to be in your path. On OS X, the makefile will work for absoft f77.)
To run: ./swe2d
See the m-files for examples of how to generate initial conditions and a bathymetry file, and how to view the output. Otherwise, feel free to ask me.
Contact Information:
Samar Khatiwala
Oceanography 201
Lamont Doherty Earth Observatory
Columbia University
Palisades, NY 10964
Email: spk@ldeo.columbia.edu
Phone: 845-365-8454
Fax: 845-365-8736
samar khatiwala
