There are mainly two types of files to convert, parameter files with readings and related parameters and binary waveform files.
|CAT_AGA:||Records the S-file header lines according to agency|
|EDRNOR:||Converts USGS monthly bulletins (EDR files) to Nordic format|
|GIINOR:||Converts from Geophysical Institute of Israel parameter format to Nordic|
|HARNOR:||Converts standard Harvard CMT solutions to Nordic format|
|HYPNOR:||Converts from Hypo71 readings files to Nordic format files|
|HINNOR:||Similar to HYPNOR for Hypoinverse files, archive format|
|HSUMNOR:||Converts from Hypo71 summary file format to SEISAN format|
|ISCNOR:||Converts from ISC 96 column format to Nordic format|
|ISCCSV2NOR:||Converts ISC csv foromat catalog to nordic|
|ISCSTA:||Converts ISC station list to SEISAN station list selecting specific stations.|
|ISFNOR:||Converts between ISF1.0 and Nordic|
|KINNOR:||Converts from Kinemetrics to NORDIC|
|NORCSV:||Converts from NORDIC to CSV format using MAGNITUDE_ORDER given|
|in SEISAN.DEF to setup priority of magnitudes.|
|NORGSE:||Converts between Nordic format and GSE parametric format|
|NORHIN:||Converts from Nordic format to Hypoinverse archive format|
|NORIMS:||Converts from Nordic to and from IMS1.0|
|NORHYP:||Converts from Nordic to HYPO71 format|
|PDENOR:||Converts a PDE bulletin file to NORDIC format|
|RSANOR:||Converts Andalucian Seismic Network data to NORDIC format|
|SEIGMT:||Converts from NORDIC file to input for GMT|
|SELMAP:||Selects out a part of a MAP file, also creates file for SeismicityViewer|
|STASEI:||USGS station file or ISC station file to SEISAN|
|USGSNOR:||USGS/NEIC CDROM catalog conversion to NORDIC format|
CAT_AGA, reordering of CAT file header lines
When plotting hypocenters or doing seismic hazard work, it is the first header line in an S-file or CAT-file that is used since it is assumed that it is the prime estimate. When making compact files it is also the first header line, which is used. However, there can be a need for resorting the many type 1 header lines for one or several events so that they are ordered according to agency. It could e.g. be needed to put priority on all the ISC solutions, which then should be the first line in the file. CAT_AGA will reorder the type 1 lines in a CAT file according to the order in which the agencies (3 character codes) are given by the user. If there are many agencies, they can be given in an input file named cat_aga.par, format is one agency per line in the first 3 columns. If the file is not present, the program will ask the user to enter the agencies manually. Optionally, also hypocenter type M lines can be used as type 1 lines. The output file cat_aga.out will contain the sorted events.
EDRNOR: USGS monthly bulletins (EDR files) to Nordic format
Program to convert USGS weekly EDR files (ftp://hazards.cr.usgs.gov/weekly/mchedr*) to Nordic format. The program is written by Mohammad Raeesi (email Mohammad.Raeesi@student.uib.no).
ISCCSV2NOR: ISC CSV catalog to Nordic format
Program converts csv ascii format to nordic format, csv is used by ISC for catalog data. Only the first 3 magnitudes are used, if a magnitude type in input is b1, output type is x.
GIINOR, Geophysical Institute of Israel to SEISAN
The input files are the bulletin type files.
HARNOR, Harvard to Nordic
The standard moment tensor solutions given by Harvard (http://www.globalcmt.org/CMTsearch.html) are converted to Nordic format. Strike, dip, rake and moment tensor solution is written out. The programs can use 3 different input formats:
Standard format: Default format on screen
Full format: Full format on screen
ndk format: File format for downloaded file
The screen formats have the disadvantage that only a limited number of events can be downloaded and captured in one screen so if a lot of events are required, the best alternative is to download a file with all events in ndk format, and then, after conversion, select the desired events.
The ndk format and the full format both give both the hypocenter calculated with arrival times (main header line) and the centroid hypocenter and origin time (MT line), while the standard format does not give the hypocenter calculated with arrival times, so the first header line gives the centroid location. For that reason, it is recommended to use the full format if a screen format is used.
The ndk format does not give Mw which is then calculated from the moment. All formats have additional information not carried over to SEISAN.
HINOR, Hypoinvers archive format to Nordic
The input files are the archive type. For details, see HYPINV program.
HYPNOR, converting HYPO71 files to Nordic files
Input is just filename of HYPO71 file. A similar program for HYPOINVERSE files is HINNOR.
HINNOR, converts from Hypoinverse to NORDIC format
This program works like HYPNOR.
HSUMNOR, HYPO71 summary file format to NORDIC format
Note that the program only converts to header lines.
ISCNOR, converting ISC bulletin file to Nordic format
This program works with the ISC fixed 96-column format as e.g. distributed on CDROM (files of type FFB). The program can select out subsets of ISC data using a latitude-longitude window, depth and prime magnitude. Any of the magnitudes Ms and mb are used. Before 1978, there was only mb on the CD's. More detailed selection can be done on the output file later with SELECT. Since the amount of data is very large it is also possible to write out only the hypocenters. Note that ISC now writes in ISF format also, which can be converted with ISFNOR.
Newer CD's have compressed data and cannot be used directly. files must be copied to disk first, decompressed and then handled as single files.
The program will first check if a file with agency codes called agency.isc is present. If so the station codes are read from this file (same format as files on CDROM). The program will also check the beginning of the data input file for a possible list of agencies and station coordinates. If present, the stations coordinates are read and converted to SEISAN format and additional codes read in. The agency codes are needed in order to identify in plain text the various agencies used.
Principles in conversion:
|Phases:||The phases out can be either the phase ID's sent to ISC or the ISC reinterpreted phases (given with a number code in the input file). If the user supplied phases are used, parenthesizes are removed, and if P/PKP etc is given, it is replaced by P.|
|Times:||If day is incremented relative to origin time day, it is carried into the hours, which can be more than 24.|
|Agency:||It is assumed that it is the same agency for hypocenter and first magnitude. Magnitude is checked for agency, if blank, it is assumed also to be the same as for hypocenter. Only first 3 characters of code is used.|
|Stations:||Only first 4 characters of code are used.|
|Depth:||If no error on depth, a depth fix flag is set.|
|First motion:||Only C or D are used, ISC codes J and B are ignored.|
|Hypocenter orders:||ISC put the best solution last, here the order is reversed, and the prime estimate is first.|
|Duration magnitude:||Change D to C for type.|
|Distance indicator:||If station furthest away is less than 1000 km indicator is L, between 1000 and 3000 km indicator is R and if more than 3000 km indicator is D. If no stations are present the type is set to D.|
In order to relocate an event and compare to ISC location, the ISC reidentified phases must be used (option 2, see below). This has the disadvantage that phases not used by ISC (mainly S-phases of local earthquakes) are weighted out in the output file. If option 3 is used, the ISC identified phases are selected if there and if no ISC identification is given, the local reported phase is used. The output file for option 2 and 3 looks the same except that for option 2, the user-defined phases are weighted out. The residuals given in the output file are always relative to the ISC identified phases.
Below is an example of a run where a latitude - longitude window has been used.
Phases selected can be: User reported phases (default=return) : 1 ISC identified phases only : 2 ISC identified phases and user reported phases when not identified by ISC : 3 3 Output: All hypocenters and phases : Return All hypocenters : 1 Prime hypocenter : 2 Latitude range, return for all 60.2 70.5 Longitude range, return for all 10, 20 Depth range, return for all Magnitude range, return for all Write selected events on screen (y/n=return) No agency.isc file present If ISC CDROM, give drive letter, else return d Give first year and month, e.g. 199501 198601 Give last year and month, e.g. 199602 198602 Now opening d:\1986\198601.FFB Number of agencies in input file 244 etc, for each month 493 events converted Output file name is iscnor.out File with stations is isc.sta
The file input can be from a CDROM as in the example above. In that case, the whole CDROM can be read or a smaller time interval can be given. The input can also be from a single file and the program will then ask for the next file when the first has been converted. If many files are to be converted, a list of file names can be made with DIRF and filenr.lis entered as an input file name. The Nordic format output file is iscnor.out and the station list is in isc.sta which has the format used by SEISAN. Optionally, output can also be in the original isc format, however that requires setting a flag in the program and recompiling, see program source code.
ISCSTA, selecting stations in the complete ISC station file
The complete station list in the ISC list is very large and it is often an advantage to use a smaller subset, although HYP can use the whole list. The program can select out subsets of stations in both SEISAN and an older ISC format. The program will read an S-file, find how many different stations there are and select those stations out of a station file, which can either be in SEISAN (=HYPO71) format or ISC format (automatically determined). The output is in SEISAN format. If no S-file is given the input station file is assumed to be in ISC format and the whole file will be converted to SEISAN format.
KINNOR, Kinemetrics to NORDIC
Converts .PCK file output of EDPPICK to file in SEISAN format. Many events are converted from one file. The program is based on program from Kinemetrics by Christopher S. Lim. For info on how conversion is made, see program source code.
ISFNOR, ISF1.0 to and from Nordic
The ISF format is used by the ISC and is an extension to the IMS format. The program is based on the routines provided by the ISC for reading and writing ISF, and the SEISAN standard routines for reading and writing Nordic data. The program converts in both directions. All possible information is converted. Information on the ISF format can be found on the ISC website (http://www.isc.ac.uk). It is recommended to use ISF format for data exchange with ISC.
NORIMS, IMS1.0 to NORDIC format.
The IMS1.0 (International Monitoring System) is a new version of the GSE format and very similar. The program can partly be used for the new ISF (IASPEI Seismic Format) which will include all of the IMS format an additional information needed by ISC and NEIC. The program and the following description is by Mario Villagrán. The program works with the IMS1.0:SHORT format (phase-readings/origin files) and the program works both ways.
The IMS1.0:SHORT format is exactly the one used at the IDC International Data Center (Vienna, Austria). In addition some features used by the ISC International Data Center and the Spanish NDC National Data Center had been added. Magnitudes in IMS format use many characters, the Nordic format allows only one; the following rule is followed:
For mb 'b'
For MS 'S'
For ML 'L'
For MD 'C'
For Ml 'l'
For MN 'N'
For mblg 'G'
For ms 's'
For MB 'B'
IMSNOR do not include coda magnitude.
The maximum likelihood magnitudes mb1, mb1mx, ms1, ms1mx, etc are
pending. IDC still does not have documentation and they may be changed.
Single measurements of magnitude/station are parsed as comment lines (type 3) starting with symbol ``$''. When importing data from IMS format, only the ``Event IDC'' number is parsed and included into a comment line (type 3) of Nordic, together with the ellipse dimensions orientation and the mb standard deviation.
All parameter values read that exceed field limits of Nordic (Amplitude, velocity, snr, etc) have been set to the maximum or minimum possible, example: if snr 999.9 then snr=999. For conversion from Nordic to IMS it is necessary to use both the hyp.out and print.out files; The reason is that IMS includes many parameters that need to be searched in both files.
When converting to IMS format, the user can specify the start numbering for the first event and phase in the file; ignoring will assume (1,1). It is optionally also possible to set the no location flag in the output header lines.
Comments in the IMS input phase data is ignored.
It seems that the program works well for converting ISC data. Just capture the output on the ISC screen and use that as input to IMSNOR. No cleaning of the file is needed before conversion.
NORCSV, NORDIC to CSV format
The program converts the date and time, hypocenter location and magnitude to a simple CSV file that can be loaded into spreadsheet software. An example of a line is
2012/01/06,06:16:26.80, 76.191, 9.870, 10.0, 2.3,CBER
where the columns are date, time, latitude, longitude, depth, magnitude and magnitude type. The magnitude with type of the highest priority is selected as given by the magnitude order in SEISAN.DEF. For example ML from your agency may be given priority over MC, or you may prefer your ML over that from another angency.
NORGSE, NORDIC from and to GSE parametric format
The program (written by Mario Villagrán) converts parametric data between Nordic and GSE2 format. It can be used interactively or by giving the options as arguments. Type norgse -help to see the options.
NORHIN, From Nordic to Hypoinverse archive format
The program is started by typing norhin input-file. The output file is norhin.out. For more details, see program HYPINV. indexHYPINV
NORHYP, From Nordic to HYPO71 format (SUN and PC)
The program is written by F. Courboulex. The program asks for the input file name and the output file name is norhyp.out.
NOTE, the program expect there to be a type 7 line for each event, else it will not output anything.
PDENOR, converting PDE bulletin file to NORDIC format
PDE distributes bulletins on e-mail, both a monthly bulletin and a weekly bulletin (different formats). The program converts one of these files to Nordic format and put the file into a standard SEISAN database called PDE for the monthly files and PDEWE for the weekly files. This database must have been created before running the program. Both CAT and S-files are made and SELECT and EEV can be used afterwards. Files can be received by email or picked up at hazards.cr.usgs.gov, dir pde for monthly files and weekly for weekly files, must be ehdf type.
Program converts between format used by ``Red Sismologica de Andalucia'' and a few others in Spain.
SEIGMT, Nordic to GMT input
The program SEIGMT reads information from Nordic or compact files and writes the parametric data to files that can be used as input for GMT(Generic Mapping Tools, http://gmt.soest.hawaii.edu/). The user can choose a scaling for the magnitudes and also select a magnitude type order. The scaling option is useful if you wish to scale the symbol size of your epicenters with magnitude. The magnitude type order defines, which magnitude should be taken in case several magnitudes have been determined for one event. If you don't give a magnitude order, the program chooses the largest magnitude.
The files written by SEIGMT are:
gmtxy.out - event locations, to be plotted with psxy
gmtxyz.out - event locations and depths, to be plotted with psxy
gmtstxy.out - station coordinates (longitude, latitude and station code)
gmtpath.out - travel path data, to be plotted with psxy
psmeca.out - fault plane solutions, to be plotted with psmeca (Aki and Richards convention)
SELMAP, selecting a subsection of a MAP file
The program can retrieve parts of a large MAP file written in SEISAN map format. On the SEISAN web site or on the SEISAN CDROM, very detailed global mapfiles are available in SEISAN format. The file originally comes from the USGS. SELMAP can select out part of a MAP file in a latitude-longitude grid. The MAP files consist of several small segments and a segment is selected if at least one point is inside the specified grid. The program also creates an output file in xyz format for the Lomax SeismiciyViewer. The program can output a smaller number of points than available in the input file by using the parameter skip.
Converts the official global station file from USGS (comma format) or ISC global station file to SEISAN station format (same as HYPO71 format with SEISAN extension for 5 letter station codes). A list of most global stations are now found on the SEISAN CD. It seesm that the USGS format is no longer used.
USGSNOR, USGS catalog to NORDIC format
The program converts USGS CDROM hypocenters to NORDIC format. Most of the information is used. If more than 3 magnitudes are available, only the 3 first are used. The number of stations is included when available. The depth is indicated as fixed in all cases where the operator has been used (A,N,G). Macroseismic information is included with max intensity. The residual standard deviation is put into rms column. Event types are set to R. Magnitude types are converted as follows:
UK is made blank
b is replaced by B
s is replaced by S
D is replaced by C
w is replaced by W
WAVEFORM CONVERSION PROGRAMS
This group of programs are mostly converting waveform files from some format to SEISAN although a few also convert from SEISAN to some other, mostly standard, formats. Most programs convert from binary to binary formats.
Many instruments come with conversion programs to some standard format like PCSUDS or MINISEED, and these have often been used to convert to SEISAN instead of writing programs reading the original files directly. Many such conversion programs work on PC so the corresponding SEISAN programs only work on PC. However, since the PC files can be read directly on Sun, this should not present a problem. Many programs have VERY LITTLE documentation, look in source codes for more information.
The number of programs are forever increasing with new recorders coming onto the market and new formats coming in use and others going out of use and it is becoming increasingly difficult to keep track of it all. For this release of SEISAN it has not been possible to test all programs on all platforms. An attempt has been made to standardize the programs. A general problem is that many seismic recorders and formats do not provide proper identification of the channels. In the worst cases, there are no station codes, only channel numbers and in very many cases, there is no room for proper component information. This is being taken care of by having a definition file, and only one format for the definition file is used, see below. This is also used with program WAVFIX.
Most programs work in the standard way with a filenr.lis file as input (made with DIRF).
The response information is seldom in the original files and in most conversion programs, the response information is taken from the CAL directory. If no response information is available, a message will be given. For each program, a comment will be given as to the status of testing and on which platforms they operate. If the channel definition file option is implemented, the array dimensions will be SEISAN standard.
The program SEIPITSA might be an easy way to convert between 1-column ASCII data and SEISAN (see below).
When converting between the major analysis format (MiniSEED, SEISAN, SAC and GSE) mostly using program WAVETOOL, only SEISAN and MiniSEED will preserve the network and location codes as well as the flag for uncertain timing since the other formats only partly have room for this information.
Conversion programs definition file
The conversion programs use a common format for the definition file for naming station and channels. The definition file is named programname.def as e.g. sudsei.def. The definition file can be in the working directory or the DAT directory. The conversion program will first look in the working directory for the file and then in DAT. The conversion of codes can take place in 2 ways (see below for details): (1) An input station and component code is converted to an output station code and component, (2) an input channel number is assigned a station and component code. The advantage of (1) is that the conversion is independent of the channel number or order, however, the user must then know the default station and component names generated by the conversion program.
Default assignment of station code and component:
This is very much dependent on the conversion program used since some data files have complete
information and others very little, see description of individual programs in manual or at start of source codes. In all cases, the conversion program will make both station and component codes based on what is available of information in the input files. IT IS THESE CODES THAT are used for input code as described below. In order to find out what they are, it is easiest to run the conversion program once (without a def file) and see what codes the program assign. Alternatively, some of the programs have documentation in the manual. Some of the station codes might be instrument serial numbers, which are not always known. Therefore, running a test might be the best way to find out.
In addition to converting channel codes, the def file can also give SEISAN waveform file header information and network code as it appears in the file name. If no network code is given, the network code will be the station code of the first channel.
Principle of conversion in order of precedence:
Header line text (29 char)... NetCd (5 chars), Comment for next line Header for REFTEK NEWNT chan stati comi stato como, In and output definitions, comment for next line 1 BO11 S Z BOM B Z BO12 S N BOM B N BO13 S E BOM B E
The first line is just a comment line, must be there in any format. Here it shows where the network code is positioned as indicated by NetCd.
The second line gives the header information for the SEISAN main header, which are the first 29 characters. The file name network code is also given and is here NEWNT. Format a29,1x,a5.
The third line is just comment to indicate the position of the columns in the following lines (max 200). A line must be there. The abbreviations are:
|chan:||Channel number, optional unless no input station and component given.|
|stati:||Input station code, 1-5 chars|
|comi:||Input component code, 4 characters|
|stato:||Output station code, 1-5 characters|
|como:||Outut component code, 4 characters. First character MUST be S, L, B, A, or I, last character MUST be Z, N or E, all upper case.|
The conversion programs are listed below
|AFADSEI:||AFAD acceleration records to SEISAN|
|AHSEI:||AH ASCII to Seisan, little tested|
|ARCSEI:||Reftek archive to Seisan, windows only|
|ASCSEI:||ASCII to SEISAN|
|BGISEI:||Beijing Geodevice Institue to SEISAN|
|CITSEI:||CityShark to SEISAN|
|CNVSSA:||Kinemetrics SSA to Kinemetrics Dataseis|
|CNVSSR:||Kinemetrics SSR to Kinemetrics Dataseis|
|DELS:||Deletes specific parts of an S-file|
|DIMASSEI:||USGS Dimas to Seisan|
|DRSEI:||Sprengnether recorders to SEISAN|
|GIISEI:||Geophysical Institute of Israel to SEISAN|
|GSRSEI:||GeoSig to SEISAN|
|GSERESP:||Conversion between GSE and SEISAN response files|
|GURSEI:||Güralp to SEISAN|
|GCF2MSD_MANY:||Run Guralp gcf2msd many times from files in a file structure|
|IRISEI:||From IRIS ASCII to SEISAN|
|ISMSEI:||ISMES to SEISAN|
|KACSEI:||Kinemetrics ASCII acceleration to SEISAN|
|KINSEI:||Kinemetrics Dataseis to SEISAN|
|K2SEI:||Kinemetrics K2 to SEISAN|
|LEESEI:||Willy Lee system to SEISAN|
|LENPCQ:||Converts from Lennartz to PCEQ to PCEQ format|
|LENSEI:||Lennarts ASCII to SEISAN|
|M88SEI:||Lennartz MARS88 to SEISAN|
|NANSEI:||Converts from Nanometrics to SEISAN|
|NEISEI:||Converts from NEIC CDROM waveform data to SEISAN|
|NRWSEI:||Geol. Survey. of Northrhine-Westphalia format to SEISAN format|
|OS9SEI:||Converts SEISLOG files to SEISAN waveform files|
|PITSA:||Conversion programs described with program PITSA|
|PCQSEI:||Converts from PCEQ to SEISAN|
|PDASEI:||Geotech Instruments PDAS to SEISAN|
|PSNSEI:||Public Seismic Networks to SEISAN|
|QNXSEI:||SEISLOG QNX to SEISAN|
|RDSEED:||IRIS program to read SEED volumes|
|RSASEI:||Conversion from Andalucian Seismic network to SEISAN|
|RT_SEIS:||Reftek Passcal format to SEISAN conversion|
|SEI2PSXY:||Convert waveform data to trace input for psxy|
|SEIM88A:||Conversion from SEISAN to MARS88 ASCII format|
|SEISAN2MSEED:||From SEISAN to MiniSEED|
|SEISAF:||SEISAN to SESAME ASCII|
|SEIPITSA:||SEISAN PITSA ASCII|
|SGRSEI:||SeisGram to SEISAN|
|SISSEI:||Sismalp format to SEISAN format|
|SILSEI:||SIL network ASCII files to SEISAN|
|SUDSEI:||PCSUDS to SEISAN|
|TERSEI:||Terra ASCII to SEISAN|
|WGSSEI:||WGSN format to SEISAN|
For each program, a summary of capabilities is mentioned: The platforms
available (all for all or specific name), channel definition file
available (chan. def. yes or no) and if the program will look for
response files in the CAL directory to insert in the headers (resp. yes or no).
If you do not find the conversion program here, look on the ORFEUS website for other programs that might convert to one of the formats used above.
AFADSEI, AFAD acceleration to SEISAN or MiniSeed
Converts from AFAD ASCII acceleration to SEISAN or MiniSeed. . Components are hardwired to HN Z etc. Assume 3 channels files only, all channels same sample rate and number of samples. Input values are in gal, the output is in 1/10 000 gal or 0.001 mm/s2, station code is taken from file (a number). There is an output file with hypocenter info from input files as well as the waveform file name. For echa channel a response file is written in SEISAN format.
AHSEI, AH ASCII to SEISAN
Converts AH ASCII files to Seisan format.
ARCSEI, Reftek archive to SEISAN (windows only)
ARCSEI is a program to automate the extraction of data from a RefTek data archive and the conversion to SEISAN format. The program works interactively and with a simple text interface that asks the user to enter the details for the data extract. Based on the user selected criteria the program (1) extracts the data from the archive in Passcal format using ARCFETCH, (2) converts the Passcal data files to SEISAN format using RT_SEIS, and (3) merges the SEISAN files if merging is activated by the user, using SEISEI. The program is written in Fortran and works on Windows only.
The ARCSEI program can be used in various ways:
The ARCFETCH and RT_SEIS programs, both part of the RefTek software package, have to be installed (see RefTek documentation) and the PATH variable set to include the directory where the programs are stored. It is assumed that the RefTek data archive exists and that the user is familiar with the content of the archive. The archive content can be shown with the command ARCINFO. To test that the program is installed correctly, open the Windows command tool (from the menu, or by selecting Start . Run . cmd) and type ARCSEI RETURN .
The definition file: arcsei.def
The purpose of the definition file is to set some parameters needed to run ARCSEI, however, the program also works without. The arcsei.def file can either be stored in the seismo/DAT directory, or the current working directory. The program first checks in the current directory. The arcsei.def file should be adjusted to the user's set-up, before ARCSEI is started.
The parameters are:
|ARCHIVE:||The path of the RefTek data archive, can also be entered manually at run time.|
|OUTPATH:||The directory in which the SEISAN files are to be stored. The default is `.\' (the current directory).|
|MERGE:||Select if SEISAN files from several stations for the same time interval should be merged (Y), or not (N).|
|NETWORK_CODE:||Network code used in case SEISAN files are merged.|
|CHANNEL:||Data channel in RefTek archive consisting of the unit, stream and channel (unit,stream,channel). The * can be used as wildcard to select all streams or channels, BUT not to select all units (since ARCFETCH is used in cooked mode, which means that the time interval extracted matches the input start- and end-time.|
Example of the arcsei.def file
KEYWORD............Comments.............Par1........................ ARCHIVE Path to archive G:\CTBTO\ARCHIVE OUTPATH Path to put C:\test converted SEISAN files MERGE Y or N for Y merging Seisan files NETWORK_CODE used for merging ARCSE CHANNEL Name of channel 8020,1,* CHANNEL Name of channel 8021,1,* CHANNEL Name of channel 8022,1,* CHANNEL Name of channel 8023,1,* CHANNEL Name of channel 8024,1,* CHANNEL Name of channel 8025,1,* CHANNEL Name of channel 8028,1,*
How to run the program:
ARCSEI is started from the Windows Command Tool (cmd) by typing
arcsei <RETURN> ARCSEI - DATA EXTRACTION FROM REFTEK ARCHIVE AND CONVERSION TO SEISAN FORMAT SEISAN FILES ARE MOVED TO DIRECTORY: C:\test ENTER ARCHIVE, OR <RETURN> TO USE DEAULT (default: G:\CTBTO\ARCHIVE) Return to accept default, which is set in the arcsei.def file, or give the archive path. AVAILABLE CHANNELS ARE: 8020,1,* 8021,1,* 8022,1,* 8023,1,* 8024,1,* 8025,1,* 8028,1,* ENTER CHANNEL SELECTION (UNIT,STREAM,CHANNEL) OR TYPE ALL TO USE ALL AVAILABLE CHANNELS
Type channel and RETURN , if defined in arcsei.def channels are listed, otherwise an example is shown. The channel is given as unit,stream,channel. Wildcards can be used for stream and channel, but not for the unit.
NEXT CHANNEL OR RETURN TO CONTINUE
Additional channels can be entered, to continue press RETURN .
ENTER START-TIME (YEAR:DAY-OF-YEAR:HOUR:MINUTE:SECOND) EXAMPLES: 2000:200:12 2000:200:12:15 2000:200:12:33:15
Type start time as year:day-of-year:hour:minute:second. Minute and second can be omitted.
NEXT START-TIME OR RETURN TO CONTINUE
Additional start times can be entered, to continue press RETURN .
ENTER END TIME USING ONE OF 3 OPTIONS: - ABSOLUTE TIME AS YYYY:DDD:HH:MM:SS (LIKE START-TIME) - +SECONDS FOR TIME INTERVAL (e.g. +300) - ++SECONDS FOR MULTIPLE INTERVALS(CONTINUOUS EXTRACT, e.g. ++300)
Specify the end time, either in the same style as for the start time (only if one start time), or in some cases more useful, the desired time window in seconds, by entering +seconds. If sequential time windows are to be extracted, use ++ seconds. The user is then asked how many time windows should be extracted. It is thus possible e.g. to extract 10 consecutive windows of 600 seconds. Only if sequential extract windows specified:
ENTER NUMBER OF CONTINUOUS WINDOWS
After the program has finished, the data in SEISAN format can be found in either the current directory (default) or in the OUTPATH directory if the variable is specified in arcsei.def. The temporarily created files are deleted automatically.
How it works
ARCSEI reads the user input that specifies what should be extracted from the RefTek archive and then calls the programs ARCFETCH, RT_SEIS and SEISEI. For temporary data storage ARCSEI creates the directory arcsei_temp under the current directory. The arcsei_temp directory is automatically deleted upon program completion.
The arcfetch program performs the data extraction from the RefTek archive. A complete list of the command line input of arcfetch can be obtained by starting the program without additional options. ARCSEI starts arcfetch in the following way:
arcfetch archive channel,start-time,end-time -o OUTPATH -c
-o OUTPATH: Specifies the output path for arcfetch, always arcsei_temp
-c: Specifies cooked mode, which means that the time interval extracted matches the input start- and end-time (this is not the case, when not running in cooked mode)
ARCHIVE 8020,1,*,2000:200:12,+10 -oarcsei_temp -c
RT_SEIS converts all files with the suffix `rt' in arcsei_temp to SEISAN format. RT_SEIS reads the RTU.INI file for station definition, if the environmental variable RTU is set to point to the RTU.INI file (see RT_SEIS section below).
SEISEI, if merge is selected, merges all SEISAN files in the arcsei_temp directory.
Finally all files (single or merged) are moved to the OUTPATH directory or the current directory if OUTPATH is not defined. In case multiple stations are selected, ARCSEI runs steps (1) and (2) in a loop, before the data is merged and moved. In case several time windows are selected, ARCSEI runs steps (1) to (4) in a loop, and in addition a second loop over multiple station (1) and (2). If sequential time windows are specified, ARCSEI computes multiple start times and works as if these time windows were user specified. All, def. File yes, resp yes
ASCSEI, ASCII to SEISAN
Converts a single column file to a one channel SEISAN file. Date, time, sample rate (number of samples per second), station and component must be entered manually. The user is ask for a scale factor, normally it is 1.0. If input numbers are smaller than 1.0, a scale factor must be used since numbers are written out as integers. The input file must contain only a column of numbers and no headers.
linux$ head data.txt -1646 -1722 -2138 -942 -734 -272 -652 -1736 -102 755 linux$ ascsei File name, # or filenr.lis for all data.txt data.txt year,month,day,hour,min,sec 2017,02,22,11,30,05.12 sample rate 20 station, max 5 chars JENS component, max 4 chars BH Z Enter scale factor, return for no scaling Number of samples in file: 144132 Outfile name: 2017-02-22-1130-05S.JENS__001_BZ JENS 1117 53 2 22 11 30 5.120 7206.550 No response file for this channel -------- JENS BH Z117 53 2 22 11 30 5.120 20.00 144132
BGISEI, Beijing GEODEVICE FORMAT (BGI) to SEISAN.
Linux, PC, chan. def. yes, resp yes
The program to convert waveform files from BGI to SEISAN format is called BGISEI. The instrument response in the original files is not used. The program has only been tested with data recorded in Cuba. The program is written by Bladimir Moreno.
CITSEI, CityShark to SEISAN
Converts from CityShark ASCII to SEISAN. Components S Z, S N, S E are assumed for the 3 channels. Assume 3 channels files only, all channels same sample rate and number of samples.
CNVSSA and CNVSSR Kinemetrics accelerometers to Kinemetrics Dataseis
The programs are supplied by Kinemetrics to convert from SSA and SSR formats to Kinemetrics Dataseis. To further convert to SEISAN, use program KINSEI. Only PC executable programs are available. The data is 16 bit.
At the moment there is no direct conversion from CSS to SEISAN. It is possible to convert CSS data to SAC or GSE using other tools like codeco, Geotool and SAC, and then convert to SEISAN format.
DELS, deletes parts of S-file
DELS, program to delete specific lines in an S-file with one or many events
There is often a need to delete particular parts of an event file like the spectral lines or the P-phases etc. Below is an example run.
C:\\seismo\\WOR>dels Give input file hyp.out Give line to delete or keep, terminate with enter 1: Lines with P-phase Not delete 2: Lines with S-phase Not delete 3: Lines with SPEC-phase Not delete 4: Lines with IAML-phase Not delete 5: Lines with IASP-AMP phase Not delete 6: Lines with any phase Not delete 7: Lines with fps solutions Not delete
Any combination of options can be selected or deselected. The output is in file dels.out.
A special option is 99. This will write out IAML amplitudes for Z, N and E if all three present. Output file is dels_amp.out. The program can also be used with arguments
First argument is file name
Second argument is optionally 'overwrite' meaning that input file will be overwritten, can only be used with a file with one event, if more program stops. Third argument is operator
The last two arguments are used with EEV, command dels in EEV.
DIMASSEI, USGS DIMAS to SEISAN
Converts Dimas files to Seisan format.
DRSEI, Sprengnether data recorders to SEISAN
Converts Sprengnether DR3024 and DR3016 to SEISAN format. These two formats are slightly different, but the program makes the adjustment. Only essential information is read in and only 4 lowest digits of serial number are used. If station codes are set up, these are used, else the serial numbers are used for station codes.
GIISEI, Geophysical Institute of Israel to SEISAN
Converts Geophysical Institute of Israel imported DAQ files to SEISAN format. The initial station codes are as defined in file, can be converted with the normal .def file. If 4.th character of station name indicates component (N or E), that is blanked out and transferred to 4.th character of component name BEFORE using the def file conversions.
GURSEI, Güralp to SEISAN
Converts Güralp GCF files to SEISAN format, only works with one channel data. Maximum number of samples as defined in seisan, at least 1 200 000, channels codes can be defined using the gursei.def definition file. If no definition file, the station name is the first 4 letters from internal station name and the component is BH Z.
Runs Güralp gcf2msd many times from files in a file structure. It is only tested for one particular structure with 15 min files, so the program has to be adopted to other structures. Requires program gcf2msd from Güralp.
GSERESP, conversion between GSE and SEISAN response files
The program provides conversion between SEISAN, GSE1 and GSE2 response files. The response can be given in frequency, amplitude and phase (FAP) triplets or in poles and zeros (PAZ). Since the number of values in the GSE format is unlimited the conversion from SEISAN to GSE only changes the format, whereas converting from GSE to SEISAN, if the number of FAP triplets is more than 30 or the number of poles and zeros larger than 37, the response in SEISAN format will be approximated by 30 FAP triplets. The output files in SEISAN format will have the default SEISAN response filenames (see RESP program and SEISAN response format). Output files in GSE format will include the station name, the component, number 1 or 2 for GSE1 and GSE2 respectively and end on `.CAL' (e.g. MOR_SHZ2.CAL (GSE2), KONO_BZ_1.CAL (GSE1).
GSRSEI, GeoSig to SEISAN
Converts from GBV recorders to SEISAN. GeoSig was earlier GeoSys. Before version 8.1, there was a bug in program so start time was wrong by the amount of the prevent time.
IRISEI, IRIS ASCII to SEISAN
The input format is the variable ASCII download format used on the GSN Quanterra stations. The format is used in connection with SEISNET. The program only works if input file has more than 1000 samples. ISMSEI, ISMES to SEISAN PC, chan. def. no, resp no ISMES is an Italian seismic recorder. This is the first version of the program made by IIEES in Iran. The program can convert one file with up to 3 channels.
KACSEI: Kinemetrics ASCII acceleration to SEISAN
Kinemetrics ASCII film record acceleration files (type *.v1) to SEISAN. It is assumed that:
KINSEI, Kinemetrics DATASEIS to SEISAN
The program takes the station code from the input files. The component codes are also taken from the input file as far as Z, N and is E is concerned, but the first letter is always set to S, like 'S Z'. The program is also used if CNVSSR or CNVSSA have been used first.
K2SEI, Kinemetrics K2 to SEISAN
Program for K2 binary files. The program works by first converting the binary files to ASCII by internally running the Kinemetrics program kw2asc (PC only). If no definition file is present, channel 1-3 will be A Z, A N and A E. If more channels they will be called A 04, A 05, etc.
LEESEI, Willy Lee binary files to SEISAN
The number of channels is fixed to 16 and the time information is not read, it must be entered when converting the file.
LENSEI, Lennartz ASCII to SEISAN
LENPCQ, converting Lennartz to PCEQ format
Only executable code for this program and only PC (made by the Royal Belgian Observatory). The format is used by an older version Lennartz tape recorder. The output files have the same names as the input files and are placed in a directory c:
qcoda, WHICH MUST BE THERE.
M88SEI, Lennartz MARS88 to SEISAN
NANSEI, Nanometrics to SEISAN
The program converts from the Y-file format to SEISAN. This is done by first making an ASCII file with Nanometrics y5dump program (done internally in NANSEI). NOTE: The y5dump program requires some special Nanometrics libraries (Solaris) or *.DLL files (PC), which are included and installed with SEISAN (see installation section). The program converts single channel files only.
NEISEI, NEIC digital data to SEISAN
NEIC earthquake digital data comes on CDROM. The data is extracted with a program coming with the data and then converted to SEISAN binary waveform data. The response information is given as poles and zeros in the SEISAN waveform file header.
OS9SEI, converting SEISLOG files to SEISAN
The program takes a SEISLOG ASCII (downloaded in CMP6 format) or binary file and converts to a SEISAN file. The input can be several files from a filenr.lis or an ASCII downloaded file either compressed or uncompressed. The program will look for the calibration file in the CAL directory and add it to the SEISAN file, or give a message if it is not there. The program will work with SEISLOG files recorded under operating system OS9 or QNX up to version 7.6. For QNX version 7.0, use program QNXSEI.
PCQSEI, converting PCEQ format to SEISAN
PCEQ format to SEISAN. Earlier used with IASPEI software libraries.
PDASEI, converting PDAS files to SEISAN
The program converts a single channel PDAS file to a single channel file in SEISAN format. Several of these files can then be merged with SEISEI. PDASEI in previous SEISAN versions (before version 6.0) only worked with PDAS in 16-bit format, so if 32 bit or gain ranged format was input, the output would have been in error. The current version of PDASEI should be able to convert all 3 types of input files. A description of the PDAS format is found in the PDASEI program.
PSNSEI, Public Seismic Networks to SEISAN
The Public Seismic Network recording system makes one file pr channel. Since component is not well defined, several files from the same recording system might get the same SEISAN file name. Do some testing when setting up the recording system. The one component files can be assembled into multichannel files with SEISEI. There might be a newer version of PSN format not supported.
QNXSEI, SEISLOG QNX version to SEISAN
This program works as OS9SEI except that it does not read the ASCII files. The program must be used with Seislog 8.0. The program is currently the only program that put in the time synchronization flag in SEISAN waveform files except for data logging programs under Seislog Windows. See format description in Appendix B. The program recalculates the sample rate based on the time in the first blocks in the file and the last blocks in the file (each block is one second long). For very long files, this might be of importance since the digitizer might not have exactly the nominal sample rate.
RSASEI, Andalucian Seismic Network to SEISAN
Conversion of network and broad band files to SEISAN format. Covers several versions of the DTS format also used by other institutions in Spain. Not tested on Linux.
RT_SEIS, Reftek Passcal to SEISAN
The RT_SEIS program converts Reftek Passcal format to SEISAN. This program is provided by Refraction Technology Inc. The program does not use the filenr.lis as input file. To see the options of RT_SEIS, start the program without any arguments. In order to make use of the RTU.INI definition file, the environmental variable RTU needs to be set to for example c:
dat (see RefTek documentation for more details). This file can be used to set station names for respective unit IDs.
Example of RTU.INI:
Converts waveform file to GMT psxy trace plotting ASCII file. The output files have one line for each sample giving the date and time and amplitude value, e.g.:
Each trace will have its own file and the output names are given for each like
Output file name: 2005-07-31-1544-10S.STOK__003_STOK__S__Z
Output file name: 2005-07-31-1544-10S.STOK__003_STOK__S__N
Output file name: 2005-07-31-1544-10S.STOK__003_STOK__S__E
To plot the trace data with psxy, use projection `-JX xsize T ysize ' and option `-R' giving time range in the same style as the data. To plot the data the gmtdefaults should be set to `gmtset INPUT_DATE_FORMAT yyyy/mm/dd INPUT_CLOCK_FORMAT hh:mm:ss.xx'. See psxy man pages for more details.
SeisGram binary to SEISAN. Only 3 component data has been tested. Channel order is assumed to be Z, N, E. The input real values have been multiplied by 100 000 before being converted to integers. Program little tested.
The Standard for Exchange of Earthquake Data (SEED) format is defined by the Federation of Digital Seismographic Networks (FDSN). The rdseed program is distributed with SEISAN to extract data from SEED volumes. RDSEED is an IRIS program to read SEED volumes. The program provides conversions to SAC (ASCII and binary), AH, CSS and miniseed. It is described in the file `rdseed.txt' in the INF directory. Updated versions of rdseed will be available at http://www.iris.washington.edu/pub/programs. A PC version (rdseed.exe) is distributed with SEISAN CD (also on home page). SEED volumes contain the complete response information, details on how to convert the SEED response to GSE response format can be found in Havskov and Alguacil .
SEIM88A, conversion from SEISAN to MARS88 ASCII format
The program converts SEISAN waveform files to Lennartz-ASCII MARS88 format. The program will write one file per channel. Output files are either mars.xxx if a single file is converted or marsxxx.yyy if the `filenr.lis' file is used as input.
The program converts from SEISAN to PITSA ASCII format and back. The ASCII format has one file per channel. The user will be asked for a name of the output file-system. If a single file is converted, the channel number will be added to the output file-system name (e.g. data.001). If the `filenr.lis' file is used the filenumber will be added to the file-system name (e.g. pitsa001.004, first file and fourth channel). The program is no longer used for conversion when PITSA is started from EEV, but might be useful, since it creates one column ASCII data and can easily be modified.
SEISAF, SEISAN to SESAME ASCII
The 3 first channels in SEISAN file are read. There is no check if from same station. It is assumed that the order in SEISAN file is Z,N,E, that all 3 channels have the same start time, number of samples and sample rate. These values are taken from the first trace.
By Chad Trabant, IRIS Data Management Center
Program developed at IRIS to convert from SEISAN to mseed, all platforms and all mseed formats. This program can be used as alternative to converting data with wavetool, advantage is that SEISAN2MSEED supports STEIM2 compression.
Source code can be found at https://seiscode.iris.washington.edu/projects/seisan2mseed
seisan2mseed [options] file1 [file2 file3 ...]
Seisan2mseed converts SeisAn waveform data files to Mini-SEED. One or more input files may be specified on the command line. If an input file name is prefixed with an '@' character or explicitly named 'filenr.lis' the file is assumed to contain a list of input data files, see LIST FILES below. The default translation of SeisAn components to SEED channel codes is as follows: a 3 character SEED channel is composed of the first, second and fourth characters of the component; furthermore if the second character is a space and the first and fourth are not spaces an 'H' is substituted for the 2nd character (i.e. 'S__Z' 'SHZ'). The default SEED location code is '00', if the third character of the SeisAn component is not a space it will be placed in the first character of the SEED location code. Other translations may be explicitly specified using the -T command line option. If the input file name is a standard SeisAn file name the default output file name will be the same with the 'S' at character 19 replaced by an 'M'. Otherwise the output file name is the input file name with a "_MSEED" suffix. The output data may be redirected to a single file or stdout using the -o option.
Print program version and exit.
Print program usage and exit.
Be more verbose. This flag can be used multiple times ("-v -v" or "-vv") for more verbosity.
Include SEED blockette 100 in each output record with the sample rate in floating point format. The basic format for storing sample rates in SEED data records is a rational approximation (numerator/denominator). Precision will be lost if a given sample rate cannot be well approximated. This option should be used in those cases.
Buffer all input data into memory before packing it into Mini-SEED records. The host computer must have enough memory to store all of the data. By default the program will flush it's data buffers after each input block is read. An output file must be specified with the -o option when using this option.
"-n Inetcode P"
Specify the SEED network code to use, if not specified the network code will be blank. It is highly recommended to specify a network code.
"-l Iloccode P"
Specify the SEED location code to use, if not specified the location code will be blank.
"-r Ibytes P"
Specify the Mini-SEED record length in Ibytes P, default is 4096.
"-e Iencoding P"
Specify the Mini-SEED data encoding format, default is 11 (Steim-2 compression). Other supported encoding formats include 10 (Steim-1 compression), 1 (16-bit 3 integers) and 3 (32-bit integers). The 16-bit integers encoding should only be used if all data samples can be represented in 16 bits.
"-b Ibyteorder P"
Specify the Mini-SEED byte order, default is 1 (big-endian or most significant byte first). The other option is 0 (little-endian or least significant byte first). It is highly recommended to always create big-endian SEED.
"-o Ioutfile P"
Write all Mini-SEED records to Ioutfile P, if Ioutfile P is a single dash (-) then all Mini-SEED output will go to stdout. All diagnostic output from the program is written to stderr and should never get mixed with data going to stdout.
"-T Icomp=chan P"
Specify an explicit SeisAn component to SEED channel mapping, this option may be used several times (e.g. "-T SBIZ=SHZ -T SBIN=SHN -T SBIE=SHE"). Spaces in components must be quoted, i.e. "-T 'S Z'=SHZ".
If an input file is prefixed with an '@' character the file is assumed to contain a list of file for input. As a special case an input file named 'filenr.lis' is always assumed to be a list file. Multiple list files can be combined with multiple input files on the command line.
The last, space separated field on each line is assumed to be the file name to be read. This accommodates both simple text, with one file per line, or the formats created by the SeisAn dirf command (filenr.lis).
An example of a simple text list:
An example of an equivalent list in the dirf (filenr.lis) format:
# 1 2003-06-20-0643-41S.EDI___003
# 2 2005-07-23-1452-04S.CER___030
Conversion from the Icelandic SIL system to SEISAN. Only conversion from ASCII files.
SISSEI, Sismalp to SEISAN
The program converts from Sismalp to SEISAN. Sismalp is a French field recording system. The input consists of 2 files pr event, a header file and a data file. It is assumed that the Sismalp ndx files have the same file name as the header file except for the file extension. It is also assumed that the file names are 12 characters long.
SUDSEI, PCSUDS to SEISAN
The program converts from PCSUDS to SEISAN. This is done by first running the program SUD2ASC (included) and then converting to SEISAN. The SUD2ASC program and test data was supplied by REFTEK through the distribution of PC-SUDS Utilities by Banfill .
TERSEI, Terra ASCII to SEISAN
Program converts from Terra Technology ASCII files to SEISAN. Only tested with 1-3 channel files
WGSSEI to SEISAN
Program converts from WGSN files to SEISAN. The format is used on IRIS stations as processing format. Little tested.