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
28
29
30
31
32
33
34
35
36
37
38
39
40
|
REM
REM GMT EXAMPLE 27
REM $Id: example_27.bat 17012 2016-08-26 21:00:21Z pwessel $
REM
REM Purpose: Illustrates how to plot Mercator img grids
REM GMT progs: makecpt, grdgradient, grdimage, pscoast
REM GMT supplement: img2grd (to read Sandwell/Smith img files)
REM DOS calls: del, grep
REM
echo GMT EXAMPLE 27
set ps=example_27.ps
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 First get gradients.
gmt grdgradient tasman_grav.nc -Nt1 -A45 -Gtasman_grav_i.nc
REM Make a suitable cpt file for mGal
gmt makecpt -T-120/120 -Crainbow > grav.cpt
REM Since this is a Mercator grid we use a linear gmt projection
gmt grdimage tasman_grav.nc=ns/0.1 -Itasman_grav_i.nc -Jx0.25i -Cgrav.cpt -P -K > %ps%
REM Then use gmt pscoast to plot land; get original -R from grid info
REM and use Mercator gmt projection with same scale as above on a spherical Earth
gmt pscoast -R145/170/-50.0163575733/-24.9698584055 -Jm0.25i -Ba10f5 -BWSne -O -K -Gblack --PROJ_ELLIPSOID=Sphere -Cwhite -Dh+ --FORMAT_GEO_MAP=dddF >> %ps%
REM Put a color legend in top-left corner of the land mask
gmt psscale -DjTL+o1c+w2i/0.15i -R -J -Cgrav.cpt -Bx50f10 -By+lmGal -F+gwhite+p1p -I -O >> %ps%
REM Clean up
del grav.cpt
del *_i.nc
del .gmt*
|
