1#!/usr/bin/env bash
2# Test grdlandmask across Greenwich
3
4ps=greenwich.ps
5gmt makecpt -Cred -T0.5,1.5 -N > t.cpt
6# ORIGINAL:
7gmt grdlandmask -R-1.05/0.05/50.95/52.05 -I3s -Glandmask.grd=ns -Df -NNaN/1/1/1/1
8gmt grdimage landmask.grd -Jx2.5id -P -Ba1g1 -BWSne -Ct.cpt -K -Xc > $ps
9gmt pscoast -R -J -O -K -Df -W0.25p >> $ps
10gmt grdlandmask -R-2.05/-0.95/50.95/52.05 -I3s -Glandmask.grd=ns -Df -NNaN/1/1/1/1
11gmt grdimage landmask.grd -J -O -Ba1g1 -BWSne -Ct.cpt -K -Y3.25i >> $ps
12gmt pscoast -R -J -O -K -Df -W0.25p >> $ps
13gmt grdlandmask -R-0.05/1.05/50.95/52.05 -I3s -Glandmask.grd=ns -Df -NNaN/1/1/1/1
14gmt grdimage landmask.grd -J -O -K -Ba1g1 -BWSne -Ct.cpt -Y3.25i >> $ps
15gmt pscoast -R -J -O -Df -W0.25p >> $ps
16rm -f landmask.grd t.cpt
17
18# FAKED TO GIVE CORRECT PS
19#grdlandmask -R-1.05/0.05/50.95/52.05 -I3s -Glandmask.grd=ns -Df -NNaN/1/1/1/1
20#psbasemap -R -Jx2.5id -P -B+gred -K -Xc > $ps
21#grdimage landmask.grd -Jx2.5id -O -Ct.cpt -K >> $ps
22#psxy -R -J -O -K -A -Gred << EOF >> $ps
23#-1.05	50.95
24#-0.95	50.95
25#-0.95	52.05
26#-1.05	52.05
27#EOF
28#pscoast -R -J -O -K -Df -W0.25p -Ba1g1 -BWSne >> $ps
29#grdlandmask -R-2.05/-0.95/50.95/52.05 -I3s -Glandmask.grd=ns -Df -NNaN/1/1/1/1
30#grdimage landmask.grd -J -O -Ba1g1 -BWSne -Ct.cpt -K -Y3.25i >> $ps
31#pscoast -R -J -O -K -Df -W0.25p >> $ps
32#grdlandmask -R-0.05/1.05/50.95/52.05 -I3s -Glandmask.grd=ns -Df -NNaN/1/1/1/1
33#grdimage landmask.grd -J -O -K -Ct.cpt -Y3.25i >> $ps
34#psxy -R -J -O -K -A -Gred << EOF >> $ps
35#-0.05	50.95
36#0.05	50.95
37#0.05	52.05
38#-0.05	52.05
39#EOF
40#pscoast -R -J -O -Df -W0.25p -Ba1g1 -BWSne >> $ps
41
42