-
Notifications
You must be signed in to change notification settings - Fork 0
/
Copy pathWWKS22_Fig_03.sh
50 lines (49 loc) · 1.27 KB
/
WWKS22_Fig_03.sh
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
#!/usr/bin/env -S bash -e
#
# Wessel, P., Watts, A. B., Kim, S.-S., and Sandwell, D. T., 2022
# Models for the Evolution of Seamount, Geophys. J. Int.
#
# Fig: Parabolic seamount model
# P. Wessel, April 2022
# Single-column figure in GJI so aim for W = 8.4 cm
# Determine if we need to specify an output directory or not
if [ "X${1}" = "X" ]; then
dir=
else
dir="${1}/"
fi
gmt begin ${dir}WWKS22_Fig_03 $1
gmt set MAP_VECTOR_SHAPE 0.5
echo "0 1" > /tmp/tmp
echo "1 1" >> /tmp/tmp
gmt math -T1/4.1/0.1 T 4 DIV 2 POW NEG 1 ADD 1 0.25 2 POW SUB DIV = >> /tmp/tmp
gmt math -T-4/-1/0.1 T 4 DIV 2 POW NEG 1 ADD 1 0.25 2 POW SUB DIV = /tmp/body
gmt math -T1/4/0.1 T 4 DIV 2 POW NEG 1 ADD 1 0.25 2 POW SUB DIV = >> /tmp/body
gmt math -T-4/4/0.1 T 4 DIV 2 POW NEG 1 ADD 1 0.25 2 POW SUB DIV = /tmp/line
gmt plot -R-5/5/0/1.5 -JX8.4c/1.6c -Glightgray /tmp/body
gmt plot -W2p /tmp/tmp
gmt plot -W0.5p,- /tmp/line
echo -0.4 1.10 | gmt plot -Ss8p -Gwhite
gmt plot -Sv0.1i+e+s -Gblack -W0.5p -N <<- EOF
-5 0 5 0
0 0 0 1.5
EOF
gmt plot -W0.25p,- <<- EOF
>
-0.5 0.2
4.3 0.2
>
1 0
1 1
>
3.60555 0
3.60555 0.3
EOF
gmt text -F+f9p,Times-Italic+j -N <<- EOF
1 -0.05 TC r@-f@- = fr@-0@-
4 -0.05 TC r@-0@-
3.60555 -0.05 TC r@-n@-
4.3 0.2 LM h@-n@-
-0.2 1 RM h@-0@-
EOF
gmt end show