mgd77listgmt - Online na nuvem

PROGRAMA:

NOME

mgd77list - Extrair dados de arquivos MGD77

SINOPSE

lista mgd77 IDs NGDC colunas[,lógica] [:testes de bit] [[+]c|d|f|m|tcódigo ] [ f|g|e ] [
A|adata de início ] [ B|bdata de parada ] [] [ astartrec ] [ bStoprec ] [ ignorar ] [[acondicionável]]
[ d|sunidade ] [ a|c|vminutos/max ] [ região ] [ aponto de partida[unidade] ] [ bponto de parada[unidade] ] [ [m|e]]
[[nível]] [ peso ] [ +|- ] [ -bo] [ -h]

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

DESCRIÇÃO

lista mgd77 lê . [mgd77 | nc] e produz uma tabela ASCII [ou binária]. o
Os arquivos. [mgd77 | nc] contêm informações da faixa, como leg-id, hora e posição,
observáveis ​​geofísicos, como gravidade, magnetismo e batimetria, e códigos de controle e
correções como Eotvos e correções diurnas. Os arquivos netCDF estendidos MGD77 + podem
também contêm colunas de usuário adicionais (para uma lista de colunas disponíveis, use mgd77info
-C, e para aprender como adicionar suas próprias colunas personalizadas, consulte mgd77gerenciar) O usuário pode
extrair qualquer combinação desses parâmetros, qualquer uma das 8 grandezas calculadas (distância,
direção, mudança de curso, velocidade, correção de Carter, correção de Eotvos e gravidade e
campos magnéticos de referência globais), subunidades de calendário de tempo (ano, mês, dia, hora,
min, sec), a id do NGDC e, finalmente, um peso predefinido (consulte -W) Uma subseção pode ser
especificado passando intervalos de tempo ou distância ao longo da trilha ou selecionando um
região geográfica. Finalmente, cada registro de saída pode ser necessário para passar qualquer número de
testes lógicos envolvendo valores de dados ou sinalizadores de bits. Se vários cruzeiros forem solicitados,
eles são separados por cabeçalhos de segmento.

É REQUERIDO ARGUMENTOS

IDs NGDC
Pode ser um ou mais dos cinco tipos de especificadores:

1. IDs NGDC de 8 caracteres, por exemplo, 01010083, JA010010etc., Etc.

2. 2 caracteres códigos que retornarão todos os cruzeiros de cada agência.

3. 4 caracteres códigos, que retornarão todos os cruzeiros daqueles
embarcações.

4. = , Onde é uma tabela com IDs NGDC, um por linha.

5. Se nada for especificado, retornamos todos os cruzeiros na base de dados.

(Ver mgd77info -L para códigos de agência e navio). As extensões ".mgd77" ou ".nc"
será automaticamente anexado, se necessário (use -I para ignorar certos tipos de arquivo).
Os arquivos do cruzeiro serão procurados primeiro no diretório atual e em segundo lugar em todos
diretórios listados em $ MGD77_HOME/mgd77_paths.txt [Se $ MGD77_HOME não está definido
será padrão para $ GMT_SHAREDIR/ mgd77].

-Fcolunas[,lógica] [:testes de bit]
As informações necessárias colunas string deve ser uma lista de parâmetros separados por vírgulas
abreviações fornecidas na ordem de saída desejada. Quaisquer parâmetros fornecidos em maiúsculas
não deve ser NaN em um registro para que ocorra a saída. A menos que especificado separadamente, o
formato de saída (se ASCII) é controlado pelo parâmetro GMT FORMAT_FLOAT_OUT. o
abreviações de coluna disponíveis para informações armazenadas nos arquivos (algumas colunas
pode ser NaN) são:

recno O contador do número do registro.

drt O tipo de registro digital, geralmente 3 ou 5 (para cruzeiros compatíveis com o Y2K).

id A string do ID da pesquisa (nome da perna).

ngdcid A string de ID de cruzeiro NGDC de 8 caracteres (geralmente o prefixo do arquivo).

