Strike Code Page: version 6.1

McNeice, G.W. and A.G. Jones, 2001, Multisite, multifrequency tensor decomposition of magnetotelluric data. Geophysics, 66, 158-173.PDF

For more information and discussion, contact Alan Jones at Strike Decomp.

It is based on Groom-Bailey tensor decomposition. It currently solves for the electric effects of the galvanic-charges only. Tests with including magnetic effects (see McNeice and Jones, 2001) show the inherent instability of trying to solve for so many parameters.

Documentation can be read here. An example input data file from a single station can be downloaded from here, and an output file for one-decade-width bands from here.

The input data must be in either J-format or EDI-format. The program will either read in a single data file (for single-site mode), or read in a list of data files (for multi-site mode). The user specifies frequency bands to operate on, and an overlap. Accordingly, mixed datasets are permitted with data at different frequencies at different sites.

IMPORTANT: In NAG Mark 19 the main minimization routine we use E04UPF was replaced by E04UNF.
According to the NAG web site, the only change that has to be made is as follows:

E04UPF
Withdrawn at Mark 19

     Old: CALL E04UPF(M,N,NCLIN,LDA,LDCJ,LDFJ,LDR,A,BL,BU,
         +            CONFUN,OBJFUN,ITER,ISTATE,C,CJAC,F,FJAC,
         +            CLAMDA,OBJF,R,X,IWORK,LIWORK,WORK,LWORK,
         +            IUSER,USER,IFAIL)
     New: CALL E04UNF(M,N,NCLIN,LDA,LDCJ,LDFJ,LDR,A,BL,BU,Y,
         +            CONFUN,OBJFUN,ITER,ISTATE,C,CJAC,F,FJAC,
         +            CLAMDA,OBJF,R,X,IWORK,LIWORK,WORK,LWORK,
         +            IUSER,USER,IFAIL)

E04UNF contains one additional parameter as follows:

     Y(M) - real array.

Note that a call to E04UPF is the same as a call to E04UNF with Y(i) = 0.0, for i = 1,2,...,M.

but I have not tested this. If you use Mark 19 and test this replacement, please let me know that it works.

We have attempted to find alternative freeware routines for the minimization search, but we have found the NAG routines to be stable and fast compared to any others. Also, the NAG routine is the only one we know of that permits the user to supply the Jacobian and the Hermetian matrices, as well as have bounds on the model parameters.
Note that we are using a routine that is no longer available in the current versions of NAG releases.
[We would be willing to work with someone who is prepared to test out other routines so we would no longer be dependent on NAG.]

We ask you to do this not to restrict use of strike in any way, but to ensure that each user acknowledges the terms for its use, that we know who is using it so you will be informed immediately of bugs and new releases, and to try to ensure that it is freely available to those who use it for free. Commercial use is certainly encouraged - please contact Alan Jones to discuss licencing terms.

Note: Each person who wished to use strike must send in a scanned signed document agreeing to these terms. That includes all students as well as their supervisors.

A list of official strike users can be viewed here. If you are using strike, and are not on that list, then you must register your usage by sending a scanned signed form to Alan Jones as described above.

Conditions for the use of the McNeice-Jones multi-site, multi-frequency magnetotelluric tensor decomposition code.

1. I accept the McNeice-Jones multi-site, multi-frequency magnetotelluric tensor decomposition code, hereinafter called "strike", on a caveat emptor basis.
2. I will use strike for non-profit purposes only.
3. I will not accept any payment for use of strike.
4. I will not give strike to any other person.
5. I will inform the authors of strike, Gary McNeice and Alan Jones, of any coding errors that I find.
6. I will inform the authors of strike, Gary McNeice and Alan Jones, of any improvements and additions that I make.
7. I will acknowledge the use of strike in any publications and presentations I write or give in which I present results based on strike-based decomposition or decomposed data.
8. I understand that the current version of strike, v6.1, uses NAG v16 routines that I have to supply myself or have to replace with newer versions.

Signed:____________________________________ Date:_________________

Name:
Address:
Institute/Dept./Group:
Organization:
Street Address:
City:
Province/State:
Postal Code:
Country:

Email address:

Authors:

---

Alan G Jones / 13 June 2014 / alan-at-cp.dias.ie