sábado, 23 de febrero de 2019

4.2. Vectores y movimientos de placas

Se explica como dibujar vectores que indican el desplazamiento relativo entre las placas de Sudamérica y de Scotia. Se utilizan dos métodos: 1) de forma personalizada y 2) mediante el sitio UNAVCO.



Script

0. Mapa previo: Mapa previo a partir de una grilla con el límite transforme dibujado de Plates Project (ver 4.1).

gmt grdimage -R -J -O -K %CUT% -C%color% -I%SHADOW% >> %OUT%
gmt psxy -R -J -O -K "transform.gmt" -A  >> %OUT% -W0.4,30 -Sf1.5/0.25+l+s -Gblack




1. Vectores personalizados 1: Los datos se ingresan manualmente.

1A. Datos: longitud, latitud (del origen del vector), acimut (en grados desde el norte), largo del vector (en cm). Los datos se guardan en el archivo temp_vector.

echo    -68,-55,45,2c >  temp_vector
echo    -68,-55,50,3c >> temp_vector



1B. Graficar vectores: con psxy se grafican con los siguientes argumentos:
  • -SV+e: dibujar cabeza de vector de 0,3 cm de tamaño.
  • -W: propiedades de la línea.
  • -G: color del relleno de la cabeza del vector.
gmt psxy -J -R -O -K >> %OUT% "temp_vector" -SV0.3+e -W0.4p,black -Gorange



2. Vectores personalizados 2: echo y psxy se pueden combinar en único comando mediante una pleca. Con este método se pueden indicar los parámetros de cada vector por separado. Se utilizan los siguientes modificadores de los vectores:
  • +bl: Mitad izquierda (left) del vector al inicio (beginning) de la línea.
  • +er: Mitad derecha  (rigth) del vector al final (end) de la línea.
echo    -66,-53.5,85,4c | gmt psxy -J -R -O -K >> %OUT% -SV1+bl -W1.5p,red -Gred
echo    -66,-55.0,85,4c | gmt psxy -J -R -O -K >> %OUT% -SV1+er -W1.5p,red -Gred




3. Vectores UNAVCO: Los datos se descargan desde el sitio web (ver Anexo 5).

3A. Datos: Utilizo dos conjuntos de datos.
3B. Graficar: se utiliza psvelo con los siguientes argumentos:
  • -W: atributos de la flecha. Por defecto línea negra, continua.
  • -Se: Selecciona el significado de las columnas de los datos de entrada. Para los de UNAVCO se utiliza e. Los valores indican la escala de los vectores (0,25), el límite de confidencia del 95% (0.95) y el tamaño de fuente para el nombre del sitio (0 para que NO anote el nombre del sitio).
  • -G: Color de la cabeza del vector.
gmt psvelo -J -R -O -K >> %OUT% "UNAVCO_MODEL.txt" -W -Se0.25/0.95/0 -Gwhite
gmt psvelo -J -R -O -K >> %OUT% "UNAVCO_EULER.txt"  -W -Se0.25/0.95/0 -Ggreen

 


Referencias:
UNAVCO Plate Motion Calculator (link).
UNAVCO Plate Tectonics and GPS 'Reference Frames' (link).

2 comentarios:

  1. Hola Federico, estuve haciendo este mapa y tengo el problema de que no me presenta las áreas en relieve (creo que es el efecto hillshade) como se ve que a tí te sale. A mí me sale plano todo. Tienes alguna sugerencia? Gracias

    ResponderBorrar
  2. Hola. Si te sale todo plano, entonces te debería faltar el dem de tu zona. Fijate de descargarlo (https://mapasgeologicos.blogspot.com/2017/01/anexo-4-descargar-grilla-gmrt.html).

    ResponderBorrar