mgd77listgmt: en línea en la nube

PROGRAMA:

NOMBRE

mgd77list - Extrae datos de archivos MGD77

SINOPSIS

Lista mgd77 ID de NGDC columnas[,lógica] [:mordiscos] [[+]c|d|f|m|tcódigo ] [ f|g|e ] [
A|afecha de inicio ] [ B|bfecha de parada ] [] [ astartrec ] [ bStoprec ] [ pasar por alto ] [[cortable]]
[ d|sunidad ] [ a|c|vmin/max ] [ región ] [ adist.inicio[unidad]] [ bparadadistante[unidad]] [[m|e]]
[[nivel]] [ peso ] [ +|- ] [ -bo] [ -h]

Nota: No se permite ningún espacio entre el indicador de opción y los argumentos asociados.

DESCRIPCIÓN

Lista mgd77 lee . [mgd77 | nc] y genera una tabla ASCII [o binaria]. los
Los archivos. [mgd77 | nc] contienen información de la pista, como la identificación del tramo, la hora y la posición,
observables geofísicos como gravedad, magnetismo y batimetría, y códigos de control y
correcciones como Eotvos y correcciones diurnas. Los archivos netCDF extendidos MGD77 + pueden
también contienen columnas de usuario adicionales (para obtener una lista de las columnas disponibles, utilice mgd77info
-Cy para aprender a agregar sus propias columnas personalizadas, consulte mgd77gestionar). El usuario puede
extraer cualquier combinación de estos parámetros, cualquiera de las 8 cantidades calculadas (distancia,
rumbo, cambio de rumbo, velocidad, corrección de Carter, corrección de Eotvos y gravedad y
campos magnéticos de referencia global), subunidades de tiempo del calendario (año, mes, día, hora,
min, seg), la identificación de NGDC y finalmente un peso preestablecido (ver -W). Una subsección puede ser
especificado pasando intervalos de tiempo o distancia a lo largo de la pista o seleccionando un
región geográfica. Finalmente, se puede requerir que cada registro de salida pase cualquier número de
pruebas lógicas que involucran valores de datos o banderas de bits. Si se solicitan varios cruceros,
están separados por encabezados de segmento.

REQUERIDO ARGUMENTOS

ID de NGDC
Puede ser uno o más de cinco tipos de especificadores:

1. ID de NGDC de 8 caracteres, por ejemplo, 01010083, JA010010etc., Etc.

2. 2 caracteres códigos que devolverán todos los cruceros de cada agencia.

3. 4 caracteres códigos, que devolverán todos los cruceros de aquellos
vasos

4. = , dónde es una tabla con ID de NGDC, uno por línea.

5. Si no se especifica nada, devolvemos todos los cruceros en la base de datos.

(Ver mgd77info -L para códigos de agencia y embarcación). Las extensiones ".mgd77" o ".nc"
se agregará automáticamente, si es necesario (use -I para ignorar ciertos tipos de archivos).
Los archivos de crucero se buscarán primero en el directorio actual y en segundo lugar en todos
directorios listados en $ MGD77_HOME/mgd77_paths.txt [Si $ MGD77_HOME no está configurado
por defecto a $ GMT_SHAREDIR/ mgd77].

-Fcolumnas[,lógica] [:mordiscos]
Lo requerido columnas la cadena debe ser una lista de parámetros separados por comas
abreviaturas dadas en el orden de salida deseado. Cualquier parámetro dado en MAYÚSCULAS
no debe ser NaN en un registro para que se produzca la salida. A menos que se especifique por separado, el
El formato de salida (si es ASCII) está controlado por el parámetro GMT FORMAT_FLOAT_OUT. los
abreviaturas de columna disponibles para la información almacenada en los archivos (algunas columnas
puede ser NaN) son:

reconocimiento El contador de números de registro.

DRT El tipo de registro digital, generalmente 3 o 5 (para cruceros que cumplen con el año 2).

id La cadena de identificación de la encuesta (nombre de la rama).

ngdcid La cadena de identificación de crucero NGDC de 8 caracteres (generalmente el prefijo del archivo).

