-
Notifications
You must be signed in to change notification settings - Fork 348
/
layout.sh
executable file
·25 lines (25 loc) · 1.25 KB
/
layout.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
#!/usr/bin/env bash
# Test script to ensure forum message 3511 has been dealt with
ps=layout.ps
function plot_it {
gmt grdtrack @Matthews_2016_subduction_subset.txt -Gdummy.nc -C500k/100/100${1} > sz_pol_left.gmt
awk '{ if ( $1 == ">") print $0 ; else if ($3 <= 0) print $0 }' sz_pol_left.gmt > Lhalfxprofile.gmt
awk '{ if ( $1 == ">") print $0 ; else if ($3 >= 0) print $0 }' sz_pol_left.gmt > Rhalfxprofile.gmt
gmt psxy -R -J -O -K -W1p @Matthews_2016_subduction_subset.txt
gmt psxy -R -J -O Lhalfxprofile.gmt -W1p,red -K
gmt psxy -R -J Rhalfxprofile.gmt -W1p,blue -O -K
gmt psxy -R -J Rhalfxprofile.gmt -Sc0.4c -Gblue -O -K
gmt psxy -R -J Lhalfxprofile.gmt -Sc0.4c -Gred -O -K
gmt pstext -R -J -O -K sz_pol_left.gmt -F+f9p,Helvetica,white+r0
gmt psxy -R -J -O -K -Ss0.25i -Gblack @Matthews_2016_subduction_subset.txt
gmt pstext -R -J -O -K @Matthews_2016_subduction_subset.txt -F+f12p,Helvetica,white+r0
echo $1 | gmt pstext -R -J -O -K -F+cTL+jTL+f18p -Dj0.2i -Gwhite
}
gmt grdmath -R135/162/42/51 -I1 X = dummy.nc
gmt psbasemap -R135/162/42/51 -JM6i -Bafg -BWSnE -P -K -Xc > $ps
plot_it +v >> $ps
gmt psbasemap -R -J -O -Bafg -BWsnE -K -Y3.1i >> $ps
plot_it +a >> $ps
gmt psbasemap -R -J -O -Bafg -BWsnE -K -Y3.1i >> $ps
plot_it >> $ps
gmt psxy -R -J -O -T >> $ps