THC_4box_With_Seaice: references and running instructions

------------------------------------------------------------------------
This box model is closely related to those used for:

E. Tziperman and H. Gildor. The stabilization of the thermohaline circulation by the temperature-precipitation feedback. J. Phys. Oceanogr., 32:2707–2714, 2002.
http://www.seas.harvard.edu/climate/eli/reprints/Tziperman-Gildor-2002a.pdf

Y. Kaspi, R. Sayag, and E. Tziperman. A ‘triple sea-ice state’ mechanism for the abrupt warming and synchronous ice sheet collapses during heinrich events. Paleoceanography, 19(3):PA3004, 10.1029/2004PA001009, 2004.
http://www.seas.harvard.edu/climate/eli/reprints/Kaspi-Sayag-Tziperman-2004.pdf

As well as

E. Tziperman, J. R. Toggweiler, Y. Feliks, and K. Bryan. Instability of the thermohaline circulation with respect to mixed boundary-conditions: Is it really a problem for realistic models. J. Phys. Oceanogr., 24(2):217–232, February 1994.
http://www.seas.harvard.edu/climate/eli/reprints/Tziperman-Toggweiler-Feliks-Bryan-1994.pdf

S. M. Griffies and E. Tziperman. A linear thermohaline oscillator driven by stochastic atmospheric forcing. J. Climate, 8(10):2440–2453, October 1995.
http://www.seas.harvard.edu/climate/eli/reprints/Griffies-Tziperman-1995.pdf

------------------------------------------------------------------------

To run:

compile using

make

Run using

Output/THC_4box_with_seaice.exe

Can compile, run and plot in a single command by using the Matlab program

> run_plot_THC_4box_with_seaice