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