equipo Elija entre el tiempo absoluto del calendario (un momento, el predeterminado) en el formato
dictado por los parámetros GMT FORMAT_DATE_OUT y FORMAT_CLOCK_OUT,
Tiempo relativo (tiempo) en el formato dictado por los parámetros GMT
FORMAT_FLOAT_OUT y TIME_SYSTEM (o TIME_EPOCH y TIME_UNIT)), o
Año fraccionario (ytiempo) en el formato dictado por FORMAT_FLOAT_OUT.

lon Longitud en el formato dictado por el parámetro GMT FORMAT_GEO_OUT.

lat Longitud en el formato dictado por el parámetro GMT FORMAT_GEO_OUT.

twt Tiempo de viaje de ida y vuelta (en s).

profundidad Batimetría corregida (en m, positivo por debajo del nivel del mar).

mtf1 Intensidad del campo magnético total del sensor 1 (en nTesla).

mtf2 Intensidad del campo magnético total del sensor 2 (en nTesla).

mag Anomalía magnética residual (en nTesla).

tragos Gravedad observada (en mGal).

faa Anomalía de la gravedad en aire libre (en mGal).

ptc Código de tipo de posición (1 = fijo, 3 = interpolado, 9 = sin especificar).

bcc Código de corrección batimétrica, que indica el procedimiento utilizado para convertir viajes
tiempo de profundidad. (01-55 = Zona de Matthews utilizada para corregir la profundidad, 59 =
Se utilizaron las correcciones de Matthews, pero las zonas no están especificadas en el registro de datos,
60 = fórmula de S. Kuwahara para TS, 61 = fórmula de Wilson para TS, 62 = Del
Fórmula de Grosso para TS, 63 = tablas de Carter, 88 = Otro, descrito en
secciones de encabezado, 99 = sin especificar).

btc Código de tipo batimétrico, que indica cómo se obtuvo el valor de batimetría (1 =
observado, 3 = interpolado, 9 = no especificado).

msens Sensor magnético que se utiliza para evaluar el campo residual (1 = 1er o líder
sensor, 2 = segundo sensor o sensor posterior, 2 = no especificado).

MSD Profundidad (o altitud) del sensor magnético (en m, positivo por debajo del nivel del mar).

diurno Corrección magnética diurna (en nTesla).

EOT Corrección de Eotvos (en mGal).

sln Cadena de número de línea sísmica.

sspn Cadena de número de punto de disparo sísmico.

nqc Código de calidad de navegación (5 = sospechoso, por institución de origen, 6 =
sospechado, por NGDC, 9 = no se identificaron problemas).

Además, se pueden solicitar las siguientes cantidades de navegación derivadas:

año El año de cada registro.

mes El mes de cada registro.

día El día del mes de cada registro.

horas. La hora de cada registro.

min Los minutos de cada registro.

AMF Los segundos decimales de cada registro.

dist Distancia a lo largo de la pista desde el inicio de la etapa. Para conocer el método de cálculo, consulte -C
[distancias de círculo máximo esférico], y para unidades de distancia, ver -N [kilómetros].

az Azimut (rumbo) del barco medido en el sentido de las agujas del reloj desde el norte (en grados).

cc Cambio de rumbo del barco (cambio de rumbo) medido en el sentido de las agujas del reloj desde el norte (en
grados).

Acelerar Velocidad del barco; ver -N para unidades [m / s].

Finalmente, se pueden solicitar las siguientes cantidades calculadas:

peso Peso asignado a este conjunto de datos (ver -W).

carretero Corrección de profundidad de Carter, si twt está presente en el archivo (en m). Signo: Corrección
se agregará a profundidades no corregidas para producir una profundidad corregida.

igf Campo de referencia geomagnético internacional (campo total) (en nTesla).

ngrav Campo de referencia de gravedad internacional ("gravedad normal") (en mGal). Campo
se selecciona en función del parámetro Código de fórmula teórica de gravedad en el
Cabecera MGD77 del crucero. Si esto no está configurado o no es válido, usamos por defecto el
IGF 1980. Alternativamente, especifique el campo directamente usando -Af (mira eso
opción para más detalles).

