Supplements modules
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)
gmtgravmag3d(cmd0::String=""; kwargs...)
Compute the gravity/magnetic anomaly of a 3-D body by the method of Okabe.
See full GMT (not the GMT.jl
one) docs at gmtgravmag3d
Parameters
C | density :: [Type => Str | GMTgrid]
Sets body density in SI. Provide either a constant density or a grid with a variable one.
F | track :: [Type => Str | Matrix | GMTdataset]
Provide locations where the anomaly will be computed. Note this option is mutually exclusive with
outgrid
.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 = gmtgravmag3d(....) form.
H | mag_params :: [Type => Number]
Sets parameters for computation of magnetic anomaly. Alternatively, provide a magnetic intensity grid.
I | inc | increment | spacing :: [Type => Str] $Arg = xinc[unit][+e|n][/yinc[unit][+e|n]]]$
x_inc [and optionally y_inc] is the grid spacing. Optionally, append an increment unit.
L | z_obs | observation_level :: [Type => Number]
Sets level of observation [Default = 0]. That is the height (z) at which anomalies are computed.
M | body :: [Type => Str | Tuple]
Create geometric bodies and compute their grav/mag effect.
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.
S | radius :: [Type => Number]
Set search radius in km (valid only in the two grids mode OR when
thickness
) [Default = 30 km].Tv | index :: [Type => Str]
Tr | raw_triang :: [Type => Str]
Ts | stl :: [Type => Str]
Gives names of a xyz and vertex (ndex="vert_file") files defining a close surface.
Z | z_level | reference_level :: [Type => Number]
Level of reference plane [Default = 0].
Example
G = gmtgravmag3d(M=(shape=:prism, params=(1,1,1,5)), inc=1.0, region="-15/15/-15/15", mag_params="10/60/10/-10/40");
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 conventionN | 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...)_
gravfft(cmd0::String="", arg1=nothing, arg2=nothing; kwargs...)
Spectral calculations of gravity, isostasy, admittance, and coherence for grids.
See full GMT (not the GMT.jl
one) docs at gravfft
Parameters
D | density :: [Type => Str | GMTgrid]
Sets body density in SI. Provide either a constant density or a grid with a variable one.
E | n_terms :: [Type => Number]
Number of terms used in Parker expansion [Default = 3].
F | geopotential :: [Type => Str | Tuple]
Specify desired geopotential field: compute geoid rather than gravity
I | admittance :: [Type => Number]
Use ingrid2 and ingrid1 (a grid with topography/bathymetry) to estimate admittance|coherence and return a GMTdataset.
N | inquire :: [Type => Str] $Arg = [a|f|m|r|s|nx/ny][+a|[+d|h|l][+e|n|m][+twidth][+v][+w[suffix]][+z[p]]$
Choose or inquire about suitable grid dimensions for FFT and set optional parameters. Control the FFT dimension:
Q | flex_topo | flexural_topography :: [Type => Bool]
Computes grid with the flexural topography.
S | subplate | subplate_load :: [Type => Bool]
Computes predicted gravity or geoid grid due to a subplate load produced by the current bathymetry and the theoretical model.
T | topo_load :: [Type => Str]
Compute the isostatic compensation from the topography load (input grid file) on an elastic plate of thickness
te
.W | z_obs | observation_level :: [Type => Number]
Set water depth (or observation height) relative to topography in meters [0]. Append k to indicate km.
Z | moho_depth :: [Type => Number]
Moho [and swell] average compensation depths (in meters positive down).
V | verbose :: [Type => Bool or Str] $Arg = [level]$
Select verbosity level, which will send progress reports to stderr.
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).
Example. Compute the gravity effect of the Gorringe bank.
G = grdcut("@earth_relief_10m", region=(-12.5,-10,35.5,37.5));
G2 = gravfft(G, density=1700, F=(faa=6,slab=4), f=:g);
imshow(G2)
grdgravmag3d(cmd0::String="", arg1=nothing, arg2=nothing; kwargs...)
Compute the gravity/magnetic anomaly of the volume contained between a surface provided by one grid and a plane, or between a top and a bottom surface provided by two grids.
See full GMT (not the GMT.jl
one) docs at grdgravmag3d
Parameters
C | density :: [Type => Str | GMTgrid]
Sets body density in SI. Provide either a constant density or a grid with a variable one.
F | track :: [Type => Str | Matrix | GMTdataset]
Provide locations where the anomaly will be computed. Note this option is mutually exclusive with
outgrid
.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 = grdgravmag3d(....) form.
E | thickness :: [Type => Number]
Provide the layer thickness in m [Default = 500 m].
H | mag_params :: [Type => Number]
Sets parameters for computation of magnetic anomaly. Alternatively, provide a magnetic intensity grid.
I | inc | increment | spacing :: [Type => Str] $Arg = xinc[unit][+e|n][/yinc[unit][+e|n]]]$
x_inc [and optionally y_inc] is the grid spacing. Optionally, append an increment unit.
L | z_obs | observation_level :: [Type => Number]
Sets level of observation [Default = 0]. That is the height (z) at which anomalies are computed.
Q | pad :: [Type => Number]
Extend the domain of computation with respect to output
region
.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.
S | radius :: [Type => Number]
Set search radius in km (valid only in the two grids mode OR when
thickness
) [Default = 30 km].Z | z_level | reference_level :: [Type => Number]
Level of reference plane [Default = 0].
V | verbose :: [Type => Bool or Str] $Arg = [level]$
Select verbosity level, which will send progress reports to stderr.
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).
x | cores | n_threads :: [Type => Str or Number] $Arg = [[-]n]$
Limit the number of cores to be used in any OpenMP-enabled multi-threaded algorithms. (http://docs.generic-mapping-tools.org/latest/gmt.html#x-full)
Example. Compute the gravity effect of the Gorringe bank.
G = grdgravmag3d("@earth_relief_10m", region=(-12.5,-10,35.5,37.5), density=2700, inc=0.05, pad=0.5, z_level=:bottom, f=:g);
imshow(G)
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)
These docs were autogenerated using GMT: v1.20.0