tempo Escolha entre o tempo absoluto do calendário (um tempo, o padrão) no formato
ditado pelos parâmetros GMT FORMAT_DATE_OUT e FORMAT_CLOCK_OUT,
Tempo relativo (hora) no formato ditado pelos parâmetros GMT
FORMAT_FLOAT_OUT e TIME_SYSTEM (ou TIME_EPOCH e TIME_UNIT)), ou
Ano fracionário (vez) no formato ditado por FORMAT_FLOAT_OUT.

lon Longitude no formato ditado pelo parâmetro GMT FORMAT_GEO_OUT.

lat Longitude no formato ditado pelo parâmetro GMT FORMAT_GEO_OUT.

twt Tempo de viagem bidirecional (em s).

profundidade Batimetria corrigida (em m, positivo abaixo do nível do mar).

mtf1 Intensidade do campo magnético total do sensor 1 (em nTesla).

mtf2 Intensidade do campo magnético total do sensor 2 (em nTesla).

mag Anomalia magnética residual (em nTesla).

gobs Gravidade observada (em mGal).

faa Anomalia de gravidade ao ar livre (em mGal).

ptc Código do tipo de posição (1 = fixo, 3 = interpolado, 9 = não especificado).

bcc Código de correção batimétrica, indicando o procedimento usado para converter a viagem
tempo para profundidade. (01-55 = zona de Matthews usada para corrigir a profundidade, 59 =
As correções de Matthews são usadas, mas as zonas não são especificadas no registro de dados,
60 = Fórmula S. Kuwahara para TS, 61 = Fórmula de Wilson para TS, 62 = Del
Fórmula de Grosso para TS, 63 = tabelas de Carter, 88 = Outros, descrito em
seções de cabeçalho, 99 = não especificado).

btc Código de tipo batimétrico, indicando como o valor da batimetria foi obtido (1 =
observado, 3 = interpolado, 9 = não especificado).

msens Sensor magnético usado para avaliar o campo residual (1 = 1º ou principal
sensor, 2 = 2º ou sensor posterior, 9 = não especificado).

MSD Profundidade (ou altitude) do sensor magnético (em m, positivo abaixo do nível do mar).

diur Correção magnética diurna (em nTesla).

et Correção de Eotvos (em mGal).

sln String do número da linha sísmica.

sspn Sequência do número do ponto de tiro sísmico.

nqc Código de qualidade de navegação (5 = suspeito, por instituição de origem, 6 =
suspeita, pelo NGDC, 9 = nenhum problema identificado).

Além disso, as seguintes quantidades de navegação derivadas podem ser solicitadas:

ano O ano de cada registro.

mês O mês de cada registro.

dia O dia do mês de cada registro.

hora A hora de cada registro.

minutos As atas de cada registro.

seca Os segundos decimais de cada registro.

dist Distância ao longo da pista desde o início da perna. Para método de cálculo, veja -C
[distâncias do grande círculo esférico], e para unidades de distância, consulte -N [km].

az Azimute do navio (rumo) medido no sentido horário do norte (em graus).

cc Mudança de curso do navio (mudança na direção) medida no sentido horário do norte (em
graus).

vel Velocidade do navio; Vejo -N para unidades [m / s].

Finalmente, as seguintes quantidades calculadas podem ser solicitadas:

peso Peso atribuído a este conjunto de dados (ver -W).

carreteiro Correção de profundidade de Carter, se twt está presente no arquivo (em m). Sinal: Correção
deve ser adicionado às profundidades não corrigidas para produzir uma profundidade corrigida.

IGRF Campo de referência geomagnético internacional (campo total) (em nTesla).

ngrav Campo de referência de gravidade internacional ("gravidade normal") (em mGal). Campo
é selecionado com base no parâmetro Código de Fórmula Teórica da Gravidade no
cabeçalho MGD77 do cruzeiro. Se isso não for definido ou for inválido, o padrão é o
IGF 1980. Alternativamente, especifique o campo diretamente usando -Af (veja isso
opção para mais detalhes).

ceot Correção de Eotvos calculada a partir da navegação, usando E = 7.5038 * V *
cos (lat) * sin (az) + 0.004154 * V ^ 2 (em mGal). Sinal: Correção deve ser
adicionado ao faa não corrigido para produzir um faa corrigido.

As seguintes bandeiras abreviadas também são reconhecidas:

mgd77 Isso resulta em todos os 27 campos MGD77 sendo escritos no MGD77 oficial
ordem.

