grdredpolgmt - Online na nuvem

PROGRAMA:

NOME

grdredpol - Calcula a redução contínua para o pólo, também conhecido como RTP diferencial.

SINOPSE

Grdredpol anom_grd rtp_grd [dec / mergulho] [ dec_grd / dipgrd] [ ] [ m | r] [] [ largura_win] [
[nível]] [ ano ] [ filtro ] [[nível]] [ -n]

Observação: Nenhum espaço é permitido entre o sinalizador de opção e os argumentos associados.

DESCRIÇÃO

Grdredpol vai levar um .nc arquivo com uma anomalia magnética e calcular a redução para o
anomalia do pólo (RTP). Esta anomalia é a que teria sido produzida se os corpos
foram magnetizados verticalmente e as anomalias foram observadas no pólo geomagnético.
O procedimento RTP padrão assume que a direção da magnetização é uniforme em todo o
corpo causador, e o campo geomagnético deve ser uniforme na direção ao longo do estudo
região. Embora essas suposições sejam razoáveis ​​para áreas pequenas, elas não se aplicam para
Grandes áreas.

No método usado aqui, os cálculos são realizados tanto na frequência quanto no espaço
domínios. A ideia é que uma grande área possa ser decomposta em janelas de pequeno porte onde ambos
o campo ambiente e o vetor de magnetização mudam muito. Dentro de cada
dessas janelas, ou bins, um conjunto de coeficientes de filtro são calculados e reconstruídos para
cada indivíduo aponta o filtro de componente usando uma expansão de série de Taylor de primeira ordem.

É REQUERIDO ARGUMENTOS

anom_grd
A grade de anomalias a ser convertida.

-Grtp_grd
é o nome do arquivo para arquivo grdfile de saída com a solução RTP

OPCIONAL ARGUMENTOS

-Cdec / mergulho
Use esta declinação (constante) e ângulos de inclinação para o campo e
magnetização. Esta opção consiste no procedimento RTP clássico.

-Edip_grd / dec_grd
Obtenha DIP e DEC de magnetização dessas grades [padrão: use IGRF]. Observe que estes
duas grades não precisam ter a mesma resolução que a grade de anomalias. Eles podem ser
mais grosseiro.

-Fm / n O tamanho da janela do filtro em termos de linhas / colunas. O valor padrão é 25x25.

-Mm | r Defina as condições de limite. m | r significa espelhar ou replicar bordas (o padrão é zero
preenchimento).

-N NÃO use a expansão Taylor.

-Roeste/leste/sul/norte
define a região dos pontos de saída. [Padrão: igual à entrada.]

-Tano Ano decimal usado pela rotina IGRF para calcular a declinação e inclinação em
cada ponto [padrão: 2000]

-Wlargura
O tamanho da janela móvel em graus [5].

-Zfiltro_grd
Grave o arquivo de filtro no disco.

-V [nível] (mais ...)
Selecione o nível de verbosidade [c].

-n [b | c | l | n] [+ a] [+ bBC] [+ c] [+ tlimiar] (mais ...)
Selecione o modo de interpolação para grades.

CONSEQUÊNCIAS OF GRID REAMOSTRAGEM

A nova amostra ou amostragem de grades usará vários algoritmos (ver -n) que pode levar a
possíveis distorções ou resultados inesperados nos valores reamostrados. Um efeito esperado
de reamostragem com splines é a tendência para os novos valores reamostrados excederem ligeiramente
os limites mín. / máx. globais da grade original. Se isso for inaceitável, você pode impor
recorte dos valores dos valores reamostrados para que não excedam os valores mín / máx de entrada por
acrescentando +c para o seu -n opção.

EXEMPLOS

Suponha que anom.grd é um arquivo com a anomalia magnética reduzida à época de 2010 e
que o dez.grd e mergulho.grd conter a declinação e inclinação de magnetização
respectivamente para uma área que inclui aquela do anom.grd, calcule o RTP usando caixas de
2 graus e um filtro de 45 coeficientes.

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

Para calcular o mesmo RTP mas agora com os vetores de campo e magnetização colineares e
calculado do IGRF:

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

REFERÊNCIA

Luis, JL e Miranda, JM (2008), Reavaliação de crons magnéticos no Atlântico Norte
entre 35N e 47N: Implicações para a formação da Tríplice Junção dos Açores e
planalto associado. JGR, vol. 113, B10105, doi: 10.1029 / 2007JB005573

DIREITOS AUTORAIS

2015, P. Wessel, WHF Smith, R. Scharroo, J. Luis e F. Wobbe

Use grdredpolgmt online usando serviços onworks.net