Rayleigh dispersion curves

- Gutenberg model

Author: pawbz

src/introduction/waves-on-string.jl

@@ -225,7 +225,7 @@ PlutoUI = "~0.7.65"
 PLUTO_MANIFEST_TOML_CONTENTS = """
 # This file is machine-generated - editing it directly is not advised
 
-julia_version = "1.11.5"
+julia_version = "1.11.6"
 manifest_format = "2.0"
 project_hash = "b9d0c54cc1aed75c67eb364c4722edbde56ef19c"
 

src/surfacewaves/Love-wave-analytic1.jl

@@ -630,7 +630,7 @@ Symbolics = "~6.12.0"
 PLUTO_MANIFEST_TOML_CONTENTS = """
 # This file is machine-generated - editing it directly is not advised
 
-julia_version = "1.11.5"
+julia_version = "1.11.6"
 manifest_format = "2.0"
 project_hash = "6148517c7c28f34361da64be253ff296fe4fb695"
 

src/surfacewaves/Love-wave-dispersion-curves.jl

@@ -46,11 +46,6 @@ Instructor: *Pawan Bharadwaj*,
 Indian Institute of Science, Bengaluru, India
 """
 
-# ╔═╡ 1af153bc-e471-4608-984b-61f1116dfa16
-md"""
-Number of layers: $(@bind n_layers Slider(1:5, default=2, show_value=true))
-"""
-
 # ╔═╡ d2c48242-92ce-11f0-01e6-cb1b8fe1fecb
 
 struct Layer
@@ -100,7 +95,7 @@ struct DispersionResult
 end
 
 # ╔═╡ 756e06e3-4f25-4f55-b242-26a45a1b7dd3
-periods = 5.0:5.0:100.0
+periods = collect(range(5.0, stop=100.0, length=64))
 
 # ╔═╡ 9c9ad5f6-debf-4376-aaf6-759e9b0c949b
 """
@@ -317,27 +312,62 @@ end
 # ╔═╡ 80086cb5-7e39-444b-a85e-ddf46a67ec42
 md"## Appendix"
 
+# ╔═╡ 5a43781a-b3e4-4bac-86d1-bd822c169804
+const GUTENBERG_MODEL = [
+    Layer(19.0, 6.14, 3.55, 2.74, 1000.0, 1000.0),
+    Layer(19.0, 6.58, 3.80, 3.00, 1000.0, 1000.0),
+    Layer(12.0, 8.20, 4.65, 3.32, 1000.0, 1000.0),
+    Layer(10.0, 8.17, 4.62, 3.34, 1000.0, 1000.0),
+    Layer(10.0, 8.14, 4.57, 3.35, 1000.0, 1000.0),
+    Layer(10.0, 8.10, 4.51, 3.36, 1000.0, 1000.0),
+    Layer(10.0, 8.07, 4.46, 3.37, 1000.0, 1000.0),
+    Layer(10.0, 8.02, 4.41, 3.38, 1000.0, 1000.0),
+    Layer(25.0, 7.93, 4.37, 3.39, 1000.0, 1000.0),
+    Layer(25.0, 7.85, 4.35, 3.41, 1000.0, 1000.0),
+    Layer(25.0, 7.89, 4.36, 3.43, 1000.0, 1000.0),
+    Layer(25.0, 7.98, 4.38, 3.46, 1000.0, 1000.0),
+    Layer(25.0, 8.10, 4.42, 3.48, 1000.0, 1000.0),
+    Layer(25.0, 8.21, 4.46, 3.50, 1000.0, 1000.0),
+    Layer(50.0, 8.38, 4.54, 3.53, 1000.0, 1000.0),
+    Layer(50.0, 8.62, 4.68, 3.58, 1000.0, 1000.0),
+    Layer(50.0, 8.87, 4.85, 3.62, 1000.0, 1000.0),
+    Layer(50.0, 9.15, 5.04, 3.69, 1000.0, 1000.0),
+    Layer(50.0, 9.45, 5.21, 3.82, 1000.0, 1000.0),
+    Layer(100.0, 9.88, 5.45, 4.01, 1000.0, 1000.0),
+    Layer(100.0, 10.30, 5.76, 4.21, 1000.0, 1000.0),
+    Layer(100.0, 10.71, 6.03, 4.40, 1000.0, 1000.0),
+    Layer(100.0, 11.10, 6.23, 4.56, 1000.0, 1000.0),
+    Layer(100.0, 11.35, 6.32, 4.63, 1000.0, 1000.0)
+]
+
+# ╔═╡ 05d910d2-527c-416f-9d63-c259fbe8a45d
+default_layers = GUTENBERG_MODEL;
+
+# ╔═╡ 1af153bc-e471-4608-984b-61f1116dfa16
+md"""
+Number of layers: $(@bind n_layers Slider(1:length(default_layers), default=length(default_layers), show_value=true))
+"""
+
 # ╔═╡ 99e9fc2c-fe13-4f94-b613-86cedb0f3653
 function layer_table_input(n_layers::Int; vp_vs_ratio=1.73)
     ui = PlutoUI.combine() do Child
         # header
         header = @htl("""
-        <tr style="background:#f0f3f8; text-align:center;">
+        <tr style="text-align:center;">
           <th style="padding:6px;">#</th>
-          <th style="padding:6px;">Thickness</th>
-          <th style="padding:6px;">Vs</th>
-          <th style="padding:6px;">Density</th>
+          <th style="padding:6px;">Thickness (km)</th>
+          <th style="padding:6px;">Vs (km/s)</th>
+          <th style="padding:6px;">Density (gm/cc)</th>
         </tr>
         """)
 
