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.