#!/usr/bin/env bash
#

#	Variables del Mapa
#	-----------------------------------------------------------------------------------------------------------
#	Titulo del mapa
	title=Movie_world_label_BUG

#	Grillas
	DEM="/home/federico/E:/Facultad/Datos_Geofisicos/Batimetria/SRTM15+v21/SRTM15+V2.1.nc"
	DEM=@earth_relief_15m_p

#	Region geografica y Dominio del grafico
	Dominio=0/1/0/360
	REGION=-67/-66/-55/-54

# 	Nombre archivo de salida
	CUT=Az_Pendientes.nc

#	-----------------------------------------------------------------------------------------------------------
#	Crear lista de fechas para la animacion: Inicio/Fin/Intervalo. o: meses. y: años
	gmt grdinfo -R$REGION -I1 -D0 > SR.txt

#	Crear Scrips pre con procesamiento previo
#	-----------------------------------------------------------------------------------------------------------
cat << EOF > pre.sh
	gmt set IO_NAN_RECORDS skip
gmt begin

#	1. Preparar Datos
#	-------------------------------------------------------------
#	Dirección de las pendientes
	gmt grdgradient $DEM -Dn -G$CUT -R$REGION

#	Extraer valores de rumbo de la grilla. Ponderar según área del nodo (-Wa). Incluir long y lat (-o0,1)
 	gmt grd2xyz $CUT > "Az_Pendientes.txt" -s -Wa -o0,1,3,2 -bof

gmt end
EOF

#	2. Scripts principal para el procesamiento
#	-----------------------------------------------------------------------------------------------------------	
cat << EOF > main.sh
gmt begin

#	Seleccionar Datos
	gmt select "Az_Pendientes.txt" \${MOVIE_TEXT} -o2,3 -fg > "temp_rumbo" -bi4f

#	Dibujar rosa.
	gmt rose "temp_rumbo" -Gorange -Wfaint -Bx0.2g0.2 -By30g30 -B+glightblue -A3r -S5cn -D -LW,E,S,N -JX10 -R$Dominio -Yc -Xc -T
  	
gmt end
EOF

#	----------------------------------------------------------------------------------------------------------
# 	3. Movie
	gmt movie main.sh -Sbpre.sh -C11.8cx11.8cx60 -TSR.txt -N$title -H2 -Fmp4 -D12 -Ml,png -Lt0+t"%s"

#	Borrar archivos temporales
	rm gmt.history SR.txt