maxrjones
/
gmt-fork
mirror of https://github.com/meghanrjones/gmt.git


  
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

	
            #!/usr/bin/env bash
# Test grdlandmask across Greenwich

ps=greenwich.ps
gmt makecpt -Cred -T0.5,1.5 -N > t.cpt
# ORIGINAL:
gmt grdlandmask -R-1.05/0.05/50.95/52.05 -I3s -Glandmask.grd=ns -Df -NNaN/1/1/1/1
gmt grdimage landmask.grd -Jx2.5id -P -Ba1g1 -BWSne -Ct.cpt -K -Xc > $ps
gmt pscoast -R -J -O -K -Df -W0.25p >> $ps
gmt grdlandmask -R-2.05/-0.95/50.95/52.05 -I3s -Glandmask.grd=ns -Df -NNaN/1/1/1/1
gmt grdimage landmask.grd -J -O -Ba1g1 -BWSne -Ct.cpt -K -Y3.25i >> $ps
gmt pscoast -R -J -O -K -Df -W0.25p >> $ps
gmt grdlandmask -R-0.05/1.05/50.95/52.05 -I3s -Glandmask.grd=ns -Df -NNaN/1/1/1/1
gmt grdimage landmask.grd -J -O -K -Ba1g1 -BWSne -Ct.cpt -Y3.25i >> $ps
gmt pscoast -R -J -O -Df -W0.25p >> $ps
rm -f landmask.grd t.cpt

# FAKED TO GIVE CORRECT PS
#grdlandmask -R-1.05/0.05/50.95/52.05 -I3s -Glandmask.grd=ns -Df -NNaN/1/1/1/1
#psbasemap -R -Jx2.5id -P -B+gred -K -Xc > $ps
#grdimage landmask.grd -Jx2.5id -O -Ct.cpt -K >> $ps
#psxy -R -J -O -K -A -Gred << EOF >> $ps
#-1.05	50.95
#-0.95	50.95
#-0.95	52.05
#-1.05	52.05
#EOF
#pscoast -R -J -O -K -Df -W0.25p -Ba1g1 -BWSne >> $ps
#grdlandmask -R-2.05/-0.95/50.95/52.05 -I3s -Glandmask.grd=ns -Df -NNaN/1/1/1/1
#grdimage landmask.grd -J -O -Ba1g1 -BWSne -Ct.cpt -K -Y3.25i >> $ps
#pscoast -R -J -O -K -Df -W0.25p >> $ps
#grdlandmask -R-0.05/1.05/50.95/52.05 -I3s -Glandmask.grd=ns -Df -NNaN/1/1/1/1
#grdimage landmask.grd -J -O -K -Ct.cpt -Y3.25i >> $ps
#psxy -R -J -O -K -A -Gred << EOF >> $ps
#-0.05	50.95
#0.05	50.95
#0.05	52.05
#-0.05	52.05
#EOF
#pscoast -R -J -O -Df -W0.25p -Ba1g1 -BWSne >> $ps