ceo Corrección de Eotvos calculada a partir de la navegación, usando E = 7.5038 * V *
cos (lat) * sin (az) + 0.004154 * V ^ 2 (en mGal). Signo: La corrección debe ser
añadido a faa sin corregir para producir un faa corregido.

También se reconocen las siguientes banderas de mano corta:

mgd77 Esto da como resultado que los 27 campos MGD77 se escriban en el MGD77 oficial
orden.

mgd77t El conjunto completo de las 26 columnas de la especificación MGD77T.

geo Esto limita la salida a 10 campos (equipo, lon, lat más los siete
observaciones geofísicas twt, profundidad, mtf1, mtf2, mag, tragosy faa). Por
adjuntando + a cualquiera de estos conjuntos también agregaremos dist, azim, cc, Acelerar,
y peso como se menciona arriba.

todos Esto devuelve todas las columnas de datos del archivo.

todo t Como mgd77t pero con elementos de hora escritos como una cadena de fecha y hora.

Como opción, se pueden agregar pruebas lógicas para cualquiera de las observaciones agregando
,lógica, que a su vez se compone de una o más instrucciones separadas por comas del
formulario porOPpropuesta de, Donde por es uno de los parámetros enumerados anteriormente, OP es una lógica
operador (<, <=, =,! =,> =,>, |), y propuesta de es una constante utilizada en la comparación.
Los parámetros de coma flotante se comparan numéricamente; los parámetros de carácter son
comparados léxicamente (después de eliminar los espacios en blanco iniciales y finales). El poco
comparación (|) significa que al menos uno de los bits en propuesta de debe estar encendido en
por. Al menos una de las pruebas debe ser verdadera para que se genere el registro, excepto para
pruebas utilizando parámetros de mayúsculas y minúsculas que deben ser verdaderos para que se produzca la salida. Nota
que especificar una prueba no implica que la columna correspondiente sea
incluido en el flujo de salida; debe estar presente en columnas para que eso ocurra.
Nota: algunos de los operadores son caracteres especiales de UNIX y se recomienda
colocar comillas alrededor de todo el argumento para -F.

Finalmente, para los archivos MGD77 +, puede agregar opcionalmente:mordiscos que es: (dos puntos)
seguido de uno o más separados por comas + -columna condiciones. Esto compara específicos
bitflags solo para cada columna de la lista. Aquí, + significa que el bit elegido debe ser 1 (ON)
mientras que - significa que debe ser 0 (APAGADO). Todas las pruebas de bits dadas deben aprobarse. Por defecto,
Archivos MGD77 + que tienen el especial MGD77_banderas la columna presente usará esas banderas,
y las observaciones asociadas con los bits ON (lo que significa que están marcados como incorrectos) serán
establecido en NaN; añadir: sin información final para desactivar este comportamiento (es decir,
no se consultarán banderas de bits).

OPCIONAL ARGUMENTOS

-A [+] c | d | f | m | tcódigo
De forma predeterminada, la profundidad corregida (profundidad), anomalía residual magnética (mag), Aire libre
anomalía de la gravedadfaa) y la cantidad derivada de corrección de profundidad de Carter (carretero)
son todas salidas como están (si se selecciona en -F); esta opción ajusta ese comportamiento. Para
En cada una de estas columnas hay de 2 a 4 formas de ajustar los datos. Adjuntar c(arteria),
d(epto), f(aa), o m(ag) y seleccione el código para el procedimiento que desea aplicar.
Puede seleccionar más de un procedimiento para una columna de datos sumando sus
códigos (1, 2, 4 y 8). P.ej, -C.A3 probaré primero el método -C.A1 para estimar un Carter
corrección pero si profundidad es NaN, lo siguiente que intentaremos -C.A2 que solo usa twt. En todo
casos, si alguno de los valores requeridos por un procedimiento de ajuste es NaN, entonces el
el resultado será NaN. Esto también es cierto si la anomalía original es NaN. Especificar -A +
para recalcular anomalías incluso si la anomalía en el archivo es NaN. Además, usted
puede utilizar -A crear tiempos falsos para cruceros que no tienen tiempo; estos se basan en
distancias y duración del crucero.

