coupe(cmd0::String="", arg1=nothing; kwargs...)

Plot cross-sections of focal mechanisms.

See full GMT (not the GMT.jl one) docs at pscoupe. Essentially the same as meca plus A. Run gmthelp(coupe) to see the list of options.

earthtide(cmd0::String=""; kwargs...)

Compute grids or time-series of solid Earth tides.

See full GMT (not the GMT.jl one) docs at earthtide

	G = earthtide();
	imshow(G)

gmtisf(cmd0::String; kwargs...)

Read seismicity data in the a ISF formated file.

Parameters

• R | region | limits :: [Type => Str or list or GMTgrid|image] $Arg = (xmin,xmax,ymin,ymax)$

  Specify the region of interest. Set to data minimum BoundinBox if not provided.

• D | date :: date="datestart[/dateend]"

  Limit the output to data >= datestart, or between datestart and dateend. <date> must be in ISO format, e.g, 2000-04-25.

• F | focal :: [Type => Bool or Str or Symbol]

  Select only events that have focal mechanisms. The default is Global CMT convention. Use focal=:a for the AKI convention

• N | notime :: [Type => Bool]

  Do NOT output time information.

• abstime or unixtime :: [Type => Integer]

  Convert the YYYY, MM, DD, HH, MM columns into a unixtime. Default puts it as first column, use abstime=2 to put it as last column.

• yx :: [Type => Str or Bool or []] $Arg = [i|o]$

  Swap 1st and 2nd column on input and/or output.

