R-SCHA software

The Revised Spherical Cap Harmonic Analysis (R-SCHA) is a regional modeling technique that can account for magnetic field measurements taken at various altitudes. It provides also a regional power spectrum that can be interpreted in the same way as a global Spherical Harmonic power spectrum. This power spectrum can be used to provide an estimate for the mean magnetization and mean magnetic thickness over the region of interest.

Here I provide the necessary matlab scripts to produce a R-SCHA model, its power spectrum and the estimates of the mean magnetization and mean magnetic thickness. The matlab scripts rely on a number of fortran functions for which I provide the linux compiled mex files.

I have developed these scripts during my PhD studies. They are heavily based on software I obtained from my PhD advisor Dr. Erwan Thébault, CNRS researcher at the Laboratoire Magma et Volcans, Université Clermont Auvergne, France.

If you are not familiar with the R-SCHA technique and are interested in using these scripts, I strongly advise you to read the bibliography suggested in each of the following sections.

A) If you have a set of magnetic field measurements, taken at various altitudes, and you want to calculate a regional R-SCHA model, click here for instructions and the suggested bibliography.

If the signal of your study averages to zero within the spherical cap, you can use the R-SCHA coefficients obtained to calculate the regional surface R-SCHA power spectrum of your model. For instructions and the suggested bibliography click here.

B) If you have a set of magnetic field measurements taken at a single altitude, and the signal under study averages to zero within the spherical cap, and you only want to estimate the R-SCHA surface power spectrum of your measurements, click here for instructions and the suggested bibliography.

C) If you want to use the R-SCHA surface power spectrum to calculate estimates of the mean magnetic thickness and mean magnetization over the spherical cap, click here for instructions and the suggested bibliography.

If you encounter any problem in the use of these scripts, please let me know.

When using any of these scripts, please cite the respective papers mentioned in the "suggested bibliography" sections.

Thank you for your interest in this work!