3.6. Diagrama Day 1977

Introducción: El diagrama de Day et al. (1977) sirve para determinar el estado de dominio de la de una muestra para magnetita. Se realiza a partir de los parámetros de un ciclo histéresis para diferenciar granos de dominios simple (SD), pseudo-simple (PSD) y múltiple (MD).
 Es ejemplo es una adaptación de la entrada 3.1. Grafico Cartesiano.

Script

Preparación de Datos: Los datos deben estar en archivos de texto simple (bloc de notas) en 2 columnas (separados por una tabulación). Las líneas del encabezado conviene que se inicien con el símbolo #.

• Jrs/Js    : Magnetización Remanente de Saturación sobre Magnetización de Saturación.
• Hrc/Hc : Coercitividad de la remanencia sobre Coercitividad.

1. Variables del gráfico, título, ejes y fondo del gráfico: Se utilizan los siguientes argumentos de psbasemap. Ver explicación en 3.1. Grafico Cartesiano.

• -B+t: Título
• -B+g: Color del fondo
• -By: Anotaciones y leyenda del eje Y
• -Bx: Anotaciones y leyenda del eje X

REM    Dominio de datos (eje X e Y)
SET    REGION=0/5.9/0/0.65
REM    Proyeccion Lineal (X). Ancho y alto (en cm)
SET    PROJ=X15c/10c

gmt psbasemap -R -J -O -K >> %OUT% -B+t"Diagrama de Day et al. 1977"
gmt psbasemap -R -J -O -K >> %OUT% -B+g200 
gmt psbasemap -R -J -O -K >> %OUT% -Bxa1f1+l"H@-RC@-/H@-C@-"
gmt psbasemap -R -J -O -K >> %OUT% -Bya0.1f1+l"J@-RS@-/J@-S@-" 

2. Campos del Diagrama Day: Primero se crea un archivo auxiliar con las líneas (1,5 y 4 en Y, y 0,05 y 0,5 en X) que definen los campos que se guardan en el archivo temporal "temp_Lineas". Luego, este archivo se grafican con psxy -W.

REM    Crear archivo auxiliar para las lineas
echo 1.5,0   >  "temp_Lineas"
echo 1.5,10  >> "temp_Lineas"
echo         >> "temp_Lineas"
echo 4,0     >> "temp_Lineas"
echo 4,10    >> "temp_Lineas"
echo         >> "temp_Lineas"
echo 0,0.05  >> "temp_Lineas"
echo 10,0.05 >> "temp_Lineas"
echo         >> "temp_Lineas"
echo 0,0.5   >> "temp_Lineas"
echo 10,0.5  >> "temp_Lineas"

REM    Graficar archivo auxiliar
gmt psxy -R -J -O -K "temp_Lineas" -Wthin >> %OUT%

3. Nombres de los campos: Del mismo modo (que en el punto 2), se crea un archivo auxiliar con los nombres de los campos (SD, PSD, MD) y su correspondiente posición. Luego, con pstext se grafican.

• -F+f: Define el tamaño del texto (en puntos tipográficos).

REM    Crear archivo auxiliar con los nombres de los campos
echo 0.75,0.575,SD  >  "temp_Campos"
echo 2.75,0.275,PSD >> "temp_Campos"
echo 5.00,0.025,MD  >> "temp_Campos"

REM    Poner Nombre de los Campos
gmt pstext -R -J -O -K >> %OUT% "temp_Campos" -F+f12

4. Graficar Datos: con psxy se dibujan los datos con símbolos con los siguientes argumentos:

• -S: forma y tamaño (en cm)
• -G: color de relleno
• -W: borde fino
• -: intercambia el orden de entrada de los datos. Por defecto, se espera que los datos estén ordenados en XY. En este ejemplo, los datos están en YX por lo tanto es necesario invertirlo.

gmt psxy -R -J -O -K "Datos.txt" >> %OUT% -Sc0.25 -Gdeepskyblue3 -Wthinnest -:

Referencias:
Day R., Fuller M., Schmidt V.A. 1977. Hysteresis properties of titanomagnetites: Grain-size and compositional dependence, Physics of the Earth and Planetary Interiors, 260-267. https://doi.org/10.1016/0031-9201(77)90108-X