NonLinLoc programs file formats

Overview - Definitions - 3D Grid Files - Phase Files - Location Files - [NonLinLoc Home]

Overview

This page describes the input and output file formats used by the various NonLinLoc programs. Some formats (such as the 3D Grid format) are used by many of the programs, while others are specific to a particular program.

Many files have an easy to read ascii format with spaces between fields. To minimize disk space usage and to speed input/output some files are binary (such as the 3D Grid "buffer" file containing the grid data.

Definitions

General - ASCII Datatypes - Binary Datatypes - [NonLinLoc Home]

• General:

  Non-GLOBAL
  formats for Non-GLOBAL mode operation (TRANS SIMPLE and TRANS LAMBERT)
  GLOBAL
  formats for GLOBAL mode operation (TRANS GLOBAL)
  FileExtension=
  required extension for file name
  required
  item must be present in file
  optional
  item is optional in file
  repeatable
  item may be present multiple times file
  ignored
  may be present in file but is not used

• ASCII Datatypes:

  integer
  ascii decimal integer (i.e. 0, 5, 285)
  float
  ascii decimal floating point number (i.e. 1.0, 34.68, -4.5)
  expFloat
  ascii exponential floating point number (i.e. -5.44e+06)
  chars
  sequence of ascii characters without spaces (i.e. NO_SAVE, abcdef, /data/bigevent.dat)
  string
  sequence of ascii characters continuing the end of line, spaces are allowed (i.e. The biggest earthquake sequence in history)

• Binary Datatypes:

  char
  1 byte character
  int
  signed integer (usually 4 bytes)
  short
  short integer (usually 2 bytes)
  long
  long integer (usually 4? bytes)
  float
  4 byte floating point number
  double
  8 byte floating point number

3D Grid Files Format

3D Grid Header file - 3D Grid Data buffer file - [NonLinLoc Home]

Each 3D Grid is stored on disk with a small, simple, ASCII header file and a (possibly very large) binary data buffer file. Both files have identical names except for the extension.

• 3D Grid Header file (ASCII, FileExtension=*.hdr)

  Line 1: (required) Specifies the size and type of the 3D velocity grid.

  xNum yNum zNum xOrig yOrig zOrig dx dy dz gridType
  Fields:
  xNum yNum zNum (integer)
  number of grid nodes in the x, y and z directions
  xOrig yOrig zOrig (float)
  Non-GLOBAL: x, y and z location of the grid origin in km relative to the geographic origin.
  GLOBAL: longitude and latitude in degrees, and z in km of the location of the south-west corner of the grid.
  dx dy dz (float)
  Non-GLOBAL: grid node spacing in kilometers along the x, y and z axes
  GLOBAL: grid node spacing in degrees along the x and y axes, and in kilometers along the z axes..
  gridType (chars)
  specifies type of grid and physical quantity stored on grid:
  VELOCITY = velocity (km/sec);
  VELOCITY_METERS = velocity (m/sec);
  SLOWNESS = slowness (sec/km);
  VEL2 = velocity**2 ((km/sec)**2);
  SLOW2 = slowness**2 ((sec/km)**2);
  SLOW2_METERS = slowness**2 (sec/m)**2;
  SLOW_LEN = slowness*length (sec);
  TIME = time (sec) 3D grid;
  TIME2D = time (sec) on 2D grid / 1D model;
  PROB_DENSITY = probability density;
  MISFIT = misfit (sec);
  ANGLE = take-off angles 3D grid;
  ANGLE2D = take-off angles on 2D grid / 1D model;
  Notes:
  1. The 3D grid horizontal dimensions are in kilometers (Non-GLOBAL) or degrees (GLOBAL) with Z always in kilometers and positive down (left-handed coordinate system).
  2. The length of the grid in the x direction is (xNum-1)*dx, and similar for y and z.
  3. The location of the grid maximum in the x direction is xOrig+(xNum-1)*dx, and similar for y and z.
  Examples:
  Non-GLOBAL: 81 81 81 22.000000 42.500000 0.000000 0.200000 0.200000 0.500000 PROB_DENSITY
  GLOBAL:8 7 101 -126.000000 35.000000 0.000000 1.000000 1.000000 1.000000 PROB_DENSITY

  Line 2: (required for gridType = TIME, TIME2D, ANGLE, ANGLE2D, ignored otherwise) Specifies the label and x, y, z location of a source/station for the grid data.

  label xSrce ySrce zSrce
  Fields:
  label (chars)
  source/station label (i.e. a station code: ABC)
  xSrce ySrce (float)
  Non-GLOBAL: x and y positions relative to geographic origin in kilometers for source.
  GLOBAL: longitude and latitude in degrees for source.
  zSrce (float)
  z grid position (depth) in kilometers for source
  Example:
  CALN -8.882080 -27.537037 -1.430000

• 3D Grid Data buffer file (Binary, FileExtension=*.buf)

  Buffer: (required) Sequence of one float for each node on the grid specifying a physical value at the node (i.e. travel-time between the node and the source), or a coded set of values (i.e. take-off angles at the node for rays to the source)

  d(N) (N = 0, xNum*yNum*zNum - 1)
  Fields:
  d(N) (float)
  value at node N = (Nx*yNum + Ny)*zNum + Nz, where xNum yNum zNum are the number of grid nodes in the x, y and z directions and Nx Ny Nz are the node indexes (starting at 0) along the x, y and z axes, respectively.
  Notes:
  1. The values are stored as a succession of planes in x, each plane consisting of a succession of vectors in y, each vector consisting of a succession of points in z. Thus, to read/write the buffer requires a loop over z, contained within a loop over y, contained within a loop over x. If the file is read directly into a memory buffer, it might be accessed in C with an array d[Nx][Ny][Nz], and in FORTRAN with an array d(Nz,Ny,Nx).
  2. For take-off angle data (gridType = ANGLE, ANGLE2D) each float value contains a coded set of 3 values - the ray take-off azimuth (16 bit integer, 0 to 3600 in tenths deg clockwise from North), the ray take-off dip (12 bit integer, 0 (down) to 1800 (up) in tenths deg) and a quality value (4 bit integer, 0 (low) to 10 (high)). The 3 values are related to the float with a C union:

  typedef union

  {
      float fval;               /* float value (dummy) */
      unsigned short ival[2];   /* unsigned short values:
                                     ival[0] bits 0-3 = quality
                                             bits 4-15 = dip
                                     ival[1] = azimuth */
  } TakeOffAngles;

  The dummy variable float fval allows accessing, reading and writing of the buffer with the same functions that operate on the other data buffer types. An analogue to this union can be formed in FORTRAN using an EQUIVALENCE statement. The unsigned short ival[2] can be coded and decoded with the folowing C functions:

  /*** function to set angle values in take-off angles union */

  TakeOffAngles SetTakeOffAngles(double azim, double dip, int iqual)
  {
          TakeOffAngles angles;
          angles.ival[1] = (unsigned short) (0.5 + 10.0 * azim);
          angles.ival[0] = (unsigned short) iqual
                  + (unsigned short) 16
                  * (unsigned short) (0.5 + 10.0 * dip);
          return(angles);
  }
  /*** function to get values in take-off angles union */
  int GetTakeOffAngles(TakeOffAngles *pangles,
                          double *pazim, double *pdip, int *piqual)
  {
          *pazim = ((double) pangles->ival[1]) / 10.0;
          *pdip = ((double) (pangles->ival[0] / (int) 16)) / 10.0;
          *piqual = (int) pangles->ival[0] % (int) 16;
          return(*piqual);
  }

Phase File Formats

NonLinLoc Format - HYPO71 Format - NEIC Format - CSEM-EMSC Alert Format - NCSN Format - SIMULPS Format - SEISAN Format - RéNaSS DEP Format - [NonLinLoc Home]

The seismic phase time-pick observations for the program NLLoc may be specified in a number of different file formats. These files may contain picks for a single event or a number of events. NLLoc can read multiple observation files using a "wild-card" specification. (see the LOCFILES statement of the NLLoc control file).

Most phase file types may contain non phase pick data before or after the pick data for each event.

There are many undocumented formats supported by NLLoc, the LOCFILES statement code for supported formats include:
ISC_IMS1.0 (http://www.isc.ac.uk/Bulletin/Info/format.htm), CSEM_GSE2.0 (http://www.emsc-csem.org/)
NLLOC_OBS HYPO71 HYPOELLIPSE NEIC CSEM_ALERT SIMULPS HYPOCENTER HYPODD SEISAN NORDIC NCSN_Y2K_5 NCEDC_UCB ETH_LOC RENASS_WWW RENASS_DEP INGV_BOLL INGV_BOLL_LOCAL INGV_ARCH
See the NLLocLib.c source file for more details.

• NonLinLoc Phase file format (ASCII, NLLoc obsFileType = NLLOC_OBS)

  The NonLinLoc Phase file format is intended to give a comprehensive phase time-pick description that is easy to write and read.

  For each event to be located, this file contains one set of records. In each set there is one "arrival-time" record for each phase at each seismic station. The final record of each set is a blank. As many events as desired can be included in one file.

  Each record has a fixed format, with a blank space between fields. A field should never be left blank - use a "?" for unused characther fields and a zero or invalid numeric value for numeric fields.

  The NonLinLoc Phase file record is identical to the first part of each phase record in the NLLoc Hypocenter-Phase file output by the program NLLoc. Thus the phase list output by NLLoc can be used without modification as time pick observations for other runs of NLLoc.

  NonLinLoc phase record:

  Fields:
  Station name (char*6)
  station name or code
  Instrument (char*4)
  instument identification for the trace for which the time pick corresponds (i.e. SP, BRB, VBB)
  Component (char*4)
  component identification for the trace for which the time pick corresponds (i.e. Z, N, E, H)
  P phase onset (char*1)
  description of P phase arrival onset; i, e
  Phase descriptor (char*6)
  Phase identification (i.e. P, S, PmP)
  First Motion (char*1)
  first motion direction of P arrival; c, C, u, U = compression; d, D = dilatation; +, -, Z, N; . or ? = not readable.
  Date (yyyymmdd) (int*6)
  year (with century), month, day
  Hour/minute (hhmm) (int*4)
  Hour, min
  Seconds (float*7.4)
  seconds of phase arrival
  Err (char*3)
  Error/uncertainty type; GAU
  ErrMag (expFloat*9.2)
  Error/uncertainty magnitude in seconds
  Coda duration (expFloat*9.2)
  coda duration reading
  Amplitude (expFloat*9.2)
  Maxumim peak-to-peak amplitude
  Period (expFloat*9.2)
  Period of amplitude reading
  PriorWt (expFloat*9.2)
  A-priori phase weight Currently can be 0 (do not use reading) or 1 (use reading). (NLL_FORMAT_VER_2 - WARNING: under development)

  Example:

  GRX    ?    ?    ? P      U 19940217 2216   44.9200 GAU  2.00e-02 -1.00e+00 -1.00e+00 -1.00e+00
  GRX    ?    ?    ? S      ? 19940217 2216   48.6900 GAU  4.00e-02 -1.00e+00 -1.00e+00 -1.00e+00
  CAD    ?    ?    ? P      D 19940217 2216   46.3500 GAU  2.00e-02 -1.00e+00 -1.00e+00 -1.00e+00
  CAD    ?    ?    ? S      ? 19940217 2216   50.4000 GAU  4.00e-02 -1.00e+00 -1.00e+00 -1.00e+00
  BMT    ?    ?    ? P      U 19940217 2216   47.3500 GAU  2.00e-02 -1.00e+00 -1.00e+00 -1.00e+00
  BMT    ?    ?    ? S      ? 19940217 2216   52.8700 GAU  4.00e-02 -1.00e+00 -1.00e+00 -1.00e+00
  ESC    ?    ?    ? P      D 19940217 2216   47.4700 GAU  2.00e-02 -1.00e+00 -1.00e+00 -1.00e+00
  ESC    ?    ?    ? S      ? 19940217 2216   52.8100 GAU  4.00e-02 -1.00e+00 -1.00e+00 -1.00e+00
  BST    ?    ?    ? P      D 19940217 2216   48.0000 GAU  1.00e+05 -1.00e+00 -1.00e+00 -1.00e+00
  BST    ?    ?    ? S      ? 19940217 2216   54.6600 GAU  4.00e-02 -1.00e+00 -1.00e+00 -1.00e+00

• HYPO71 Phase file format (ASCII, NLLoc obsFileType = HYPO71)

  The HYPO71 Phase file format allows direct reading by NLLoc of HYPO71 (Lee and Lahr, 1972) and of HYPOELLIPSE (Lahr, 1989) phase data files.

  The following description of the HYPO71 Phase file format is modified from Lahr, 1989. Some HYPO71 fields are not used by NLLoc.

  For each event to be located, this file contains one set of records. In each set there is one "arrival-time" record for each seismic station with the arrival times of P and S phases, the maximum amplitude and the period for amplitude determination, and a code (0 through 4) that indicates the precision of the arrival times. The final record of each set is a blank. As many events as desired can be included in one file.

  HYPO71 phase record:

  Fields:
  Station name (Columns 1-4, char*4)
  station name or code
  P phase onset (Columns 5, char*1)
  description of P phase arrival onset; i, e
  Phase descriptor (Columns 6, char*1)
  P phase type P, N or E
  First Motion (Columns 7, char*1)
  first motion direction of P arrival; c, C, u, U = compression; d, D = dilatation; +, -, Z, N; . or blank = not readable.
  P weight (Columns 8, int*1)
  P phase weight code based on estimated standard deviation, 0, 1, 2, 3, 4
  Date/time (yymmddhhmm) (Columns 10-19, int*10)
  year (without century), month, day, hour, min
  P seconds (Columns 20-24, float*5.2)
  seconds of P arrival
  S seconds (Columns 32-36, float*5.2)
  seconds of S arrival
  S remark (Columns 37-39, char*3)
  S phase type
  S weight (Columns 40, int*1)
  S phase weight code based on estimated standard deviation, 0, 1, 2, 3, 4
  Amplitude (Columns 44-47, float*4)
  Maxumim peak-to-peak amplitude
  Period (Columns 48-50, float*3)
  period of maxumim amplitude in 1/100 sec
  Tme correction (Columns 66-70, float*5)
  time correction sec (not used by NLLoc)
  F-P time (Columns 71-75, float*5)
  F-P time interval in sec for FMAG calculation
  Example:

  ANNMiPc0 961217114029.07
  BRUMeP 1 961217114035.97
  CANM P 0 961217114029.57
  STV iPd0 961217114027.40       45.64 S 1
  ROB  P 0 961217114031.59       52.80 S 2
  FIN  P 1 961217114032.23       53.65 S 2

• NEIC "Phase (Arrival Time) Data" file format ( ASCII, NLLoc obsFileType = NEIC)

  This format reads the ASCII information in the NEIC HTML page "Phase (Arrival Time) Data" available in many of the event pages under the NEIC Current Earthquake Information page. To use this format, you must select, copy and paste to a file all of the text in the "Phase Data" frame. For some events, NLLoc may successfully read the phase data directly from the HTML source of the "Phase Data" page.

  Example:
  see http://neic.usgs.gov/neis/phase_data/phase_sample.html

• European-Mediterranean Seismological Centre (CSEM-EMSC) "EMSC ALERT - INFORMATION MESSAGE" file format ( ASCII, NLLoc obsFileType = CSEM_ALERT)

  This format reads the ASCII information in the CSEM-EMSC HTML page "EMSC Relocations of Epicentres" available for many of the events under the CSEM-EMSC "Real Time Seismicity" page. To use this format, you must save the HTML source of the "EMSC Relocations of Epicentres" frame to a file. NLLoc will read the phase data directly from the HTML source.

  Example:
  <header text>

Sta Phase Date Time Res Dist Azm Net
---------------------------------------------------------------
MELI a P 2004/02/24 02:28:02.7 -0.9 0.88 85 GFZ
IFR m i P 2004/02/24 02:28:20.0 0.3 1.94 208 SPGM
SFS mde P 2004/02/24 02:28:24.9 2.2* 2.17 305 NEIC
ELUQ m P 2004/02/24 02:28:27.2 1.9 2.34 355 MAD
EQES m P 2004/02/24 02:28:32.0 1.8 2.68 16 MAD
EHOR m P 2004/02/24 02:28:31.8 0.4 2.77 339 MAD
EHUE m P 2004/02/24 02:28:33.5 1.2 2.83 23 MAD
EADA m P 2004/02/24 02:28:35.2 0.9 2.97 351 MAD
EMIN m P 2004/02/24 02:28:39.1 -0.1 3.32 320 MAD
CART a P 2004/02/24 02:28:38.0 -2.1* 3.39 45 GFZ
AVE m i P 2004/02/24 02:28:40.2 -0.2 3.41 236 SPGM
EGRO m P 2004/02/24 02:28:42.5 -1.0 3.62 310 MAD
ETOB m P 2004/02/24 02:28:48.5 0.4 3.94 29 MAD
PBEJ m e Pn 2004/02/24 02:28:51.0 0.2 4.17 313 IMP
PBEJ m e Sn 2004/02/24 02:29:38.0 -1.1 4.17 313 IMP
ESLA mde P 2004/02/24 02:28:54.8 -0.4 4.44 0 NEIC
PTEO m e Pn 2004/02/24 02:28:54.9 0.5 4.45 302 IMP
PTEO m e Sn 2004/02/24 02:29:44.7 -0.6 4.45 302 IMP
TIO m i P 2004/02/24 02:29:03.0 -1.4 5.08 213 SPGM
PCBR m e Pn 2004/02/24 02:29:07.8 1.4 5.36 330 IMP
...

• Northern California Seismic Netowrk (NCSN) "Raw phase data (Y2K compliant)" file format ( ASCII, NLLoc obsFileType = NCSN_Y2K_5)

  The Northern California Seismic Netowrk (NCSN) "Raw phase data (Y2K compliant)" file format allows direct reading by NLLoc of Northern California Earthquake Catalog Search phase data files (available from the Northern California Earthquake Data Center (http://quake.geo.berkeley.edu/ncedc)).

  The documentation for the NCSN "Raw phase data (Y2K compliant)" file format is at ftp://quake.geo.berkeley.edu/pub/doc/cat5/ncsn.phase.y2k.5

  Note: NLLoc ignores all NCSN phases for which the P phase onset is "X", "Y", or "Z". This filter should remove all RTP picks.

  Example:

  198910180009294237  910121 5738  196340  5240 19   334579 73214810 244 00LOM  59 -  1 240 719  8  13   0  5  0LOMT2 2D  9X                10090510X342  13
  JLX  NC VVHZ EPD3198910180009 3047 -49  0    0   0   0     00 0  0  -8   0  62 6504  0    0334 00 00   0   02
  JBC  NC VVHZ EPU3198910180009 3041 -68  0    0   0   0     00 0  0   2   0  67 8104  0    0283 00 00   0   02
  JBL  NC EVLE EPU3198910180009 3274 -40  0 3604ES 3   0     00 0 52 -37 -66 191 6902  0   17262 00 00   0 9992
  JSF  NC ZVLZ IPD1198910180009 3638   1258    0   0   0  15200 0  0  33   0 341 5634 25M   0326 00340 959   02 X
  CAL  NC ZVLZ EPU3198910180009 3705   3 52    0   0   0   4700 0  0  66   0 359 5524 32M   0 22 00330 997   02 X
  JMP  NC FVFZ EPD2198910180009 3694  -7103    0   0   0   1200 0  0  21   0 385 5424 14P   0332 00350 512   02
  CAO  NC ZVLZ EPD3198910180009 3924 109  0    0   0   0   4900 0  0  55   0 433 5224 35M   0 59 00350   0   02 X
  HMO  NC VVHZ EPD4198910180009 3871-176  0    0   0   0     00 0  0 -30   0 613 5504  0    0176 00 00   0   02
  CMC  NC VVHZ EPU3198910180009 4281  -4 35    0   0   0     00 0  0  58   0 725 4904  0    0345 00 00 530   02
                                                                  10090510

• SIMULPS file format ( ASCII, NLLoc obsFileType = SIMULPS)

  The SIMULPS file format allows direct reading by NLLoc of SIMULPS "Traveltime Data for Earthquakes" format.).

  Example:

  84 3 6  6 7 58.58 46N19.78   7E26.85   3.77   2.00
  SIE_IP-1  1.45DIX2IP-1  4.94DIX_IP-1  4.95EMS_IP-1  8.14EMV_IP-1  8.54MMK2IP-1  8.58
  MMK_IP-1  8.67STG_IP-1 16.48SLE_IP-1 27.54WIL_IP-1 28.14SAX_IP-1 30.78
  84 3 7  4 7 56.92 46N20.89   7E26.59   6.09   1.50
  ZZB_IP-1  1.03ZZE_IP-1  1.06ZZA_IP-1  1.07ZZC_IP-1  1.11ZZD_IP-1  1.11ZZF_IP-1  1.11

• SEISAN Phase file format (ASCII, NLLoc obsFileType = SEISAN)

  The SEISAN Phase file format allows direct reading by NLLoc of (SEISAN/HYPOCENTER) phase data files.

  The description of the SEISAN/HYPOCENTER Phase file format below is taken from the HYPOCENTER 3.2 reference manual (hypocent.ps) (Lienert, 1994, and HYPOCENTER download). Some SEISAN/HYPOCENTER fields are not used by NLLoc.

  For each event to be located, this file contains one set of records. In each set there is one "phase header line" with general date and hyocenter information followed by optional information lines followed by a "phase data line" for each arrival times of a single P, S or other phase. The final record of each set is a blank. As many events as desired can be included in one file.

  description of the SEISAN/HYPOCENTER Phase file format

  Example:

   1997 1014 2347 54.5 L  43.069  12.817  8.0  ITA  5 0.1 0.4CITA
                  1.23       7.2     6.5  9.5 -0.2402E+01 -0.7041E+01 -0.1050E+02E
   ACTION:UPD 98-05-07 09:49 OP:fran STATUS:               ID:19971014234754     I
   9710-14-2347-11S.ITA_018                                                      6
   STAT SP IPHASW D HRMM SECON CODA AMPLIT PERI AZIMU VELO SNR AR TRES W  DIS CAZ7
   ANNI    IP     D 2347 56.36                                    -0.110    3 110
   ANNI    IS   4   2347 57.80                                    -0.2 0    3 110
   COLF    IP     D 2347 56.97   18                                0.110    8 135
   COLF    IS   1   2347 58.96                                     0.2 7    8 135
   FRAN    IP   1   2347 57.09                                     0.1 7    9 173
   FRAN    IS   1   2347 58.85                                    -0.1 7    9 173
   SVIT    IP     C 2347 57.18   22                                0.110   10 216
   SVIT    IS   1   2347 59.15                                     0.0 7   10 216
   FEMA    IP     D 2347 59.06                                    -0.1 9   21 135
   FEMA    IS   1   2348  2.70                                     0.0 7   21 135

• RéNaSS DEP Phase file format (ASCII, NLLoc obsFileType = RENASS_DEP)

  The RéNaSS DEP Phase file format allows direct reading by NLLoc of the DEP ("depouillement") format for phase picks available on the internet from the Réseau National de Surveillance Sismique, Ecole et Observatoire des Sciences de la Terre de Strasbourg.

  Some RéNaSS DEP fields are not used by NLLoc.

  This file contains one set of records for a single event. The file name must follow the RéNaSS DEP specification as this name is used by NLLoc to determine the date of the event.

  description of the RéNaSS DEP Phase file format (see section III - Les fichier DEP ("depouillement"))

Event Location Files Format

NLLoc Hypocenter-Phase file - Phase Statistics file - Scatter file - Confidence Level file - Quasi-HYPOELLIPSE file - HYPO71 Hypocenter/Station file - HypoInverse Archive file - [NonLinLoc Home]

The location results from a run of the program NLLoc are stored in individual event files and in summary files for the run. A number of different file types are output, some of these are created and saved only if specifically requested in the NLLoc control file.

The individual event files have names of the form:

path-Name.date.time.gridN.loc.FileExtension

and the summary files have names of the form:

path-Name.sum.gridN.loc.FileExtension

where path-Name is the output file path and root name specified in the LOCFILES statement of the NLLoc control file, date.time is an automatically generated data and time for the event of the form yyyymmdd.hhmmss, N is the grid index starting from 0 for the initial search grid, and FileExtension is the required file extension for each output file type as specified below.

• NLLoc Hypocenter-Phase file (ASCII, FileExtension=*.hyp)

  The NLLoc Hypocenter-Phase file is an easy to read ASCII file output by the program NLLoc. This file contains a description of the search results and the optimal/maximum-likelihood hypocenter, it may also include phase information for this hypocenter. The summary NLLoc Hypocenter-Phase file contains hypocenter information for all events from a single run of the program NLLoc but does not contain phase information. The event NLLoc Hypocenter-Phase files contain hypocenter and phase information for a single events. The NLLoc Hypocenter-Phase file output by program LocSum also include a Scatter block (SCATTER . . . END_SCATTER) containing samples of the location PDF.

  NLLOC Line: (required) Specifies the beginning of a NLLoc Hypocenter-Phase description block and gives the event file name. The end of the block is denoted by an END_NLLOC line (see below).

  NLLOC "eventFileName"
  Fields:
  eventFileName (string)
  path and root name of individual event files corresponding to this hypocenter

  SIGNATURE Line: (required) Signature text and program run stamp.

  SIGNATURE "signature verRunTime"
  Fields:
  signature (string)
  is the signature text specified in the LOCSIG statement of the NLLoc control file
  verRunTime (string)
  is an automatically generated text giving the NLLoc program version and the date/time of the run

  COMMENT Line: (required) Comment text.

  COMMENT "commentText"
  Fields:
  commentText (string)
  is the comment text specified in the LOCCOM statement of the NLLoc control file

  GRID Line: (required) 3D Grid description.

  GRID gridDescription
  Fields:
  gridDescription (string)
  is a standard grid description - see Line 1 of 3D Grid Header file above.

  SEARCH Line: (required) Search type description.

  SEARCH GRID numSamplesDraw

  or
  SEARCH METROPOLIS nSamp numSamples nAcc numAccepted nSave numSaved nClip numClipped Dstep0 stepInit Dstep stepFinal
  or
  SEARCH OCTREE nInitial nInitial nEvaluated nEvaluated smallestNodeSide sideX/sideY/sideZ

  Fields: (GRID)
  numSamplesDraw (integer)
  specifies the number of samples drawn from a PDF grid
  Fields: (METROPOLIS)
  Fields: (METROPOLIS)
  numSamples (integer)
  total number of accepted samples obtained
  numSaved (integer)
  number of accepted samples saved
  numClipped (integer)
  number of samples rejected because they were outside of the search grid (a high value indicates that event maximum likelihood location may be outside search grid)
  stepInit (integer)
  initial step size in km used for the learning stage
  stepFinal (integer)
  step size in km used for the saving stage
  Fields: (OCTREE)
  nInitial (integer)
  total initial number of octtree cells
  nEvaluated (integer)
  total number of octtree cells evaluated
  sideX/sideY/sideZ (integers)
  smallest octtree node side length (km) evaluated in x/y/z directions
  
  HYPOCENTER Line: (required) Maximum likelihood hypocenter - xyz coordinates.
  HYPOCENTER x xLoc y yLoc z zLoc OT originSec ix iX iy iY iz iZ
  Fields:
  xLoc yLoc (float)
  Non-GLOBAL: x and y positions of hypocenter relative to geographic origin in kilometers
  GLOBAL: longitude and latitude coordinates of hypocenter in degrees
  zLoc (float)
  z position (depth) of hypocenter in kilometers
  originSec (float)
  seconds of origin time of hypocenter
  iX iY iZ (integer)
  x, y and z index of grid node location of hypocenter

  GEOGRAPHIC Line: (required) Maximum likelihood hypocenter - Geographic coordinates.

  GEOGRAPHIC OT year month day hour minute second Lat latitude Long longitude Depth depth
  Fields:
  year month day (integer)
  date of origin time
  hour minute (float)
  hour and minute of origin time
  second (float)
  seconds of origin time
  latitude longitude depth (float)
  latitude and longitude in decimal degrees and depth in kilometers of hypocenter

  STAT_GEOG Line: (optional) Gaussian Expectation (Statistical) hypocenter - Geographic coordinates.

  STAT_GEOG ExpectLat latitude Long longitude Depth depth
  Fields:
  latitude longitude depth (float)
  latitude and longitude in decimal degrees and depth in kilometers of expectation hypocenter

  QUALITY Line: (required) Maximum likelihood hypocenter - Geographic coordinates.

  QUALITY Pmax probMax MFmin misfitMin MFmax misfitMax RMS rms Nphs nPhases Gap maxGap Dist minStaDist Mamp magAmp nMagAmp Mdur magDur nMagDur
  Fields:
  probMax (float)
  maximum probability on grid
  misfitMin misfitMax (float)
  maximum and minimum weighted misfit on grid
  rms (float)
  root-mean-square of residuals at maximum likelihood hypocenter
  nPhases (integer)
  number of readings used for location
  maxGap (integer)
  maximum azimuth gap in degrees between stations used for location
  minStaDist (float)
  distance of closest station to maxium likelihood hypocnter (Non-GLOBAL: kilometers, GLOBAL: degrees)
  magAmp (float)
  amplitude (i.e. ML) magnitude for maximum likelihood hypocenter
  nMagAmp (integer)
  number of readings used for amplitude magnitude
  magDur (float)
  duration magnitude for maximum likelihood hypocenter
  nMagDur (integer)
  number of readings used for duration magnitude

  VPVSRATIO Line: (required) Estimated Vp/Vs ratio. (ver 2.0)

  VPVSRATIO VpVsRatio ratio Npair numPair
  Fields:
  ratio (float)
  Vp/Vs ratio estimated following the iterative minimization technique of Lahr (1989; ch. 5)
  Npair (integer)
  number of P-S pairs used to calculate Vp/Vs ratio

  STATISTICS Line: (required) "Traditional" Gaussian (normal) statistics of PDF (evaluated for PROB_DENSITY grids only).

  STATISTICS ExpectX expX Y expY Z expZ CovXX covXX XY covXY XZ covXZ YY covYY YZ covYZ ZZ covZZ EllAz1 Azimuth1 Dip1 Dip1 Len1 StdErr1 Az2 Azimuth2 Dip2 Dip2 Len2 StdErr2 Len3 StdErr3
  Fields:
  expX expY expZ (expFloat)
  Non-GLOBAL: expectation values of PDF for x, y and z relative to geographic origin in kilometers
  GLOBAL: expectation values of PDF for longitude and latitude in degrees and depth in kilometers
  covXX covXY ... (expFloat)
  covariance values of PDF
  Azimuth1 (float)
  azimuth of axis 1 of confidence ellipsoid
  Dip1 (float)
  dip of axis 1 of confidence ellipsoid
  StdErr1 (expFloat)
  length of semi-axis 1 of confidence ellipsoid
  Azimuth2 (float)
  azimuth of axis 2 of confidence ellipsoid
  Dip2 (float)
  dip of axis 2 of confidence ellipsoid
  StdErr2 (expFloat)
  length of semi-axis 2 of confidence ellipsoid
  StdErr3 (expFloat)
  length of semi-axis 3 of confidence ellipsoid
  Notes:
  1. See NLLoc Program Inversion Approach for more information on the calculation of the expectation, covariance and confidence ellipsoid "traditional" Gaussian statistics.
  2. The confidence ellipsoid represents the 68% confidence region for 3 degrees of freedom (x, y and z).

  TRANS Line: (required) Geographic to rectangular transformation parameters.

  TRANS SIMPLE latOrig latOrig longOrig longOrig rotAngle rotAngle

  or TRANS LAMBERT refEllipsoid refEllipsoid latOrig latOriglongOrig longOrigfirstStdParal firstStdParalsecondStdParal secondStdParal
  or TRANS GLOBAL
  Specifies geographic to rectangular transformation parameters.
  Fields: (SIMPLE)
  latOrig (float)
  latitude in decimal degrees of the rectangular coordinates origin
  longOrig (float)
  longitude in decimal degrees of the rectangular coordinates origin
  rotAngle (float)
  rotation angle of the rectangular coordinates system Y-axis in degrees counter-clockwise (?) from geographic north (not implemented)
  Fields: (LAMBERT)
  refEllipsoid (chars)
  reference ellipsoid (one of: WGS-84, GRS-80, WGS-72, Australian, Krasovsky, International, Hayford-1909, Clarke-1880, Clarke-1866, Airy, Bessel, Hayford-1830, Sphere)
  latOrig (float)
  latitude in decimal degrees of the rectangular coordinates origin
  longOrig (float)
  longitude in decimal degrees of the rectangular coordinates origin
  firstStdParal secondStdParal (float)
  first and second standard parallels (meridians) in decimal degrees

  FOCALMECH Line: (optional) Focal mechanism description line. Specifies the hypocenter, focal mechanism, and mechanisn statistics

  FOCALMECH dlat dlong depth Mech dipDir dipAng rake mf misfit nObs nObs
  Fields:
  dlat dlong (floats)
  latitude and logitude in decimal degrees of the hypocenter
  depth (float)
  depth in kilometers of the hypocenter
  dipDir dipAng rake (floats)
  dip strike and rake of double couple mechanism
  mf (float)
  misfit of mechanism determination
  nObs (float)
  number of first motion observations used in mechanism determination

  PHASE Line: (optional) Phase format description line. Specifies the beginning of a block with a list of phase records. The end of the block is denoted by an END_PHASE line (see below)

  PHASE ID Ins Cmp On Pha FM Date HrMn Sec Err ErrMag Coda Amp Per > TTpred Res Weight StaLoc(X Y Z) SDist SAzim RAz RDip RQual Tcorr
  Fields:
  None

  next N Lines: (optional) Phase data for N phases.

  ID Ins Cmp On Pha FM Q Date HrMn Sec Coda Amp Per > Err ErrMag TTpred Res Weight StaLoc(X Y Z) SDist SAzim RAz RDip RQual TcorrTcorr
  Fields:
  ID (char*6)
  station name or code
  Ins (char*4)
  instrument identification for the trace for which the time pick corresponds (i.e. SP, BRB, VBB)
  Cmp (char*4)
  component identification for the trace for which the time pick corresponds (i.e. Z, N, E, H)
  On (char*1)
  description of P phase arrival onset; i, e
  Pha (char*6)
  Phase identification (i.e. P, S, PmP)
  FM (char*1)
  first motion direction of P arrival; c, C, u, U = compression; d, D = dilatation; +, -, Z, N; . or ? = not readable.
  Date (yyyymmdd) (int*6)
  year (with century), month, day
  HrMn (hhmm) (int*4)
  Hour, min
  Sec (float*7.4)
  seconds of phase arrival
  Err (char*3)
  Error/uncertainty type; GAU
  ErrMag (expFloat*9.2)
  Error/uncertainty magnitude in seconds
  Coda (expFloat*9.2)
  coda duration reading
  Amp (expFloat*9.2)
  Maxumim peak-to-peak amplitude
  Per (expFloat*9.2)
  Period of amplitude reading
  > (char*1)
  Required separator between first part (observations) and second part (calculated values) of phase record.
  TTpred (float*9.4)
  Predicted travel time
  Res (float*9.4)
  Residual (observed - predicted arrival time)
  Weight (float*9.4)
  Phase weight
  StaLoc(X Y Z) (3 * float*9.4)
  Non-GLOBAL: x, y, z location of station in transformed, rectangular coordinates
  GLOBAL: longitude, latitude, z location of station
  SDist (float*9.4)
  Maximum likelihood hypocenter to station epicentral distance in kilometers
  SAzim (float*6.2)
  Maximum likelihood hypocenter to station epicentral azimuth in degrees CW from North
  RAz (float*5.1)
  Ray take-off azimuth at maximum likelihood hypocenter in degrees CW from North
  RDip (float*5.1)
  Ray take-off dip at maximum likelihood hypocenter in degrees upwards from vertical down (0 = down, 180 = up)
  RQual (float*5.1)
  Quality of take-off angle estimation (0 = unreliable, 10 = best)
  Tcorr (float*9.4)
  Time correction (station delay) used for location

  PHASE Line (NLL_FORMAT_VER_2 - WARNING: under development): (optional) Phase format description line. Specifies the beginning of a block with a list of phase records. The end of the block is denoted by an END_PHASE line (see below)

  PHASE ID Ins Cmp On Pha FM Date HrMn Sec Err ErrMag Coda Amp Per PriorWt > TTpred Res Weight StaLoc(X Y Z) SDist SAzim RAz RDip RQual Tcorr TTerr
  Fields:
  None

  next N Lines: (optional) Phase data for N phases.

  ID Ins Cmp On Pha FM Q Date HrMn Sec Coda Amp Per PriorWt > Err ErrMag TTpred Res Weight StaLoc(X Y Z) SDist SAzim RAz RDip RQual Tcorr TTerrTcorr
  Fields:
  ID (char*6)
  station name or code
  Ins (char*4)
  instrument identification for the trace for which the time pick corresponds (i.e. SP, BRB, VBB)
  Cmp (char*4)
  component identification for the trace for which the time pick corresponds (i.e. Z, N, E, H)
  On (char*1)
  description of P phase arrival onset; i, e
  Pha (char*6)
  Phase identification (i.e. P, S, PmP)
  FM (char*1)
  first motion direction of P arrival; c, C, u, U = compression; d, D = dilatation; +, -, Z, N; . or ? = not readable.
  Date (yyyymmdd) (int*6)
  year (with century), month, day
  HrMn (hhmm) (int*4)
  Hour, min
  Sec (float*7.4)
  seconds of phase arrival
  Err (char*3)
  Error/uncertainty type; GAU
  ErrMag (expFloat*9.2)
  Error/uncertainty magnitude in seconds
  Coda (expFloat*9.2)
  coda duration reading
  Amp (expFloat*9.2)
  Maxumim peak-to-peak amplitude
  Per (expFloat*9.2)
  Period of amplitude reading
  PriorWt (expFloat*9.2)
  A-priori phase weight
  > (char*1)
  Required separator between first part (observations) and second part (calculated values) of phase record.
  TTpred (float*9.4)
  Predicted travel time
  Res (float*9.4)
  Residual (observed - predicted arrival time)
  Weight (float*9.4)
  Phase weight (covariance matrix weight for LOCMETH GAU_ANALYTIC, posterior weight for LOCMETH EDT EDT_OT_WT)
  StaLoc(X Y Z) (3 * float*9.4)
  Non-GLOBAL: x, y, z location of station in transformed, rectangular coordinates
  GLOBAL: longitude, latitude, z location of station
  SDist (float*9.4)
  Maximum likelihood hypocenter to station epicentral distance in kilometers
  SAzim (float*6.2)
  Maximum likelihood hypocenter to station epicentral azimuth in degrees CW from North
  RAz (float*5.1)
  Ray take-off azimuth at maximum likelihood hypocenter in degrees CW from North
  RDip (float*5.1)
  Ray take-off dip at maximum likelihood hypocenter in degrees upwards from vertical down (0 = down, 180 = up)
  RQual (float*5.1)
  Quality of take-off angle estimation (0 = unreliable, 10 = best)
  Tcorr (float*9.4)
  Time correction (station delay) used for location
  TTerr (expFloat*9.2)
  Traveltime error used for location

  END_PHASE Line: (required if Phase format description line is present) Specifies the end of a block with a list of phase records.

  END_PHASE
  Fields:
  None

  END_NLLOC Line: (required) Specifies the end of a NLLoc Hypocenter-Phase description block.

  END_NLLOC
  Fields:
  None

  SCATTER Line: (optional) Specifies the begining of a NLLoc Scatter block containing samples of the location PDF.

  SCATTER Nsamples nSamples
  Fields:
  Nsamples (integer)
  number of scatter samples

  next Nsamples Lines: (optional) PDF samples.

  x y z prob
  Fields:
  x y z (float)
  Non-GLOBAL: x, y and z location of the sample in kilometers relative to the geographic origin.
  GLOBAL: longitude and latitude in degrees, and z location in kilometers of the sample.
  prob (float)
  relative probability density value at sample (for Grid Search, the probability density values are normalized so that the volume integral over the corresponding search grid of the PDF = 1.0)

  END_SCATTER Line: (required if SCATTER line is present) Specifies the end of a NLLoc Scatter block.

  END_SCATTER
  Fields:
  None

  Examples:

  Non-GLOBAL

  NLLOC "/temp/nlloc_tmp/durance/loc_test/dur_OCT.19990103.212657.22.grid0" "LOCATED" "Location completed."
  SIGNATURE "  IRSN - Fontenay-aux-Roses   NLLoc:v3.10.5 02May2005 10h46m48"
  COMMENT "       IRSN Reseau Durance (Oct-tree search / vox 3D Model)"
  GRID  256 332 88  0.000000 0.000000 -2.100000  0.250000 0.250000 0.250000 PROB_DENSITY
  SEARCH OCTREE nInitial 1600 nEvaluated 50000 smallestNodeSide 0.062256/0.064648/0.067969
  HYPOCENTER  x 24.497681 y 32.291895 z 2.827734  OT 56.341531  ix -1 iy -1 iz -1
  GEOGRAPHIC  OT 1999 01 03  21 26 56.341531  Lat 43.711240 Long 5.665205 Depth 2.827734
  QUALITY  Pmax 6.929e-01 MFmin 0.349723 MFmax 94617.715955 RMS 0.070806 Nphs 6 Gap 167 Dist 2.498662 Mamp  0.57 3 Mdur -9.90 0
  VPVSRATIO  VpVsRatio 5.000e-01  Npair 2  Diff 1.027e+00
  STATISTICS  ExpectX 2.4282e+01 Y 3.2226e+01 Z 2.6602e+00  CovXX 6.90e-01 XY 1.05e-01 XZ 2.66e-01 YY 2.43e-01 YZ 3.10e-01 ZZ 1.81e+00 EllAz1  353.4 Dip1  -9.6 Len1  7.92e-01 Az2  85.8 Dip2  -13.8 Len2  1.49e+00 Len3  2.61e+00
  STAT_GEOG  ExpectLat 43.710727 Long 5.662503 Depth 2.660174
  TRANSFORM  LAMBERT RefEllipsoid Clarke-1880  LatOrig 43.430100  LongOrig 5.346580  FirstStdParal 45.898900  SecondStdParal 47.696000  RotCW -2.190000
  FOCALMECH  Hyp  43.711240 5.665205 2.827734 Mech  0.0 0.0 0.0 mf  0.00 nObs 0
  PHASE ID Ins Cmp On Pha  FM Date     HrMn   Sec     Err  ErrMag    Coda      Amp       Per      >   TTpred    Res       Weight    StaLoc(X  Y         Z)        SDist    SAzim  RAz   RDip  RQual    Tcorr
  ESC4   ?    ?    ? P      ? 19990103 2126   57.2160 GAU  2.00e-02 -1.00e+00 -1.00e+00 -1.00e+00 >    0.8397    0.0348    1.0146   25.0170   34.7360   -0.4050    2.4987  14.19  17.0 139.2  7     0.0000
  ESC4   ?    ?    ? S      ? 19990103 2126   57.8090 GAU  4.00e-02 -1.00e+00  4.61e-08  2.47e-01 >    1.4947   -0.0272    0.9854   25.0170   34.7360   -0.4050    2.4987  14.19  17.0 139.2  7     0.0000
  JOU1   ?    ?    ? P      ? 19990103 2126   57.9910 GAU  2.00e-02 -1.00e+00 -1.00e+00 -1.00e+00 >    1.7869   -0.1374    1.0146   24.3315   23.6050   -0.3000    8.6884 183.29 183.2  99.7  7     0.0000
  JOU1   ?    ?    ? S      ? 19990103 2126   59.6110 GAU  4.00e-02 -1.00e+00  6.27e-09  4.77e-01 >    3.1807    0.0888    0.9854   24.3315   23.6050   -0.3000    8.6884 183.29 183.2  99.7  7     0.0000
  CAD1   ?    ?    ? S      ? 19990103 2126   59.7670 GAU  4.00e-02 -1.00e+00  1.80e-08  3.26e-01 >    3.4149    0.0105    0.9854   32.9072   28.7836   -0.3700    9.1120 114.83 112.2  97.2  9     0.0000
  BST2   ?    ?    ? P      ? 19990103 2126   58.6070 GAU  2.00e-02 -1.00e+00 -1.00e+00 -1.00e+00 >    2.2328    0.0327    1.0146   22.4312   42.3289   -0.5000   10.2476 350.56 341.4 101.9  9     0.0000
  END_PHASE
  END_NLLOC

  GLOBAL

  NLLOC "/temp/nlloc_tmp/20041226/loc/AUTO_90_all_RMS.20041226.011057.grid0" "LOCATED" "Location completed."
  SIGNATURE "Anthony Lomax (www.alomax.net)   NLLoc:v3.10.5 18Mar2005 08h16m37"
  COMMENT "Coda location"
  GRID  361 181 601  -180.000000 -90.000000 0.000000  1.000000 1.000000 1.000000 PROB_DENSITY
  SEARCH OCTREE nInitial 3888 nEvaluated 50016 smallestNodeSide 8.680556/8.680556/0.781250
  HYPOCENTER  x 92.695312 y 11.679688 z 0.390625  OT 6.530129  ix -1 iy -1 iz -1
  GEOGRAPHIC  OT 2004 12 26  01 07  6.530129  Lat 11.679688 Long 92.695312 Depth 0.390625
  QUALITY  Pmax 1.323e-79 MFmin 1.851182 MFmax 3.439646 RMS 49.635872 Nphs 106 Gap 35 Dist 12.629659 Mamp -9.90 0 Mdur -9.90 0
  VPVSRATIO  VpVsRatio -1.000e+00  Npair 0  Diff -2.000e+30
  STATISTICS  ExpectX 9.2724e+01 Y 1.1697e+01 Z 1.2764e+02  CovXX 5.24e+03 XY -3.11e+02 XZ 5.37e+02 YY 3.53e+03 YZ 1.97e+02 ZZ 1.07e+04 EllAz1  10.6 Dip1  -2.3 Len1  1.11e+02 Az2  100.8 Dip2  -5.2 Len2  1.36e+02 Len3  1.95e+02
  STAT_GEOG  ExpectLat 11.696777 Long 92.723993 Depth 127.640604
  TRANSFORM  GLOBAL
  FOCALMECH  Hyp  11.679688 92.695312 0.390625 Mech  0.0 0.0 0.0 mf  0.00 nObs 0
  PHASE ID Ins Cmp On Pha  FM Date     HrMn   Sec     Err  ErrMag    Coda      Amp       Per      >   TTpred    Res       Weight    StaLoc(X  Y         Z)        SDist    SAzim  RAz   RDip  RQual    Tcorr
  PALK   (cou Z    ? P      ? 20041226 0110   57.9600 GAU  3.44e+01 -1.00e+00  1.15e+12 -1.00e+00 >  180.8273   50.6026    0.6321   80.7022    7.2728   -0.4600   12.6297 250.63   0.0   0.0  0     0.0000
  LSA    (cou Z    ? P      ? 20041226 0113    7.2350 GAU  4.15e+01 -1.00e+00  8.07e+10 -1.00e+00 >  252.6406  108.0643    0.2560   91.1500   29.7000   -3.7890   18.0912 355.67   0.0   0.0  0     0.0000
  COCO   (cou Z    ? P      ? 20041226 0111   19.3200 GAU  2.41e+01 -1.00e+00  2.09e+11 -1.00e+00 >  318.4248  -65.6349    3.2995   96.8349  -12.1901    0.0690   24.2392 170.10   0.0   0.0  0     0.0000
  ENH    (cou Z    ? P      ? 20041226 0112   32.2600 GAU  2.33e+01 -1.00e+00  1.69e+11 -1.00e+00 >  318.7041    7.0258    1.0282  109.4868   30.2718   -0.4870   24.2698  37.40   0.0   0.0  0     0.0000
  TPUB   (cou Z    ? P      ? 20041226 0113   56.6700 GAU  3.65e+01 -1.00e+00  2.59e+09 -1.00e+00 >  361.4501   48.6898    0.5882  120.6300   23.3000    0.0000   29.0119  62.59   0.0   0.0  0     0.0000
  TWGB   (cou Z    ? P      ? 20041226 0113   18.8000 GAU  1.96e+01 -1.00e+00  5.21e+09 -1.00e+00 >  363.6115    8.6584    2.0749  121.0800   22.8200   -0.1270   29.2555  63.80   0.0   0.0  0     0.0000
  SSLB   (cou Z    ? P      ? 20041226 0113   26.1300 GAU  1.22e+01 -1.00e+00  4.06e+09 -1.00e+00 >  365.3290   14.2709    5.1465  120.9500   23.7900    0.0000   29.4496  61.84   0.0   0.0  0     0.0000
  NACB   (cou Z    ? P      ? 20041226 0113   22.9200 GAU  4.14e+01 -1.00e+00  2.57e+09 -1.00e+00 >  371.3581    5.0318    0.4544  121.5900   24.1800    0.0000   30.1311  61.49   0.0   0.0  0     0.0000
  TATO   (cou Z    ? P      ? 20041226 0112   45.7700 GAU  6.18e+01 -1.00e+00  3.78e+09 -1.00e+00 >  373.0050  -33.7651    0.2000  121.4881   24.9754   -0.0530   30.3172  59.94   0.0   0.0  0     0.0000
  ...
  LAST   (cou Z    ? P      ? 20041226 0117   25.0300 GAU  4.31e+01 -1.00e+00  1.71e+09 -1.00e+00 >    0.0000    0.0000    0.0000    0.0000    0.0000    0.0000    0.0000   0.00   0.0   0.0  0     0.0000
  ATTU   (cou Y    ? P      ? 20041226 0118   30.3000 GAU  1.20e+01 -1.00e+00  3.04e+12 -1.00e+00 >    0.0000    0.0000    0.0000    0.0000    0.0000    0.0000    0.0000   0.00   0.0   0.0  0     0.0000
  IL31   (cou Y    ? P      ? 20041226 0122   18.4500 GAU  4.50e+01 -1.00e+00  1.88e+11 -1.00e+00 >    0.0000    0.0000    0.0000    0.0000    0.0000    0.0000    0.0000   0.00   0.0   0.0  0     0.0000
  END_PHASE
  END_NLLOC

• Phase Statistics file (ASCII, FileExtension=*.stat, *.stat_totcorr)

  The Phase Statistics file (*.stat) contains the average residuals (denoted by Average Phase Residuals) for P and S phases at each station. The average residuals are calculated for all station/phase combinations that are present in the observations files and for which the corresponding time grids exist, regardless of the weight of the station/phase reading used for location.

  In the Phase Statistics file there are 2 sets of phase statistics. One is based on the residuals for the maximum likelihood hypocenter. The other is based on a PDF-weighted average of the residuals for locations throughout the search region, where the weight is given by the PDF function likelihood values (Grid Search) or the density of samples (Metroplolis).

  The Phase Statistics file also lists a cumulative delay (denoted by Total Phase Corrections) for P and S phases at each station. The cumulative delay is given by the sum of the average residuals and the input delay specified in LOCDELAY statements in the NLLoc Statements section of the Input Control File. The cumulative delay is listed for the maximum likelihood hypocenter residuals and for the PDF-weighted residuals.

  The body of a Phase Statistics file can be used directly as a set of time delay LOCDELAY statements in the Input Control File. Thus the average phase residuals from a run of NLLoc can be used as the station corrections for later runs of NLLoc.

  The list of cumulative delay (denoted by Total Phase Corrections) for P and S phases at each station is also written to a second file with the extension *.stat_totcorr. This simplifies iterative development of station corrections: the *.stat_totcorr file for a run on NLLoc can be included (using the INCLUDE statement) in the control file for the next run, which will then use these total corrections as station delays.

  Line 1: (required) Title and maximum values for P and S average residual calculation.

  Average Phase Residuals (CalcResidual) P_Residual_Max: P_residualMax S_Residual_Max: S_residualMax
  or
  Total Phase Corrections (CalcResidual + InputDelay) P_Residual_Max: P_residualMax S_Residual_Max: S_residualMax
  or
  Average Phase Residuals PDF (CalcPDFResidual) P_Residual_Max: P_residualMax S_Residual_Max: S_residualMax
  or
  Total Phase Corrections PDF (CalcPDFResidual + InputDelay) P_Residual_Max: P_residualMax S_Residual_Max: S_residualMax
  Fields:
  P_residualMax S_residualMax (float)
  maximum (cutoff) value for P and S average residual calculation (set in Input Control File statement LOCPHSTAT). Residuals larger than the cutoff are not used in the calculation.

  Line 2: (required) Column headings.

  ID Phase Nres AveRes

  Lines 3-N: (required) Average residuals.

  LOCDELAY label phaseType nResiduals residual
  Fields:
  label (chars)
  source/station label
  phaseType (chars)
  phase type P or S
  nResiduals (int)
  number of residuals used to calculate average residual
  residual (float)
  average residual

  Example:

  Average Phase Residuals (CalcResidual)  P_Residual_Max: 1.000000  S_Residual_Max: 1.000000
            ID        Phase     Nres    AveRes
  LOCDELAY  BST     P       19       0.033835
  LOCDELAY  BST     S       17       0.032845
  LOCDELAY  CAD     P       22       0.155235
  LOCDELAY  CAD     S       24       0.070664
  Total Phase Corrections (CalcResidual + InputDelay)  P_Residual_Max: 1.000000  S_Residual_Max: 1.000000
            ID        Phase     Nres    TotCorr
  LOCDELAY  BST     P       19       0.033835
  LOCDELAY  BST     S       17       0.032845
  LOCDELAY  CAD     P       22       0.155235
  LOCDELAY  CAD     S       24       0.070664
  Average PDF Phase Residuals (CalcPDFResidual)  P_Residual_Max: 1.000000  S_Residual_Max: 1.000000
            ID        Phase     Nres    AveRes
  LOCDELAY  BST     P       19       0.051917
  LOCDELAY  BST     S       17       0.026063
  LOCDELAY  CAD     P       22       0.151105
  LOCDELAY  CAD     S       24       0.058676
  Total PDF Phase Corrections (CalcPDFResidual + InputDelay)  P_Residual_Max: 1.000000  S_Residual_Max: 1.000000
            ID        Phase     Nres    TotCorr
  LOCDELAY  BST     P       19       0.051917
  LOCDELAY  BST     S       17       0.026063
  LOCDELAY  CAD     P       22       0.151105
  LOCDELAY  CAD     S       24       0.058676
          

• Scatter file (Binary, FileExtension=*.scat)

  The Scatter file contains the x,y,z locations and PDF value of each sample of the location PDF. The number of samples to save is specified in the LOCSEARCH statement in the NLLoc Statements section of the Input Control File.

  Header: (required) one integer and 3 float values

  nSamples dummy dummy dummy
  Fields:
  nSamples (integer)
  number of PDF samples in the following buffer
  dummy (float)
  unused

  Buffer: (required) Sequence of four float values for each PDF sample

  x(N), y(N), z(N), pdf(N) (N = 0, nSamples - 1)
  Fields:
  x(N), y(N), z(N) (float)
  Non-GLOBAL: x, y and z location of the sample in kilometers relative to the geographic origin.
  GLOBAL: longitude and latitude in degrees, and z location in kilometers of the sample.
  pdf(N) (float)
  PDF value of the sample, normalized so that the volume integral over the corresponding search grid of the PDF = 1.0

• Confidence Level file (ASCII, FileExtension=*.conf)

  The Confidence Level file lists the values of the location PDF corresponding to confidence levels from 0.1 to 1.0 in increments of 0.1. Each PDF value denotes a constant-PDF surface in the 3D, x,y,z, PDF function; the integral of the PDF over the volume within this surface will evaluate to the corresponding confidence level. The algorithm to determine the PDF value is described in Moser, van Eck and Nolet (1992).

  Line 1-11: (required) PDF value and confidence level.

  PDFvalue C confidenceLevel
  Fields:
  PDFvalue (float)
  PDF value of the constant PDF surface
  confidenceLevel (float)
  confidence level corresponding to the constant PDF surface defined by PDFvalue
  Notes:
  1. The character C is required in each line but is not currently used.

  Example:

  0.000000 C 1.00
  0.011893 C 0.90
  0.023787 C 0.80
  0.040438 C 0.70
  0.054710 C 0.60
  0.073740 C 0.50
  0.095148 C 0.40
  0.118935 C 0.30
  0.142722 C 0.20
  0.183160 C 0.10
  0.237870 C 0.00

• Quasi-HYPOELLIPSE hypocenter file (ASCII, FileExtension=*.hypo_ell)

  The Quasi-HYPOELLIPSE hypocenter format is an expanded version of the hypocenter summary record output by HYPOELLIPSE (Lahr, 1989). In the original HYPOELLIPSE format there are no spaces between fields and all decimal points are removed; the Quasi-HYPOELLIPSE hypocenter format includes spaces and decimal points, but does not include all fields of the original HYPOELLIPSE format.

  Line 1: (required) Field description.

  DATE ORIGIN LAT LONG DEPTH MAG NO GAP D1 RMS AZ1 DIP1 SE1 AZ2 DIP2 SE2 SE3

  Line 2-N+1: (required) Hypocenter summary record for N events.

  Fields:
  Origin date (Columns 1-8, int*8) (yyyymmdd)
  origin year (including century), month, day
  Origin hour/minute (Columns 10-13, int*4) (hhmm)
  origin hour, minute
  Origin seconds (Columns 15-19, float*5.2) (dd.dd)
  origin seconds
  LatitudeDeg (Columns 21-23, int*3) (ddd)
  latitude degrees
  LatitudeNS (Columns 25, char*1) (A)
  latitude direction N or S
  LatitudeMin (Columns 27-31, float*5.2) (dd.dd)
  latitude decimal minutes
  LongitudeDeg (Columns 33-36, int*4) (dddd)
  longitude degrees
  LongitudeEW (Columns 38, char*1) (A)
  longitude direction W or E
  LongitudeMin (Columns 40-44, float*5.2) (dd.dd)
  longitude decimal minutes
  Depth (Columns 46-52, float*7.3) (ddd.ddd)
  depth in kilometers
  Magnitude (Columns 54-57, float*4.2) (d.dd)
  preferred magnitude
  Number Obs (Columns 59-61, int*3) (ddd)
  number of observations used in solution
  Gap (Columns 63-65, int*3) (ddd)
  largest azimuth separation in degrees between stations
  Distance (Columns 67-72, float*6.2) (ddd.dd)
  distance in kilometers to closest station used in solution
  RMS (Columns 74-78, float*5.2) (dd.dd)
  RMS of residuals of observations used in solution
  Azimuth 1 (Columns 80-83, int*4) (dddd)
  azimuth of axis 1 of confidence ellipsoid
  Dip 1 (Columns 85-88, int*4) (dddd)
  dip of axis 1 of confidence ellipsoid
  StdErr 1 (Columns 90-95, float*6.2) (ddd.dd)
  length of semi-axis 1 of confidence ellipsoid
  Azimuth 2 (Columns 97-100, int*4) (dddd)
  azimuth of axis 2 of confidence ellipsoid
  Dip 2 (Columns 102-105, int*4) (dddd)
  dip of axis 2 of confidence ellipsoid
  StdErr 2 (Columns 107-112, float*6.2) (ddd.dd)
  length of semi-axis 2 of confidence ellipsoid
  StdErr 3 (Columns 114-119, float*6.2) (ddd.dd)
  length of semi-axis 3 of confidence ellipsoid
  Example:

  DATE     ORIGIN     LAT         LONG         DEPTH   MAG  NO  GAP D1     RMS   AZ1  DIP1 SE1    AZ2  DIP2 SE2    SE3
  19961019 1136 22.48  44 N 19.60    7 E 15.63   4.500 0.00  16  78   1.24  0.13    0    0   0.00    0    0   0.00   0.00

• HYPO71 Hypocenter/Station file (ASCII, FileExtension=*.hypo_71)

  The HYPO71 Hypocenter format reproduces "hypocenter output" and "station output" produced by the program HYPO71 (Lee and Lahr, 1972). This format includes all of the fields of the original HYPO71 output, but many of the fields are not used and have dummy values such as 0.0.

  Example:

    DATE    ORIGIN    LAT      LONG      DEPTH    MAG NO DM GAP M  RMS  ERH  ERZ Q SQD  ADJ IN NR  AVR  AAR NM AVXM SDXM NF AVFM SDFM I
   940217 1945 10.28 43 49.04   5 47.40   4.54   0.00 10  9 235 0 0.11  1.7  0.7 C A D 0.00  0  0-0.00 0.00  0  0.0  0.0  0  0.0  0.0 0
    STN  DIST AZM AIN PRMK HRMN P-SEC TPOBS TPCAL DLY/H1 P-RES P-WT AMX PRX CALX K XMAG RMK FMP FMAG SRMK S-SEC TSOBS S-RES  S-WT    DT
   GRX    8.8 138 105  PD4 1945 12.43  2.15  2.19  0.00 -0.04  1.31 0.0 0.0 0.00 0 0.00 000 00.0 0.00 ??4 00.00 00.00 00.00   0.0
   GRX    8.8 138 105  S?4 1945 14.28  4.00  4.05  0.00 -0.05  1.28 0.0 0.0 0.00 0 0.00 000 00.0 0.00 ??4 00.00 00.00 00.00   0.0
   BST   11.8 262  88  PD4 1945 13.00  2.72  2.74  0.00 -0.02  0.00 0.0 0.0 0.00 0 0.00 000 00.0 0.00 ??4 00.00 00.00 00.00   0.0
   BST   11.8 262  87  S?4 1945 15.20  4.92  5.06  0.00 -0.15  1.28 0.0 0.0 0.00 0 0.00 000 00.0 0.00 ??4 00.00 00.00 00.00   0.0
   BMT   12.6 221  81  PD4 1945 13.34  3.06  2.87  0.00 0.19   0.90 0.0 0.0 0.00 0 0.00 000 00.0 0.00 ??4 00.00 00.00 00.00   0.0
   BMT   12.6 221  81  S?4 1945 15.74  5.46  5.31  0.00 0.15   0.88 0.0 0.0 0.00 0 0.00 000 00.0 0.00 ??4 00.00 00.00 00.00   0.0
   ESC   13.3 225  77  PD4 1945 13.34  3.06  2.98  0.00 0.08   0.90 0.0 0.0 0.00 0 0.00 000 00.0 0.00 ??4 00.00 00.00 00.00   0.0
   ESC   13.3 225  77  S?4 1945 15.73  5.45  5.51  0.00 -0.07  0.88 0.0 0.0 0.00 0 0.00 000 00.0 0.00 ??4 00.00 00.00 00.00   0.0
   CAD   15.7 187  72  PD4 1945 13.79  3.51  3.40  0.00 0.11   1.31 0.0 0.0 0.00 0 0.00 000 00.0 0.00 ??4 00.00 00.00 00.00   0.0
   CAD   15.7 187  72  S?4 1945 16.45  6.17  6.29  0.00 -0.12  1.28 0.0 0.0 0.00 0 0.00 000 00.0 0.00 ??4 00.00 00.00 00.00   0.0

  Notes:
  1. The AZM field contains the ray take-off azuimuth, not the source to station aziumth. These two values may differ if 3D velocity models are used.
  2. The quality value in the PRMK field contains the integer quality value coressponding to the phase time uncertainty in seconds. This correspondence is obtained by applying the inverse of the mapping specified by the LOCQUAL2ERR statement of the NLLoc control file.

• HypoInverse Archive file (ASCII, FileExtension=*.hypo_inv)

  The HypoInverse Archive file contains a hypocenter summary record followed by optional phase records. The results for multiple events may present in a single file.

  HYPOINVERSE format, described below in a document from Fred Klein:

  -----------------------------------------------------------------------------

                                                 February 13, 1991
  HYPOINVERSE SUMMARY CARD FORMAT
  ------- ------ ------
  Cols.  Format    Data
  _____  ______    ____
  1-10   5I2       Year, month, day, hour and minute.
  11-14  F4.2      Origin time seconds.
  15-16  F2.0      Latitude (deg).
  17     A1        S for south, blank otherwise.
  18-21  F4.2      Latitude (min).
  22-24  F3.0      Longitude (deg).
  25     A1        E for east, blank otherwise.
  26-29  F4.2      Longitude (min).
  30-34  F5.2      Depth (km).
  35-36  F2.1      Primary amplitude magnitude XMAG1.
  37-39  I3        Number of P & S times with final weights greater than 0.1.
  40-42  I3        Maximum azimuthal gap.
  43-45  F3.0      Distance to nearest station (km).
  46-49  F4.2      RMS travel time residual.
  50-52  F3.0      Azimuth of smallest principal error (deg E of N).
  53-54  F2.0      Dip of smallest principal error (deg).
  55-58  F4.2      Magnitude of smallest principal error (km).
  59-61  F3.0      Azimuth of intermediate principal error.
  62-63  F2.0      Dip of intermediate principal error.
  64-67  F4.2      Magnitude of intermediate principal error (km).
  68-69  F2.1      Primary coda magnitude FMAG1.
  70-72  A3        Event location remark. (See table 7 below).
  73-76  F4.2      Magnitude of largest principal error (km).
  77-78  2A1       Auxiliary remarks (See note below).
  79-80  I2        Number of S times with weights greater than 0.1.
  81-84  F4.2      Horizontal error (km).
  85-88  F4.2      Vertical error (km).
  89-90  I2        Number of P first motions.
  91-93  F3.1      Total of amplitude magnitude weights.
  94-96  F3.1      Total of duration magnitude weights.
  97-99  F3.2      Median-absolute-difference of amplitude magnitudes.
  100-02 F3.2      Median-absolute-difference of duration magnitudes.
  103-05 A3        3-letter code of crust and delay model. (See table 8 below).
  106    A1        Crust model type code (H or T).
  107    A1        Most common P & S data source code. (See table 1 below).
  108    A1        Most common FMAG data source code.
  109    A1        Most common XMAG data source code.
  110    A1        Primary coda magnitude type code (from FC1 command).
  111-13 I3        Number of valid P & S readings (assigned weight > 0).
  114    A1        Primary amplitude magnitude type code (from XC1 command).
                 The following are  written  only  if secondary magnitudes are
                 present.  Secondary magnitudes may appear in either position
                 and use the label codes given by the FC2 and XC2 commands.
                 The code L is for ML calculated by the USGS from Berkeley amp-
                 litudes, and B is for ML from the Berkeley catalog.
  115    A1        Secondary magnitude label or type code.
  116-18 F3.2      Secondary magnitude.
  119-21 F3.1      Total of the secondary magnitude weights.
  122    A1        Secondary magnitude label or type code.
  123-25 F3.2      Secondary magnitude.
  126-28 F3.1      Total of the secondary magnitude weights.

  MAGNITUDES

  The magnitudes and their label codes are:

  E   Primary coda magnitude. FMAG of Eaton (1991). Uses all components.
  X   Primary amplitude magnitude. Jerry Eaton's XMAG. Uses all components.
                 Secondary magnitudes:
  L   Local magnitude computed from UC Berkeley Wood Anderson amplitudes.
  B   Local magnitude from UC Berkeley's catalog.
                 Magnitudes no longer used in CALNET catalog:
  T   Lapse time coda magnitude of Michaelson (1990) from high gain verticals.
  Z   Low gain (Z component) magnitude of Hirshorn and Lindh (1989?).

  AUXILIARY EVENT REMARKS (Summary card cols. 77-78)

  Assigned by analyst (col. 77):

  Q  Quarry blast (or NTS explosion)
  N  NTS blast
  F  Felt

  Assigned by HYPOINVERSE (col. 78):

  *  Location had convergence problems such as maximum number of iterations or
     failure to reach a minimum RMS.
  -  Depth was poorly constrained and was held fixed at its current value.

  Example:

  9608011344195144 2727  7E2285 4000 0  6317 45  14  0 0   0  0 0   0 0XXX   0   0   0   0 0
  SURF P?0 9608011344 2857 -19100 3524 S 0  18  0100   0   0 451  000  0    0  0228 0   0   0
  JAUF P?0 9608011344 2929   2100 3581 S 0 -10  0100   0   0 501  000  0    0  0228 0   0   0
  OG30 P?0 9608011344 3069  16100 3795 S 0  -7  0100   0   0 614  000  0    0  0228 0   0   0
  9608020434148944 2619  7E1611  500 0  6212 15   6  0 0   0  0 0   0 0XXX   0   0   0   0 0
   PZZ P?0 9608020434 1800   4100 1999 S 0  -6  0100   0   0 153  000  0    0  0161 0   0   0
   STV P?0 9608020434 1906   7100 2179 S 0   2  0100   0   0 218  000  0    0  0161 0   0   0
   ENR P?0 9608020434 1974   2100 2293 S 0  -8  0100   0   0 263  000  0    0  0161 0   0   0

Back to the NonLinLoc site Home page.

Anthony Lomax