Sheila,
Yes there is a correction. Assuming you have the source code, you can
look at src/ucf/distaz.c. Here is a clip from the comments:
* @note Calculations are based upon the reference spheroid of 1968 and
* are defined by the major radius (RAD) and the flattening (FL).
*
* @note Distance is computed using Rudoe's forumula given in GEODESY,
* secontion 2.15(b).
* (There is some numerical problem with the following formulae.
* If the station is in the southern hemisphere and the event in
* in the northern, these equations give the longer, not the
* shorter distance between the two locations. Since the
* equations are fairly messy, the simplist solution is to reverse
* the meanings of the two locations for this case.)
*
* @note Geodesy, G. Bomford, Carlendon Press, 4th ed, 1980, Section
2.15(b)
* page 121
*
* @bug Documentation for computing distance and azimumth should be added
*
* @bug Variables here are declared static, which is really not what they
* should be declared. The actual values are difficult to identify
* what they are referring to as well.
*
This method is also an option within GMT.
Arthur
On Wed, 20 Mar 2013, Sheila Peacock wrote:
Dear all,
can someone tell me whether there is a correction for earth
ellipticity when the epicentre-station distance and
azimuth are worked out by SAC for various purposes?
Thanks.
Sheila Peacock
Blacknest.
Sheila Peacock
Arthur Snoke
George Helffrich
