Program: MKLSGAPMAP
mklsgapmap.......compute azimuthal gap for each point of surface
call: [map,lat,lon]=mklsgapmap(slon,slat);
[map,lat,lon]=mklsgapmap(slon,slat,elon,elat);
slat: vector of station latitudes [deg N]
slon: vector of station longitudes [deg E]
elat: event latitudes [deg N]
elon: event longitudes [deg E]
result: map: array containing the azimuthal gap for each point of surface
lat: latitudes for plotting of MAP
lon: longitudes for plotting of MAP
plot result via surf(lon,lat,map); or similar commands.
The azimuthal gap is the largest angular gap between azimtuhally adjacent
stations, measured at the point of the epicenter. The larger tha gap, the larger
the errors in horizontal and vertical event localization
longitues should be in the range 0...360
longitudes between 180 and 360 will be corrected. others will not be corrected.
latitudes should be in the range -90...90
Latitudes outside this range are in the sole responsibility of the user.
(This program is derived from an earlier one which used a coordinate
convention different of the LOCSMITH conventions. sorry for the
confusion)
If only one station is given, the program terminates without plotting a
map.
Martin Knapmeyer, 14.08.2002, 02.06.2003, 08.03.2007, 11.12.2007
Read M-File Source Code
eof.