-        # rows
-        rows = Any[]
+		 rows = Any[]
         for i in 1:n_layers
-            # no thickness slider for last layer (half-space)
-            tw = i < n_layers ? Child("thickness_$i", Slider(0.5:0.5:200, default=35, show_value=true)) : "∞"
-
-            vs = Child("vs_$i", Slider(2.5:0.01:5.5, default=(i == 1) ? 3.4 : 4.6, show_value=true))
-            rw = Child("rho_$i", Slider(1.0:0.01:6.0, default=(i == 1) ? 2.6 : 3.4, show_value=true))
+            # Use Gutenberg model values if available, else fallback
+            layer = i <= length(default_layers) ? default_layers[i] : default_layers[end]
+            tw = i < n_layers ? Child("thickness_$i", Slider(0.5:0.5:200, default=layer.thickness, show_value=true)) : "∞"
+            vs = Child("vs_$i", Slider(2.5:0.01:8.0, default=layer.vs, show_value=true))
+            rw = Child("rho_$i", Slider(1.0:0.01:6.0, default=layer.rho, show_value=true))
 
             push!(rows, @htl("""
                 <tr>
@@ -349,6 +379,7 @@ function layer_table_input(n_layers::Int; vp_vs_ratio=1.73)
             """))
         end
 
+
         tbl = @htl("""
         <table style="border-collapse:collapse; border:1px solid #ddd; width:100%;">
           <thead>$header</thead>
@@ -363,8 +394,7 @@ function layer_table_input(n_layers::Int; vp_vs_ratio=1.73)
         Number of layers: **$n_layers**
 
         $tbl
-        _Notes:_ Thickness in km.  
-        Constraint enforced: **Vp = $(vp_vs_ratio) × Vs > Vs**  
+        _Notes:_  
         The last layer is treated as a half-space (no thickness).
         """
     end
@@ -453,7 +483,7 @@ Roots = "f2b01f46-fcfa-551c-844a-d8ac1e96c665"
 
 [compat]
 HypertextLiteral = "~0.9.5"
-PlutoPlotly = "~0.3.4"
+PlutoPlotly = "~0.6.4"
 PlutoUI = "~0.7.71"
 Roots = "~2.2.6"
 """
@@ -464,7 +494,7 @@ PLUTO_MANIFEST_TOML_CONTENTS = """
 
 julia_version = "1.11.6"
 manifest_format = "2.0"
-project_hash = "9038c24234f0293a4bad7fb0208acaec306d02e2"
+project_hash = "0264277e532391f8f3beb74860e3dc79469a39a9"
 
 [[deps.AbstractPlutoDingetjes]]
 deps = ["Pkg"]
@@ -516,19 +546,15 @@ version = "3.30.0"
 
 [[deps.ColorTypes]]
 deps = ["FixedPointNumbers", "Random"]
-git-tree-sha1 = "67e11ee83a43eb71ddc950302c53bf33f0690dfe"
+git-tree-sha1 = "b10d0b65641d57b8b4d5e234446582de5047050d"
 uuid = "3da002f7-5984-5a60-b8a6-cbb66c0b333f"
-version = "0.12.1"
-weakdeps = ["StyledStrings"]
-
-    [deps.ColorTypes.extensions]
-    StyledStringsExt = "StyledStrings"
+version = "0.11.5"
 
 [[deps.ColorVectorSpace]]
 deps = ["ColorTypes", "FixedPointNumbers", "LinearAlgebra", "Requires", "Statistics", "TensorCore"]
-git-tree-sha1 = "8b3b6f87ce8f65a2b4f857528fd8d70086cd72b1"
+git-tree-sha1 = "a1f44953f2382ebb937d60dafbe2deea4bd23249"
 uuid = "c3611d14-8923-5661-9e6a-0046d554d3a4"
-version = "0.11.0"
+version = "0.10.0"
 
     [deps.ColorVectorSpace.extensions]
     SpecialFunctionsExt = "SpecialFunctions"
@@ -538,9 +564,9 @@ version = "0.11.0"
 
 [[deps.Colors]]
 deps = ["ColorTypes", "FixedPointNumbers", "Reexport"]
-git-tree-sha1 = "37ea44092930b1811e666c3bc38065d7d87fcc74"
+git-tree-sha1 = "362a287c3aa50601b0bc359053d5c2468f0e7ce0"
 uuid = "5ae59095-9a9b-59fe-a467-6f913c188581"
-version = "0.13.1"
+version = "0.12.11"
 
 [[deps.CommonSolve]]
 git-tree-sha1 = "0eee5eb66b1cf62cd6ad1b460238e60e4b09400c"
@@ -607,6 +633,11 @@ git-tree-sha1 = "05882d6995ae5c12bb5f36dd2ed3f61c98cbb172"
 uuid = "53c48c17-4a7d-5ca2-90c5-79b7896eea93"
 version = "0.8.5"
 
+[[deps.HashArrayMappedTries]]
+git-tree-sha1 = "2eaa69a7cab70a52b9687c8bf950a5a93ec895ae"
+uuid = "076d061b-32b6-4027-95e0-9a2c6f6d7e74"
+version = "0.2.0"
+
 [[deps.Hyperscript]]
 deps = ["Test"]
 git-tree-sha1 = "179267cfa5e712760cd43dcae385d7ea90cc25a4"
@@ -771,10 +802,18 @@ version = "0.8.21"
     JSON3 = "0f8b85d8-7281-11e9-16c2-39a750bddbf1"
 
 [[deps.PlutoPlotly]]
-deps = ["AbstractPlutoDingetjes", "Dates", "HypertextLiteral", "InteractiveUtils", "LaTeXStrings", "Markdown", "PlotlyBase", "PlutoUI", "Reexport"]
-git-tree-sha1 = "b470931aa2a8112c8b08e66ea096c6c62c60571e"
+deps = ["AbstractPlutoDingetjes", "Artifacts", "ColorSchemes", "Colors", "Dates", "Downloads", "HypertextLiteral", "InteractiveUtils", "LaTeXStrings", "Markdown", "Pkg", "PlotlyBase", "PrecompileTools", "Reexport", "ScopedValues", "Scratch", "TOML"]
+git-tree-sha1 = "232630fee92e588c11c2b260741b4fa70784b4c5"
 uuid = "8e989ff0-3d88-8e9f-f020-2b208a939ff0"
-version = "0.3.4"
+version = "0.6.4"
+
+    [deps.PlutoPlotly.extensions]
+    PlotlyKaleidoExt = "PlotlyKaleido"
+    UnitfulExt = "Unitful"
+
+    [deps.PlutoPlotly.weakdeps]
+    PlotlyKaleido = "f2990250-8cf9-495f-b13a-cce12b45703c"
+    Unitful = "1986cc42-f94f-5a68-af5c-568840ba703d"
 
 [[deps.PlutoUI]]
 deps = ["AbstractPlutoDingetjes", "Base64", "ColorTypes", "Dates", "Downloads", "FixedPointNumbers", "Hyperscript", "HypertextLiteral", "IOCapture", "InteractiveUtils", "JSON", "Logging", "MIMEs", "Markdown", "Random", "Reexport", "URIs", "UUIDs"]
@@ -844,6 +883,18 @@ version = "2.2.6"
 uuid = "ea8e919c-243c-51af-8825-aaa63cd721ce"
 version = "0.7.0"
 
+[[deps.ScopedValues]]
+deps = ["HashArrayMappedTries", "Logging"]
+git-tree-sha1 = "c3b2323466378a2ba15bea4b2f73b081e022f473"
+uuid = "7e506255-f358-4e82-b7e4-beb19740aa63"
+version = "1.5.0"
+
+[[deps.Scratch]]
+deps = ["Dates"]
+git-tree-sha1 = "9b81b8393e50b7d4e6d0a9f14e192294d3b7c109"
+uuid = "6c6a2e73-6563-6170-7368-637461726353"
+version = "1.3.0"
+
 [[deps.Serialization]]
 uuid = "9e88b42a-f829-5b0c-bbe9-9e923198166b"
 version = "1.11.0"
@@ -958,6 +1009,8 @@ version = "17.4.0+2"
 # ╠═29afa36b-391f-4861-aead-d62918daf3c6
 # ╠═8d5d9594-2197-4ccc-a7a4-0f0d54a2370a
 # ╠═99e9fc2c-fe13-4f94-b613-86cedb0f3653
+# ╠═5a43781a-b3e4-4bac-86d1-bd822c169804
+# ╠═05d910d2-527c-416f-9d63-c259fbe8a45d
 # ╟─fc1fd90c-020a-4f2d-aedd-66115ac6a287
 # ╟─00000000-0000-0000-0000-000000000001
 # ╟─00000000-0000-0000-0000-000000000002

src/surfacewaves/Rayleigh-function.jl

@@ -1,5 +1,5 @@
 ### A Pluto.jl notebook ###
-# v0.20.5
+# v0.20.13
 
 #> [frontmatter]
 #> title = "Rayleigh Function"
@@ -418,7 +418,7 @@ Symbolics = "~6.12.0"
 PLUTO_MANIFEST_TOML_CONTENTS = """
 # This file is machine-generated - editing it directly is not advised
 
-julia_version = "1.11.4"
+julia_version = "1.11.6"
 manifest_format = "2.0"
 project_hash = "a4c9aa60b01da02b7637d7404f9276db98797133"
 
@@ -1044,7 +1044,7 @@ version = "0.3.27+1"
 [[deps.OpenLibm_jll]]
 deps = ["Artifacts", "Libdl"]
 uuid = "05823500-19ac-5b8b-9628-191a04bc5112"
-version = "0.8.1+4"
+version = "0.8.5+0"
 
 [[deps.OpenSpecFun_jll]]
 deps = ["Artifacts", "CompilerSupportLibraries_jll", "JLLWrappers", "Libdl"]