-C.A Determina cómo carretero Se calcula el plazo de corrección. Abajo, C (twt) representa
la profundidad corregida por Carter (también depende de lon, lat), U (twt, v) Es la
profundidad sin corregir (= twt * v / 2) usar como v el parámetro "Velocidad del sonido supuesta"
en el encabezado MGD77 (si es una velocidad válida, de lo contrario, el valor predeterminado es 1500 m / s);
alternativamente, agregue su velocidad preferida v en m / s, TU (profundidad, v) es el bidireccional
tiempo de viaje estimado a partir del (presumiblemente) no corregido profundidady TC (profundidad) Es la
El tiempo de viaje de 2 vías obtenido al invertir el (presumiblemente) corregido profundidad usando el
Fórmula de corrección de Carter. Seleccionar de

-Ac1[,v] devuelve la diferencia entre U (twt, v) y profundidad [Defecto].

-Ac2[,v] devuelve la diferencia entre U (twt, v) y Carter (twt).

-Ac4[,v] devuelve la diferencia entre (asumido sin corregir) profundidad y Carter
(TU (profundidad)).

-Ac8[,v] devuelve la diferencia entre U (TC (profundidad), v) y profundidad.

-Anuncio Determina cómo profundidad se obtiene la salida de la columna:

-Anuncio1 devoluciones profundidad tal como está almacenado en el conjunto de datos [Predeterminado].

-Anuncio2[,v] devuelve la profundidad calculada sin corregir U (twt, v).

-Anuncio4 devuelve la profundidad corregida calculada C (twt).

-Af Determina cómo faa se obtiene la salida de columna. Si ngrav (es decir, la Internacional
Se requiere campo de referencia de gravedad (IGF), o "gravedad normal") si se selecciona
basado en el parámetro de encabezado MGD77 "Código de fórmula de gravedad teórica"; si esto
el código no está presente o no es válido, el valor predeterminado es 4. Alternativamente, agregue el
preferido campo (1-4) para seleccionar 1 (Heiskanen 1924), 2 (IGF 1930), 3 (IGF 1967) o 4
(IGF 1980). Seleccionar de

-Af1[,campo] devoluciones faa tal como está almacenado en el conjunto de datos [Predeterminado]. Opcionalmente, establece el
IGF campo para usar si también ha solicitado ngrav como columna de salida en -F.

-Af2[,campo] devuelve la diferencia entre tragos y ngrav (con opcional campo
directiva).

-Af4[,campo] devuelve la combinación de tragos + EOT - ngrav (con opcional campo
directiva).

-Af8[,campo] devuelve la combinación de tragos + pred_eot - ngrav (con opcional
campo directiva).

-Soy Determina cómo mag se obtiene la salida de columna. Puede haber uno o dos en total
mediciones de campo en el archivo (mtf1 y mtf2) y la columna msens puede decir
cuál es el sensor principal (1 o 2; también puede estar indefinido). Seleccionar de

-Am1 devoluciones mag tal como está almacenado en el conjunto de datos [Predeterminado].

-Am2 devuelve la diferencia entre mgfx y igf, Donde x es el sensor principal1
or 2) indicado por el msens campo de datos (predeterminado en 1 si no se especifica).

-Am4 devuelve la diferencia entre mgfx y igf, Donde x es el sensor2 or 1)
no indicado por el msens campo de datos (predeterminado en 2 si no se especifica).

-Amc [unidad] Aplicar una corrección que intente compensar el hecho de que el
El campo magnético no se adquirió en la misma posición que la referencia de navegación.
(posición del barco). Esto se logra re-interpolando el total magnético
campo a lo que habría tenido si se hubiera medido en la posición del barco
(recuerde, probablemente se midió metros por detrás). Debido a esto
paso de interpolación, mala navegación, es decir, demasiados puntos repetidos, pueden causar
nubes. Se toman medidas para minimizar este efecto pero no están al 100%
profe. El método de interpolación está controlado por el GMT predeterminado GMT_INTERPOLANT.
Adjuntar e por metro, f para los pies, k por km, m por millas, n para millas náuticas, o u
para pies topográficos [El valor predeterminado es e (metros)].

