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

	
            #!/usr/bin/env bash
#
# Check if we can wrap global grids over longitude
#

ps=lonshift.ps

# Make a data set (z values)
$AWK 'BEGIN{n=12;m=0;for (j=0;j<6*n;j++) {if (j%n==0) m--;print sin(m*3.14159265/6);m++}}' > tmp.txt
# Build gridline-registered grid
gmt xyz2grd -I30 -Gtmp.nc -ZTLa -fg -R15/345/-75/75 tmp.txt
gmt makecpt -Crainbow -T-1/1/0.1 > tmp.cpt
# Plot orig grid and force actual region -R15/345/-75/75 despite a global grid
gmt grdimage -Ctmp.cpt tmp.nc -JX6i/3i -Bx60f10 -By30f10 -BWeSn+t"Plotted as 15/345/-75/75" --MAP_FRAME_TYPE=plain --FONT_ANNOT_PRIMARY=10p --FORMAT_GEO_MAP=DF -R15/345/-75/75 -K -Y6.5i -Xc -P > $ps
# plot location of grid nodes
gmt grd2xyz tmp.nc | gmt psxy -R15/345/-75/75 -JX6i/3i -O -K -Sc0.1i -Gblack -N >> $ps
# Plot orig grid but shifted by giving -R0/360/-75/75
gmt grdimage -Ctmp.cpt tmp.nc -J -Bx60f10 -By30f10 -BWeSn+t"Plotted as 0/360/-75/75" --MAP_FRAME_TYPE=plain --FONT_ANNOT_PRIMARY=10p --FORMAT_GEO_MAP=DF -R0/360/-75/75 -Y-5i -O -K >> $ps
# plot location of grid nodes
gmt grd2xyz tmp.nc | gmt psxy -R0/360/-75/75 -JX6i/3i -O -Sc0.1i -Gblack -N >> $ps