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
|
#!/bin/bash
# $Id: GMT_atan.sh 11490 2013-05-16 06:26:21Z pwessel $
#
gmt grdgradient -A45 ../tutorial/us.nc -N -fg -Gtt.t.nc
gmt grd2xyz -Z tt.t.nc > tt.d
gmt pshistogram tt.d -R-0.75/0.75/0/20 -JX1.5/1 -Bx0.5 -By5f5 -BWSne -W0.01 -P -K -Gblack -Z1 > GMT_atan.ps
gmt pstext -R -J -O -K -F+f9p+jLB << EOF >> GMT_atan.ps
-0.7 17 Raw
-0.7 15 slopes
EOF
gmt gmtmath -T-5/5/0.01 T ATAN PI DIV 2 MUL = > tt.d
gmt psxy tt.d -R-5/5/-1/1 -JX1.5/1 -Bx2f1g1 -By1f0.5g0.25 -BWSne -O -K -Wthick -X1.85 >> GMT_atan.ps
gmt psxy -R -J -O -K -Sv5p+e -Gblack -W0.5p << EOF >> GMT_atan.ps
3 0.8 180 0.45
3 0.8 -90 0.4
EOF
gmt psxy -R -J -O -K -Wthinnest << EOF >> GMT_atan.ps
>
-5 0
5 0
>
0 -1
0 1
EOF
gmt grdgradient -A45 ../tutorial/us.nc -Nt -fg -Gtt.tt.nc
gmt grd2xyz -Z tt.tt.nc > tt.d
gmt pshistogram tt.d -R-0.75/0.75/0/5 -JX1.5/1 -Bx0.5 -By2f1 -BWSne -W0.01 -O -K -Gblack -X1.85 -Z1 >> GMT_atan.ps
gmt pstext -R -J -O -F+f9p+jLB << EOF >> GMT_atan.ps
-0.7 4.3 tan@+-1@+
-0.7 3.7 transformed
EOF
|
