3.3. Gráfico Tierras Raras.

Este ejemplo muestra como realizar un gráfico de tierras raras normalizado. Se explica como agregar anotaciones personalizadas en un eje y como hacer cálculos básicos (normalizar los datos).

Preparación de los Datos:
Para hacer el gráfico se necesitan 3 archivos de tipo de tablas de texto simple (ASCII).
A. Eje Tierras Raras: El primero archivo (REE) será utilizado para personalizar las anotaciones (a) del eje X y contiene 3 columnas. La primera con los números del 1 al 15 y la tercera con los símbolos químicos de las tierras raras (La, Ce, Pr, etc). La "a" de la segunda columna indica que las anotaciones del 1 al 15 serán reemplazadas por los símbolos químicos.


#Tierras Raras
1 a La
2 a Ce
3 a Pr
4 a Nd
5 a Pm
6 a Sm
7 a Eu
8 a Gd
9 a Tb
10 a Dy
11 a Ho
12 a Er
13 a Tm
14 a Yb
15 a Lu


B. Condrita: El segundo archivo (Condrita.txt) contiene los valores de concentración (columna 2) de las tierras raras de la condrita que se utilizará para normalizar los datos. La primera columna tiene los números correspondiente a cada Tierra Rara según el archivo REE.
Nota: En este archivo se omite la fila 5, correspondiente al Prometio (Pr) dado que es un elemento que no se encuentra en la naturaleza.
Nota 2: Los datos de la Condrita son sólo para hacer el gráfico. Cada uno debería actualizar los datos según la fuente que quiera utilizar.

1    0.237
2    0.613
3    0.0928
4    0.457
6    0.148
7    0.0563
8    0.199
9    0.0361
10    0.246
11    0.0546
12    0.16
13    0.0247
14    0.161
15    0.0246

C. Muestra: El tercer archivo es similar al B (AB1.txt en este ejemplo) pero con los datos de la muestra en la columna 2. Es importante que estos dos últimos archivos tengan la misma cantidad de filas (14) para poder normalizar los datos (punto 6).



Descripción del Script (Link):

Variables de la Figura: 


1. Región: Los valores del eje Y (que será logarítmica) son 1 a 1000 para este ejemplo. Eventualmente estos valores se puede editar para incluir datos que queden fuera del gráfico o para ajustarlo. En cambio, para el eje X elegimos los valores de 0.5 a 15.5. Los valores 1 a 15 serán asignados (punto 5) a los distintos elementos de las tierras raras.

rem    Region Geografica/Figura (WESN)
    SET    REGION=0.5/15.5/1/1000


2. Proyección: Utilizamos la proyección no geográfica (X) con de 7.5 cm ancho y 5 de alto. Para el eje Y utilizamos escala logarítmica.


REM    Proyeccion No Geografica. Linear, Logaritmica, Exponencial. JXwidth[/height]
    SET    PROJ=X7.5c/5cl


Hacer Figura:

3. Comando Previos: Agregamos el fondo y el título del gráfico según lo explicado previamente.

REM    Definir color de fondo del grafico (+g) 
    gmt psbasemap -R -J -O -K >> %OUT% -B+g200

REM    Tìtulo de la figura
    gmt psbasemap -R -J -O -K >> %OUT% -B+t"Gr\341fico Tierras Raras"

4. Marco: En este ejemplo, a diferencia de los anteriores, queremos que los ejes superior y derecho (N y E) estén dibujados pero sin anotaciones (a) ni marco (f). Para dibujar un recuadro sin anotaciones, utilizamos el siguiente comando:

    gmt psbasemap -R -J -O -K >> %OUT% -B0wesn

5. Ejes: Para incluir los ejes X e Y primero configuramos que se dibujen los ejes inferior (S) e izquierdo (W) con la herramienta gmtset. Para el eje Y utilizamos el comando explicado previamente. Para personalizar al eje X con los valores personalizados (custom; c) del archivo REE (ver Preparación de Datos) con el comando -BxcREE.

    gmtset    MAP_FRAME_AXES WS
    gmt psbasemap -R -J -O -K >> %OUT% -BxcREE+l"Tierras Raras"
    gmt psbasemap -R -J -O -K >> %OUT% -Bya1f3+l"Roca/N-MORB"

6. Cálculos: Para graficar los datos, primero hay que normalizarlos en relación a la Condrita. Utilizamos gmtmath para dividir (DIV)  los datos de los archivos AB1.txt y Condrita.txt. Los datos se guardan en el archivo Muestra que serán graficados en el próximo paso. GMT utiliza la notación polaca inversa para hacer cálculos.

rem    Calcular. Normalizar valores de las tierras raras
    gmt math AB1.txt Condrita.txt DIV = Muestra

7. Graficar Datos: Los datos los graficamos 2 veces con la psxy explicado previamente. Primero, con -W, dibujamos una línea (thin,red) que una todos los puntos. Luego, dibujamos círculos (azules de 0,15 cm)  sobre la línea.

rem    Graficar Datos

    gmt psxy -R -J -O -K Muestra -Wthin,red >> %OUT%
    gmt psxy -R -J -O -K Muestra -Sc0.15 -Gblue >> %OUT%