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
|
#!/bin/bash
#
# $Id$
# Makes the data for GMT_App_O_[1-9].sh
#
gmt grdcut ../examples/ex01/osu91a1f_16.nc -R50/160/-15/15 -Ggeoid.nc
# Complex line algorithm points
cat << EOF > cross.d
> 1st branch
59 -12
62 -7
66 -3
71 -1
77 3
> 2nd branch
94 -11
100 -10.5
105 -9.5
109 -8.6
114 -6.5
119 -4
126 2
> 3rd branch
148 3
158 13
EOF
info=`gmt grdinfo -M -C geoid.nc`
x0=`echo $info | cut -f12 -d ' '`
y0=`echo $info | cut -f13 -d ' '`
x1=`echo $info | cut -f14 -d ' '`
y1=`echo $info | cut -f15 -d ' '`
gmt project -C$x0/$y0 -E$x1/$y1 -G10 -Q > tt.d
dist=`gmt gmtconvert tt.d --FORMAT_FLOAT_OUT=%.0lf -El -o2`
R=`gmt info -I1 tt.d`
echo "# Geoid Extrema Separation is $dist km" > transect.d
gmt grdtrack tt.d -Ggeoid.nc | gmt grdtrack -GGMT_App_O.nc >> transect.d
rm -f tt.d gmt.history
|
