Bash script to plot fracture zone cross-profiles


fzprofiler [ -Hheight[c|i|p] ] [ -Iprofile ] [ -Linc ] [ -Mncols ] [ -Tprefix ] [ -V[level] ] [ -Wwidth,[c|i|p] ]

Note: No space is allowed between the option flag and the associated arguments.


fzprofiler is a script developed as part of the Global Seafloor Fabric and Magnetic Lineation Project [see GSFML for a full description of the project]. It plots the cross-profiles and the information produced by fzanalyzer. All plots are combined into one plot using whatever custom page size is needed to fit the individual plots given their specifications.

Optional Arguments


Sets the plot height of an individual profile. Unless the measure unit is appended the unit is assumed to be whatever the GMT default PROJ_LENGTH_UNIT is currently set to.


By default we plot all the cross-profiles in one stack. To select a single profile only, append the running number of the profile, where 0 is the first profile.


Determines which profiles to plot. By default we plot every profile (inc = 1). Use -L to plot every inc profile instead [1].


Spread profiles across ncols columns. If ncols = 1 then all profiles are stacked vertically in one long panel; if ncols = 2 then we split the profiles evenly between two columns, etc.


Select verbosity level [w]. (See full description) (See cookbook information).


Sets the file name prefix used for all input files as produced by fzanalyzer [default is fztrack]. The files are prefix_cross.txt, prefix_analysis.txt, prefix_par.txt as well as the resampled output from grdtrack which should be called prefix_resampled.txt.


Similarly sets the plot width of an individual profile. For units, see -H.

Plot Features

fzprofiler packs much information into each plot by using different symbols and colors. The cross-profile VGG data are plotted as red circles connected by a faint red dotted line. The central corridor set with fzanalyzer -D is shown in light blue. Black line is the best symmetric trough model; inverted triangle indicates the best location of the FZ (trough) and its error bars reflect the half-width. Blue line is the best asymmetric blend model; circle shows the best location of the maximum VGG slope. Finally, the inverted red triangle and half-width error bar shows the best empirical trough location. Dashed orange line shows the crustal age variations. Profile orientation is indicated by the W-E or S-N letters, with vertical scale shown on the left (VGG) and right (age). Finally, each panel prints the best blend parameters.


To look at the cross-profiles and the best-fit models and optimal FZ locations on a per-profiles basis, with the prefix used previously as “traces”, and placing the profiles into two columns, each 6 inches wide with individual plots 2 inches tall, try:

fzprofiler -Ttraces -W6i -H2i -N2 -V

The final plot will be named prefix_cross.pdf. To see the same profiles in map view, use fzmapper. To plot a synthetic profile, see fzmodeler.

See Also

gmt fzanalyzer, fzblender, mlconverter, fzinformer, fzmodeler, fzmapper, grdtrack


Wessel, P., Matthews, K. J., Müller, R. D., Mazzoni, A., Whittaker, J. M., Myhill, R., Chandler, M. T., 2015, “Semiautomatic fracture zone tracking”, Geochem. Geophys. Geosyst., 16 (7), 2462–2472. https://doi.org/10.1002/2015GC005853.