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

maxrjones
/
gmt-fork
mirror of https://github.com/meghanrjones/gmt.git

Watch 3
Issues   Discussions 0

fields.sh 4.3 KB
Permalink History Raw

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
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93 
  1. #!/usr/bin/env bash
    2. 

  2. #
    3. 

  3. # Compute FAA, VGG, and geoid over synthetic seamount
    4. 

  4. ps=fields.ps
    5. 

  5. order=2
    6. 

  6. dx=1
    7. 

  7. gmt set GMT_FFT kiss
    8. 

  8. function Ugly {
    9. 

  9. 	# sub-function for integral I (z, b, c)
    10. 

  10. 	z=$1
    11. 

  11. 	b=$2
    12. 

  12. 	c=$3
    13. 

  13. 	beta2=$(gmt math -Q 1.0 $b $b MUL ADD =)
    14. 

  14. 	beta=$(gmt math -Q $beta2 SQRT =)
    15. 

  15. 	beta3=$(gmt math -Q $beta2 $beta MUL =)
    16. 

  16. 	P1=$(gmt math -Q $z $z MUL 2 $b MUL $c MUL $z MUL $c $c MUL ADD $beta2 DIV ADD SQRT $beta DIV =)
    17. 

  17. 	P2=$(gmt math -Q $z $b $c MUL $beta2 ADD DIV $z $z MUL 2 $b MUL $c MUL $z MUL $c $c MUL ADD $beta2 DIV ADD SQRT ADD LOG =)
    18. 

  18. 	gmt math -Q $z $P1 SUB $b $c MUL $P2 MUL $beta3 DIV ADD =
    19. 

  19. #	v = z - sqrt (z * z + (2.0 * b * c * z  + c * c) / beta2) / beta + ...
    20. 

  20. #		b * c * log (z + b * c / beta2 + sqrt (z * z + (2.0 * b * c * z  + c * c) / beta2)) / beta3;
    21. 

  21. #	  = z - P1 + b *c * P2 / beta3
    22. 

  22. }
    23. 

  23. 

  24. # 2 panels of topo and grav, with top profile of admittance & coherence
    25. 

  25. # NOT FINISHED
    26. 

  26. # 1. Create a bathymetry data set with one circular truncated seamount
    27. 

  27. #    as in Fig 3. of Marks & Smith, 2007 [GRL], with R_base = 35 km,
    28. 

  28. # R_top = 10 km, height = 3751 m, depth = -5084 m, density d_rho = 2800-1030
    29. 

  29. # = 1670 kg/m^3, so the flattening is 10/25 = 0.4.
    30. 

  30. echo "0	0	25	3751" | gmt grdseamount -R-256/256/-256/256 -I$dx -r -C -Gsmt.nc -F0.4 -Z-5084
    31. 

  31. # BL Plot the bathymetry
    32. 

  32. gmt makecpt -Crainbow -T-5100/-1000 > t.cpt
    33. 

  33. gmt grdimage smt.nc -R-100/100/-100/100 -JX3i -P -Bag -BWSne -Ct.cpt -K > $ps
    34. 

  34. gmt grdtrack -Gsmt.nc -ELM/RM+d > smt.trk
    35. 

  35. gmt psxy -R -J -O -K -W5p,white smt.trk >> $ps
    36. 

  36. gmt psxy -R -J -O -K -W1p smt.trk >> $ps
    37. 

  37. echo "-100 100 BATHYMETRY" | gmt pstext -R -J -O -K -F+jTL+f14p -Dj0.1i -Gwhite -C+tO >> $ps
    38. 

  38. # 2. Compute the VGG anomaly
    39. 

  39. gmt gravfft smt.nc+uk -D1670 -Nf+a -Fv -E$order -Gvgg.nc
    40. 

  40. # BR plot the VGG anomaly
    41. 

  41. gmt makecpt -Crainbow -T-50/250 > t.cpt
    42. 

  42. gmt grdimage vgg.nc -R-100/100/-100/100 -JX3i -O -Bag -BwSne -Ct.cpt -K -X3.5i >> $ps
    43. 

  43. gmt grdtrack -Gvgg.nc -ELM/RM+d > vgg.trk
    44. 

  44. gmt psxy -R -J -O -K -W5p,white vgg.trk >> $ps
    45. 

  45. gmt psxy -R -J -O -K -W1p,blue vgg.trk >> $ps
    46. 

  46. echo "-100 100 VGG" | gmt pstext -R -J -O -K -F+jTL+f14p -Dj0.1i -Gwhite -C+tO >> $ps
    47. 

  47. # 3. Compute the FAA anomaly
    48. 

  48. gmt gravfft smt.nc+uk -D1670 -Nf+a -Ff -E$order -Gfaa.nc
    49. 

  49. # Compute the exact analytical result for peak amplitude at center
    50. 

  50. r0=10000
    51. 

  51. z0=$(gmt math -Q 5084 3751 SUB =)
    52. 

  52. r1=35000
    53. 

  53. z1=5084
    54. 

  54. rho=$(gmt math -Q 2800 1030 SUB =)
    55. 

  55. b=$(gmt math -Q $r1 $r0 SUB $z1 $z0 SUB DIV =)
    56. 

  56. c=$(gmt math -Q $r0 $b $z0 MUL SUB =)
    57. 

  57. I1=$(Ugly $z1 $b $c)
    58. 

  58. I2=$(Ugly $z0 $b $c)
    59. 

  59. gmax=$(gmt math -Q 2 PI MUL 6.673e-6 MUL $rho MUL $I1 $I2 SUB MUL =)
    60. 

  60. echo "Max FAA should be $gmax mGal"
    61. 

  61. # ML plot the FAA anomaly
    62. 

  62. gmt makecpt -Crainbow -T-50/250 > t.cpt
    63. 

  63. gmt grdimage faa.nc -R-100/100/-100/100 -JX3i -O -Bag -BWsne -Ct.cpt -K -X-3.5i -Y3.25i >> $ps
    64. 

  64. gmt grdtrack -Gfaa.nc -ELM/RM+d > faa.trk
    65. 

  65. gmt psxy -R -J -O -K -W5p,white faa.trk >> $ps
    66. 

  66. gmt psxy -R -J -O -K -W1p,red faa.trk >> $ps
    67. 

  67. echo "-100 100 FAA" | gmt pstext -R -J -O -K -F+jTL+f14p -Dj0.1i -Gwhite -C+tO >> $ps
    68. 

  68. # 4. Compute the geoid anomaly
    69. 

  69. gmt gravfft smt.nc+uk -D1670 -Nf+a -Fg -E$order -Ggeoid.nc
    70. 

  70. # MR plot the VGG anomaly
    71. 

  71. gmt makecpt -Crainbow -T0/5 > t.cpt
    72. 

  72. gmt grdimage geoid.nc -R-100/100/-100/100 -JX3i -O -Bag -Bwsne -Ct.cpt -K -X3.5i >> $ps
    73. 

  73. gmt grdtrack -Ggeoid.nc -ELM/RM+d > geoid.trk
    74. 

  74. gmt psxy -R -J -O -K -W5p,white geoid.trk >> $ps
    75. 

  75. gmt psxy -R -J -O -K -W1p,orange geoid.trk >> $ps
    76. 

  76. echo "-100 100 GEOID" | gmt pstext -R -J -O -K -F+jTL+f14p -Dj0.1i -Gwhite -C+tO >> $ps
    77. 

  77. # 5 Plot crossections of bathy and faa crossections
    78. 

  78. # TL plot the bathy and faa canomaly
    79. 

  79. gmt psxy -R-100/100/-5100/1000 -JX3i/2.5i -O -K -W1p -i0,3 smt.trk -Baf -BWsN -X-3.5i -Y3.2i >> $ps
    80. 

  80. echo "-100 1000 TOPO" | gmt pstext -R -J -O -K -F+jTL+f12p -Dj0.1i/0.15i >> $ps
    81. 

  81. gmt psxy -R-100/100/-50/250 -J -O -K -W1p,red -i0,3 faa.trk -Bafg1000 -BENs >> $ps
    82. 

  82. gmt psxy -R -J -O -K -W0.5p,- << EOF >> $ps
    83. 

  83. -100	$gmax
    84. 

  84. +100	$gmax
    85. 

  85. EOF
    86. 

  86. echo "100 250 FAA" | gmt pstext -R -J -O -K -F+jTR+f12p,Helvetica,red -Dj0.1i/0.15i >> $ps
    87. 

  87. # Add VGG and geoid crossections
    88. 

  88. # TRL plot the VGG and geoid anomaly
    89. 

  89. gmt psxy -R-100/100/-50/250 -JX3i/2.5i -O -K -W1p,blue -i0,3 vgg.trk -Bafg1000 -BwsN -X3.5i >> $ps
    90. 

  90. echo "-100 250 VGG" | gmt pstext -R -J -O -K -F+jTL+f12p,Helvetica,blue -Dj0.1i >> $ps
    91. 

  91. gmt psxy -R-100/100/0/4 -J -O -K -W1p,orange -i0,3 geoid.trk -Baf -BE >> $ps
    92. 

  92. echo "100 4 GEOID" | gmt pstext -R -J -O -K -F+jTR+f12p,Helvetica,orange -Dj0.1i >> $ps
    93. 

  93. gmt psxy -R -J -O -T >> $ps
    94.
