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

	
            REM		GMT EXAMPLE 27
REM
REM Purpose:	Illustrates how to plot Mercator img grids
REM GMT modules:	makecpt, grdimage, grdinfo, coast, colorbar
REM
gmt begin ex27
	REM Gravity in tasman_grav.nc is in 0.1 mGal increments and the grid
	REM is already in projected Mercator x/y units.

	REM Make a suitable cpt file for mGal
	gmt makecpt -T-120/120 -Crainbow

	REM Since this is a Mercator grid we use a linear projection
	gmt grdimage @tasman_grav.nc=ns+s0.1 -I+d -Jx0.6c

	REM Then use gmt coast to plot land; get original -R from grid img remark
	REM and use Mercator gmt projection with same scale as above on a spherical Earth
	gmt grdinfo @tasman_grav.nc -Ii > R.txt
	set /p R=<R.txt
	gmt coast %R% -Jm0.6c -B -BWSne -Gblack --PROJ_ELLIPSOID=Sphere -Cwhite -Dh+ --FORMAT_GEO_MAP=dddF

	REM Put a color legend in top-left corner of the land mask
	gmt colorbar -DjTL+o1c+w5c/0.4c -Bx -By+lmGal -I -F+gwhite+p1p
gmt end show