-Cf | g | e
Agregue un código de una letra para seleccionar el procedimiento para la distancia a lo largo de la trayectoria
cálculo (ver -N para seleccionar unidades):

f Distancias terrestres planas.

g Grandes distancias circulares [predeterminado].

e Distancias geodésicas en elipsoide GMT actual.

-dafecha de inicio
No enumere los datos recopilados antes fecha de inicio (aaaa-mm-ddBD (T) [hh: mm: ss]) [El valor predeterminado es
inicio del crucero]. Usar -DA para excluir registros cuyo tiempo no está definido (es decir, NaN).
[Los informes predeterminados de esos registros].

-Dbfecha de parada
No enumere los datos recopilados en o después fecha de parada (aaaa-mm-ddBD (T) [hh: mm: ss]).
[El valor predeterminado es fin de crucero]. Usar -DB para excluir registros cuyo tiempo no está definido
(es decir, NaN). [Los informes predeterminados de esos registros].

-E Coincidencia exacta: solo los registros de salida que coinciden con todas las columnas geofísicas solicitadas
[Registros de salida predeterminados que coinciden con al menos una de las columnas observadas].

-Georgiastartrec
No enumere registros antes startrec [El valor predeterminado es 0, el primer registro].

-GBStoprec
No enumere los datos después Stoprec. [El valor predeterminado es el último registro].

-Ipasar por alto
No tenga en cuenta ciertos formatos de archivo de datos. Adjuntar a | c | t ignorar MGD77
ASCII, MGD77 + netCDF o archivos de tabla ASCII simples separados por tabulaciones, respectivamente. los
La opción puede repetirse para ignorar más de un formato. [El valor predeterminado ignora ninguno].

-L [cortable]
Aplique correcciones óptimas a las columnas donde dichas correcciones estén disponibles. Adjuntar
la tabla de corrección a utilizar [El valor predeterminado utiliza la tabla de corrección
mgd77_corrections.txt en el $ MGD77_HOME directorio]. Para el formato de este archivo,
ver CORRECCIONES a continuación.

-n Emita un registro de encabezado de segmento con ID de crucero para cada crucero.

-Nd | sunidad
Adjuntar d por distancia o s para velocidad, luego dale el unidad as e (metro o
Sra), f (pies o pies / s), k (km o km / h), m (millas o millas / h), n (millas náuticas
o nudos), o u (encuesta pies o pies / s). [El valor predeterminado es -Ndk -Nse (km y m / s)].

-Qamin/max
Especifique un rango aceptado (min/max) de acimutes. Registros cuyo azimut de trayectoria cae
fuera de este rango se ignoran [0/360].

-Qcmin/max
Especifique un rango aceptado (min/max) por supuesto cambios. Registros cuyo curso de pista
los cambios que caen fuera de este rango se ignoran [-360 / + 360]. Usar -control de calidad para asumir el
valor absoluto del cambio de curso antes de la prueba [El valor predeterminado utiliza un curso firmado
cambios].

-Qvmin/max
Especifique un rango aceptado (min/max; o solo min si no hay límite superior) de
velocidades. Los registros cuya velocidad de pista cae fuera de este rango se ignoran
[0 / infinito].

