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.

New 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: [] [] 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 [].

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 also 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. []. FuncLab is also available with topography and coastline data for use with the Mapping Toolbox [].

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: 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: []. Please note that I’ve stopped working on this code as the original French group have been making strong strides in improving SplitLab.


44 thoughts on “Software”

  1. Hi Rob. Would you like us to put bug reports here?

    • Hi Shahar. This will work for now. Make sure to include the red error information if possible. This medium allows myself and other users to discuss the codes.

  2. SplitLab: I suggest to add the following line to SL_changeServiceDataCenter:

    case ‘’
    config.serviceDataCenter = ‘’;

    And perhaps even change “” to “” in the calling function, since does not exist anymore.

  3. Your current Matlab version (R2015a)
    Funclab 1.7.7

    Adding FuncLab1.7.7 receiver function seismic tools to your MATLAB paths

    You have the Mapping, Signal Processing Toolbox, and Parallel computing toolboxes, which provide full
    access to all the functionality of the FuncLab package.
    Error using textscan
    Invalid file identifier. Use fopen to generate a valid file

    Error in setup_funclab (line 181)
    javaClassPathCell = textscan(javaClasspathFid,’%s’);

  4. Known bug with SplitLab1.2.1 and Matlab versions 2014b and newer: after associating events, if you have the mapping toolbox, the statistics plotting tool will fail due to a change in plotm and its interaction with matlab plotting primitives. To work around, force the lines which test if the mapping toolbox is available to evaluate to false.

    • I’m sorry, how do I “force the lines which test if the mapping toolbox is available to evaluate to false.” I’m not any better with MatLab than I used to be.

  5. Thank you, quite useful.
    I was poking around the SplitLab code to try and make it work, but my extensive knowledge of disp(‘Hello world’) wasn’t enough.

  6. I am a student and have a very basic experience with Linux, C and coding. But i wanna use either PQLX or your station_analysis_tools software for PSD and PDF. Could you help me to make a step by step installation guide so i may be able to use your station_analysis_tools software package in CentOS 6.3 32 bit version.

    Warm regards,
    Abid Ali

    • Biggest thing is to make sure you have FFTW installed. Its free and easy and I highly recommend it. To use the tools, just try them out on the command line by going to the directory and trying a command like:
      > ./compute_psd
      and it’ll give you the usage. Basically you’re going to need SAC formatted files and the SAC Polezero response file.

      • I am quite sorry to say that i tried using command line by ‘compute_psd’ but i did not fully comprehend how to give the input .SAC file see the command below
        >compute_psd [ KS.ADO2..BGE.D.2014.001.000000.SAC]
        but there was no output file i got . If you have some ‘getting started ‘ manual or you can show how to do that job it would be highly appreciable.

        Thanks for your cooperation !

  7. Check out the README file or the seiscode wiki page. But simply, the way it works is you need to create a new text file with specific parameters. Then your first argument on the command line will be the file name of that text file.

  8. Ju Changhui said:

    Hi Rob. I have downloaded SplitLab 1.0.5 from website you offered and found that the ATD example project data has something difficult. The theoretical travel time caculated by taup toolbox doesnot match the actual SKS phase. I thought it may be my fault in the installing and setting berore the process. But I am not sure about it.Could you please do me a favour to solve it ?

    • Hi Ju,I’ve found that sometimes there is a problem with the depth definition – ie some events are in meters and taup assumes km for depth. This could cause a noticeable offset if the event is deep.

      • Changhui Ju said:

        Rob, I checked the event data head file and confirmed that the event depth is written in km style. But still the problem exists. However, the difference between taup and SKS phase of ATD data is not so big but not regular. So I don’t know what causes the difference.

      • How many seconds off are the prediction and the phase? If small, it may just be heterogeneity along the raypath. It may also just be that the event has a weak SKS arrival. However, I can’t be sure as I haven’t yet checked the example data you’ve mentioned.

      • Changhui Ju said:

        The time difference may be about 10 seconds. And sometimes the “A” and “F” given by the author is before the therotical SKS travel time and sometimes after it. I don’t know the reason. So I thought maybe it’s my wrong during the installing or setting before processing.

  9. Sorry I don’t have more chance to look into this right now. You might want to double check the SAC manual to find out the meaning of those header variables and their relation to timing of the seismogram. Try adjusting them either in sac or matlab and then testing out Splitlab’s and Taup’s capabilities.

  10. CHIH-WEI CHENG said:

    Hi Rob. I am a student and I am trying to work out funclab fetch traces and create project. I’ve tried several computers and different matlab version but that’s not the solution. Below is the bug in matlab:
    Matlab: Unable to create directory C:\Program Files\MATLAB\R2009b\toolbox\lib\
    Please add matTaup.jar and IRIS-WS jar file to your local path and update
    your classpath.txt file.
    funclab:??? Error using ==> taupTime at 117
    Java exception occurred! Please check model name.

    Error in ==> funclab>FL_new_fetch_project_setup_Callback at 3138
    tt = taupTime(‘iasp91′,eventDepth,traceStartPhase,’deg’,gcarc);

    ??? Error while evaluating uicontrol Callback
    Could you please do me a favour to solve it ?

  11. Michael said:

    Hey Rob,
    during searching for updates regarding SplitLab I found your extended version.
    I think for all users it would be nice to add the fixed determination of the degrees of freedom for the SC error estimation published by Walsh et al. (2013).

    Walsh, E., R. Arnold and M. K. Savage (2013), “Silver and Chan Revisited”, JGR


    • Additionally I found a potential error in function preSplit.m, lines 143-145, where the original T component is used to calculate the errors instead of the resulting T component after linearisation. Can you please cross-check this?


      • Thanks for the comments Michael. I haven’t had much chance to look into this, but will try to look into soon. I’ve got some other ideas to improve SplitLab too, just been busy recently with other projects.

  12. Kasemsak Saetang said:

    Dear Rab,
    I am using Funlab 1.7.7 for H-k and CCP stackings beneath Thailand (Lat 8-20, Long 98-105). I am not sure that DNA-13.mat could be used in CCP stacking for an area beneath Thailand.
    Kasemsak Saetang

  13. Kasemsak Saetang said:

    Dear Rob,
    I am using Funclab 1.7.7 for H-k and CCP stacking in Thailand (Lat 8-20, Long 98-105). I am not sure that DNA-13 tomographic model could be used for Thailand.
    Kasemsak Saetang

    • Hi, I’d suggest either finding and adding a model more appropriate for Thailand, or use a global model, such as GyPSuM which is included.

      • Kasemsak Saetang said:

        I just completely created a model beneath Thailand by using CCPV2 with the GyPSuM tomographic model. I don’t know how can I find a method used in CCPV2 for writing a scientific paper. Please give me more information about CCPV2.
        Kasemsak Saetang

  14. Hey Rob, I am using Funclab 1.7.10 for H-k stacking in Caspian Sea. The Matlab version that I use is 2016a. I get the error when I am trying to create the project:
    Error using taupTime (line 117)
    Java exception occurred! Please check model name.

    Error in funclab>FL_new_fetch_project_setup_Callback (line 3244)
    tt = taupTime(‘iasp91′,eventDepth,traceStartPhase,’deg’,gcarc);

    Could you please send a version that works with 2016a or maybe an older version where this is not an issue?

    • Hi Elshan,
      The problem is with the install of MatTaup. As far as I’ve found, there is not an easy way to install java packages with an automated installer.

  15. Yusuf Haidar Ali said:

    Hello Rob

    Is this Funclab working under matlab for windows ? I study about receiver function in Java Island, Indonesia, right now. Could you send me the manual user for your Funclab version ? Thank you


    • Hi Yusuf, the code should work under Windows, however a few file system errors may occur. It was designed on mac and linux. There is no manual yet, but I’m working on making it a bit more self-documenting. For instance, version 1.8.0 (just uploaded) has a menu with seminal citations for reference.

  16. Hello Rob,

    We are a group of students who are now studying geophysics (especially seismology) at Academia Sinica Taiwan. We are also trying to learn/practice writing program for our work. Could you have a visit in our website?

    Thank you and wish you a nice day.

  17. Martin Cárdenas said:

    Dear Rob,

    Could you facilitate me a gmt script to plot probability density spectra functions (general_pdf output)? Many thanks

    • Say the output is tmp.prob and you’re using GMT5, then:
      gmt makecpt -Cgray -I -D -Z > tmp.cpt
      gmt psxy tmp.prob -R${xmin}/${xmax}/${ymin}/${ymax} -JX7i/7i -Bxa${xannot}f${xannot} -Bya${yannot}f${ytick} -BnSeW -Ss0.1i -Ctmp.cpt > tmp.prob.eps

  18. Hi Rob,

    I am trying to use CCP Stacking V2 using my own 3D model. However, when I run the Ray Tracing I get the following error:

    Error using griddedInterpolant
    Interpolation requires at least two sample points in each dimension.

    Error in interp1 (line 161)
    F = griddedInterpolant(X,V,method);

    Error in ccpv2_do_ray_tracing_serial_PtoS (line 396)
    interpEposition = interp1(deps{1},epos{1},InterpDepthVector,’linear’,’extrap’);

    Error in ccpv2_raytrace>saveraytrace_params_and_go_Callback (line 287)

    The size of deps{1} and epos{1} vectors varies between ~1000 to 1 for different stations. As expected, I get the error above when their size is 1. I couldn’t figure out what deps{1} and epos{1} are. The code in “ccpv2_do_ray_tracing_serial_PtoS” is a bit cryptic. I was hoping you might know.

    Thanks for the help!

  19. Hello Rob,

    Great program! I’ve found it very useful.

    Out of curiosity, why does it save the data twice? Once in /RAWTRACES, and once in /RAWDATA/RAWTRACES. It seems like this is a waste of computing time and storage space.

    Just curious!

    • The original version of FuncLab makes a copy of the data into RAWDATA/ from wherever you have the data already stored. That subroutine is called again at the end of the fetching, so that copy command is executed. More importantly, I’ve decided to keep this system to make sure the user has a clean copy of the data either as backup or for other usages. It shouldn’t cost much computation time, but the extra storage can be an issue.

      • I’m backing everything up over Google Drive, and it makes a considerable difference if I have 5,000 files as opposed to 2,500 files. Do you know if there’s an easy way to disable the double back up? If not I’m sure I can find a workaround. It would just be nice to be able to do it by altering the source code.

  20. Shahar B. said:

    Funclab 1.8.0. – feature suggestion
    Currently, the auto trace editing options can be applied either to a selected table or on all the tables. In the latter option, it seems to compare traces from different stations. It would be useful to have an option to do (in one click) all stations separately.

    Rob, thanks again for the good work you are putting into Funclab!

  21. Hey Rob,
    Like your program.
    I am running it Funclab1.8 on Windows machine and tried to make a new project “from exisiting RFs, but run into following error:

    Error using dataread
    Delimiter has bad \ constant.

    Error in strread (line 48)
    [varargout{1:nlhs}]=dataread(‘string’,varargin{:}); %#ok

    Error in funclab>new_project_setup_Callback (line 1547)
    DirName = strread(DataDirectory,’%s’,’delimiter’,filesep);

    Error while evaluating UIControl Callback

    I have preprocessed the files according to the Funclab Manual and used your file/folder name convention. Maybe you might know.

    Thanks for your help!

  22. Hello Rob,
    In order to use funclab functions, I added the paths to the files in the MATLAB session. When I run the setup.m in Matlab, it showing:
    Warning: Unable to save path to file ‘/usr/local/MATLAB/R2015b/toolbox/local/pathdef.m’. You can save
    your path to a different location by calling SAVEPATH with an input argument that specifies the full
    path. For MATLAB to use that path in future sessions, save the path to ‘pathdef.m’ in your MATLAB
    startup folder.
    > In savepath (line 169)
    In setup_funclab (line 378)

    Please help me to set up the Funclab as I am a beginner.

Leave a Reply

Fill in your details below or click an icon to log in: Logo

You are commenting using your account. Log Out / Change )

Twitter picture

You are commenting using your Twitter account. Log Out / Change )

Facebook photo

You are commenting using your Facebook account. Log Out / Change )

Google+ photo

You are commenting using your Google+ account. Log Out / Change )

Connecting to %s