Using the TTD toolbox:

Note 1:

If you have version 7.5 and higher of Matlab, you are all set. If you don't, you 
will need to download the following vectorized quadrature routine:
http://faculty.smu.edu/shampine/quadva.m
(email me if you cannot find it.)

Note 2:

The functions provided here allow you to convolve known surface boundary 
conditions with the 'standard' 1-d TTDs (parameterized by a mean age (gamma) 
and width (delta)). However, to apply them to real oceanographic data, you 
need to go a bit further. For example, if you have CFC measurements and 
would like to fit a TTD to the data, you can proceed as follows:

1) Create a lookup table between tracer age and gamma/delta. You do this by sweeping 
through a range of gamma (delta) values, for each of which you compute the corresponding
(theoretical) pCFC or (equivalently) age by convolving the pCFC boundary condition with
the TTD(gamma,delta). Once you have this, you create a reverse lookup table that predicts
gamma (delta) given age.

See the script 'make_tracer_age_ttd_lookup_table.m' for an example of how to do this.

2) Use the above lookup table on actual data. Every time you have a CFC concentration, you
convert it to an age (assumptions re. saturation are applied here, hence my preference for an
age-gamma lookup table). Then you lookup the corresponding gamma (delta) value.

See the script 'make_tracer_age_ttd_lookup_table.m' for an example of how to do this.

3) Use tracer-constrained gamma and delta to predict the concentration of any other tracer.