Next: 20. Signal processing programs Up: 19. File conversion and modification Previous: 19.1.1 Nordic to QuakeML converter Contents Index
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. |
Format i5,1x,a5,2x,a4,1x,a5,2x,a4
The conversion programs are listed below
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.
(http://orfeus.knmi.nl/other.services/conversion.shtml).
AFADSEI, AFAD acceleration to SEISAN or MiniSeed
all, chan. def. yes, resp yes
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 all, chan. def. yes, resp yes
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
Where:
-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)
Example:
arcfetch G:
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 all, chan. def. yes, resp yes
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.
Example:
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 all, chan. def. yes, resp yes
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 PC
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.
CSS
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 AM, AT and AS phase Not delete 7: Lines with END phase Not delete 8: Lines with BAZ phase Not delete 9: Lines with any phase Not delete 10: Lines with fp solutions Not delete 11: Phases with given stations Not delete
Any combination of options can be selected or deselected. The output is in file dels.out. If used withing EEV, it is used with current event S-file. Note: Option 11 can use several stations and all phases for those stations are deleted, however not SPEC phases.
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 all, chan. def. yes, resp yes
Converts Dimas files to Seisan format.
DRSEI, Sprengnether data recorders to SEISAN
all, chan. def. yes, resp yes
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
all, chan. def. yes, resp yes
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 all, chan. def. yes, resp yes
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.
If a file name indicates that the file is thought to be from AFAD,
no def file is used since it is assumed that intern codes are correct.
.
GCF2MSD_MANY
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 all
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 all, chan. def. yes, resp yes
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 all, chan. def. no, resp yes
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 all, chan. def. yes, resp yes
Kinemetrics ASCII film record acceleration files (type *.v1) to SEISAN. It is assumed that:
| - | channel 1 is N, 2 is Z and 3 is E |
| - | there are always 3 channels in file |
| - | input values are in 1/10 g, the output is in 1/1 000 000 g |
| - | station code is taken from file name as given in first line of input file |
| - | the 3 channels can have different number of samples, however it is assumed that they all start at the same time |
KINSEI, Kinemetrics DATASEIS to SEISAN PC, chan. def. yes, resp yes
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 PC,Linux, chan. def. yes, resp yes
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 PC, chan. def. no, resp no
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 all, chan. def. yes, resp yes
LENPCQ, converting Lennartz to PCEQ format PC
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 all, chan. def. yes, resp yes
NANSEI, Nanometrics to SEISAN or mSEED PC, Sun, chan. def. yes, resp no
The program converts from the Y-file format to SEISAN waveform format or MSEED. 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
(Solarisi, only available in older versions of SEISAN) or *.DLL files (PC), which are included and installed with
SEISAN (see installation section). These libraries are not available un Linux. The program converts single channel files only.
NEISEI, NEIC digital data to SEISAN PC, chan. def. no, resp no
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 PC, SUN, chan. def. no, resp yes
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 PC, chan. def. yes, resp no
PCEQ format to SEISAN. Earlier used with IASPEI software libraries.
PDASEI, converting PDAS files to SEISAN
all, chan. def. yes, resp yes
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 all, chan. def. yes, resp yes
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 all, chan. def. no, resp yes
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 all, chan. def. yes, resp yes
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 PC, chan. def. no, resp no
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:
seismo
dat (see RefTek documentation for more details). This file can be used to set station names for respective unit IDs.
Example of RTU.INI:
[8020]
Station=SB00
Network=CTBTO
CH1Band=
CH1Type=
CH1Axis=a
CH1Loc =
CH2Band=
CH2Type=
CH2Axis=b
CH2Loc=
CH3Band=
CH3Type=
CH3Axis=c
CH3Loc =
[8021]
Station=SB01
Network=CTBTO
CH1Band=
...
SEI2PSXY
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.:
2005/06/16T00:59:59.51 -40.0000
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.
SGRSEI PC, chan. def. yes, resp yes
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.
SEED
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 (2004).
SEIM88A, conversion from SEISAN to MARS88 ASCII format
all, chan. def. no, resp no
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.
SEIPITSA all, chan. def. yes, resp yes
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
all, chan. def. no, resp no
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.
SEISAN2MSEED All chan.def. no resp no
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
SYNOPSIS
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.
OPTIONS
"-V "
Print program version and exit.
"-h "
Print program usage and exit.
"-v "
Be more verbose. This flag can be used multiple times ("-v -v" or "-vv") for more verbosity.
"-S "
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.
"-B "
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".
LIST FILES
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:
2003-06-20-0643-41S.EDI___003
2005-07-23-1452-04S.CER___030
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
SILSEI all, chan. def. no, resp no
Conversion from the Icelandic SIL system to SEISAN. Only conversion from ASCII files.
SISSEI, Sismalp to SEISAN all, chan. def. yes, resp yes
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 PC, chan. def. yes, resp yes
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 (1996).
TERSEI, Terra ASCII to SEISAN all, chan. def. yes, resp yes
Program converts from Terra Technology ASCII files to SEISAN. Only tested with 1-3 channel files
An example of the input format is seen below:
191213 043 55.85 41.4825N 19.6017E 26.19 -99.00 98 0.09 AB18P1 3.99AB26P0 4.07AB06P0 4.08AB30P0 4.22AB10P0 4.74AB10S0 9.01 AB17P0 4.78AB09P0 4.85AB09S1 9.21AB03P1 5.95AB03S2 11.12AB02P0 6.02 AB01P1 6.27AB13P1 7.12 191213 131 28.93 41.4204N 19.5849E 10.39 -99.00 93 0.08 AB26P1 1.91AB18P3 2.04AB30P0 2.05AB06P0 2.18AB09P1 2.22AB10P0 2.26 AB17P0 2.54AB02P1 3.66AB03P1 3.68AB01P1 3.93AB13P2 4.54
WGSSEI to SEISAN all, chan. def. yes, resp yes
Program converts from WGSN files to SEISAN. The format is used on IRIS stations as processing format. Little tested.
WAVASC, ASCII file from any WAV file
Program to convert any seisan waveform file to a 2 column ASCII file. One file will be made per channel as well as one file for all the converted data. The start of the converted file looks as follows:
COCO BH E 10 II 2007 928 1345 31.440 42953 40.000
31.440 1874.00
31.465 1787.00
31.490 2159.00
31.515 3002.00
where the header has station, component, location, network, time, number of samples and sample rate. The following times are referred to the minute in the header.
the channel output filename is like
2007-09-28-1345-31A.COCO_BH_E
where a stands for ASCII
The input can be an s-file where all waveform files in s-file will be used or it can be a waveform file or a filenr.lis of waveform files.