11.2. Perfiles Multiples.

Extraer y crear perfiles topográficos, de pendientes, magnetométricos y gravimétricos a partir de las respectivas grillas y agregar coordenadas al perfil.

El script permite crear 5 perfiles E-W, sobre el paralelo 32 °S entre 76° y 66 °W. Los perfiles incluyen topografía, pendientes, magnetometría, anomalías de Aire Libre, Gradiente Vertical de la Gravedad (VGG). Esta entrada es una adaptación de la anterior (11.1) para graficar múltiples perfiles en lugar de uno sólo. Se incluye como agregar las coordenadas geográficas que indican la orientación del perfil.

0. Grillas: Se definen 4 variables con la ruta absoluta (según mi sistema de archivos) de las grillas GEBCO (DEM), EMAG2 (MAG), Anomalías de Aire Libre de Sandwell y Smith v. 29.1 (GRA) y de gradiente Vertical (VGG). Ver Base de Datos para más información. También se definen otras 4 variables para las grillas que necesitan ser procesadas.

REM    Base de datos de GRILLAS
    SET    DEM="E:\Facultad\Datos_Geofisicos\Batimetria\GEBCO\GEBCO_2014_1D.nc"
    SET    MAG="E:\Facultad\Datos_Geofisicos\Magnetometria\Geomag\EMAG2_2009\EMAG2_V2.grd"
    SET    GRA="E:\Facultad\Datos_Geofisicos\Gravimetria\Sandwell-Smith\Datos\grav.29.1.img"
    SET    VGG="E:\Facultad\Datos_Geofisicos\Gravimetria\Sandwell-Smith\Datos\curv.29.1.img"

REM    Grillas Temporales
    SET    TPEN=temp_Pend.grd
    SET    TGRA=temp_Grav.grd
    SET    TMAG=temp_Mag.grd
    SET    TVGG=temp_VGG.grd

1. Procesar Grillas: Es necesario procesar las grillas para extraer los datos. Se utilizan los comandos utilizados en entradas (6.1, 7.1, 7.2, 8.1) previas para preparar las grillas (salvo la topográfica). Las grillas definen una franja que coincide con la extensión del perfil y con un ancho de 0.2° (-R-76/-66/-32.1/-31.9) que luego se utilizarán para extraer los datos.
   
REM    Grilla de Pendientes  
    gmt grdgradient %DEM% -D -Stemp_mag.grd -R-76/-66/-32.1/-31.9 -fg
    gmt grdmath temp_mag.grd ATAN R2D = %TPEN%  

REM    Grilla Magnetometria
    gmt grdconvert %MAG%=nf %TMAG%=nf/1/0/-99999 -R

REM    Grilla Aire Libre
    gmt img2grd %GRA% -R -G%TGRA% -T1 -I1 -E -S0.1

REM    Grilla Gradiente Vertical de la Gravedad (VGG)
    gmt img2grd %VGG% -R -G%TVGG% -T1 -I1 -E -S0.02


2. Extraer datos perfil: Los pasos a seguir son los mismos de la entrada anterior (11.1). Sólo es necesario adaptar el siguiente grdtrack, agregando las 5 grillas (temporales) de las cuales se extraen los datos a graficar (variables %DEM%, %TPEN%, %TMAG%, %TGRA%, %TVGG%).

    gmt grdtrack temp_track -G%DEM% -G%TPEN% -G%TMAG% -G%TGRA% -G%TVGG% > temp_data

3. Información Perfil: gmtinfo nos brinda información de cada columna del archivo temp_data (ver figura). Los datos de las columnas 3 a 8 son utilizados luego para crear los distintos perfiles.

El comando pause permite ajustar las variables del script (punto 4) antes de que se siga ejecutando.
 

4. Definir Variables Figuras: Definir variables que se utilizarán para hacer los gráficos: Longitud (L=15 cm), Altura (H= 4 cm). La longitud representada en los perfiles (945 km), y los rangos de valores de los perfiles (Topo, Pend, Magn, Grav, VGra) que se extraen del comando gmtinfo (punto 3).

