You are free to use this example for academic purposes.
If you use it in a publication,
please reference the following article, which contains a detailed
description of the problem:
We used this test problem to illustrate the flexibility of the
IDR(\(s\))
code
idrs.m,
by combining it with subdomain deflation.
The corresponding codes are included.
The file idrs_ocean_example.tgz
contains the data files needed to generate the test problem
as well as a routine to read the data files.
Uncompressing and unzipping this file
(for example, by entering the command
tar -xzf idrs_ocean_example.tgz)
will create a directory (folder)
IDRS_OCEAN_EXAMPLE
containing:
Data Files
In the following,
"?"
can be either
1,
2,
3,
or
4,
and specifies grid spacing in degrees.
stommel?.mtx:
the system matrices.
stommel?_b.mtx:
the right-hand-side vectors. Each file contains twelve vectors, one for each
month.
stommel?_P.mtx:
map from the system's degrees of freedom to Earth's degrees of freedom
needed to plot the solution.
bathymetry.mtx:
depth information needed to plot the solution.
Matlab Scripts
ocean.m:
running this program opens a GUI that allows you to choose the grid
and set the solver parameters. It solves the systems for the 12
right-hand-sides and shows the resulting ocean circulation patterns as a movie.
© Martin van Gijzen
21 August 2020