-R [unidad]oeste/este/sur/norte[/zmín/zmáx] [r]
oeste, este, sury norte especificar la región de interés, y puede especificar
en grados decimales o en formato [+ -] dd: mm [: ss.xxx] [W | E | S | N]. Adjuntar r si es menor
Las coordenadas del mapa izquierda y superior derecha se dan en lugar de w / e / s / n. Los dos
taquigrafía -Rg y -Rd significa dominio global (0/360 y -180 / + 180 de longitud
respectivamente, con -90 / + 90 en latitud). Alternativamente, para la creación de cuadrículas, dé
Rcodelón/lat/nx/ny, Donde código es una combinación de 2 caracteres de L, C, R (para la izquierda,
centro o derecha) y T, M, B para arriba, medio o abajo. por ejemplo, BL para la parte inferior izquierda.
Esto indica en qué punto de una región rectangular lon/lat coordinar se refiere
a, y las dimensiones de la cuadrícula nx y ny con espaciamiento de cuadrícula vía -I se usa para crear
la región correspondiente. Alternativamente, especifique el nombre de un archivo de cuadrícula existente
y del -R la configuración (y el espaciado de la cuadrícula, si corresponde) se copian de la cuadrícula.
Usar -Runidad espera coordenadas proyectadas (cartesianas) compatibles con las elegidas -J
y proyectamos inversamente para determinar la región geográfica rectangular real. Para
vista de perspectiva (-p), opcionalmente anexar /zmín/zmáx. En caso de vista en perspectiva
(-p), un rango z (zmín, zmáx) se puede agregar para indicar la tercera dimensión. Esta
debe hacerse solo cuando se utiliza el -Jz opción, no cuando se usa solo la -p .
En el último caso, se traza una vista en perspectiva del plano, sin una tercera
dimensión.

-sadist.inicio[unidad]
No enumere datos que sean inferiores a dist.inicio metro a lo largo de la vía desde el puerto de
salida. Adjuntar e por metro, f para los pies, k por km, m por millas, n para náutica
millas, o u para pies topográficos [El valor predeterminado es 0e (metros)].

-Sbparadadistante[unidad]
No enumere los datos que son paradadistante o más metros a lo largo de la vía desde el puerto de
salida. Adjuntar e por metro, f para los pies, k por km, m por millas, n para náutica
millas, o u para pies topográficos [El valor predeterminado es el final de la pista].

-T [m | e]
Desactiva el ajuste automático de valores basado en términos de corrección.
que se almacenan en el archivo MGD77 + y se utilizan para contrarrestar cosas como errores
unidades utilizadas por la institución de origen al crear el archivo MGD77 original a partir de
que deriva el archivo MGD77 + (la opción no tiene ningún efecto en MGD77 ASCII simple
archivos). Adjuntar m or e para limitar la opción al MGD77 o al conjunto de columnas extendidas
solo [El valor predeterminado se aplica a ambos].

-V [nivel] (más ...)
Seleccione el nivel de verbosidad [c].

-Wpeso
Establezca el peso de estos datos. La opción de salida de peso debe establecerse en -F. Es
útil si los datos se van a procesar con las técnicas de promediado ponderado
ofrecido por bloquear, medio bloquey Modo bloque [1].

-Z+|- Añade el letrero que quieras profundidad, carreteroy MSD valores por debajo del nivel del mar (Z-
da batimetría negativa) [El valor predeterminado es positivo hacia abajo].

-bo [ncoles][tipo] (más ...)
Seleccione salida binaria nativa. ignorado si -bo está seleccionado. Del mismo modo, los campos de cadena
no se puede seleccionar. Tenga en cuenta que si el tiempo es una de las columnas de salida binaria,
almacenarse como tiempo Unix (segundos desde 1970). Para leer esta información en GMT para
obtener el tiempo absoluto del calendario requerirá que utilice --TIME_SYSTEM = 1.

-h [i | o] [n] [+ c] [+ d] [+ rcomentario] [+ rtítulo] (más ...)
Omitir o producir registros de encabezado.

-^ or solo -
Imprima un mensaje corto sobre la sintaxis del comando, luego sale (NOTA: en Windows
usar solo -).

-+ or solo +
Imprima un mensaje de uso extenso (ayuda), incluida la explicación de cualquier
opción específica del módulo (pero no las opciones comunes de GMT), luego sale.

-? or no argumentos
Imprima un mensaje de uso completo (ayuda), incluida la explicación de las opciones, luego
salidas

--versión
Imprime la versión GMT y sal.

--show-datadir
Imprima la ruta completa al directorio compartido GMT y salga.

EJEMPLOS

Para obtener una lista (distancia, rumbo, gravedad, batimetría) de 01010047.mgd77, comenzando en
3 de junio de 1971 20:45 y terminando a una distancia = 5000 km, use el siguiente comando:

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

