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

	
            #!/bin/bash
#
#	$Id$

ps=imgmap.ps
IMG=topo.15.1.img
if [ -f /tmp/gmt/$IMG ]; then # use one placed/linked from /tmp
	IMG=/tmp/gmt/topo.15.1.img
else
	OK=`gmt gmtwhich -C $IMG`
	if [ $OK = N ]; then
		echo "File $IMG not available for testing imgmap.sh"
		exit
	fi
fi
# Get merc grid
gmt img2grd $IMG -R180/200/-5/5 -I1m -T1 -D -S1 -Gimg.nc -M
gmt makecpt -Crainbow -T-8000/0/500 -Z > t.cpt
gmt grdimage img.nc -Jx0.25i -Ct.cpt -P -K -Xc > $ps
gmt psbasemap -R -Jm0.25i -Ba -BWSne -O -K >> $ps
# Get geo grid
gmt img2grd $IMG -R -I1m -T1 -D -S1 -Gimg.nc
gmt grdimage img.nc -Jm -Ct.cpt -O -K -Ba -BWSne -Y3.25i >> $ps
# Get resampled geo grid
gmt img2grd $IMG -R -I1m -T1 -D -S1 -Gimg.nc -E
gmt grdimage img.nc -Jm -Ct.cpt -O -K -Ba -BWSne -Y3.25i >> $ps
gmt psxy -R -J -O -T >> $ps