REM    Dimensiones del Grafico: Longitud (L), Altura (H).
    SET    L=15C
    SET    H=4C
REM    Datos del perfil segun gmtinfo
    SET    KM=945
    SET    Topo=-6200/5300
    SET    Pend=-1/27
    SET    Magn=-270/172
    SET    Grav=-169/450
    SET    VGra=-62/90

 
5.  Gráfico y Eje X: Crear el gráfico cartesiano con las dimensiones de los gráficos (-JX%L%/%H%) y la región, que en este contexto, corresponde al rango de valores de los perfiles (km en eje X).

REM    Crear Grafico
    gmt psxy -JX%L%/%H% -R0/%KM%/%Topo% -T -K -P > %OUT%

REM    Dibujar Eje X (Inferior: S)
    gmt psbasemap -R -J -O -K -Bxaf+l"Distancia (km)" >> %OUT%  --MAP_FRAME_AXES=S


6. Dibujar Ejes Y: Se configura MAP_FRAME_AXES para que en los gráficos se escriban las anotaciones sólo en el eje derecho (wEsn, que correspondería al Este si fuera un mapa). Luego, con psxy se grafican los 5 perfiles siempre a partir del mismo archivo (temp_data). La opción -By permite especificar el intervalo de anotación (a), marcas (f) y lineas de grillado (g) y la leyenda del eje (+l""). El argumento -i permite indicar en que columna está cada dato (que dependerá del orden en que los datos fueron extraídos de las grillas en el punto 2). -Y%H% desplaza hacia arriba los subsiguientes gráficos según la variable de la altura de los gráficos (H). El argumento -W permite especificar las propiedades de las líneas (ancho, color y estilo).

REM    Dibujar Eje Y y datos de columnas 3 y 4 (-i2,3)
    gmtset MAP_FRAME_AXES wEsn
    gmt psxy -R0/%KM%/%Topo% -J -O -K "temp_data" -i2,3 -W1,blue   -Bya2500f500g2500+l"Topograf\355a (m)" >> %OUT%
    gmt psxy -R0/%KM%/%Pend% -J -O -K "temp_data" -i2,4 -W1,red    -Byafg+l"Pendiente (\232)"      >> %OUT% -Y%H%
    gmt psxy -R0/%KM%/%Magn% -J -O -K "temp_data" -i2,5 -W1,orange -Byafg+l"Magnetometr\355a (nT)" >> %OUT% -Y%H%
    gmt psxy -R0/%KM%/%Grav% -J -O -K "temp_data" -i2,6 -W1,green  -Byafg+l"Anomal\355a Aire Libre (mGal)" >> %OUT% -Y%H%
    gmt psxy -R0/%KM%/%Vgra% -J -O -K "temp_data" -i2,7 -W1,violet -Byafg+l"Grad. Vert. Grav. (E\366tv\366s)" >> %OUT% -Y%H

7. Orientación Perfil: Para agregar las coordenadas geográficas que indiquen la orientación del perfil (E y O en este caso) se utilizan los comandos echo y pstext vinculados a través de una pleca o barra vertical ("|") (ver pipe en GMT o tubería). De esta manera el texto de la salida del primer comando ("E") es ingresado en el segundo (pstext) el cual nos permite definir las características para graficarlo: -F permite indica la fuente a utilizar (+f14p; tamaño 14 puntos) y su ubicación dentro del gráfico (+cTR; Top Right para "E" y +cTL; Top Left para "W"). -W dibuja un rectángulo debajo del texto y con -G se indica el color de fondo.

REM    Coordenadas Perfil (E, O)
    echo E | gmt pstext -R -J -O -K -F+cTR+f14p -Gwhite -W1 >> %OUT%
    echo O | gmt pstext -R -J -O -K -F+cTL+f14p -Gwhite -W1 >> %OUT%