Para hacer entrada para bloquear y superficie utilizando anomalías al aire libre de todos los cruceros
enumerados en el archivo cruises.lis, pero solo los datos que están dentro del área especificada, y
hacer que la salida sea binaria:

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

Extraer las ubicaciones de profundidades superiores a 9000 metros que no fueron interpoladas (btc !=
1) de todos los cruceros enumerados en el archivo cruises.lis:

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

Para extraer dist, faa y grav12_2 de registros cuyas profundidades sean inferiores a 3 km y
donde ninguno de los campos solicitados es NaN, de todos los archivos MGD77 + netCDF cuyo crucero
Los identificadores se enumeran en el archivo cruises.lis, intentamos

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

Para extraer dist, faa y grav12_2 de todos los archivos MGD77 + netCDF cuyos ID de crucero son
enumerados en el archivo cruises.lis, pero solo recuperan registros cuyo bitflag para faa indique
MALOS valores, intentamos

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

Para generar lon, lat, mag y faa de todos los cruceros enumerados en el archivo cruises.lis, pero
recalcular los dos residuales en función de los últimos campos de referencia, intente:

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

RECALCULADO ANOMALIAS

Cuando se solicitan anomalías recalculadas (ya sea explícitamente a través del -A Opción o
implícitamente a través de metadatos E77 en el archivo MGD77 +) solo lo hacemos para los registros cuyo
La anomalía original no era un NaN. Esta restricción se implementa debido a que muchas columnas de anomalías
contiene correcciones, generalmente en forma de cambios editados a mano, que no se pueden
duplicado de la observación correspondiente.

IGRF

Los cálculos de IGRF se basan en un programa de Fortran escrito por Susan Macmillan, británica
Geological Survey, traducido a C vía f2c por Joaquim Luis, U Algarve, y adaptado a
Estilo GMT de Paul Wessel.

IGF

Las ecuaciones utilizadas se reproducen aquí utilizando coeficientes extraídos directamente de la
código fuente (avísenos si encuentra errores):

(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))]

CORRECCIONES

La tabla de corrección es un archivo ASCII con coeficientes y parámetros necesarios para realizar
correcciones. Se permiten los registros de comentarios que comienzan con #. Todos los registros de corrección son de
la forma

ID de crucero observación elimina necesidad de pilares angulados

donde ID de crucero es un prefijo NGDC, observación es una de las abreviaturas de geofísica
observaciones enumeradas en -F arriba, y elimina necesidad de pilares angulados consta de uno o más términoes eso
ser resumido y luego restado de la observación antes de la salida. Cada término debe tener
esta sintaxis exacta:

factor[* [función] ([escala](abrev[-natural])) [^industria ]]

donde los términos entre paréntesis son opcionales (los corchetes en sí no se utilizan, pero son regulares
los paréntesis deben usarse como se indica). No se permiten espacios excepto entre términos. La
factor es la amplitud de la función base, mientras que el opcional función puede ser uno de
pecado, cos o exp. El opcional escala y natural se puede utilizar para traducir el argumento
(antes de dárselo a la función opcional). El argumento abrev es una de las
abreviaturas de las observaciones enumeradas anteriormente. Si natural se da como T significa que nosotros
debería reemplazarlo con el valor de abrev para el primer registro en el archivo (esto es
generalmente solo se hace para equipo). Si la primera entrada de registro es NaN revertimos natural a cero.
Opcionalmente, eleve la expresión completa al valor dado industria , antes de multiplicar por el
amplitud. El siguiente es un ejemplo de correcciones ficticias al crucero 99999999,
implicando el profundidad debería eliminarse la corrección de Carter, faa debería tener un lineal
tendencia eliminada, la anomalía magnética (mag) debe corregirse por una extraña dependencia de
rumbo y latitud del barco, y trozos necesita tener 10 mGal agregados (por lo tanto, dado como -10):

99999999 profundidad 1.0 * ((carter))

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

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

99999999 gotas -10

Use mgd77listgmt en línea usando los servicios de onworks.net