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
|
#!/usr/bin/env bash
# Show the seamount density model offered in grdseamount
gmt begin GMT_seamount_density
echo 0 0 50 6000 | gmt grdseamount -H6000/2500/3000+p1+d0 -Kmodel.grd -Cc -F0.2
gmt makecpt -Cbilbao -T2500/3000 --COLOR_NAN=white
# Plot density reference model above
gmt grdimage model.grd -R0/1.1/0/1.3 -JX15c/4c -Bxaf+l"Normalized radial distance, @%6%r@%%" -Byafg1+l"Normalized height, @%6%h(r)@%%"
printf "0 1.25\n1 0\n" | gmt plot -W0.25p,-
printf "0 1\n 0.2 1\n1 0\n" | gmt plot -W1p
z=0.5
y=$(gmt math -Q 1.25 $z ADD 2 DIV =)
echo 0.6 0.15 | gmt plot -Sc4p -Gwhite -W0.25p
gmt plot -Sv12p+b+e+s -W1p -Gblack -N <<- EOF
0.6 0.15 0.6 ${z}
0.6 1.3 0.6 ${z}
1.1 0 1.1 1.3
EOF
printf "0.6 %s\n1 %s\n" $z $z | gmt plot -W0.25p,-
printf "0.6 %s\n0.95 %s\n" $y $y | gmt plot -W0.25p
printf "0.6 0.3\n0.95 0.3\n" | gmt plot -W0.25p
printf "0.6 0\n0.6 0.15\n0.95 0.15\n" | gmt plot -W0.25p
echo 0.0 1.25 "@~r@~@-l@-" | gmt text -F+f10p+jTL -Dj2p/4p -N
echo 0 0 "@~r@~@-h@-" | gmt text -F+f10p,white+jBL -Dj2p/4p -N
gmt text -F+f12p+jLM -Dj3p <<- EOF
0.95 $y @[h_r - h(r)@[
0.95 0.3 @[h(r) - z@[
0.95 0.15 @[(r, z)@[
1.05 0.65 @[h_r@[
EOF
gmt text -F+f9p+j -Dj3p <<- EOF
0.75 $z BL Outside seamount
0.75 $z TL Inside seamount
EOF
gmt colorbar -DJBC+w80%/0.25c -Bx -By+l"kg m@+-3@+"
rm -f model.grd
gmt end show
|