mgd77t O conjunto completo de todas as 26 colunas na especificação MGD77T.

geo Isso limita a saída a 10 campos (tempo, lon, lat mais os sete
observações geofísicas twt, profundidade, mtf1, mtf2, mag, gobs e faa) Por
acrescentando + a qualquer um desses conjuntos também anexaremos dist, Azim, cc, vel,
e peso conforme listado acima.

todos os Isso retorna todas as colunas de dados no arquivo.

tudo T Como mgd77t, mas com itens de tempo escritos como uma string de data-hora.

Como uma opção, testes lógicos podem ser adicionados para qualquer uma das observações anexando
,lógica, que é composto por uma ou mais instruções separadas por vírgulas do
formulário pelaOPvalor, Onde pela é um dos parâmetros listados acima, OP é um lógico
operador (<, <=, =,! =,> =,>, |) e valor é uma constante usada na comparação.
Os parâmetros de ponto flutuante são comparados numericamente; parâmetros de personagem são
comparados lexicamente (após a remoção dos espaços em branco à esquerda e à direita). A parte
comparação (|) significa que pelo menos um dos bits em valor deve ser ligado em
pela. Pelo menos um dos testes deve ser verdadeiro para que o registro seja gerado, exceto para
testa usando parâmetros de maiúsculas e minúsculas, todos devem ser verdadeiros para que a saída ocorra. Observação
que especificar um teste não implica que a coluna correspondente será
incluído no fluxo de saída; deve estar presente em colunas para que isso ocorra.
Observação: alguns dos operadores são caracteres especiais do UNIX e você é aconselhado a
coloque aspas em todo o argumento para -F.

Finalmente, para arquivos MGD77 +, você pode opcionalmente anexar:testes de bit que é: (dois pontos)
seguido por um ou mais separados por vírgula + -col termos. Isso se compara a
bitflags apenas para cada coluna listada. Aqui, + significa que o bit escolhido deve ser 1 (ON)
enquanto - significa que deve ser 0 (OFF). Todos os testes de bits fornecidos devem ser aprovados. Por padrão,
Arquivos MGD77 + que têm o especial MGD77_flags coluna presente usará esses sinalizadores,
e as observações associadas aos bits ON (o que significa que são marcados como ruins) serão
definido como NaN; append: sem informações finais para desligar esse comportamento (ou seja,
nenhum sinalizador de bit será consultado).

OPCIONAL ARGUMENTOS

-A [+] c | d | f | m | tcódigo
Por padrão, profundidade corrigida (profundidade), anomalia residual magnética (mag), ao ar livre
anomalia de gravidade (faa), e a quantidade derivada de correção de profundidade de Carter (carreteiro)
são todas as saídas como estão (se selecionado em -F); esta opção ajusta esse comportamento. Para
cada uma dessas colunas existem 2-4 maneiras de ajustar os dados. Acrescentar c(artéria),
d(epta), f(aa), ou m(ag) e selecione o código para o procedimento que você deseja aplicar.
Você pode selecionar mais de um procedimento para uma coluna de dados somando seus
códigos (1, 2, 4 e 8). Por exemplo, -Ar condicionado3 tentará primeiro o método -Ar condicionado1 para estimar um Carter
correção, mas se profundidade é NaN, vamos tentar na próxima -Ar condicionado2 que usa apenas twt. Em tudo
casos, se algum dos valores exigidos por um procedimento de ajuste for NaN, então o
o resultado será NaN. Isso também é verdadeiro se a anomalia original for NaN. Especificamos -A +
para recalcular anomalias, mesmo se a anomalia no arquivo for NaN. Além disso, você
pode usar - Ato criar horários falsos para cruzeiros que não têm horário; estes são baseados em
distâncias e duração do cruzeiro.

-Ar condicionado Determina como o carreteiro o termo de correção é calculado. Abaixo, C (twt) apoia
a profundidade corrigida por Carter (também depende de lon, lat), VOCÊ(twt, v) é o
profundidade não corrigida (= twt * v / 2) usando como v o parâmetro "Velocidade do som presumida"
no cabeçalho MGD77 (se for uma velocidade válida, caso contrário, o padrão é 1500 m / s);
alternativamente, acrescente sua velocidade preferida v em m / s, TU (profundidade, v) é a via dupla
tempo de viagem estimado a partir do (presumivelmente) não corrigido profundidade, e TC (profundidade) é o
Tempo de viagem de ida e volta obtido invertendo o (presumivelmente) corrigido profundidade usando o
Fórmula de correção de Carter. Selecione de

