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

transfer.sh 2.9 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
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
  1. #!/usr/bin/env bash
  2. # Testing firmoviscous flexure code by dumping all the various transfer functions
  3. gmt set GMT_FFT kiss FONT_TAG 7p FONT_ANNOT_PRIMARY 10p
  4. gmt begin transfer ps
  5. gmt subplot begin 6x3 -Fs5.5c/3c -Srl -Scb -JX5.5cl/3.25c -R1/5000/0/1.1 -Bxaf+ukm -Byafg1 -Bwesn -M3p -Y2c -T"Transfer Functions" -A
  6. # Elastic
  7. gmt grdflexure -D3300/2700/2700/1000 -Q -E
  8. gmt subplot set 0 -A"T@-e@- = 1 km"
  9. gmt plot grdflexure_transfer_function_te_001_km.txt -i0,2 -W1p
  10. gmt subplot set 1 -A"T@-e@- = 10 km"
  11. gmt plot grdflexure_transfer_function_te_010_km.txt -i0,2 -W1p
  12. gmt subplot set 2 -A"T@-e@- = 100 km"
  13. gmt plot grdflexure_transfer_function_te_100_km.txt -i0,2 -W1p
  14. # Viscoelastic
  15. gmt grdflexure -D3300/2700/2700/1000 -Q -E10 -M20k
  16. gmt subplot set 3 -A"M@-t@- = 20 ky E = 10 km t = 5 ky"
  17. gmt plot grdflexure_transfer_function_te_010_km.txt -i0,4 -W1p
  18. gmt subplot set 4 -A"M@-t@- = 20 ky E = 10 km t = 100 ky"
  19. gmt plot grdflexure_transfer_function_te_010_km.txt -i0,8 -W1p
  20. gmt subplot set 5 -A"M@-t@- = 20 ky E = 10 km t = 1 My"
  21. gmt plot grdflexure_transfer_function_te_010_km.txt -i0,11 -W1p
  22. # Firmoviscous 1
  23. gmt grdflexure -D3300/2700/2700/1000 -Q -E10 -F2e20
  24. gmt subplot set 6 -A"@~h@~ = 2e20 E = 5 km t = 5 ky"
  25. gmt plot grdflexure_transfer_function_te_005_km.txt -i0,4 -W1p
  26. gmt subplot set 7 -A"@~h@~ = 2e20 E = 10 km t = 100 ky"
  27. gmt plot grdflexure_transfer_function_te_010_km.txt -i0,6 -W1p
  28. gmt subplot set 8 -A"@~h@~ = 2e20 E = 20 km t = 1 My"
  29. gmt plot grdflexure_transfer_function_te_020_km.txt -i0,8 -W1p
  30. # Firmoviscous 2
  31. gmt grdflexure -D3300/2700/2700/1000 -Q -E10 -F2e20/100k/1e21
  32. gmt subplot set 9 -A"@~h@~ = 2e20/1e21 E = 5 km A = 100 km t = 5 ky"
  33. gmt plot grdflexure_transfer_function_te_005_km.txt -i0,4 -W1p
  34. gmt subplot set 10 -A"@~h@~ = 2e20/1e21 E = 10 km A = 100 km t = 100 ky"
  35. gmt plot grdflexure_transfer_function_te_010_km.txt -i0,6 -W1p
  36. gmt subplot set 11 -A"@~h@~ = 2e20/1e21 E = 20 km A = 100 km t = 1 My"
  37. gmt plot grdflexure_transfer_function_te_020_km.txt -i0,8 -W1p
  38. # Viscous 1
  39. gmt grdflexure -D3300/2700/2700/1000 -Q -E0 -F2e20
  40. gmt subplot set 12 -A"@~h@~ = 2e20 E = 0 t = 5 ky"
  41. gmt plot grdflexure_transfer_function_te_000_km.txt -i0,4 -W1p
  42. gmt subplot set 13 -A"@~h@~ = 2e20 E = 0 t = 100 ky"
  43. gmt plot grdflexure_transfer_function_te_000_km.txt -i0,6 -W1p
  44. gmt subplot set 14 -A"@~h@~ = 2e20 E = 0 t = 1 My"
  45. gmt plot grdflexure_transfer_function_te_000_km.txt -i0,8 -W1p
  46. # Viscous 2
  47. gmt grdflexure -D3300/2700/2700/1000 -Q -E0 -F2e20/100k/1e21
  48. gmt subplot set 15 -A"@~h@~ = 2e20/1e21 A = 100 km E = 0 t = 5 ky"
  49. gmt plot grdflexure_transfer_function_te_000_km.txt -i0,4 -W1p
  50. gmt subplot set 16 -A"@~h@~ = 2e20/1e21 A = 100 km E = 0 t = 100 ky"
  51. gmt plot grdflexure_transfer_function_te_000_km.txt -i0,6 -W1p
  52. gmt subplot set 17 -A"@~h@~ = 2e20/1e21 A = 100 km E = 0 t = 1 My"
  53. gmt plot grdflexure_transfer_function_te_000_km.txt -i0,8 -W1p
  54. gmt subplot end
  55. gmt end show
Tip!

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

Comments

Loading...