LiMS Downlodable Codes Page: version 5.1
Version 5.1 released July 2004.
This version is Y2K compliant. The major change is in the date format,
which was "yyddmmhhmmss" (char*12) but is now "yyyy-dd-mm hh:mm:ss"
(char*19). A simple code is provided to update your windows files
from the old format to the new one (wv3v4tov5).
Additional features have been added to control the processing
through the process configuration file (process.cfg), including:
REMOTE : E or H whether to use the remote E or H as remote
COHMIN : 20*0.5 minimum coherence to use when restacking
A number of Unix scripts have been added, including ones for
multiple remote-reference processing (lims_process_all) and for
choosing to process a list of times (lims_process_times_all).
Latest version build date: 20 JULY 2004
The LIMS processing codes consist of two main programs, and other ancillary
codes for conversion, plotting, etc. The codes read in data in "ts" format.
This is a general format (I hope), so you can reformat your time series files
into ts format and use these codes.
To see the time series format for input to LiMS v4.x and higher, please read
time series format documentation
The time series for each window from a site
are read into tscascade, which performs cascade
decimation, and writes the decimate estimates of spectra out (to a .cas file).
The cascade decimates files from each window are then read by tsrestack,
and computations performed of MT impedance and GDS estimates using the
Jones/Joedicke jackknife selection scheme with a choice of weights.
For each ensemble average estimate, if there are more than four individual
estimates that are acceptible, then the errors are computed using a jackknife
method. Estimates and error estimates for MT and GDS are performed
independently, and so the resulting estimates cannot be obtained from a
single spectral matrix. Spectra for both MT and GDS estimates are written
As well as the ts file, the programs need some other files to run.
One is the windows file, which contains the parameters of the
windows from each site.
The other file is a hardware file, that contains information
regarding the systems used for calibration (filter coefficients etc.).
More information can be obtained by reading the
Local Operating System.....
There should only be one file that needs to be edited for your local
operating system. This is the include file "op_sys.inc". It sets UNIX
to .TRUE. or .FALSE. depending on whether you are on a Unix box or
In order to ensure the same binaries on both DOS and UNIX platforms,
I have made use of preprocessor options. These are "DOS" and "SUN"
respectively. For the PC you must compile with DOS set as your preprocessor
keyword. Virtual Fortran makefiles are available on request.
Most programs use unix system call "fdate" to get the local time and date.
Program x2edi uses unix system call "idate" to get the current date and
You can either replace these with an appropriate call to your
operating system, or comment out those lines.
Downloadable source codes.....
If you download source codes, please inform Alan Jones so that you
will be placed on a list for automatic notification of updates and
bug reports, etc.
The two files contain all the source codes necessary.
For UNIX users, there is an "install" command which installs the libraries,
executables, and scripts in the relevant directories.
For now, PC users must build the executables manually. If you use
Visual Fortran, the .dsw file is available.
If the download will take too long, contact Alan Jones for a floppy disk by mail.
Version date: 20 JULY 2004
The compressed files contain the following routines...
- tscascade: Derives cascade decimates from a time series for a single
window, and performs initial processing. Outputs a
- tsrestack: Stacks together cascade decimates read in as various
- tsplot: Plots a TS file. Uses either UNIRAS, or routines to write
PostScript, or plots to a PC screen.
- tsactivity: Determines activity level of the data and writes out in
- ts2asc: Converts a BINARY TS file into an ASCII file
- ts2birrp: Converts a TS file (or files for remote-reference
processing) into input files for Chave's BIRRP code
- birrp2j: Converts BIRRP outfile files to J-format
- tscombine: Combines components from two TS files
- tsextract: Extracts a portion of a TS file
- tsresamp: Resamples a TS file from one digitizing rate to a lower
one (set for LiMS standard rates of 64 Hz, 16 Hz, 4 Hz,
1 Hz/1 s, 5 s, 30 s)
- tssplice: Splices two TS files together
- casinfo: Gives information about the contents of a cascade file
- casread: Writes out the cascade decimate spectral
estimates (optionally headers only)
- cas2tab: Computes individual estimates of impedance tensors at a
given period from all available cascade spectral estimates.
Writes estimates out in table format for plotting with e.g.
- mp2ts:Converts a LiMS binary file (in .1mp format) to .ts format
- x2edi:Converts a spectral file output from tsrestack into an EDI
format file for input into e.g. Geotools
It is recommended that j2edi be used to convert the
responses in J-format to EDI-format instead of x2edi
- Utility routines for calculating the time length of a window etc.
- Routines for converting MARIA data and Canadian Geomagnetic Observatory
into TS format (maria2ts & obs2ts respectively)
- Unix C-shell cripts for running all of these codes from a menu-based
- Conversion codes to convert Canadian Observatory data (obs2ts)
and CANOPUS MARIA (maria2ts) data into TS format.