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
  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. echo "100 75 50 5000 1 0" > t.txt
  6. 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
  7. 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
  8. gmt subplot begin 5x2 -SCb+tc -SRl -A -R40/160/0/5.010 -Fs3i/0.5i -M8p
  9. gmt set FONT_TAG 16p,Times-Italic,black
  10. gmt subplot set 0,0 -A"t = 4"
  11. gmt grdtrack -Gsmtc_00.00.nc -E40000/75000/160000/75000 -o0,2 | gmt plot -W1p -i0+s0.001,1+s0.001 -L+y0 -Gblack
  12. gmt subplot set 1,0 -A"t = 3"
  13. gmt grdtrack -Gsmtc_00.20.nc -E40000/75000/160000/75000 -o0,2 | gmt plot -W1p -i0+s0.001,1+s0.001 -L+y0 -Gblack
  14. gmt subplot set 2,0 -A"t = 2"
  15. gmt grdtrack -Gsmtc_00.40.nc -E40000/75000/160000/75000 -o0,2 | gmt plot -W1p -i0+s0.001,1+s0.001 -L+y0 -Gblack
  16. gmt subplot set 3,0 -A"t = 1"
  17. gmt grdtrack -Gsmtc_00.60.nc -E40000/75000/160000/75000 -o0,2 | gmt plot -W1p -i0+s0.001,1+s0.001 -L+y0 -Gblack
  18. gmt subplot set 4,0 -A"t = 1"
  19. gmt grdtrack -Gsmtc_00.80.nc -E40000/75000/160000/75000 -o0,2 | gmt plot -W1p -i0+s0.001,1+s0.001 -L+y0 -Gblack
  20. #
  21. gmt subplot set 0,1 -A"t = 4"
  22. gmt grdtrack -Gsmti_00.00.nc -E40000/75000/160000/75000 -o0,2 | gmt plot -W1p -i0+s0.001,1+s0.001 -L+y0 -Gblack
  23. gmt subplot set 1,1 -A"t = 3"
  24. gmt grdtrack -Gsmti_00.20.nc -E40000/75000/160000/75000 -o0,2 | gmt plot -W1p -i0+s0.001,1+s0.001 -L+y0 -Gblack
  25. gmt subplot set 2,1 -A"t = 2"
  26. gmt grdtrack -Gsmti_00.40.nc -E40000/75000/160000/75000 -o0,2 | gmt plot -W1p -i0+s0.001,1+s0.001 -L+y0 -Gblack
  27. gmt subplot set 3,1 -A"t = 1"
  28. gmt grdtrack -Gsmti_00.60.nc -E40000/75000/160000/75000 -o0,2 | gmt plot -W1p -i0+s0.001,1+s0.001 -L+y0 -Gblack
  29. gmt subplot set 4,1 -A"t = 0"
  30. gmt grdtrack -Gsmti_00.80.nc -E40000/75000/160000/75000 -o0,2 | gmt plot -W1p -i0+s0.001,1+s0.001 -L+y0 -Gblack
  31. gmt subplot end
  32. gmt end show
  33. 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...