.\" Man page generated from reStructuredText. . .TH "PSWIGGLE" "1gmt" "Nov 05, 2016" "5.3.1" "GMT" .SH NAME pswiggle \- Plot z = f(x,y) anomalies along tracks . .nr rst2man-indent-level 0 . .de1 rstReportMargin \\$1 \\n[an-margin] level \\n[rst2man-indent-level] level margin: \\n[rst2man-indent\\n[rst2man-indent-level]] - \\n[rst2man-indent0] \\n[rst2man-indent1] \\n[rst2man-indent2] .. .de1 INDENT .\" .rstReportMargin pre: . RS \\$1 . nr rst2man-indent\\n[rst2man-indent-level] \\n[an-margin] . nr rst2man-indent-level +1 .\" .rstReportMargin post: .. .de UNINDENT . RE .\" indent \\n[an-margin] .\" old: \\n[rst2man-indent\\n[rst2man-indent-level]] .nr rst2man-indent-level -1 .\" new: \\n[rst2man-indent\\n[rst2man-indent-level]] .in \\n[rst2man-indent\\n[rst2man-indent-level]]u .. .SH SYNOPSIS .sp \fBpswiggle\fP [ \fItable\fP ] \fB\-J\fP\fIparameters\fP \fB\-R\fP\fIwest\fP/\fIeast\fP/\fIsouth\fP/\fInorth\fP[/\fIzmin\fP/\fIzmax\fP][\fBr\fP] \fB\-Z\fP\fIscale\fP [ \fB\-A\fP\fIazimuth\fP ] [ \fB\-B\fP[\fBp\fP|\fBs\fP]\fIparameters\fP ] [ \fB\-C\fP\fIcenter\fP ] [ \fB\-G\fP[\fB+\fP|\fB\-\fP|\fB=\fP]\fIfill\fP ] [ \fB\-I\fP\fIfix_az\fP ] [ \fB\-K\fP ] [ \fB\-O\fP ] [ \fB\-P\fP ] [ \fB\-S\fP[\fBx\fP]\fIlon0\fP/\fIlat0\fP/\fIlength\fP[/\fIunits\fP] ] [ \fB\-T\fP\fIpen\fP ] [ \fB\-U\fP[\fIstamp\fP] ] [ \fB\-V\fP[\fIlevel\fP] ] [ \fB\-W\fP\fIpen\fP ] [ \fB\-X\fP\fIx_offset\fP ] [ \fB\-Y\fP\fIy_offset\fP ] [ \fB\-bi\fPbinary ] [ \fB\-c\fP\fIcopies\fP ] [ \fB\-di\fPnodata ] [ \fB\-f\fPflags ] [ \fB\-g\fPgaps ] [ \fB\-h\fPheaders ] [ \fB\-i\fPflags ] [ \fB\-p\fPflags ] [ \fB\-t\fPtransp ] [ \fB\-:\fP[\fBi\fP|\fBo\fP] ] .sp \fBNote:\fP No space is allowed between the option flag and the associated arguments. .SH DESCRIPTION .sp \fBpswiggle\fP reads (\fIx\fP,\fIy\fP,\fIz\fP) triplets from files [or standard input] and plots z as a function of distance along track. This means that two consecutive (\fIx\fP,\fIy\fP) points define the local distance axis, and the local \fIz\fP axis is then perpendicular to the distance axis. The user may set a preferred positive anomaly plot direction, and if the positive normal is outside the plus/minus 90 degree window around the preferred direction, then 180 degrees are added to the direction. Either the positive or the negative wiggle may be shaded. The resulting PostScript code is written to standard output. .SH REQUIRED ARGUMENTS .INDENT 0.0 .TP \fB\-J\fP\fIparameters\fP (more ...) Select map projection. .UNINDENT .INDENT 0.0 .TP \fB\-R\fP[\fIunit\fP]\fIxmin\fP/\fIxmax\fP/\fIymin\fP/\fIymax\fP[\fBr\fP] (more ...) Specify the region of interest. .UNINDENT .sp For perspective view \fBp\fP, optionally append /\fIzmin\fP/\fIzmax\fP\&. (more ...) .INDENT 0.0 .TP \fB\-Z\fP\fIscale\fP Gives anomaly scale in data\-units/distance\-unit. .UNINDENT .SH OPTIONAL ARGUMENTS .INDENT 0.0 .TP .B \fItable\fP One or more ASCII (or binary, see \fB\-bi\fP[\fIncols\fP][\fItype\fP]) data table file(s) holding a number of data columns. If no tables are given then we read from standard input. .UNINDENT .INDENT 0.0 .TP \fB\-A\fP\fIazimuth\fP Sets the preferred positive azimuth. Positive wiggles will "gravitate" towards that direction. .UNINDENT .INDENT 0.0 .TP \fB\-B\fP[\fBp\fP|\fBs\fP]\fIparameters\fP (more ...) Set map boundary intervals. .UNINDENT .INDENT 0.0 .TP \fB\-C\fP\fIcenter\fP Subtract \fIcenter\fP from the data set before plotting [0]. .UNINDENT .INDENT 0.0 .TP \fB\-G\fP[\fB+\fP|\fB\-\fP|\fB=\fP]\fIfill\fP Set fill shade, color or pattern for positive and/or negative wiggles [Default is no fill]. Optionally, prepend \fB+\fP to fill positive areas (this is the default behavior). Prepend \fB\-\fP to fill negative areas. Prepend \fB=\fP to fill both positive and negative areas with the same fill. .UNINDENT .INDENT 0.0 .TP \fB\-I\fP\fIfix_az\fP Set a fixed azimuth projection for wiggles [Default uses track azimuth, but see \fB\-A\fP]. .UNINDENT .INDENT 0.0 .TP \fB\-Jz\fP|\fBZ\fP\fIparameters\fP (more ...) Set z\-axis scaling; same syntax as \fB\-Jx\fP\&. .UNINDENT .INDENT 0.0 .TP \fB\-K\fP (more ...) Do not finalize the PostScript plot. .UNINDENT .INDENT 0.0 .TP \fB\-O\fP (more ...) Append to existing PostScript plot. .UNINDENT .INDENT 0.0 .TP \fB\-P\fP (more ...) Select "Portrait" plot orientation. .UNINDENT .INDENT 0.0 .TP \fB\-S\fP[\fBx\fP]\fIlon0\fP/\fIlat0\fP/\fIlength\fP[/\fIunits\fP] Draws a simple vertical scale centered on \fIlon0/lat0\fP\&. Use \fB\-Sx\fP to specify cartesian coordinates instead. \fIlength\fP is in z units, append unit name for labeling. \fBFONT_ANNOT_PRIMARY\fP is used as font. .UNINDENT .INDENT 0.0 .TP \fB\-T\fP\fIpen\fP Draw track [Default is no track]. Append pen attributes to use [Defaults: width = 0.25p, color = black, style = solid]. .UNINDENT .INDENT 0.0 .TP \fB\-U\fP[\fIjust\fP/\fIdx\fP/\fIdy\fP/][\fBc\fP|\fIlabel\fP] (more ...) Draw GMT time stamp logo on plot. .UNINDENT .INDENT 0.0 .TP \fB\-V\fP[\fIlevel\fP] (more ...) Select verbosity level [c]. .UNINDENT .sp \fB\-W\fP\fIpen\fP .INDENT 0.0 .TP \fB\-bi\fP[\fIncols\fP][\fBt\fP] (more ...) Select native binary input. [Default is 3 input columns]. .UNINDENT .INDENT 0.0 .TP \fB\-c\fP\fIcopies\fP (more ...) Specify number of plot copies [Default is 1]. .UNINDENT .INDENT 0.0 .TP \fB\-di\fP\fInodata\fP (more ...) Replace input columns that equal \fInodata\fP with NaN. .UNINDENT .INDENT 0.0 .TP \fB\-f\fP[\fBi\fP|\fBo\fP]\fIcolinfo\fP (more ...) Specify data types of input and/or output columns. .UNINDENT .INDENT 0.0 .TP \fB\-g\fP[\fBa\fP]\fBx\fP|\fBy\fP|\fBd\fP|\fBX\fP|\fBY\fP|\fBD\fP|[\fIcol\fP]\fBz\fP[+|\-]\fIgap\fP[\fBu\fP] (more ...) Determine data gaps and line breaks. .UNINDENT .INDENT 0.0 .TP \fB\-h\fP[\fBi\fP|\fBo\fP][\fIn\fP][\fB+c\fP][\fB+d\fP][\fB+r\fP\fIremark\fP][\fB+r\fP\fItitle\fP] (more ...) Skip or produce header record(s). .UNINDENT .INDENT 0.0 .TP \fB\-i\fP\fIcols\fP[\fBl\fP][\fBs\fP\fIscale\fP][\fBo\fP\fIoffset\fP][,\fI\&...\fP] (more ...) Select input columns (0 is first column). .UNINDENT .INDENT 0.0 .TP \fB\-p\fP[\fBx\fP|\fBy\fP|\fBz\fP]\fIazim\fP/\fIelev\fP[/\fIzlevel\fP][\fB+w\fP\fIlon0\fP/\fIlat0\fP[/\fIz0\fP]][\fB+v\fP\fIx0\fP/\fIy0\fP] (more ...) Select perspective view. .UNINDENT .INDENT 0.0 .TP \fB\-t\fP[\fItransp\fP] (more ...) Set PDF transparency level in percent. .UNINDENT .INDENT 0.0 .TP \fB\-:\fP[\fBi\fP|\fBo\fP] (more ...) Swap 1st and 2nd column on input and/or output. .UNINDENT .INDENT 0.0 .TP \fB\-^\fP or just \fB\-\fP Print a short message about the syntax of the command, then exits (NOTE: on Windows use just \fB\-\fP). .TP \fB\-+\fP or just \fB+\fP Print an extensive usage (help) message, including the explanation of any module\-specific option (but not the GMT common options), then exits. .TP \fB\-?\fP or no arguments Print a complete usage (help) message, including the explanation of options, then exits. .UNINDENT .SH EXAMPLES .sp To plot the magnetic anomaly stored in the file track.xym along track @ 1000 nTesla/cm (after removing a mean value of 32000 nTesla), using a 15\-cm\-wide Polar Stereographic map ticked every 5 degrees in Portrait mode, with positive anomalies in red on a blue track of width 0.25 points, use .INDENT 0.0 .INDENT 3.5 .INDENT 0.0 .INDENT 3.5 .sp .nf .ft C gmt pswiggle track.xym \-R\-20/10/\-80/\-60 \-JS0/90/15c \-Z1000 \-B5 \e \-C32000 \-P \-Gred \-T0.25p,blue \-S1000 \-V > track_xym.ps .ft P .fi .UNINDENT .UNINDENT .UNINDENT .UNINDENT .SH BUGS .sp Sometimes the (x,y) coordinates are not printed with enough significant digits, so the local perpendicular to the track swings around a lot. To see if this is the problem, you should do this: .INDENT 0.0 .INDENT 3.5 .INDENT 0.0 .INDENT 3.5 .sp .nf .ft C awk \(aq{ if (NR > 1) print atan2(y\-$1, x\-$2); y=$1; x=$2; }\(aq yourdata.xyz | more .ft P .fi .UNINDENT .UNINDENT .UNINDENT .UNINDENT .sp (note that output is in radians; on some machines you need "nawk" to do this). Then if these numbers jump around a lot, you may do this: .INDENT 0.0 .INDENT 3.5 .INDENT 0.0 .INDENT 3.5 .sp .nf .ft C awk \(aq{ print NR, $0 }\(aq yourdata.xyz | filter1d \-Fb5 \-N4/0 \e \-\-FORMAT_FLOAT_OUT=%.12g > smoothed.xyz .ft P .fi .UNINDENT .UNINDENT .UNINDENT .UNINDENT .sp and plot this data set instead. .SH SEE ALSO .sp gmt, gmtcolors, filter1d, psbasemap, splitxyz .SH COPYRIGHT 2016, P. Wessel, W. H. F. Smith, R. Scharroo, J. Luis, and F. Wobbe .\" Generated by docutils manpage writer. .