-Ac1[,v] retorna a diferença entre U (twt, v) e profundidade [Predefinição].

-Ac2[,v] retorna a diferença entre U (twt, v) e Carter (twt).

-Ac4[,v] retorna a diferença entre (assumido não corrigido) profundidade e Carter
(TU (profundidade)).

-Ac8[,v] retorna a diferença entre U (TC (profundidade), v) e profundidade.

-De Anúncios Determina como o profundidade a saída da coluna é obtida:

-Anúncio1 Retorna profundidade conforme armazenado no conjunto de dados [Padrão].

-Anúncio2[,v] retorna a profundidade calculada não corrigida U (twt, v).

-Anúncio4 retorna a profundidade corrigida calculada C (twt).

-Af Determina como o faa a saída da coluna é obtida. Se ngrav (ou seja, o Internacional
Campo de referência de gravidade (IGF), ou "gravidade normal") é necessário que seja selecionado
com base no parâmetro de cabeçalho MGD77 "Código da fórmula de gravidade teórica"; se este
código não está presente ou é inválido, o padrão é 4. Como alternativa, anexe o
preferido campo (1-4) para selecionar 1 (Heiskanen 1924), 2 (IGF 1930), 3 (IGF 1967) ou 4
(IGF 1980). Selecione de

-Af1[,campo] retorna faa conforme armazenado no conjunto de dados [Padrão]. Opcionalmente, define o
IGF campo para usar se você também solicitou ngrav como uma coluna de saída em -F.

-Af2[,campo] retorna a diferença entre gobs e ngrav (com opcional campo
diretiva).

-Af4[,campo] retorna a combinação de gobs + et - ngrav (com opcional campo
diretiva).

-Af8[,campo] retorna a combinação de gobs + pred_eot - ngrav (com opcional
campo diretiva).

-Sou Determina como o mag a saída da coluna é obtida. Pode haver um ou dois no total
medições de campo no arquivo (mtf1 e mtf2), e a coluna msens pode declarar
qual é o sensor principal (1 ou 2; também pode ser indefinido). Selecione de

-Am1 Retorna mag conforme armazenado no conjunto de dados [Padrão].

-Am2 retorna a diferença entre mgfx e IGRF, Onde x é o sensor principal (1
or 2) indicado pelo msens campo de dados (o padrão é 1 se não for especificado).

-Am4 retorna a diferença entre mgfx e IGRF, Onde x é o sensor (2 or 1)
não indicado pelo msens campo de dados (o padrão é 2 se não for especificado).

-Amc [unidade] Aplique uma correção que tenta compensar o fato de que o
campo magnético não foi adquirido na mesma posição que a referência de navegação
(posição do navio). Isso é feito reinterpolando o valor magnético total
campo para o que teria se fosse medido na posição do navio
(lembre-se, provavelmente foi medido metros atrás). Devido a isso
etapa de interpolação, navegação ruim, ou seja, muitos pontos repetidos, podem causar
problemas. Medidas são tomadas para minimizar este efeito, mas elas não estão 100% completas
prof. O método de interpolação é controlado pelo GMT GMT_INTERPOLANT padrão.
Acrescentar e para metro, f para os pés, k por km, m por milhas, n para milhas náuticas, ou u
para pés de levantamento [o padrão é e (metros)].

-Cf | g | e
Acrescente um código de uma letra para selecionar o procedimento para a distância ao longo da trilha
cálculo (ver -N para selecionar unidades):

f Distâncias planas da Terra.

g Grandes distâncias do círculo [padrão].

e Distâncias geodésicas no elipsóide GMT atual.

-Simdata de início
Não liste os dados coletados antes data de início (aaaa-mm-ddBD (T) [hh: mm: ss]) [O padrão é
início do cruzeiro]. Usar -DA para excluir registros cujo tempo é indefinido (ou seja, NaN).
[O padrão relata esses registros].

