Program: MKLSGRIDEVALINCOMP
mklsgridevalincomp........test degree of data compatibility of grid
call: [incomp,incomplist]=mklsgridevalincomp(mygrid);
mygrid: (structure array)
LOCSMITH search grid as returned by MKLSGRIDNEW.
result: incomp: (number) [counter]
Degree of incompatibility of data.
0: data fully compatible at some nodes.
>0: best fitting node is not compatible with INCOMP data
incomplist: (numeric array) [counter]
Degree of incompatibility of the individual layers.
incomplist(i,2) is the incompatiblity degree of layer
with ID incomplist(i,1).
This can be used to determine which layers may be
removed.
A node that is acceptable as earthquake hypocenter should be compatible
with all observed data. If there are observations no node is compatible
with, these data might be outliers, or there is another reason that makes
a fully consitens location impossible (e.g. layer not in possible hypocenter depth,
or, more seriusly, wrong velocity model).
This routine determines how many of such "outlier" data exist.
Martin Knapmeyer, 16.03.2005, 26.04.2007
Read M-File Source Code
eof.