#!/bin/bash
################################################################################
gmt begin GMT6_Canaries_Madere

# Location of the GMT binaries: 
#gmtbin=/Applications/GMT-6.1.1.app/Contents/Resources

gmt set FONT_ANNOT_PRIMARY 10p FONT_ANNOT_SECONDARY 8p FONT_TITLE 16p,AvantGarde-Book FONT_LABEL 8p FORMAT_GEO_MAP ddd:mm:ssF
gmt set PS_CHAR_ENCODING ISOLatin1+
# gmt set PS_MEDIA A4
#
LANG=C
######Type de projection (Mercator) et centre de la carte
projection=-JM-16/35/12c
######Région cartographiée
region=-R-27/-5/25/45
# Resolution des données raster
raster_inc=-I1m

# main directory for data, mainly raster (grd) type: 
rasterpath=$HOME/wrk/data
grid=$rasterpath/etopo1/ETOPO1_Ice_g_gmt4.grd

#Ombrage
gmt grdgradient $grid $region -A0/30 -GCanary_intens.grd -Ne0.6

# Appel des données topo
gmt grdimage $grid $region $projection -Cgeo -ICanary_intens.grd -V

# tracé de la carte
gmt coast $region $projection -Di -N1/thinner,sienna -Ir/thin,steelblue -Wthin,black -V

# Age des fonds océaniques (anomalies magnétiques)
gmt grdcontour $rasterpath/seafloor_age/age.3.6.grd $region $projection -C10 -A10+f8p,Helvetica,white -Wthin,white

#Echelle d'altitude
gmt colorbar -C -Dx9c/13c+w4c/0.25c+h+e+jCT -F+gghostwhite+r -Bxaf -By+l"m" -V

# Noms det âges des îles des Canaries 
gmt text $region $projection -Dj0.4vthick,black -Wklimegreen \
-Gwhite -F+f8pt,Helvetica,forestgreen+a+j -V << CANYDAT 
-14.019 28.181 0 LM Fuerteventura, 24 Ma
-13.618	28.9 0 LB Lanzarote, 24 Ma
-15.601	27.817 0 LT Grande Canarie, 15 Ma
-16.606 28.2 0 CB Tenerife, 12 Ma
-17.133 27.948 0 RM Gomera, 10 Ma
-17.869 28.519 0 RM La Palma, 2 Ma
-18.051 27.584 0 CT Hierro, 1 Ma
-16.00 30.074 0 RB Selvagen, 30 Ma
-12.85 30.78 0 LM Anika, 55 Ma
-13.18 32.8 0 LM Lars, 68 Ma
CANYDAT

# Position des îles des Canaries
gmt plot $region $projection -Wthick,black -Glimegreen -Sc0.2c << CANY
-14.019	28.181
-13.618	28.9
-15.601	27.817
-16.606 28.2
-17.133 27.948
-17.869 28.519
-18.051 27.584
-16.00 30.074
-12.85 30.78
-13.18 32.8
CANY

# Noms et âges des îles de Madère 
gmt text $region $projection -Wkred -Dj0.3vthick,red \
-Glightgoldenrodyellow -F+f8pt,Helvetica,red+a+j -V << MADRDAT
-17.001 32.605 0 RM Mad\350re, 5 Ma
-16.359 32.919 0 RB Porto Santo, 15 Ma
-14.39 33.76 0 RB Seine, 22 Ma
-14.5 34.63 0 RB Unicorn, 27 Ma
-12.0612 34.896 0 LM Amp\350re, 31 Ma
-11.38 36.6 0 CB Ormonde, 65-67 Ma
MADRDAT

# Position des îles de l'archipel de Madère
gmt plot $region $projection -Wthick,black -Gred -Sc0.2c << MADR
-17.001 32.605
-16.359 32.919
-14.39 33.76
-14.5 34.63
-12.0612 34.9
-11.38 36.6
MADR

# 
gmt basemap $region $projection -Bxaf -Byaf -B+t"Archipels des Canaries et de Mad\350re"

# Conversion du resultat en fichiers PDF [option -Tf] ou PNG [option -Tg] ou JPG [-Tj].
gmt psconvert -Tj -A -P -FGMT6_Canaries_Madere

# Nettoyage des fichiers intermédiaires :
rm -f .gmt* *_intens.grd

# fin du script
gmt end show