Program: MKLSPOINTCHK
mklspointchk.........determine degree of compatibility of a single point
call: incomp=mklspointchk(lsopt,statcat,seisdat,hlon,hlat,hdepth);
lsopt: (structure)
LOCSMITH options structure as returned by MKLSOPT
statcat: (structure)
Optional LSD structure containing Station descriptions
This struct may contain station descriptions
(coordinates) to be used when the SEISDAT structure does
not contain all station information (or no station info
at all)
This might be empty [], then only the station
information from the SEISDAT structure is used.
IMPORTANT:
It is not checked if STATCAT and SEISDAT contain
contradictory informations, e.g. different coordinates
for the same station. If such contradictory data is
given, it is undefined which information is actually
used, and phantastic things may happen!
seisdat: (structure)
LSD structure containing seismic arrival data for a
single event. The structure may optionally contain
station information data.
NOTE: _some_ station information must be given! If
neither STATCAT nor SEISDAT contain station info, the
location cannot be carried out!
hlon: (number) [deg]
candidate hypocenter longitude
hlat: (number) [deg]
candidate hypocenter latitude
hdepth: (number) [km]
candidate hypocenter depth below surface
results:
incomp: (numeric array) [counter]
final incompatiblity status as determined by
MKLSGRIDEVALINCOMP.
incomp(1): This value shows how many of the data
are not compatible with the returned solution.
0: all data compatible. congratulations, you have located
the quake!
>0: a number of INCOMP data values is not compatible with
the observations. You should look into MYGRID and
into the data to find out the reason (outlier data,
wrong velocity model or other problems)
Running MKLOCSMITH at higher verbosity (level 4 or higher)
will give hints which data are not compatible.
incomp(2): This is just the total number of data against
which the nodes were checked. You need this to
find out how severe the incompatibility really is.
This routine takes LSD input data as MKLOCSMITH itself, but tests only
one single point in the planet for compatibility with the given data.
The routine constructs a dummy grid which contains the desired point and
then uses the same subroutines for evaluation as MKLOCSMITH does. To do
so, it creates a dummy LSOPT structure which is partly derived from the
LSOPT structure given by the user. However, it is assumed that a
tringular grid is used.
Martin Knapmeyer, 27.06.2007
Read M-File Source Code
eof.