Prev
Next
Tip!

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

Comments

Write Preview
Loading...
Integrate AWS S3
Use S3 remote
Select bucket
Access key
Finish

Use AWS S3 as storage!

Browsing data directories saved to S3 is possible with DAGsHub. Let's configure your repository to easily display your data in the context of any commit!

Specify your S3 bucket

Access key (If needed)

Congratulations!

gmt-fork is now integrated with AWS S3!

Delete Storage Key

Are you sure you want to delete this access key?

No
Yes
Integrate Google Cloud Storage
Use Google Storage
Select bucket
Upload key
Finish

Use Google Cloud Storage!

Browsing data directories saved to Google Cloud Storage is possible with DAGsHub. Let's configure your repository to easily display your data in the context of any commit!

Specify your Google Storage bucket

Service Account Key

Congratulations!

gmt-fork is now integrated with Google Cloud Storage!

Delete Storage Key

Are you sure you want to delete this access key?

No
Yes
Integrate Azure Cloud Storage
Use Azure Storage
Select bucket
Set key
Finish

Use Azure Cloud Storage!

Browsing data directories saved to Azure Cloud Storage is possible with DAGsHub. Let's configure your repository to easily display your data in the context of any commit!

Specify your Azure Storage bucket

Access key (If needed)

Congratulations!

gmt-fork is now integrated with Azure Cloud Storage!

Delete Storage Key

Are you sure you want to delete this access key?

No
Yes
Integrate S3 compatible storage
Use S3 like remote
Select bucket
Access key
Finish

Use any S3 compatible storage!

Browsing data directories saved to S3 compatible storage is possible with DAGsHub. Let's configure your repository to easily display your data in the context of any commit!

Specify your S3 bucket

Bucket name cannot be the same as the repository name. Please change one of them.

Access key (If needed)

Congratulations!

gmt-fork is now integrated with your S3 compatible storage!

Delete Storage Key

Are you sure you want to delete this access key?

No
Yes