Register
Login
Resources
Docs Blog Datasets Glossary Case Studies Tutorials & Webinars
Product
Data Engine LLMs Platform Enterprise
Pricing Explore
Connect to our Discord channel

GMT_seamount_cum_inc.sh 2.0 KB

You have to be logged in to leave a comment. Sign In
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
  1. #!/usr/bin/env bash
  2. # Plot 5 life-stages in growing volcanoes in crossection
  3. # Using un-truncated Gaussian shapes with linear flux
  4. gmt begin GMT_seamount_cum_inc ps
  5. gmt set GMT_THEME cookbook
  6. echo "100 75 50 5000 1 0" > t.txt
  7. gmt grdseamount -R40/160/74/76+uk -I100 -Gsmtc_%05.2f.nc t.txt -T0.8/0/0.2 -Qc/l -Dk -Cg -Mc.lis
  8. gmt grdseamount -R40/160/74/76+uk -I100 -Gsmti_%05.2f.nc t.txt -T0.8/0/0.2 -Qi/l -Dk -Cg -Mi.lis
  9. gmt subplot begin 5x2 -Scb+tc -Srl -A -R40/160/0/5.010 -Fs3i/0.5i -M8p
  10. gmt set FONT_TAG 16p,Times-Italic,black
  11. gmt subplot set 0,0 -A"t = 4"
  12. gmt grdtrack -Gsmtc_00.00.nc -E40000/75000/160000/75000 -o0,2 | gmt plot -W1p -i0+s0.001,1+s0.001 -L+y0 -Gblack
  13. gmt subplot set 1,0 -A"t = 3"
  14. gmt grdtrack -Gsmtc_00.20.nc -E40000/75000/160000/75000 -o0,2 | gmt plot -W1p -i0+s0.001,1+s0.001 -L+y0 -Gblack
  15. gmt subplot set 2,0 -A"t = 2"
  16. gmt grdtrack -Gsmtc_00.40.nc -E40000/75000/160000/75000 -o0,2 | gmt plot -W1p -i0+s0.001,1+s0.001 -L+y0 -Gblack
  17. gmt subplot set 3,0 -A"t = 1"
  18. gmt grdtrack -Gsmtc_00.60.nc -E40000/75000/160000/75000 -o0,2 | gmt plot -W1p -i0+s0.001,1+s0.001 -L+y0 -Gblack
  19. gmt subplot set 4,0 -A"t = 1"
  20. gmt grdtrack -Gsmtc_00.80.nc -E40000/75000/160000/75000 -o0,2 | gmt plot -W1p -i0+s0.001,1+s0.001 -L+y0 -Gblack
  21. #
  22. gmt subplot set 0,1 -A"t = 4"
  23. gmt grdtrack -Gsmti_00.00.nc -E40000/75000/160000/75000 -o0,2 | gmt plot -W1p -i0+s0.001,1+s0.001 -L+y0 -Gblack
  24. gmt subplot set 1,1 -A"t = 3"
  25. gmt grdtrack -Gsmti_00.20.nc -E40000/75000/160000/75000 -o0,2 | gmt plot -W1p -i0+s0.001,1+s0.001 -L+y0 -Gblack
  26. gmt subplot set 2,1 -A"t = 2"
  27. gmt grdtrack -Gsmti_00.40.nc -E40000/75000/160000/75000 -o0,2 | gmt plot -W1p -i0+s0.001,1+s0.001 -L+y0 -Gblack
  28. gmt subplot set 3,1 -A"t = 1"
  29. gmt grdtrack -Gsmti_00.60.nc -E40000/75000/160000/75000 -o0,2 | gmt plot -W1p -i0+s0.001,1+s0.001 -L+y0 -Gblack
  30. gmt subplot set 4,1 -A"t = 0"
  31. gmt grdtrack -Gsmti_00.80.nc -E40000/75000/160000/75000 -o0,2 | gmt plot -W1p -i0+s0.001,1+s0.001 -L+y0 -Gblack
  32. gmt subplot end
  33. gmt end show
  34. rm -f smt[ci]*.nc ?.lis t.txt
Tip!

Press p or to see the previous file or, n or to see the next file

Comments

Loading...