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.