This module can also be called via gmtread. I.,e. `gmtread("file.isf", opts...)_

grdrotater(cmd0::String="", arg1=nothing; kwargs...)

Takes a geographical grid and reconstructs it given total reconstruction rotations.

See full GMT (not the GMT.jl one) docs at grdrotater

Parameters

• A | rot_region :: [Type => Str | Tuple | Vec]

  Specify directly the region of the rotated grid.

• D | rot_outline :: [Type => Bool or Str] $Arg = true | filename$

  Name of the grid polygon outline file. This represents the outline of the grid reconstructed to the specified time.

• F | rot_polyg | rot_polygon :: [Type => Str | GMTdaset | Mx2 array] $Arg = filename | dataset)$

  Specify a multisegment closed polygon file that describes the inside area of the grid that should be rotated.

• G | save | outgrid | outfile :: [Type => Str]

  Output grid file name. Note that this is optional and to be used only when saving the result directly on disk. Otherwise, just use the G = grdrotater(....) form.

• R | region | limits :: [Type => Str or list or GMTgrid|image] $Arg = (xmin,xmax,ymin,ymax)$

  Specify the region of interest. Set to data minimum BoundinBox if not provided.

• T | ages :: [Type => Str | Tuple]

  Sets the desired reconstruction times. For a single time append the desired time. (http://docs.generic-mapping-tools.org/latest/grdrotater.html#t)

• V | verbose :: [Type => Bool or Str] $Arg = [level]$

  Select verbosity level, which will send progress reports to stderr.

• b | binary :: [Type => Str]

• d | nodata :: [Type => Str or Number] $Arg = [i|o]nodata$

  Control how user-coded missing data values are translated to official NaN values in GMT.

• f | geog | colinfo | coltypes | coltype :: [Type => Str] $Arg = [i|o]colinfo$

  Specify the data types of input and/or output columns (time or geographical data).

• g | gap :: [Type => Str] $Arg = [a]x|y|d|X|Y|D|[col]z[+|-]gap[u]$

  Examine the spacing between consecutive data points in order to impose breaks in the line.

• h | header :: [Type => Str] $Arg = [i|o][n][+c][+d][+rremark][+ttitle]$

  Primary input file(s) has header record(s).

• n | interp | interpolation :: [Type => Str] $Arg = [b|c|l|n][+a][+bBC][+c][+tthreshold]$

  Select grid interpolation mode by adding b for B-spline smoothing, c for bicubic interpolation, l for bilinear interpolation, or n for nearest-neighbor value.

• o | outcols | outcol :: [Type => Str] $Arg = cols[,…]$

  Select specific data columns for primary output, in arbitrary order.

Example

G = grdmath("-R-5/5/-5/5 -I0.1 -fg X Y HYPOT");
tri = [-2.411 -1.629; -0.124 2.601; 2.201 -1.629; -2.410 -1.629];
Gr, tri_rot = grdrotater(G, rotation="-40.8/32.8/-12.9", rot_outline=true, rot_polygon=tri);
imshow(Gr, plot=(data=tri_rot,))

meca(cmd0::String="", arg1=nothing; kwargs...)

Plot focal mechanisms.

See full GMT (not the GMT.jl one) docs at meca

Parameters

• J | proj | projection :: [Type => String]

  Select map projection. Defaults to 15x10 cm with linear (non-projected) maps.

• R | region | limits :: [Type => Str or list or GMTgrid|image] $Arg = (xmin,xmax,ymin,ymax)$

  Specify the region of interest. Set to data minimum BoundinBox if not provided.

• B | frame | axis | xaxis yaxis:: [Type => Str]

  Set map boundary frame and axes attributes.

• A | offset :: [Type => Bool | Str | GMTcpt]

  Offsets focal mechanisms to the longitude, latitude specified in the last two columns of the input

• C | color | cmap :: [Type => Number | Str | GMTcpt]

  Give a CPT and let compressive part color be determined by the z-value in the third column.

• D | depth_limits :: [Type => Str | Tuple]

  Plots events between depmin and depmax.

• E | fill_extensive | extensionfill :: [Type => Str | Number]

  Selects filling of extensive quadrants. [Default is white].

• Fa | Fe | Fg | Fo | Fp | Fr | Ft | Fz :: [Type => ]

  Sets one or more attributes.

• G | fill | compressionfill :: [Type => Str | Number]

  Selects shade, color or pattern for filling the sectors [Default is no fill].

• P | portrait :: [Type => Bool or []]

  Tell GMT to NOT draw in portrait mode (that is, make a Landscape plot)

• L | outline_pen | pen_outline :: [Type => Str | Number | Tuple]

  Draws the “beach ball” outline with pen attributes instead of with the default pen set by pen

• M | same_size | samesize :: [Type => Bool]

  Use the same size for any magnitude. Size is given with S

• N | no_clip | noclip :: [Type => Str | []]

  Do NOT skip symbols that fall outside frame boundary.

• Sc|aki | Sc|CMT|gcmt | Sm|mt|moment_tensor | ... :: [Type => Str]

  Selects the meaning of the columns in the input data.

• **convention=:Sa|:aki|:Sc|:CMT|:gcmt|:Sm|:mt|:momenttensor|:Sd|:mtclosest|:momentclosest|:Sz|:mtdeviatoric :momentdeviatoric|:Sp :partial|:Sx|:principal|:principalaxis|:Sy|:principalclosest|:St|:principaldeviatoric

  Alternative way of selecting the meaning of the columns in the input data.

• T | nodal :: [Type => Number | Str]

  Plots the nodal planes and outlines the bubble which is transparent.

• W | pen :: [Type => Str | Tuple]

  Set pen attributes for all lines and the outline of symbols.

• U | time_stamp | timestamp :: [Type => Str or Bool or []] $Arg = [[just]/dx/dy/][c|label]$

  Draw GMT time stamp logo on plot.

• V | verbose :: [Type => Bool or Str] $Arg = [level]$

  Select verbosity level, which will send progress reports to stderr.

• X | x_offset | xshift :: [Type => Str] $Arg = [a|c|f|r][x-shift[u]]$

• Y | y_offset | yshift :: [Type => Str] $Arg = [a|c|f|r][y-shift[u]]$

  Shift plot origin relative to the current origin by (x-shift,y-shift) and optionally append the length unit (c, i, or p).

• di | nodata_in :: [Type => Str or Number] $Arg = nodata$

  Examine all input columns and if any item equals nodata we interpret this value as a missing data item and substitute the value NaN.

• e | pattern | find :: [Type => Str] $Arg = [~]”pattern” | -e[~]/regexp/[i]$

  Only accept ASCII data records that contains the specified pattern.

• h | header :: [Type => Str] $Arg = [i|o][n][+c][+d][+rremark][+ttitle]$

  Primary input file(s) has header record(s).

• i | incols | incol :: [Type => Str] $Arg = cols[+l][+sscale][+ooffset][,…]$

  Select specific data columns for primary input, in arbitrary order.

• p | view | perspective :: [Type => Str or List] $Arg = [x|y|z]azim[/elev[/zlevel]][+wlon0/lat0[/z0]][+vx0/y0]$

  Selects perspective view and sets the azimuth and elevation of the viewpoint [180/90].

• t | alpha | transparency :: [Type => Str] $Arg = transp$

  Set PDF transparency level for an overlay, in (0-100] percent range. [Default is 0, i.e., opaque].

• yx :: [Type => Str or Bool or []] $Arg = [i|o]$

  Swap 1st and 2nd column on input and/or output.

Example: Plot a focal mechanism using the Aki & Richards convention

    psmeca([0.0 3.0 0.0 0 45 90 5 0 0], aki=true, fill=:black, region=(-1,4,0,6), proj=:Merc, show=1)

The same but add a Label

    psmeca(mat2ds([0.0 3.0 0.0 0 45 90 5 0 0], ["Thrust"]), aki=true, fill=:black, region=(-1,4,0,6), proj=:Merc, show=1)

velo(cmd0::String="", arg1=nothing; kwargs...)

Plot velocity vectors, crosses, and wedges.

See full GMT (not the GMT.jl one) docs at velo

    velo(mat2ds([0. -8 0 0 4 6 0.5; -8 5 3 3 0 0 0.5], ["4x6", "3x3"]), pen=(0.6,:red), fill_wedges=:green, outlines=true, Se="0.2/0.39/18", arrow="0.3c+p1p+e+gred", region=(-15,10,-10,10), show=1)