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
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
|
REM GMT EXAMPLE 49
REM $Id$
REM
REM Purpose: Illustrate data analysis using the seafloor depth/age relationship
REM GMT progs: blockmode, gmtmath, grdcontour, grdimage, grdsample, makecpt,
REM psbasemap, pslegend, psscale, psxy, xyz2grd
echo GMT EXAMPLE 49
set ps=example_49.ps
REM Convert coarser age grid to pixel registration to match bathymetry grid
gmt grdsample age_gridline.nc -T -Gage_pixel.nc
REM Image depths with color-coded age contours
gmt makecpt -Cabyss -T-7000/0 > depth.cpt
gmt grdimage depth_pixel.nc -Cdepth.cpt -JM6i -P -Baf -BWSne -X1.5i -K --FORMAT_GEO_MAP=dddF > %ps%
gmt psxy -Rdepth_pixel.nc -J -O -K -W1p ridge.gmt >> %ps%
gmt makecpt -Chot -T0/100/10 > age.cpt
gmt grdcontour age_pixel.nc -J -O -K -A+f14p -Cage.cpt -Wa0.1p+c -GL30W/22S/5E/13S >> %ps%
gmt psscale -Rdepth_pixel.nc -J -DjTR+w2i/0.15i+h+o0.3i/0.15i -Cdepth.cpt -Baf+u" km" -W0.001 -F+p1p+gbeige -O -K >> %ps%
REM Obtain depth, age pairs by dumping grids and pasting results
gmt grd2xyz age_pixel.nc -bof > age.bin
gmt grd2xyz depth_pixel.nc -bof > depth.bin
gmt convert -A age.bin depth.bin -bi3f -o2,5,5 -bo3f > depth-age.bin
REM Create and map density grid of (age,depth) distribution
gmt xyz2grd -R0/100/-6500/0 -I0.25/25 -r depth-age.bin -bi3f -An -Gdensity.nc
gmt makecpt -Chot -T0/100 > density.cpt
gmt grdimage density.nc -JX6i/4i -Q -O -K -Cdensity.cpt -Y4.8i >> %ps%
REM Obtain modal depths every ~5 Myr
gmt blockmode -R0/100/-10000/0 -I5/10000 -r -E depth-age.bin -bi3f -o0,2,3 > modal.txt
REM Compute Parsons & Sclater [1977] depth-age curve
REM depth(t) = 350 * sqrt(t) + 2500, t < 70 Myr
REM = 6400 - 3200 exp (-t/62.8), t > 70 Myr
gmt math -T0/100/0.1 T SQRT 350 MUL 2500 ADD T 70 LE MUL 6400 T 62.8 DIV NEG EXP 3200 MUL SUB T 70 GT MUL ADD NEG = ps.txt
gmt psxy -Rdensity.nc -J -O -K ps.txt -W4p,green >> %ps%
gmt psxy -R -J -O -K ps.txt -W1p >> %ps%
REM Compute Stein & Stein [1992] depth-age curve
REM depth(t) = 365 * sqrt(t) + 2600, t < 20 Myr
REM = 5651 - 2473 * exp (-0.0278*t), t > 20 Myr
gmt math -T0/100/0.1 T SQRT 365 MUL 2600 ADD T 20 LE MUL 5651 T -0.0278 MUL EXP 2473 MUL SUB T 20 GT MUL ADD NEG = ss.txt
REM Plot curves and place the legend
gmt psxy -R -J -O -K ss.txt -W4p,white >> %ps%
gmt psxy -R -J -O -K ss.txt -W1p >> %ps%
gmt psxy -R -J -Ss0.4c -O -K -Gblue modal.txt -Ey/1p,blue >> %ps%
gmt psxy -R -J -Ss0.1c -O -K -Gwhite modal.txt >> %ps%
gmt psbasemap -R0/100/0/6.5 -JX6i/-4i -Bxaf+u" Myr" -Byaf+u" km" -BWsNe -O -K >> %ps%
echo S 0.2i - 0.35i - 4p,green 0.5i Parsons & Sclater (1977) > L
echo S 0.2i - 0.35i - 4p,white 0.5i Stein & Stein (1992) >> L
echo S 0.2i s 0.15i blue - 0.5i Modal depth estimates >> L
gmt pslegend -R -J -O -K -DjRT+w2.5i+o0.1i -F+p1p+gbeige+s L >> %ps%
echo S 0.2i - 0.35i - 1p 0.3i > L
echo S 0.2i - 0.35i - 1p 0.3i >> L
echo S 0.2i s 0.1c white - 0.3i >> L
gmt pslegend -R -J -O -K -DjRT+w2.5i+o0.1i L >> %ps%
EOF
gmt psxy -R -J -O -T >> %ps%
del age_pixel.nc
del *.bin
del *.txt
del *.cpt
del density.nc
|
