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
#
# Paste overlapping pixel registered grids along X & Y

ps=paste_pixreg_overlap.ps

# The final grid is just f(x,y) = x*y
gmt grdmath -R-15.25/15.25/-15.25/15.25 -I0.5 X Y MUL -r = lixo.nc
# The 4 pieces to assemble into final grid
gmt grdmath -R-15.25/0.25/-15.25/0.25 -I0.5 X Y MUL -r = lixo_x1.nc
gmt grdmath -R0.25/15.25/-15.25/0.25  -I0.5 X Y MUL -r = lixo_x2.nc
gmt grdmath -R-15.25/0.25/0.25/15.25  -I0.5 X Y MUL -r = lixo_y1.nc
gmt grdmath -R0.25/15.25/0.25/15.25   -I0.5 X Y MUL -r = lixo_y2.nc
gmt grdpaste lixo_x1.nc lixo_x2.nc -Glixo_x.nc
gmt grdpaste lixo_y1.nc lixo_y2.nc -Glixo_y.nc
gmt grdpaste lixo_x.nc lixo_y.nc -Glixo_xy.nc
# Top is single source, bottom is assembled
gmt grdcontour lixo_xy.nc -JX10c -C10 -B5 -P -K -Xc > $ps
gmt grdcontour lixo.nc -J -C10 -B5 -O -Y12c >> $ps