KPP-Model implementation into GOTM

It is the aim to include the KPP model by Large et al. [1994] into GOTM. This is in order to compare it to other turbulence closure schemes and to eventually extract specific features from the KPP model and integrate them into other models.

However, the current version does not represent the results of an idealistic numerical experiment as specified  by Large et al. [1994].

The experiment is as follows:

The initial temperature is linear and stable wih a surface temperature of 22 deg C and a decrease of 1 deg per 10 m, which is equivalent to N = 0.016 1/s. Salinity has a constant value of 35 psu. The surface forcing is a constant heat loss of 100 W/m^2, surface stress is zero.

The experiment is run for 3 days and should reproduce the results as shown in figure 1 in Large et al. [1994]. However, the results of the current GOTM impementaion do look significantly different. The most striking difference is a boundary layer depth which is 12 m instead of 13.6 m and buoyancy profiles which are instable throughout the whole boundary layer, instead of being neutrally or stably stratified in the lower part of the boundary layer.

The version which can be downloaded here includes the source code and the scenario for the above described experiment and a gnuplot script which shows major time series and profiles of results.

Everybody is invited to have a look at the implementation and to comment on it. Those who find the error (or figure out an eventuel misunderstanding of the Large et al. [1994]) are invited to a beer during the next common conference or work shop.

Please return your comments to Hans Burchard (burchard@ifm.uni-hamburg.de).

To run the KPP implementation under UNIX (or LINUX), do the following steps:

1. Download the gzipped tar file which contains everything and copy it to an empty directory.
2. Give the command 'gunzip gotm4kpp.tar.gz'
3. Give the command 'tar -xvf gotm4kpp.tar'
4. Change directory to gotm4kpp and give the command 'make'
5. Change directory to ../Large and give the command 'ln -s ../gotm4kpp/gotm gotm'
6. Give the command 'gotm', such that GOTM runs for the idealised test case.
7. Give the command 'gnuplot plot' in order to view the results.