Hi, Here are some Matlab scripts that I have written over the years. If you use them please give me credit. If you modify them or find errors, please let me know ciochett@phys.ocean.dal.ca or david@kikaz.com (this one should last longer). Thanks for your interest. 1. Calculate Length Scales for inversions (Thorpe Sorting modified by Galbraith and Kelley) It is not in final release form so the documentation is not completely coherent. The output of one file leads to the input to the next file. There are three steps. Look at ot_example.m to see how to use the files. The files should be run in order as done in the example. The original UNIX awk scripts are included as is their example. These are in reorange-V1.1. Files: ot_example.m study this to learn to run the overturn toolbox runlength.m ot_cutoff.m overturns.m inversions.m test_runlength.m test the runlength statistic Zoverturn.m Create a fake Z overturn in linear stratification 2. Mixed Layer Depth Calculations (5 methods - 1 restricted) A script to use 1 of 5 methods to calculate mixed layer depth. The files of Thompson and Fine were sent to me from Thomson so before I can release them, I need permission from him. The matlab scripts that they sent me were incorrect and I fixed them so you should ge tthe files from me. Main program is mld. Read the help on it for instruction. Files: mld.m Calculate mixed layer depth pbin02.m bin data into equal pressure depths, return the mean or median, std and # samppm.m number of samples per meter kara_mld_1.m Part of the method of Kara et al. 2000 3. Power Spectra To do power spectra, scaled correctly, in Matlab with labels and all that jazz, I made an interface to the psd.m built in function. It will add the 95% confidence intervals and use as many data ensembles as you feed it in columns. Do help multispec to see how to use it. Files: multispec.m estimate the power spectrum multispectest.m Generate N ensembles of data of set frequency and noise level 4. General purpose files Scripts that I wrote or obtained from the internet that I find useful and have used in these scripts. paper.m Set up the plot and printing paper characteristics. ooz.m Make a string with leading zeros for incremental counting on plot labels denan.m remove nans from a matrix