#!/bin/bash

ps=mapa-exa.ps

# Vamos definir aqui o valor da exageração
EXAG=3

# Faça um mapa base com a exageracao desejada (neste caso estamos utilizando 1:6)
psbasemap -Jx$EXAG/1 -R0/2.5/-3/0 -Ba.5f.1:"Distance along profile (km)":/a1f.5:"Depth (km)":WSne -K -P -Xc -Yc > $ps

# Para cada um dos angulos a seguir .... 
for angulo in 0 30 60 90
do
    # Obtenha o angulo aparente na exageracao utilizando a fórmula ( Ap = Atan( tan(angulo) / exageracao )
    ap=$(echo "$angulo" | gmtmath STDIN TAND = | gmtmath STDIN $EXAG DIV = | gmtmath STDIN ATAN = | gmtmath STDIN R2D = | gmtmath STDIN 90 ADD =)
    echo "O angulo $angulo tem um angulo aparente na exageração escholhida de $(echo $ap|gmtmath STDIN 90 SUB =)"
    psxy -R -J -W1,black -Gblack -SV0/0/0 -O -K << END >>$ps 
0.5 -1.0 $ap 1
END
done

#Fecha o script do GMT
psxy /dev/null -R -J -O >> $ps

gv $ps