-Dbdata de parada
Não liste os dados coletados em ou depois data de parada (aaaa-mm-ddBD (T) [hh: mm: ss]).
[O padrão é fim do cruzeiro]. Usar -DB para excluir registros cujo tempo é indefinido
(ou seja, NaN). [O padrão relata esses registros].

-E Correspondência exata: apenas registros de saída que correspondam a todas as colunas geofísicas solicitadas
[Padrão gera registros que correspondem a pelo menos uma das colunas observadas].

-Gástartrec
Não liste registros antes startrec [O padrão é 0, o primeiro registro].

-GBStoprec
Não liste dados após Stoprec. [O padrão é o último registro].

-Iignorar
Ignore certos formatos de arquivo de dados. Acrescentar a | c | t ignorar MGD77
ASCII, MGD77 + netCDF ou arquivos de tabela ASCII separados por tabulação, respectivamente. o
a opção pode ser repetida para ignorar mais de um formato. [O padrão não ignora nenhum].

-EU[acondicionável]
Aplique as correções ideais às colunas onde tais correções estão disponíveis. Acrescentar
a tabela de correção a ser usada [Padrão usa a tabela de correção
mgd77_corrections.txt no $ MGD77_HOME diretório]. Para o formato deste arquivo,
veja CORREÇÕES abaixo.

-n Emita um registro de cabeçalho de segmento com ID de cruzeiro para cada cruzeiro.

-Nd | sunidade
Acrescentar d para distância ou s para velocidade, então dê o desejado unidade as e (medidor ou
em), f (pés ou pés / s), k (km ou km / h), m (milhas ou milhas / hora), n (milhas náuticas
ou nós), ou u (levantamento de pés ou pés / pés). [O padrão é -Ndk -Nse (km e m / s)].

-Qaminutos/max
Especifique um intervalo aceito (minutos/max) de azimutes. Registros cuja trilha azimute cai
fora deste intervalo são ignorados [0/360].

-Qcminutos/max
Especifique um intervalo aceito (minutos/max) é claro que muda. Registra cujo percurso
as alterações fora deste intervalo são ignoradas [-360 / + 360]. Usar -CQ tomar o
valor absoluto da mudança de curso antes do teste [Padrão usa curso assinado
alterar].

-Qvminutos/max
Especifique um intervalo aceito (minutos/max; ou apenas minutos se não houver limite superior) de
velocidades. Registros cuja velocidade da trilha esteja fora desta faixa são ignorados
[0 / infinito].

-R [unidade]oeste/leste/sul/norte[/zmin/zmax] [r]
oeste, leste, sul e norte especifique a região de interesse, e você pode especificar
em graus decimais ou no formato [+ -] dd: mm [: ss.xxx] [W | E | S | N]. Acrescentar r se mais baixo
as coordenadas do mapa à esquerda e à direita superior são fornecidas em vez de w / e / s / n. Os dois
taquigrafias -Rg e -Rd significa domínio global (0/360 e -180 / + 180 de longitude
respectivamente, com -90 / + 90 na latitude). Alternativamente para a criação da grade, dê
Rcódigo/lat/nx/ny, Onde código é uma combinação de 2 caracteres de L, C, R (para a esquerda,
centro ou direita) e T, M, B para superior, meio ou inferior. por exemplo, BL para canto esquerdo inferior.
Isso indica qual ponto em uma região retangular o lon/lat coordenar refere-se
para, e as dimensões da grade nx e ny com espaçamentos de grade via -I é usado para criar
a região correspondente. Alternativamente, especifique o nome de um arquivo de grade existente
e os votos de -R as configurações (e espaçamento da grade, se aplicável) são copiadas da grade.
utilização -Runidade espera coordenadas projetadas (cartesianas) compatíveis com as escolhidas -J
e nós projetamos inversamente para determinar a região geográfica retangular real. Para
Visão de perspectiva (-p), opcionalmente anexar /zmin/zmax. Em caso de visão em perspectiva
(-p), um intervalo z (zmin, zmax) pode ser anexado para indicar a terceira dimensão. Esse
precisa ser feito apenas ao usar o -Jz opção, não ao usar apenas a -p opção.
Neste último caso, uma vista em perspectiva do plano é traçada, sem terceiros
dimensão.

