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 SIMPLEand 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:

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).

The length of the grid in the x direction is (xNum-1)*dx, and similar for y and z.

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-1)*yNum + (Ny-1))*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:

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).

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:
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 NLLoc.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
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 nEvaluatedsmallestNodeSide 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 rms (float)

maximum and minimum weighted misfit on grid

rms (float)

root-mean-square of residuals at maximum likelihood hypocenter

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:

See NLLoc Program Inversion Approach for more information on the calculation of the expectation, covariance and confidence ellipsoid "traditional" Gaussian statistics.

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 latOrig longOrig longOrig firstStdParal firstStdParal secondStdParal 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

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

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 CCW from North

RAz (float*5.1)

Ray take-off azimuth at maximum likelihood hypocenter in degrees CCW 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)

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

Example: Non-GLOBAL
```
NLLOC "./loc/dur.19940217.221644.grid0" "LOCATED" "Location completed."
SIGNATURE "Anthony Lomax - Geosciences Azur    NLLoc:v2.00 04Aug1998 16h40m19"
COMMENT "IPSN Reseau Durance 94-97 Grid/Layer"
GRID  51 51 43  -50.000000 -50.000000 -1.500000  2.000000 2.000000 1.000000 PROB_DENSITY
SEARCH METROPOLIS nSamp 13829 nAcc 10000 nSave 1000 nClip 0 Dstep0 1.774978 Dstep 0.789846
HYPOCENTER  x 24.705017 y -8.665524 z 1.592826  OT 40.998183  ix -1 iy -1 iz -1
GEOGRAPHIC  OT 1994 02 17  22 16 40.998183  Lat 43.671585 Long 6.056350 Depth 1.592826
QUALITY  Pmax 2.370e-02 MFmin 0.865120 MFmax 5.701586 RMS 0.157467 Nphs 10 Gap 329 Dist 18.346664 Mamp -9.90 0 Mdur -9.90 0
STATISTICS  ExpectX 2.6546e+01 Y -7.5794e+00 Z 5.1371e+00  CovXX 1.69e+00 XY 3.78e-01 XZ 9.87e-01 YY 2.08e+00 YZ 1.08e+00 ZZ 9.78e+00 EllAz1  296.7 Dip1  2.7 Len1  2.26e+00 Az2  206.2 Dip2  10.3 Len2  2.69e+00 Len3  5.96e+00
TRANSFORM  LAMBERT RefEllipsoid Clarke-1880  LatOrig 43.750000  LongOrig 5.750000  FirstStdParal 43.199300  SecondStdParal 44.996100
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
GRX    ?    ?    ? P      U 19940217 2216   44.9200 GAU  2.00e-02 -1.00e+00 -1.00e+00 -1.00e+00 >    4.1083   -0.1865    1.3136    9.1262    1.0246   -0.3350   18.3467 301.88 301.9  90.9  7
GRX    ?    ?    ? S      ? 19940217 2216   48.6900 GAU  4.00e-02 -1.00e+00 -1.00e+00 -1.00e+00 >    7.6074    0.0844    1.2757    9.1262    1.0246   -0.3350   18.3467 301.88 301.9  90.5  6
CAD    ?    ?    ? P      D 19940217 2216   46.3500 GAU  2.00e-02 -1.00e+00 -1.00e+00 -1.00e+00 >    5.1083    0.2435    1.3136    1.3977   -8.1653   -0.3700   23.3127 271.23 271.2  59.3  8
CAD    ?    ?    ? S      ? 19940217 2216   50.4000 GAU  4.00e-02 -1.00e+00 -1.00e+00 -1.00e+00 >    9.4494   -0.0476    1.2757    1.3977   -8.1653   -0.3700   23.3127 271.23 271.2  58.4  8
BMT    ?    ?    ? P      U 19940217 2216   47.3500 GAU  2.00e-02 -1.00e+00 -1.00e+00 -1.00e+00 >    6.3464    0.0054    0.8953   -5.0083   -2.0904   -0.4300   30.4321 282.48 282.5  50.9 10
BMT    ?    ?    ? S      ? 19940217 2216   52.8700 GAU  4.00e-02 -1.00e+00 -1.00e+00 -1.00e+00 >   11.7356    0.1362    0.8776   -5.0083   -2.0904   -0.4300   30.4321 282.48 282.5  50.7 10
ESC    ?    ?    ? P      D 19940217 2216   47.4700 GAU  2.00e-02 -1.00e+00 -1.00e+00 -1.00e+00 >    6.5539   -0.0821    0.8953   -6.2166   -1.8857   -0.4050   31.6561 282.37 282.4  50.9 10
ESC    ?    ?    ? S      ? 19940217 2216   52.8100 GAU  4.00e-02 -1.00e+00 -1.00e+00 -1.00e+00 >   12.1189   -0.3071    0.8776   -6.2166   -1.8857   -0.4050   31.6561 282.37 282.4  50.7 10
BST    ?    ?    ? P      D 19940217 2216   48.0000 GAU  1.00e+05 -1.00e+00 -1.00e+00 -1.00e+00 >    7.3481   -0.3463    0.0000   -8.4624    5.8008   -0.5000   36.1850 293.56 293.6  50.9 10
BST    ?    ?    ? S      ? 19940217 2216   54.6600 GAU  4.00e-02 -1.00e+00 -1.00e+00 -1.00e+00 >   13.5860    0.0758    1.2757   -8.4624    5.8008   -0.5000   36.1850 293.56 293.6  50.7 10
END_PHASE
END_NLLOC
```
Example: GLOBAL
```
NLLOC "/home/anthony/nlloc_tmp/neic/loc/EDT_0.20040317.052139.grid0" "LOCATED" "Location completed."
SIGNATURE "Anthony Lomax   NLLoc:v3.03 17Mar2004 13h39m17"
COMMENT "NEIC eq locations (IASPEI91 model)"
GRID  361 181 1001  -180.000000 -90.000000 0.000000  1.000000 1.000000 1.000000 PROB_DENSITY
SEARCH OCTREE nInitial 6480 nEvaluated 25648 smallestNodeSide 1.085069/1.085069/0.097656
HYPOCENTER  x 23.530273 y 34.536133 z 0.048828  OT 57.293902  ix -1 iy -1 iz -1
GEOGRAPHIC  OT 2004 03 17  05 20 57.293902  Lat 34.536133 Long 23.530273 Depth 0.048828
QUALITY  Pmax 1.686e-11 MFmin 0.498056 MFmax 1.581059 RMS 1.934640 Nphs 100 Gap 70 Dist 2.417298 Mamp -9.90 0 Mdur -9.90 0
VPVSRATIO  VpVsRatio -1.000e+00  Npair 0  Diff -2.000e+30
STATISTICS  ExpectX 2.3533e+01 Y 3.4567e+01 Z 6.4279e+00  CovXX 6.08e+01 XY 5.60e-02 XZ 9.58e-01 YY 1.24e+02 YZ 9.48e+00 ZZ 3.54e+01 EllAz1  198.8 Dip1  83.6 Len1  1.10e+01 Az2  90.4 Dip2  2.0 Len2  1.47e+01 Len3  2.10e+01
STAT_GEOG  ExpectLat 34.567250 Long 23.532505 Depth 6.427921
TRANSFORM  GLOBAL
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
SANT   ?    ?    e Pn     c 20040317 0521   39.9900 GAU  2.00e+00 -1.00e+00 -1.00e+00 -1.00e+00 >   40.7592    1.9369    0.8475   25.4590   36.3710    0.0000    2.4173  40.02   0.0   0.0  0
CSS    ?    ?    e Pn     ? 20040317 0522   55.1300 GAU  2.00e+00 -1.00e+00 -1.00e+00 -1.00e+00 >  118.3789   -0.5428    0.8475   33.3306   34.9622   -0.3960    8.0666  84.18   0.0   0.0  0
CII    ?    ?    e Pn     d 20040317 0523   23.9900 GAU  2.00e+00 -1.00e+00 -1.00e+00 -1.00e+00 >  147.7544   -1.0583    0.8475   14.3050   41.7230   -0.9100   10.2092 317.50   0.0   0.0  0
MALT   ?    ?    e P      c 20040317 0523   53.8200 GAU  2.00e+00 -1.00e+00 -1.00e+00 -1.00e+00 >  179.9217   -3.3956    0.8475   38.4273   38.3134   -1.1200   12.5604  68.16   0.0   0.0  0
PSZ    ?    ?    e P      ? 20040317 0524    9.0000 GAU  2.00e+00 -1.00e+00 -1.00e+00 -1.00e+00 >  195.0132   -3.3071    0.8475   19.8944   47.9184   -0.9400   13.6649 349.63   0.0   0.0  0
VLC    ?    ?    e P      d 20040317 0524   14.5200 GAU  2.00e+00 -1.00e+00 -1.00e+00 -1.00e+00 >  199.2437   -2.0176    0.8475   10.3864   44.1594   -0.5550   13.9749 317.47   0.0   0.0  0
KWP    ?    ?    e P      c 20040317 0524   30.4200 GAU  2.00e+00 -1.00e+00 -1.00e+00 -1.00e+00 >  214.8237   -1.6977    0.8475   22.7075   49.6314   -0.4480   15.1188 357.95   0.0   0.0  0
END_PHASE
END_NLLOC
```
Phase Statistics file (ASCII, FileExtension=*.stat)
The Phase Statistics file 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.

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) four float values

nSamples dummy dummy dummy

Fields:

nSamples (float)

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:

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:

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.

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