Release Notes for TTBOX patch 26.01.2010

This document is divided into the following sections:

  1. Installation
  2. Incredible New Features
  3. New/Modified Documentation
  4. New Routines
  5. Modified Routines
  6. Modified Data Structures
  7. Deprecated Routines
Installation of this patch

The installation of this patch is very simple: just copy all the m-files into your TTBOX directory, and overwrite the old versions that currently reside there. The new versions will be available in MatLab immediately.

If you encounter any problems, contact me: martin.knapmeyer(that-strange-knotty-character)dlr.de

Incredible New Features
New/Modified Documentation
New Routines
mkchkmodelintegrity
Tests velocity models (in MODEL structures as returned by MKREADND) for several bugs that make them impossible to use, and gives diagnostic informations.
mkreduceto180
Reduces epicentral distances to the interval 0°...180°, as common in seismology.
Modified Routines
mkpsampler
In case MKSMARTTAKEOFF returns an empty list of angles, MKPSAMPLER now quits instead of throwing an error.
mksmarttakeoff
An invalid assumption about sources in liquids caused MKSMARTTAKEOFF to return empty for all sources in liquids, although only S wave sources are impossible. Now only S wave sources cause empty.
mkttcurves
The travel time curve plot (figure 3 of the MKTTCURVES output) is now underlain with standard MatLab axes withour tick labels to enable the grid. The grid is switched ON by default.
mkwritend
When the model publication year was given as a number (as usual), the routine produced a warning and didn't write the year into the file. fixed.
mkfindextremedists
The "insunfficient information" error ischanged into a warning, and some actions are executed now. This is experimental.
mkraydepth
A small bug in this routine supressed most of the SKSac phase. Fixed.
mkreduceto180
Now this routine returns empty if the input is empty. (instead of crashing)
mksynopsis
The returned SYN structure now contains a field .tp in which the greatest depth reached by the described rays is contained.
mkchkmodelintegrity
The routine now tests if the proper units are used for velocities and densities (the proper units are km/s and g/ccm), and casts a warning if values are suspicious. Values are considered suspicious if they are larger than 1000, since this points to use of m/s or kg/m3
mksmarttakeoff
The output parameters minangle and maxangle were not initialized properly. This has been repaired.
mkpolarframe
The labeling of the angular scale is modified, some minor issues concerning the versatility of the routine were changed.
mkreadnd
After constructing the model structure, the new routine mkchkmodelintegrity is carried out.
mkinterpmodel
Before interpolating the model structure, mkchkmodelintegrity is carried out.
mkshootray
The previous version did not find ray solutions with epicentral distances outside 0°...180°. It did not recognize that distances Δ 360-Δ, n*360*±Delta are all the same.
mkfzeronest
This routine may now be called without input parameters. It then returns the epsilon used to control the ietration convergence. This is used internally by MKSYNOPSIS.
mkimprovemodel
An "index exceeds matrix dimensions" error was observed if a LVZ starts at the surface, i.e. if v(z>0)<v(0). This has been repaired.
mksynopsis
Now uses MKREDUCETO180 to reduce all distances into the interval 0°...180°. It also supresses all solutions that deviate from the desired distance by more than the epsilon returned by MKFZERONEST
New/Modified Data Structures
SYN structure
A new field .tp has been added to the SYN structure. This field contains the turning point depth (either from bending or from reflecting) of the ray.
Deprecated Routines

eof.