1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
|
#!/bin/bash
# $Id: GMT_mapscale.sh 15178 2015-11-06 10:45:03Z fwobbe $
#
ps=GMT_mapscale.ps
gmt psbasemap -R0/40/50/56 -JM5i -Baf -P -K -LjML+c53+w1000k+f+l"Scale at 53\\232N" -F+glightcyan+c0+p > $ps
gmt psbasemap -R -J -O -K -LjBR+c53+w1000k+l+f -F+p1p+i+gwhite+c0.1i >> $ps
h=`gmt mapproject -R -J -Wh -Di`
h=`gmt math -Q $h 2 DIV =`
lat=`echo 0 $h | gmt mapproject -R -J -I -Di -o1`
gmt psxy -R -J -O -K -Wfaint -A -N << EOF >> $ps
>
0 $lat
40 $lat
>
0 50
40 50
EOF
gmt psxy -R -J -O -T >> $ps
|