-Saponto de partida[unidade]
Não liste dados inferiores a ponto de partida metro ao longo da trilha do porto de
partida. Acrescentar e para metro, f para os pés, k por km, m por milhas, n para náutico
milhas, ou u para pés de levantamento [o padrão é 0e (metros)].

-Sbponto de parada[unidade]
Não liste dados que são ponto de parada ou mais metros ao longo da trilha do porto de
partida. Acrescentar e para metro, f para os pés, k por km, m por milhas, n para náutico
milhas, ou u para pés de levantamento [o padrão é o fim da trilha].

-T [m | e]
Desativa o ajuste de outra forma automático de valores com base em termos de correção
que são armazenados no arquivo MGD77 + e usados ​​para neutralizar coisas como errado
unidades usadas pela instituição de origem ao criar o arquivo MGD77 original a partir de
que o arquivo MGD77 + deriva (a opção não tem efeito no MGD77 ASCII simples
arquivos). Acrescentar m or e para limitar a opção ao MGD77 ou conjunto de colunas estendidas
apenas [o padrão se aplica a ambos].

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

-Wpeso
Defina o peso para esses dados. A opção de saída de peso deve ser definida em -F. Isto é
útil se os dados forem processados ​​com as técnicas de média ponderada
oferecido por meio de bloco, bloqueado e modo de bloqueio [1].

-Z+|- Anexe o sinal que você deseja profundidade, carreteiro e MSD valores abaixo do nível do mar (-Z-
dá batimetria negativa) [O padrão é positivo para baixo].

-bo [notas][tipo] (mais ...)
Selecione a saída binária nativa. ignorado se -bo é selecionado. Da mesma forma, campos de string
não pode ser selecionado. Observe que se o tempo for uma das colunas de saída binária, ele
ser armazenado como tempo Unix (segundos desde 1970). Para ler esta informação no GMT para
obter o tempo absoluto do calendário exigirá que você use --TIME_SYSTEM = 1.

-h [i | o] [n] [+ c] [+ d] [+ rcomentário] [+ rtítulo] (mais ...)
Pule ou produza registro (s) de cabeçalho.

-^ or apenas por -
Imprime uma mensagem curta sobre a sintaxe do comando e sai (NOTA: no Windows
use apenas -).

-+ or apenas por +
Imprima uma mensagem de uso extensivo (ajuda), incluindo a explicação de qualquer
opção específica do módulo (mas não as opções comuns GMT) e, em seguida, sai.

-? or não argumentos
Imprima uma mensagem completa de uso (ajuda), incluindo a explicação das opções e, em seguida,
saídas.

--versão
Imprima a versão GMT e saia.

--show-datadir
Imprima o caminho completo para o diretório de compartilhamento GMT e saia.

EXEMPLOS

Para obter uma listagem de (distância, direção, gravidade, batimetria) de 01010047.mgd77, começando em
3 de junho de 1971 20:45 e terminando na distância = 5000 km, use o seguinte comando:

gmt mgd77list 01010047 -Da1971-06-03T20: 45 -Sb5000 -Fdist, azim, faa, depth> myfile.d

Para fazer contribuições para meio de bloco e à superfície dos talhos, usando anomalias de ar livre de todos os cruzeiros
listados no arquivo cruises.lis, mas apenas os dados que estão dentro da área especificada, e
faça a saída binária:

gmt mgd77list `cat cruises.lis` -Flon, lat, faa -R-40 / -30 / 25/35 -bo> allgrav.b

Para extrair os locais de profundidades superiores a 9000 metros que não foram interpolados (btc !=
1) de todos os cruzeiros listados no arquivo cruises.lis:

gmt mgd77list `cat cruises.lis` -F" depth, DEPTH> 9000, BTC! = 1 "> really_deep.d

Extrair dist, faa e grav12_2 de registros cujas profundidades sejam inferiores a 3 km e
onde nenhum dos campos solicitados é NaN, de todos os arquivos MGD77 + netCDF cujo cruzeiro
ids estão listados no arquivo cruises.lis, tentamos

gmt mgd77list `cat cruises.lis` -E -Ia -F" dist, faa, grav12_2, profundidade <3000 "> \
raso_grav.d

