8.2. Mapa Magnetométrico. EMAG2 v3

Elaborar mapa magnetométrico del modelo EMAG2 v3 (Earth Magnetic Anomaly Grid).

Mapa Magnetométrico (EMAG2 v3) con anomalías sobre el nivel del Mar en océanos y a 4 km de altitud sobre continentes.

Introducción: Se explica como graficar la nueva versión del modelo de anomalías magnetométricas EMAG2 (v3, Meyer et al. 2016; en 8.1 se graficó la versión 2). Esta versión está conformada por 2 grillas:

1. Datos sobre el nivel del mar en océanos y a 4 km de altitud en continentes
2. Datos continuos a 4 km de altitud.

Link Script


Datos: Los datos de la versión 3 del EMAG2 están en el archivo EMAG2_V3_20160623.csv (comprimido en un zip de 1.5 Gb en este link) y constan de 8 columnas (readme):

Column 1: i ; grid column/longitude index
Column 2: j ; grid row/latitude index
Column 3: LON ; Geographic Longitude WGS84 (decimal degrees)
Column 4: LAT ; Geographic Latitude WGS84 (decimal degrees)
Column 5: SeaLevel ; Magnetic Anomaly Value at Sea Level over oceans / 4km altitude over continents (nT)
Column 6: UpCont ; Magnetic Anomaly Value at continuous 4km altitude (nT)
Column 7: Code ; Data Source Code (see table below)
Column 8: Error ; Error estimate (nT)

0. Comandos Previos: Para graficar el mapa se utilizan los siguientes comandos (explicados en detalle en 8.1) con la excepción de la preparación de los datos (ver 1).

gmt grd2cpt -Crainbow -Z %CUT% -Di -L-200/200 >%color%
gmt grdgradient %CUT% -A0/270 -G%SHADOW% -Ne0.5 
gmt grdimage -R -J -O -K %CUT% -C%color% -I%SHADOW% >> %OUT%
gmt psscale -O -K -C%color% -Baf+l"Magnetometr\355a (nT)" -I >> %OUT%-Dx15.7/1.5+w11.9/0.618c+n"Sin datos"+e
gmt psxy -R -J -O -K "E:\Facultad\Datos_Geofisicos\IGN\1_GMT\limites_politico_administrativos.gmt" -Wthinner >> %OUT% 
gmt psxy -R -J -O -K "E:\Facultad\Datos_Geofisicos\IGN\1_GMT\LIMITE_COPLA.gmt" -Wthinner,.- -Am >> %OUT% 
gmt pscoast -R -J -O -K -Df -N1/thinner >> %OUT%
gmt pscoast -R -J -O -K -Df -W1/thinner >> %OUT%
gmt psbasemap -R -J -O -K -Baf >> %OUT%
gmt psbasemap -R -J -O -K -Ln0.11/0.075+c-32:00+w800k+f+l -F+gwhite+p+i    >> %OUT%



1. Preparar Datos: xyz2grd permite crear una grilla a partir del archivo "EMAG2_V3_20160623.csv". Se utilizan los siguientes argumentos:

• -I: resolución de la grilla (2 minutos de arco).
• -G: nombre de la grilla de salida.
• -i: selecciona las columnas de datos según que se desea graficar.
• -hi1: indica que el encabezado tiene 1 línea.
• -di9999: define el valor correspondiente a NaN (sin datos).

Nota: Las columnas (con -i) se empiezan a contar desde 0.


1A. Nivel del Mar: Datos en columnas 3, 4 y 5 (-i2,3,4).

gmt xyz2grd -R "EMAG2_V3_20160623.csv" -I2m -G%CUT% -fg -hi1 -r -di99999 -i2,3,4 

1B. 4 km altitud: Datos en columnas 3, 4 y 6 (-i2,3,5).

gmt xyz2grd -R "E:\Facultad\Datos_Geofisicos\Magnetometria\Geomag\EMAG2_V3\EMAG2_V3_20160623.csv" -I2m -G%CUT% -fg -hi1 -r -di99999 -i2,3,5

Referencias:
Brian Meyer, Richard Saltus, Arnaud Chulliat (2016): EMAG2: Earth Magnetic Anomaly Grid (2-arc-minute resolution) Version 3. National Centers for Environmental Information, NOAA. Model. doi:10.7289/V5H70CVX