19. Distancia y dirección a La Meca

Dibujar curvas de equidistancia y vectores de dirección hacia La Meca y la terminación día-noche.

Introducción: Se explica como calcular distancias a un punto geográfico (La Meca en este caso), dibujar curvas de equidistancia de color, vectores de dirección hacia ese punto y el terminador día-noche para el solsticio de junio de 2018.

Script (link)

1. Mapa base: mapamundi base con áreas secas y húmedas pintadas, línea de costa y límite de países (ver figura 1.1).

2. Grilla distancia: grdmath calcula la distancia (en km) a un punto geográfico (SDIST) definido por las variables "long" y "lat". El resultado se guarda en la grilla "temp_dist.nc" que tiene una resolución de 1 minuto de arco (-I1).

SET long=39.826167
SET lat=21.4225
gmt grdmath -Rg -I1 %long% %lat% SDIST = "temp_dist.nc"

3. Curvas de nivel de Color: Se crea una CPT para colorear las curvas de nivel cada 1000 km. Luego con grdcontour se dibujan a partir de la grilla "temp_dist.nc".

• -C%color%: Utiliza los colores y los intervalos de la CPT para colorear las curvas.
• -A: define el texto de las anotaciones con la unidad (+u" km"), un borde redondeado (+p+o), tamaño de fuente 6 (+f6) y fondo blanco con transparencia del 50% (+gwhite@50).
• -W: define el tipo de línea (thin) para las curvas con anotaciones y toma los colores de la CPT para dibujar las líneas (+cl).

gmt makecpt -T0/20000/1000 > %color% -Crainbow -I
gmt grdcontour "temp_dist.nc" -R -J -O -K >> %OUT% -C%color% -A+u" km"+f6+p+o+gwhite@50 -Wthin+cl

4. Vectores dirección: los vectores de dirección se dibujan a partir de la grilla de distancia. Primero es necesario calcular las derivadas direccionales en X y en Y (DDX, DDY con grdmath). Luego se dibujan con grdvector con los siguientes argumentos:

• -Gblack: color de relleno de los vectores.
• -W: estilo de los vectores.
• -S: longitud de los vectores.
• -Ix20: ubicación de los vectores cada 20 nodos de la grilla (20 x 1 minuto).
• -Q: atributos de la cabeza de los vectores.
• MAP_VECTOR_SHAPE=2: Tipo de cabeza de vector.

gmt grdmath "temp_dist.nc" DDX -1 MUL = "temp_dzdx.nc"
gmt grdmath "temp_dist.nc" DDY -1 MUL = "temp_dzdy.nc"
gmt grdvector -R -J -O -K "temp_dzdx.nc" "temp_dzdy.nc" >> %OUT% -Gblack -W0.5p -S0.2c -Ix20 -Q0.1c+e+jc --MAP_VECTOR_SHAPE=2.0

5. Día-Noche: con pssolar se agrega el terminador día-noche.

• -W1p: propiedades de la línea terminador.
• -G: pinta de negro la zona de noche.
• -t50: transparencia (del 50%) de toda la capa.
• -Td: Dibuja el tipo de terminador día-noche (d). Opcionalmente se puede indicar en que fecha calcular el terminador (+d). En este ejemplo se calcula para el solsticio de invierno en el hemisferio sur (21 de junio de 2018 a las 10:07 hs).

gmt pssolar -R -J -O -K >> %OUT% -W1p -Gblack -t50 -Td+d2018-06-21T10:07:00

6. Proyección ortográfica: Se puede rehacer la figura con proyección ortográfica centrada en La Meca con intervalos cada 500 km, modificando PROJ y makecpt, respectivamente.

SET PROJ=G39.826167/21.4225/90/15c

gmt makecpt -T0/20000/500 > %color% -Crainbow -I