Para extrair dist, faa e grav12_2 de todos os arquivos MGD77 + netCDF cujos IDs de cruzeiro são
listado no arquivo cruises.lis, mas recupera apenas os registros cujo bitflag para faa indica
Valores RUINS, nós tentamos

gmt mgd77list `cat cruises.lis` -E -Ia -F" dist, faa, grav12_2: + faa "> bad_grav.d

Para produzir lon, lat, mag e faa de todos os cruzeiros listados no arquivo cruises.lis, mas
recalcule os dois resíduos com base nos campos de referência mais recentes, tente:

gmt mgd77list `cat cruises.lis` -Flon, lat, mag, faa -Af2,4 -Am2> data.d

RECALCULADO ANOMALIAS

Quando anomalias recalculadas são solicitadas (explicitamente por meio do -A opção ou
implicitamente por meio de metadados E77 no arquivo MGD77 +), fazemos isso apenas para os registros cujos
anomalia original não era um NaN. Esta restrição é implementada uma vez que muitas colunas de anomalias
contém correções, geralmente na forma de alterações editadas manualmente, que não podem ser
duplicado da observação correspondente.

IGRF

Os cálculos do IGRF são baseados em um programa Fortran escrito por Susan Macmillan, britânica
Geological Survey, traduzido para C via f2c por Joaquim Luis, U Algarve, e adaptado para
Estilo GMT por Paul Wessel.

IGF

As equações usadas são reproduzidas aqui usando coeficientes extraídos diretamente do
código-fonte (informe-nos se encontrar erros):

(1) g = 978052.0 * [1 + 0.005285 * sin ^ 2 (lat) - 7e-6 * sin ^ 2 (2 * lat) + 27e-6 * cos ^ 2 (lat) *
cos ^ 2 (lon-18)]

(2) g = 978049.0 * [1 + 0.0052884 * sin ^ 2 (lat) - 0.0000059 * sin ^ 2 (2 * lat)]

(3) g = 978031.846 * [1 + 0.0053024 * sin ^ 2 (lat) - 0.0000058 * sin ^ 2 (2 * lat)]

(4) g = 978032.67714 * [(1 + 0.00193185138639 * sin ^ 2 (lat)) / sqrt (1 - 0.00669437999013 *
sin ^ 2 (lat))]

CORRECÇÕES

A tabela de correção é um arquivo ASCII com coeficientes e parâmetros necessários para realizar
correções. Registros de comentários começando com # são permitidos. Todos os registros de correção são de
a forma

CruiseID observação correção

onde CruiseID é um prefixo NGDC, observação é uma das abreviações de geofísico
observações listadas em -F acima, e correção consiste em um ou mais prazoé que vai
ser resumido e então subtraído da observação antes da saída. Cada prazo deve ter
esta sintaxe exata:

fator[* [função] ([escada](abrev[-origem])) [^poder]]

onde os termos entre colchetes são opcionais (os próprios colchetes não são usados, mas são regulares
os parênteses devem ser usados ​​conforme indicado). Não são permitidos espaços, exceto entre prazos. O
fator é a amplitude da função de base, enquanto o opcional função pode ser um de
sin, cos ou exp. O opcional escada e origem pode ser usado para traduzir o argumento
(antes de atribuí-lo à função opcional). O argumento abrev é um dos
abreviações para as observações listadas acima. Se origem é dado como T isso significa que nós
deve substituí-lo pelo valor de abrev para o primeiro registro no arquivo (este é
normalmente só é feito para tempo) Se a primeira entrada de registro for NaN, nós revertemos origem para zero.
Opcionalmente, aumente a expressão inteira para o dado poder, antes de multiplicar pelo
amplitude. A seguir está um exemplo de correções fictícias para o cruzeiro 99999999,
implicando o profundidade deve ter a correção Carter removida, faa deve ter um linear
tendência removida, a anomalia magnética (mag) deve ser corrigido por uma estranha dependência de
direção e latitude do navio, e gobs precisa ter 10 mGal adicionado (portanto, dado como -10):

99999999 profundidade 1.0 * ((carter))

99999999 faa 14.1 1e-5 * ((tempo-T))

99999999 mag 0.5*cos(0.5*(azim-19))^2 1.0*exp(-1e-3(lat))^1.5

99999999 toneladas -10

Use mgd77listgmt online usando serviços onworks.net