Throughout my PhD work, I’ve built several pieces of research code, mostly in c. I feel writing software is an integral part of any science PhD, because if the necessary software already exists, then the problem is solved and you should work on something else. But also it really helps you think through the problem and better understand why things work or don’t work.
Matlab R2014b has implemented a new graphics program. This has led to a few bits to not work quite properly in SplitLab and FuncLab. Download these versions to avoid these problems: [SplitLab1.2.1.zip] [Funclab1.8.2.zip] Note that these have not been fully tested, so bugs may still exist. FuncLab is also available with topography and coastline data for use with the Mapping Toolbox [Funclab1.8.2_with_map_data.zip].
Along these lines, Matlab R2014b comes with a new default color palette called “Parula”. This is a really fantastic color palette, and if you want to use it in GMT you can download it here: parula.cpt
FuncLab I am working on updating FuncLab to improve functionality and added customization. Interface updates include color palette swapping and new plot tools to look at stacked Receiver Functions. Work in the background include improved version tracking within a given project, depth mapping based on a 1D velocity model, station and event maps without the Matlab Mapping Toolbox, Receiver Function computations within FuncLab, and much more. For a video tutorial of the new FuncLab, check out this YouTube video from the IRIS EPO. The manual has recently been updated for version 1.8.0 and is available here [funclab_manual_1.8.0.pdf]. This project is not yet complete, so some of the new functions are not yet available. [FuncLab1.8.2.zip]. FuncLab is also available with topography and coastline data for use with the Mapping Toolbox [Funclab1.8.2_with_map_data.zip].
If you use FuncLab, please cite our paper: Porritt, R. W. and Miller, M. S., (2018), Updates to FuncLab, a Matlab based GUI for handling receiver functions. Computers and Geoscience, 111, 260-271, doi:10.1016/j.cageo.2017.11.022
Michigan receiver function data example: [here]
Station Analysis Tools One useful pack I wrote is a set of c routines for computation of power spectral densities, coherence, probability density functions, and a handful of other tools for monitoring the health of a station. These are effectively like PQLX or IRIS’s Quack, but are designed for very simple inline processing with sac data files and associated pole-zero response files. I’ve packaged these tools with the standard unix method, so just download the tarball, untar, configure, and make and you should be ready to go. Includes example data files, but requires the proper installation of fftw3 and pkgconfig. [station_analysis_tools]
[Old] SplitLab1.2.1 SplitLab is the standard shear-wave splitting analysis tool built in Matlab. It was written in 2006 by Andreas Wüstefeld and was due for an update. The original program and a splitting analysis database can be found here: http://splitting.gm.univ-montp2.fr/. I’ve updated it to utilize irisFetch.m for obtaining event and waveform data, removed the dependence on the mapping toolbox, and created a new output called ‘station.mat’ which contains the full misfit space for each splitting measurement. Download here: [SplitLab1.2.1.zip]. Please note that I’ve stopped working on this code as the original French group have been making strong strides in improving SplitLab.