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
|
REM GMT EXAMPLE 16
REM
REM $Id$
REM
REM Purpose: Illustrates interpolation methods using same data as Example 12.
REM GMT progs: gmtset, grdview, grdfilter, pscontour, psscale, pstext, surface, triangulate
REM DOS calls: echo, del
REM
REM First make a cpt file as in example 12:
REM
echo GMT EXAMPLE 16
set ps=example_16.ps
REM
REM Illustrate various means of contouring, using triangulate and surface.
REM
gmt set FONT_ANNOT_PRIMARY 9p
REM
gmt pscontour -R0/6.5/-0.2/6.5 -Jx0.45i -P -K -Y5.5i -Ba2f1 -BWSne table_5.11 -Cex16.cpt -I > %ps%
echo 3.25 7 pscontour (triangulate) | gmt pstext -R -J -O -K -N -F+f18p,Times-Roman+jCB >> %ps%
REM
gmt surface table_5.11 -R -I0.2 -Graws0.nc
gmt grdview raws0.nc -R -J -B -Cex16.cpt -Qs -O -K -X3.5i >> %ps%
echo 3.25 7 surface (tension = 0) | gmt pstext -R -J -O -K -N -F+f18p,Times-Roman+jCB >> %ps%
REM
gmt surface table_5.11 -R -I0.2 -Graws5.nc -T0.5
gmt grdview raws5.nc -R -J -B -Cex16.cpt -Qs -O -K -Y-3.75i -X-3.5i >> %ps%
echo 3.25 7 surface (tension = 0.5) | gmt pstext -R -J -O -K -N -F+f18p,Times-Roman+jCB >> %ps%
REM
gmt triangulate table_5.11 -Grawt.nc -R -I0.2
gmt grdfilter rawt.nc -Gfiltered.nc -D0 -Fc1
gmt grdview filtered.nc -R -J -B -Cex16.cpt -Qs -O -K -X3.5i >> %ps%
echo 3.25 7 triangulate @~\256@~ gmt grdfilter | gmt pstext -R -J -O -K -N -F+f18p,Times-Roman+jCB >> %ps%
echo 3.2125 7.5 Gridding of Data | gmt pstext -R0/10/0/10 -Jx1i -O -K -N -F+f32p,Times-Roman+jCB -X-3.5i >> %ps%
gmt psscale -Dx3.25i/-0.4i+jTC+w5i/0.25i+h -Cex16.cpt -O >> %ps%
del *.nc
del .gmt*
del gmt.conf
|
