OnWorks favicon

grdredpolgmt - Online in the Cloud

Run grdredpolgmt in OnWorks free hosting provider over Ubuntu Online, Fedora Online, Windows online emulator or MAC OS online emulator

This is the command grdredpolgmt that can be run in the OnWorks free hosting provider using one of our multiple free online workstations such as Ubuntu Online, Fedora Online, Windows online emulator or MAC OS online emulator



grdredpol - Compute the Continuous Reduction To the Pole, AKA differential RTP.


grdredpol anom_grd rtp_grd [dec/dip] [ dec_grd/dipgrd] [ <m/n>] [ m|r] [ ] [ win_width] [
[level] ] [ year ] [ filtergrd ] [ [level] ] [ -n<flags> ]

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


grdredpol will take a .nc file with a magnetic anomaly and compute the reduction to the
pole (RTP) anomaly. This anomaly is the one that would have been produce if the bodies
were magnetized vertically and the anomalies were observed at the geomagnetic pole.
Standard RTP procedure assumes the direction of magnetization to be uniform throughout the
causative body, and the geomagnetic field to be uniform in direction throughout the study
region. Although these assumptions are reasonable for small areas, they do not hold for
large areas.

In the method used here computations are carried out in both the frequency and the space
domains. The idea is that a large area may be decomposed in small size windows where both
the ambient field and the magnetization vector change by a very small amount. Inside each
of those windows, or bins, a set of filter coefficients are calculate and reconstruct for
each individual point the component filter using a first order Taylor series expansion.


The anomaly grid to be converted.

is the filename for output grdfile with the RTP solution


Use this (constant) declination and inclination angles for both field and
magnetization. This option consists in the classical RTP procedure.

Get magnetization DIP & DEC from these grids [default: use IGRF]. Note that these
two grids do not need to have the same resolution as the anomaly grid. They can be

-Fm/n The filter window size in terms of row/columns. The default value is 25x25.

-Mm|r Set boundary conditions. m|r stands for mirror or replicate edges (Default is zero

-N Do NOT use Taylor expansion.

defines the Region of the output points. [Default: Same as input.]

-Tyear Decimal year used by the IGRF routine to compute the declination and inclination at
each point [default: 2000]

The size of the moving window in degrees [5].

Write the filter file to disk.

-V[level] (more ...)
Select verbosity level [c].

-n[b|c|l|n][+a][+bBC][+c][+tthreshold] (more ...)
Select interpolation mode for grids.


Resample or sampling of grids will use various algorithms (see -n) that may lead to
possible distortions or unexpected results in the resampled values. One expected effect
of resampling with splines is the tendency for the new resampled values to slightly exceed
the global min/max limits of the original grid. If this is unacceptable, you can impose
clipping of the resampled values values so they do not exceed the input min/max values by
adding +c to your -n option.


Suppose that anom.grd is a file with the magnetic anomaly reduced to the 2010 epoch and
that the dec.grd and dip.grd contain the magnetization declination and inclination
respectively for an area that encloses that of the anom.grd, compute the RTP using bins of
2 degrees and a filter of 45 coefficients.

gmt grdredpol anom.grd -Grtp.grd -W2 -F45/45 -T2010 -Edec.grd/dip.grd -V

To compute the same RTP but now with the field and magnetization vectors collinear and
computed from IGRF :

gmt grdredpol anom.grd -Grtp.grd -W2 -F45/45 -T2010 -V


Luis, J.L. and Miranda, J.M. (2008), Reevaluation of magnetic chrons in the North Atlantic
between 35N and 47N: Implications for the formation of the Azores Triple Junction and
associated plateau. JGR, VOL. 113, B10105, doi:10.1029/2007JB005573


2015, P. Wessel, W. H. F. Smith, R. Scharroo, J. Luis, and F. Wobbe

Use grdredpolgmt online using onworks.net services

Free Servers & Workstations

Download Windows & Linux apps

  • 1
    Strikr Free Software project. Artifacts
    released under a 'intent based'
    dual license: AGPLv3 (community) and
    CC-BY-NC-ND 4.0 international
    Download strikr
  • 3
    giflib is a library for reading and
    writing gif images. It is API and ABI
    compatible with libungif which was in
    wide use while the LZW compression
    algorithm was...
    Download GIFLIB
  • 4
    Alt-F provides a free and open source
    alternative firmware for the DLINK
    DNS-320/320L/321/323/325/327L and
    DNR-322L. Alt-F has Samba and NFS;
    supports ext2/3/4...
    Download Alt-F
  • 5
    Usm is a unified slackware package
    manager that handles automatic
    dependency resolution. It unifies
    various package repositories including
    slackware, slacky, p...
    Download usm
  • 6
    Chart.js is a Javascript library that
    allows designers and developers to draw
    all kinds of charts using the HTML5
    canvas element. Chart js offers a great
    array ...
    Download Chart.js
  • More »

Linux commands