Khatiwala, S. (2007). A Computational Framework for Simulation of Biogeochemical Tracers
in the Ocean. Glob. Biogeochem. Cycles, 21,GB3001, doi:10.1029/2007GB002923. (PDF)

Khatiwala, S., M. Visbeck, and M. A. Cane (2005). Accelerated simulation of passive tracers in
ocean circulation models. Ocean Modelling, 9, 51-69. (PDF)

For more details on how to use the TMM (beyond the sketchy information provided below), feel free to email me (spk@ldeo.columbia.edu).

Samar Khatiwala

Last updated: 5/11/08

---

To start, download the following:
http://www.ldeo.columbia.edu/~spk/Research/TMM/MIT_Matrix_Global_2.8deg.tar.gz

This contains the annual and monthly mean TM's for a coarse resolution configuration of the MITgcm, and some example code. This is my 'workhorse' configuration. If you want TM's from a higher resolution model, just email me. All the example codes included are for relatively simple problems with linear source/sink/boundary conditions, where the goal is to find a steady state solution. To use the TMM for more complex, nonlinear problems, including time-dependent ones, email me for details.

Next, download the following Matlab scripts (unarchive and place them in your Matlab path):
http://www.ldeo.columbia.edu/~spk/Code/Matlab/Matrix.tar.gz
http://www.ldeo.columbia.edu/~spk/Code/Matlab/Misc.tar.gz
http://www.ldeo.columbia.edu/~spk/Code/Matlab/MITgcm.tar.gz

Optional download:
http://www.ldeo.columbia.edu/~spk/OceanCarbon.tar.gz

---

Now untar and uncompress MIT_Matrix_Global_2.8deg.tar.gz. You should get the following:

MIT_Matrix_Global_2.8deg/
    Matrix5/
    MIT/
    PaTh_Matrix5_1/
    Ar39_Matrix5_1/
    N2O_Matrix5_1/
    Phosphate_Matrix5_1/
    grid.mat
    basepath.mat
    atl_ind_pac_basin_mask.mat
    Readme
    make_grid_mit_xyz.m

I highly recommend maintaining this directory structure. It will facilitate exchange of scripts and new TM's.

- The directory Matrix5/ contains TM's and other related files. Key files are:
    Matrix5/Scripts/matrix_nocorrection_annualmean.mat contains the annual mean TM's
    Matrix5/Scripts/matrix_nocorrection.mat contains the monthly mean TM's
    Matrix5/Data/boxes.mat contains various geometry-related arrays (e.g., the volume of
    grid box, etc.)
- The directory MIT/ contains data files used to spin-up the GCM and compute the TM's. Key files are:
   GCM Forcing data (all 32 bit, direct access binary files)
    MIT/input/lev_monthly_temp.bin  monthly mean SST for restoring
    MIT/input/lev_monthly_salt.bin  monthly mean SSS for restoring
    MIT/input/shi_qnet.bin  monthly mean heat flux

   GCM output (all 32 bit, direct access binary files)
    MIT/SeasonalDST3/Ttave*  monthly mean temperature
    MIT/SeasonalDST3/Stave*  monthly mean salinity

---

Additional notes:

1) MIT_Matrix_Global_2.8deg is the 'toplevel' directory. It represents a specific configuration of the MIT GCM.

2) Create any new experiment within this toplevel directory as a new directory. As examples, I provide the following:
    PaTh_Matrix5_1/
    Ar39_Matrix5_1/
    N2O_Matrix5_1/
    Phosphate_Matrix5_1/

I tag each experiment (e.g., PaTh_Matrix5_1) with the TM (e.g., Matrix5 here) and experiment version (e.g., "1").

3) The file "basepath.mat" contains a single string variable "base_path" that is the absolute path (on your machine) of the toplevel directory. (You will of course have to change it.) All my scripts load this file and use the variable "base_path" to reference all other directories/files. So by simply changing basepath.mat your scripts should work on any machine without change. (Within the experiment directory, I make a softlink to basepath.mat: ln -s ../basepath.mat .  )

As an example of how to exploit this (and more generally how to setup an experiment), see the file: PaTh_Matrix5_1/calc_steady_state_PaTh.m. If you follow this procedure, it will greatly facilitate the exchange of scripts and data between us (and others who use the TMM).