Release Notes for TTBOX patch 26.01.2010
This document is divided into the following sections:
- Installation
- Incredible New Features
- New/Modified Documentation
- New Routines
- Modified Routines
- Modified Data Structures
- Deprecated Routines
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
- 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.
- 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
- 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.
eof.