diff --git a/.gitignore b/.gitignore
index 51e0643..a71ce7f 100644
--- a/.gitignore
+++ b/.gitignore
@@ -181,3 +181,5 @@ data
 
 # DS_STORE...
 .DS_STORE
+
+*.pdf
diff --git a/Manifest.toml b/Manifest.toml
deleted file mode 100644
index 57bab97..0000000
--- a/Manifest.toml
+++ /dev/null
@@ -1,1311 +0,0 @@
-# This file is machine-generated - editing it directly is not advised
-
-julia_version = "1.10.0"
-manifest_format = "2.0"
-project_hash = "31bf1a950631602335b6f5022756aca63671050c"
-
-[[deps.AbstractFFTs]]
-deps = ["LinearAlgebra"]
-git-tree-sha1 = "d92ad398961a3ed262d8bf04a1a2b8340f915fef"
-uuid = "621f4979-c628-5d54-868e-fcf4e3e8185c"
-version = "1.5.0"
-weakdeps = ["ChainRulesCore", "Test"]
-
-    [deps.AbstractFFTs.extensions]
-    AbstractFFTsChainRulesCoreExt = "ChainRulesCore"
-    AbstractFFTsTestExt = "Test"
-
-[[deps.Adapt]]
-deps = ["LinearAlgebra", "Requires"]
-git-tree-sha1 = "6a55b747d1812e699320963ffde36f1ebdda4099"
-uuid = "79e6a3ab-5dfb-504d-930d-738a2a938a0e"
-version = "4.0.4"
-weakdeps = ["StaticArrays"]
-
-    [deps.Adapt.extensions]
-    AdaptStaticArraysExt = "StaticArrays"
-
-[[deps.ArgParse]]
-deps = ["Logging", "TextWrap"]
-git-tree-sha1 = "22cf435ac22956a7b45b0168abbc871176e7eecc"
-uuid = "c7e460c6-2fb9-53a9-8c5b-16f535851c63"
-version = "1.2.0"
-
-[[deps.ArgTools]]
-uuid = "0dad84c5-d112-42e6-8d28-ef12dabb789f"
-version = "1.1.1"
-
-[[deps.ArrayInterface]]
-deps = ["Adapt", "LinearAlgebra", "SparseArrays", "SuiteSparse"]
-git-tree-sha1 = "133a240faec6e074e07c31ee75619c90544179cf"
-uuid = "4fba245c-0d91-5ea0-9b3e-6abc04ee57a9"
-version = "7.10.0"
-
-    [deps.ArrayInterface.extensions]
-    ArrayInterfaceBandedMatricesExt = "BandedMatrices"
-    ArrayInterfaceBlockBandedMatricesExt = "BlockBandedMatrices"
-    ArrayInterfaceCUDAExt = "CUDA"
-    ArrayInterfaceCUDSSExt = "CUDSS"
-    ArrayInterfaceChainRulesExt = "ChainRules"
-    ArrayInterfaceGPUArraysCoreExt = "GPUArraysCore"
-    ArrayInterfaceReverseDiffExt = "ReverseDiff"
-    ArrayInterfaceStaticArraysCoreExt = "StaticArraysCore"
-    ArrayInterfaceTrackerExt = "Tracker"
-
-    [deps.ArrayInterface.weakdeps]
-    BandedMatrices = "aae01518-5342-5314-be14-df237901396f"
-    BlockBandedMatrices = "ffab5731-97b5-5995-9138-79e8c1846df0"
-    CUDA = "052768ef-5323-5732-b1bb-66c8b64840ba"
-    CUDSS = "45b445bb-4962-46a0-9369-b4df9d0f772e"
-    ChainRules = "082447d4-558c-5d27-93f4-14fc19e9eca2"
-    GPUArraysCore = "46192b85-c4d5-4398-a991-12ede77f4527"
-    ReverseDiff = "37e2e3b7-166d-5795-8a7a-e32c996b4267"
-    StaticArraysCore = "1e83bf80-4336-4d27-bf5d-d5a4f845583c"
-    Tracker = "9f7883ad-71c0-57eb-9f7f-b5c9e6d3789c"
-
-[[deps.Artifacts]]
-uuid = "56f22d72-fd6d-98f1-02f0-08ddc0907c33"
-
-[[deps.Atomix]]
-deps = ["UnsafeAtomics"]
-git-tree-sha1 = "c06a868224ecba914baa6942988e2f2aade419be"
-uuid = "a9b6321e-bd34-4604-b9c9-b65b8de01458"
-version = "0.1.0"
-
-[[deps.AxisAlgorithms]]
-deps = ["LinearAlgebra", "Random", "SparseArrays", "WoodburyMatrices"]
-git-tree-sha1 = "01b8ccb13d68535d73d2b0c23e39bd23155fb712"
-uuid = "13072b0f-2c55-5437-9ae7-d433b7a33950"
-version = "1.1.0"
-
-[[deps.BFloat16s]]
-deps = ["LinearAlgebra", "Printf", "Random", "Test"]
-git-tree-sha1 = "2c7cc21e8678eff479978a0a2ef5ce2f51b63dff"
-uuid = "ab4f0b2a-ad5b-11e8-123f-65d77653426b"
-version = "0.5.0"
-
-[[deps.Base64]]
-uuid = "2a0f44e3-6c83-55bd-87e4-b1978d98bd5f"
-
-[[deps.Blosc_jll]]
-deps = ["Artifacts", "JLLWrappers", "Libdl", "Lz4_jll", "Zlib_jll", "Zstd_jll"]
-git-tree-sha1 = "19b98ee7e3db3b4eff74c5c9c72bf32144e24f10"
-uuid = "0b7ba130-8d10-5ba8-a3d6-c5182647fed9"
-version = "1.21.5+0"
-
-[[deps.Bzip2_jll]]
-deps = ["Artifacts", "JLLWrappers", "Libdl", "Pkg"]
-git-tree-sha1 = "9e2a6b69137e6969bab0152632dcb3bc108c8bdd"
-uuid = "6e34b625-4abd-537c-b88f-471c36dfa7a0"
-version = "1.0.8+1"
-
-[[deps.CEnum]]
-git-tree-sha1 = "389ad5c84de1ae7cf0e28e381131c98ea87d54fc"
-uuid = "fa961155-64e5-5f13-b03f-caf6b980ea82"
-version = "0.5.0"
-
-[[deps.CFTime]]
-deps = ["Dates", "Printf"]
-git-tree-sha1 = "5afb5c5ba2688ca43a9ad2e5a91cbb93921ccfa1"
-uuid = "179af706-886a-5703-950a-314cd64e0468"
-version = "0.1.3"
-
-[[deps.CUDA]]
-deps = ["AbstractFFTs", "Adapt", "BFloat16s", "CEnum", "CUDA_Driver_jll", "CUDA_Runtime_Discovery", "CUDA_Runtime_jll", "Crayons", "DataFrames", "ExprTools", "GPUArrays", "GPUCompiler", "KernelAbstractions", "LLVM", "LLVMLoopInfo", "LazyArtifacts", "Libdl", "LinearAlgebra", "Logging", "NVTX", "Preferences", "PrettyTables", "Printf", "Random", "Random123", "RandomNumbers", "Reexport", "Requires", "SparseArrays", "StaticArrays", "Statistics"]
-git-tree-sha1 = "fe61a257e94621e25471071ca58d29ea45eef13b"
-uuid = "052768ef-5323-5732-b1bb-66c8b64840ba"
-version = "5.3.4"
-
-    [deps.CUDA.extensions]
-    ChainRulesCoreExt = "ChainRulesCore"
-    EnzymeCoreExt = "EnzymeCore"
-    SpecialFunctionsExt = "SpecialFunctions"
-
-    [deps.CUDA.weakdeps]
-    ChainRulesCore = "d360d2e6-b24c-11e9-a2a3-2a2ae2dbcce4"
-    EnzymeCore = "f151be2c-9106-41f4-ab19-57ee4f262869"
-    SpecialFunctions = "276daf66-3868-5448-9aa4-cd146d93841b"
-
-[[deps.CUDA_Driver_jll]]
-deps = ["Artifacts", "JLLWrappers", "LazyArtifacts", "Libdl", "Pkg"]
-git-tree-sha1 = "dc172b558adbf17952001e15cf0d6364e6d78c2f"
-uuid = "4ee394cb-3365-5eb0-8335-949819d2adfc"
-version = "0.8.1+0"
-
-[[deps.CUDA_Runtime_Discovery]]
-deps = ["Libdl"]
-git-tree-sha1 = "38f830504358e9972d2a0c3e5d51cb865e0733df"
-uuid = "1af6417a-86b4-443c-805f-a4643ffb695f"
-version = "0.2.4"
-
-[[deps.CUDA_Runtime_jll]]
-deps = ["Artifacts", "CUDA_Driver_jll", "JLLWrappers", "LazyArtifacts", "Libdl", "TOML"]
-git-tree-sha1 = "4ca7d6d92075906c2ce871ea8bba971fff20d00c"
-uuid = "76a88914-d11a-5bdc-97e0-2f5a05c973a2"
-version = "0.12.1+0"
-
-[[deps.ChainRulesCore]]
-deps = ["Compat", "LinearAlgebra"]
-git-tree-sha1 = "575cd02e080939a33b6df6c5853d14924c08e35b"
-uuid = "d360d2e6-b24c-11e9-a2a3-2a2ae2dbcce4"
-version = "1.23.0"
-weakdeps = ["SparseArrays"]
-
-    [deps.ChainRulesCore.extensions]
-    ChainRulesCoreSparseArraysExt = "SparseArrays"
-
-[[deps.ColorTypes]]
-deps = ["FixedPointNumbers", "Random"]
-git-tree-sha1 = "b10d0b65641d57b8b4d5e234446582de5047050d"
-uuid = "3da002f7-5984-5a60-b8a6-cbb66c0b333f"
-version = "0.11.5"
-
-[[deps.Colors]]
-deps = ["ColorTypes", "FixedPointNumbers", "Reexport"]
-git-tree-sha1 = "362a287c3aa50601b0bc359053d5c2468f0e7ce0"
-uuid = "5ae59095-9a9b-59fe-a467-6f913c188581"
-version = "0.12.11"
-
-[[deps.CommonDataModel]]
-deps = ["CFTime", "DataStructures", "Dates", "Preferences", "Printf", "Statistics"]
-git-tree-sha1 = "d6fb5bf939a2753c74984b11434ea25d6c397a58"
-uuid = "1fbeeb36-5f17-413c-809b-666fb144f157"
-version = "0.3.6"
-
-[[deps.CommonSubexpressions]]
-deps = ["MacroTools", "Test"]
-git-tree-sha1 = "7b8a93dba8af7e3b42fecabf646260105ac373f7"
-uuid = "bbf7d656-a473-5ed7-a52c-81e309532950"
-version = "0.3.0"
-
-[[deps.Compat]]
-deps = ["TOML", "UUIDs"]
-git-tree-sha1 = "b1c55339b7c6c350ee89f2c1604299660525b248"
-uuid = "34da2185-b29b-5c13-b0c7-acf172513d20"
-version = "4.15.0"
-weakdeps = ["Dates", "LinearAlgebra"]
-
-    [deps.Compat.extensions]
-    CompatLinearAlgebraExt = "LinearAlgebra"
-
-[[deps.CompilerSupportLibraries_jll]]
-deps = ["Artifacts", "Libdl"]
-uuid = "e66e0078-7015-5450-92f7-15fbd957f2ae"
-version = "1.0.5+1"
-
-[[deps.ConstructionBase]]
-deps = ["LinearAlgebra"]
-git-tree-sha1 = "260fd2400ed2dab602a7c15cf10c1933c59930a2"
-uuid = "187b0558-2788-49d3-abe0-74a17ed4e7c9"
-version = "1.5.5"
-
-    [deps.ConstructionBase.extensions]
-    ConstructionBaseIntervalSetsExt = "IntervalSets"
-    ConstructionBaseStaticArraysExt = "StaticArrays"
-
-    [deps.ConstructionBase.weakdeps]
-    IntervalSets = "8197267c-284f-5f27-9208-e0e47529a953"
-    StaticArrays = "90137ffa-7385-5640-81b9-e52037218182"
-
-[[deps.Crayons]]
-git-tree-sha1 = "249fe38abf76d48563e2f4556bebd215aa317e15"
-uuid = "a8cc5b0e-0ffa-5ad4-8c14-923d3ee1735f"
-version = "4.1.1"
-
-[[deps.CubedSphere]]
-deps = ["Elliptic", "FFTW", "Printf", "ProgressBars", "SpecialFunctions", "TaylorSeries", "Test"]
-git-tree-sha1 = "10134667d7d3569b191a65801514271b8a93b292"
-uuid = "7445602f-e544-4518-8976-18f8e8ae6cdb"
-version = "0.2.5"
-
-[[deps.DataAPI]]
-git-tree-sha1 = "abe83f3a2f1b857aac70ef8b269080af17764bbe"
-uuid = "9a962f9c-6df0-11e9-0e5d-c546b8b5ee8a"
-version = "1.16.0"
-
-[[deps.DataFrames]]
-deps = ["Compat", "DataAPI", "DataStructures", "Future", "InlineStrings", "InvertedIndices", "IteratorInterfaceExtensions", "LinearAlgebra", "Markdown", "Missings", "PooledArrays", "PrecompileTools", "PrettyTables", "Printf", "REPL", "Random", "Reexport", "SentinelArrays", "SortingAlgorithms", "Statistics", "TableTraits", "Tables", "Unicode"]
-git-tree-sha1 = "04c738083f29f86e62c8afc341f0967d8717bdb8"
-uuid = "a93c6f00-e57d-5684-b7b6-d8193f3e46c0"
-version = "1.6.1"
-
-[[deps.DataStructures]]
-deps = ["Compat", "InteractiveUtils", "OrderedCollections"]
-git-tree-sha1 = "1d0a14036acb104d9e89698bd408f63ab58cdc82"
-uuid = "864edb3b-99cc-5e75-8d2d-829cb0a9cfe8"
-version = "0.18.20"
-
-[[deps.DataValueInterfaces]]
-git-tree-sha1 = "bfc1187b79289637fa0ef6d4436ebdfe6905cbd6"
-uuid = "e2d170a0-9d28-54be-80f0-106bbe20a464"
-version = "1.0.0"
-
-[[deps.Dates]]
-deps = ["Printf"]
-uuid = "ade2ca70-3891-5945-98fb-dc099432e06a"
-
-[[deps.DiffResults]]
-deps = ["StaticArraysCore"]
-git-tree-sha1 = "782dd5f4561f5d267313f23853baaaa4c52ea621"
-uuid = "163ba53b-c6d8-5494-b064-1a9d43ac40c5"
-version = "1.1.0"
-
-[[deps.DiffRules]]
-deps = ["IrrationalConstants", "LogExpFunctions", "NaNMath", "Random", "SpecialFunctions"]
-git-tree-sha1 = "23163d55f885173722d1e4cf0f6110cdbaf7e272"
-uuid = "b552c78f-8df3-52c6-915a-8e097449b14b"
-version = "1.15.1"
-
-[[deps.DiskArrays]]
-deps = ["LRUCache", "OffsetArrays"]
-git-tree-sha1 = "ef25c513cad08d7ebbed158c91768ae32f308336"
-uuid = "3c3547ce-8d99-4f5e-a174-61eb10b00ae3"
-version = "0.3.23"
-
-[[deps.Distances]]
-deps = ["LinearAlgebra", "Statistics", "StatsAPI"]
-git-tree-sha1 = "66c4c81f259586e8f002eacebc177e1fb06363b0"
-uuid = "b4f34e82-e78d-54a5-968a-f98e89d6e8f7"
-version = "0.10.11"
-weakdeps = ["ChainRulesCore", "SparseArrays"]
-
-    [deps.Distances.extensions]
-    DistancesChainRulesCoreExt = "ChainRulesCore"
-    DistancesSparseArraysExt = "SparseArrays"
-
-[[deps.Distributed]]
-deps = ["Random", "Serialization", "Sockets"]
-uuid = "8ba89e20-285c-5b6f-9357-94700520ee1b"
-
-[[deps.DocStringExtensions]]
-deps = ["LibGit2"]
-git-tree-sha1 = "2fb1e02f2b635d0845df5d7c167fec4dd739b00d"
-uuid = "ffbed154-4ef7-542d-bbb7-c09d3a79fcae"
-version = "0.9.3"
-
-[[deps.Downloads]]
-deps = ["ArgTools", "FileWatching", "LibCURL", "NetworkOptions"]
-uuid = "f43a241f-c20a-4ad4-852c-f6b1247861c6"
-version = "1.6.0"
-
-[[deps.Elliptic]]
-git-tree-sha1 = "71c79e77221ab3a29918aaf6db4f217b89138608"
-uuid = "b305315f-e792-5b7a-8f41-49f472929428"
-version = "1.0.1"
-
-[[deps.ExprTools]]
-git-tree-sha1 = "27415f162e6028e81c72b82ef756bf321213b6ec"
-uuid = "e2ba6199-217a-4e67-a87a-7c52f15ade04"
-version = "0.1.10"
-
-[[deps.FFTW]]
-deps = ["AbstractFFTs", "FFTW_jll", "LinearAlgebra", "MKL_jll", "Preferences", "Reexport"]
-git-tree-sha1 = "4820348781ae578893311153d69049a93d05f39d"
-uuid = "7a1cc6ca-52ef-59f5-83cd-3a7055c09341"
-version = "1.8.0"
-
-[[deps.FFTW_jll]]
-deps = ["Artifacts", "JLLWrappers", "Libdl", "Pkg"]
-git-tree-sha1 = "c6033cc3892d0ef5bb9cd29b7f2f0331ea5184ea"
-uuid = "f5851436-0d7a-5f13-b9de-f02708fd171a"
-version = "3.3.10+0"
-
-[[deps.FileIO]]
-deps = ["Pkg", "Requires", "UUIDs"]
-git-tree-sha1 = "82d8afa92ecf4b52d78d869f038ebfb881267322"
-uuid = "5789e2e9-d7fb-5bc7-8068-2c6fae9b9549"
-version = "1.16.3"
-
-[[deps.FileWatching]]
-uuid = "7b1f6079-737a-58dc-b8bc-7a2ca5c1b5ee"
-
-[[deps.FiniteDiff]]
-deps = ["ArrayInterface", "LinearAlgebra", "Requires", "Setfield", "SparseArrays"]
-git-tree-sha1 = "2de436b72c3422940cbe1367611d137008af7ec3"
-uuid = "6a86dc24-6348-571c-b903-95158fe2bd41"
-version = "2.23.1"
-
-    [deps.FiniteDiff.extensions]
-    FiniteDiffBandedMatricesExt = "BandedMatrices"
-    FiniteDiffBlockBandedMatricesExt = "BlockBandedMatrices"
-    FiniteDiffStaticArraysExt = "StaticArrays"
-
-    [deps.FiniteDiff.weakdeps]
-    BandedMatrices = "aae01518-5342-5314-be14-df237901396f"
-    BlockBandedMatrices = "ffab5731-97b5-5995-9138-79e8c1846df0"
-    StaticArrays = "90137ffa-7385-5640-81b9-e52037218182"
-
-[[deps.FixedPointNumbers]]
-deps = ["Statistics"]
-git-tree-sha1 = "05882d6995ae5c12bb5f36dd2ed3f61c98cbb172"
-uuid = "53c48c17-4a7d-5ca2-90c5-79b7896eea93"
-version = "0.8.5"
-
-[[deps.ForwardDiff]]
-deps = ["CommonSubexpressions", "DiffResults", "DiffRules", "LinearAlgebra", "LogExpFunctions", "NaNMath", "Preferences", "Printf", "Random", "SpecialFunctions"]
-git-tree-sha1 = "cf0fe81336da9fb90944683b8c41984b08793dad"
-uuid = "f6369f11-7733-5829-9624-2563aa707210"
-version = "0.10.36"
-weakdeps = ["StaticArrays"]
-
-    [deps.ForwardDiff.extensions]
-    ForwardDiffStaticArraysExt = "StaticArrays"
-
-[[deps.Future]]
-deps = ["Random"]
-uuid = "9fa8497b-333b-5362-9e8d-4d0656e87820"
-
-[[deps.GMP_jll]]
-deps = ["Artifacts", "Libdl"]
-uuid = "781609d7-10c4-51f6-84f2-b8444358ff6d"
-version = "6.2.1+6"
-
-[[deps.GPUArrays]]
-deps = ["Adapt", "GPUArraysCore", "LLVM", "LinearAlgebra", "Printf", "Random", "Reexport", "Serialization", "Statistics"]
-git-tree-sha1 = "38cb19b8a3e600e509dc36a6396ac74266d108c1"
-uuid = "0c68f7d7-f131-5f86-a1c3-88cf8149b2d7"
-version = "10.1.1"
-
-[[deps.GPUArraysCore]]
-deps = ["Adapt"]
-git-tree-sha1 = "ec632f177c0d990e64d955ccc1b8c04c485a0950"
-uuid = "46192b85-c4d5-4398-a991-12ede77f4527"
-version = "0.1.6"
-
-[[deps.GPUCompiler]]
-deps = ["ExprTools", "InteractiveUtils", "LLVM", "Libdl", "Logging", "Scratch", "TimerOutputs", "UUIDs"]
-git-tree-sha1 = "1600477fba37c9fc067b9be21f5e8101f24a8865"
-uuid = "61eb1bfa-7361-4325-ad38-22787b887f55"
-version = "0.26.4"
-
-[[deps.Glob]]
-git-tree-sha1 = "97285bbd5230dd766e9ef6749b80fc617126d496"
-uuid = "c27321d9-0574-5035-807b-f59d2c89b15c"
-version = "1.3.1"
-
-[[deps.GnuTLS_jll]]
-deps = ["Artifacts", "GMP_jll", "JLLWrappers", "Libdl", "Nettle_jll", "P11Kit_jll", "Zlib_jll"]
-git-tree-sha1 = "383db7d3f900f4c1f47a8a04115b053c095e48d3"
-uuid = "0951126a-58fd-58f1-b5b3-b08c7c4a876d"
-version = "3.8.4+0"
-
-[[deps.HDF5_jll]]
-deps = ["Artifacts", "CompilerSupportLibraries_jll", "JLLWrappers", "LLVMOpenMP_jll", "LazyArtifacts", "LibCURL_jll", "Libdl", "MPICH_jll", "MPIPreferences", "MPItrampoline_jll", "MicrosoftMPI_jll", "OpenMPI_jll", "OpenSSL_jll", "TOML", "Zlib_jll", "libaec_jll"]
-git-tree-sha1 = "38c8874692d48d5440d5752d6c74b0c6b0b60739"
-uuid = "0234f1f7-429e-5d53-9886-15a909be8d59"
-version = "1.14.2+1"
-
-[[deps.Hwloc_jll]]
-deps = ["Artifacts", "JLLWrappers", "Libdl"]
-git-tree-sha1 = "ca0f6bf568b4bfc807e7537f081c81e35ceca114"
-uuid = "e33a78d0-f292-5ffc-b300-72abe9b543c8"
-version = "2.10.0+0"
-
-[[deps.IfElse]]
-git-tree-sha1 = "debdd00ffef04665ccbb3e150747a77560e8fad1"
-uuid = "615f187c-cbe4-4ef1-ba3b-2fcf58d6d173"
-version = "0.1.1"
-
-[[deps.IncompleteLU]]
-deps = ["LinearAlgebra", "SparseArrays"]
-git-tree-sha1 = "6c676e79f98abb6d33fa28122cad099f1e464afe"
-uuid = "40713840-3770-5561-ab4c-a76e7d0d7895"
-version = "0.2.1"
-
-[[deps.InlineStrings]]
-deps = ["Parsers"]
-git-tree-sha1 = "9cc2baf75c6d09f9da536ddf58eb2f29dedaf461"
-uuid = "842dd82b-1e85-43dc-bf29-5d0ee9dffc48"
-version = "1.4.0"
-
-[[deps.IntelOpenMP_jll]]
-deps = ["Artifacts", "JLLWrappers", "Libdl"]
-git-tree-sha1 = "be50fe8df3acbffa0274a744f1a99d29c45a57f4"
-uuid = "1d5cc7b8-4909-519e-a0f8-d0f5ad9712d0"
-version = "2024.1.0+0"
-
-[[deps.InteractiveUtils]]
-deps = ["Markdown"]
-uuid = "b77e0a4c-d291-57a0-90e8-8db25a27a240"
-
-[[deps.Interpolations]]
-deps = ["Adapt", "AxisAlgorithms", "ChainRulesCore", "LinearAlgebra", "OffsetArrays", "Random", "Ratios", "Requires", "SharedArrays", "SparseArrays", "StaticArrays", "WoodburyMatrices"]
-git-tree-sha1 = "88a101217d7cb38a7b481ccd50d21876e1d1b0e0"
-uuid = "a98d9a8b-a2ab-59e6-89dd-64a1c18fca59"
-version = "0.15.1"
-
-    [deps.Interpolations.extensions]
-    InterpolationsUnitfulExt = "Unitful"
-
-    [deps.Interpolations.weakdeps]
-    Unitful = "1986cc42-f94f-5a68-af5c-568840ba703d"
-
-[[deps.InvertedIndices]]
-git-tree-sha1 = "0dc7b50b8d436461be01300fd8cd45aa0274b038"
-uuid = "41ab1584-1d38-5bbf-9106-f11c6c58b48f"
-version = "1.3.0"
-
-[[deps.IrrationalConstants]]
-git-tree-sha1 = "630b497eafcc20001bba38a4651b327dcfc491d2"
-uuid = "92d709cd-6900-40b7-9082-c6be49f344b6"
-version = "0.2.2"
-
-[[deps.IterativeSolvers]]
-deps = ["LinearAlgebra", "Printf", "Random", "RecipesBase", "SparseArrays"]
-git-tree-sha1 = "59545b0a2b27208b0650df0a46b8e3019f85055b"
-uuid = "42fd0dbc-a981-5370-80f2-aaf504508153"
-version = "0.9.4"
-
-[[deps.IteratorInterfaceExtensions]]
-git-tree-sha1 = "a3f24677c21f5bbe9d2a714f95dcd58337fb2856"
-uuid = "82899510-4779-5014-852e-03e436cf321d"
-version = "1.0.0"
-
-[[deps.JLD2]]
-deps = ["FileIO", "MacroTools", "Mmap", "OrderedCollections", "Pkg", "PrecompileTools", "Reexport", "Requires", "TranscodingStreams", "UUIDs", "Unicode"]
-git-tree-sha1 = "bdbe8222d2f5703ad6a7019277d149ec6d78c301"
-uuid = "033835bb-8acc-5ee8-8aae-3f567f8a3819"
-version = "0.4.48"
-
-[[deps.JLLWrappers]]
-deps = ["Artifacts", "Preferences"]
-git-tree-sha1 = "7e5d6779a1e09a36db2a7b6cff50942a0a7d0fca"
-uuid = "692b3bcd-3c85-4b1f-b108-f13ce0eb3210"
-version = "1.5.0"
-
-[[deps.JSON3]]
-deps = ["Dates", "Mmap", "Parsers", "PrecompileTools", "StructTypes", "UUIDs"]
-git-tree-sha1 = "eb3edce0ed4fa32f75a0a11217433c31d56bd48b"
-uuid = "0f8b85d8-7281-11e9-16c2-39a750bddbf1"
-version = "1.14.0"
-
-    [deps.JSON3.extensions]
-    JSON3ArrowExt = ["ArrowTypes"]
-
-    [deps.JSON3.weakdeps]
-    ArrowTypes = "31f734f8-188a-4ce0-8406-c8a06bd891cd"
-
-[[deps.JuliaNVTXCallbacks_jll]]
-deps = ["Artifacts", "JLLWrappers", "Libdl", "Pkg"]
-git-tree-sha1 = "af433a10f3942e882d3c671aacb203e006a5808f"
-uuid = "9c1d0b0a-7046-5b2e-a33f-ea22f176ac7e"
-version = "0.2.1+0"
-
-[[deps.KernelAbstractions]]
-deps = ["Adapt", "Atomix", "InteractiveUtils", "LinearAlgebra", "MacroTools", "PrecompileTools", "Requires", "SparseArrays", "StaticArrays", "UUIDs", "UnsafeAtomics", "UnsafeAtomicsLLVM"]
-git-tree-sha1 = "db02395e4c374030c53dc28f3c1d33dec35f7272"
-uuid = "63c18a36-062a-441e-b654-da1e3ab1ce7c"
-version = "0.9.19"
-
-    [deps.KernelAbstractions.extensions]
-    EnzymeExt = "EnzymeCore"
-
-    [deps.KernelAbstractions.weakdeps]
-    EnzymeCore = "f151be2c-9106-41f4-ab19-57ee4f262869"
-
-[[deps.LLVM]]
-deps = ["CEnum", "LLVMExtra_jll", "Libdl", "Preferences", "Printf", "Requires", "Unicode"]
-git-tree-sha1 = "839c82932db86740ae729779e610f07a1640be9a"
-uuid = "929cbde3-209d-540e-8aea-75f648917ca0"
-version = "6.6.3"
-weakdeps = ["BFloat16s"]
-
-    [deps.LLVM.extensions]
-    BFloat16sExt = "BFloat16s"
-
-[[deps.LLVMExtra_jll]]
-deps = ["Artifacts", "JLLWrappers", "LazyArtifacts", "Libdl", "TOML"]
-git-tree-sha1 = "88b916503aac4fb7f701bb625cd84ca5dd1677bc"
-uuid = "dad2f222-ce93-54a1-a47d-0025e8a3acab"
-version = "0.0.29+0"
-
-[[deps.LLVMLoopInfo]]
-git-tree-sha1 = "2e5c102cfc41f48ae4740c7eca7743cc7e7b75ea"
-uuid = "8b046642-f1f6-4319-8d3c-209ddc03c586"
-version = "1.0.0"
-
-[[deps.LLVMOpenMP_jll]]
-deps = ["Artifacts", "JLLWrappers", "Libdl"]
-git-tree-sha1 = "d986ce2d884d49126836ea94ed5bfb0f12679713"
-uuid = "1d63c593-3942-5779-bab2-d838dc0a180e"
-version = "15.0.7+0"
-
-[[deps.LRUCache]]
-git-tree-sha1 = "b3cc6698599b10e652832c2f23db3cab99d51b59"
-uuid = "8ac3fa9e-de4c-5943-b1dc-09c6b5f20637"
-version = "1.6.1"
-weakdeps = ["Serialization"]
-
-    [deps.LRUCache.extensions]
-    SerializationExt = ["Serialization"]
-
-[[deps.LaTeXStrings]]
-git-tree-sha1 = "50901ebc375ed41dbf8058da26f9de442febbbec"
-uuid = "b964fa9f-0449-5b57-a5c2-d3ea65f4040f"
-version = "1.3.1"
-
-[[deps.LazyArtifacts]]
-deps = ["Artifacts", "Pkg"]
-uuid = "4af54fe1-eca0-43a8-85a7-787d91b784e3"
-
-[[deps.LibCURL]]
-deps = ["LibCURL_jll", "MozillaCACerts_jll"]
-uuid = "b27032c2-a3e7-50c8-80cd-2d36dbcbfd21"
-version = "0.6.4"
-
-[[deps.LibCURL_jll]]
-deps = ["Artifacts", "LibSSH2_jll", "Libdl", "MbedTLS_jll", "Zlib_jll", "nghttp2_jll"]
-uuid = "deac9b47-8bc7-5906-a0fe-35ac56dc84c0"
-version = "8.4.0+0"
-
-[[deps.LibGit2]]
-deps = ["Base64", "LibGit2_jll", "NetworkOptions", "Printf", "SHA"]
-uuid = "76f85450-5226-5b5a-8eaa-529ad045b433"
-
-[[deps.LibGit2_jll]]
-deps = ["Artifacts", "LibSSH2_jll", "Libdl", "MbedTLS_jll"]
-uuid = "e37daf67-58a4-590a-8e99-b0245dd2ffc5"
-version = "1.6.4+0"
-
-[[deps.LibSSH2_jll]]
-deps = ["Artifacts", "Libdl", "MbedTLS_jll"]
-uuid = "29816b5a-b9ab-546f-933c-edad1886dfa8"
-version = "1.11.0+1"
-
-[[deps.Libdl]]
-uuid = "8f399da3-3557-5675-b5ff-fb832c97cbdb"
-
-[[deps.Libiconv_jll]]
-deps = ["Artifacts", "JLLWrappers", "Libdl"]
-git-tree-sha1 = "f9557a255370125b405568f9767d6d195822a175"
-uuid = "94ce4f54-9a6c-5748-9c1c-f9c7231a4531"
-version = "1.17.0+0"
-
-[[deps.LineSearches]]
-deps = ["LinearAlgebra", "NLSolversBase", "NaNMath", "Parameters", "Printf"]
-git-tree-sha1 = "7bbea35cec17305fc70a0e5b4641477dc0789d9d"
-uuid = "d3d80556-e9d4-5f37-9878-2ab0fcc64255"
-version = "7.2.0"
-
-[[deps.LinearAlgebra]]
-deps = ["Libdl", "OpenBLAS_jll", "libblastrampoline_jll"]
-uuid = "37e2e46d-f89d-539d-b4ee-838fcccc9c8e"
-
-[[deps.LogExpFunctions]]
-deps = ["DocStringExtensions", "IrrationalConstants", "LinearAlgebra"]
-git-tree-sha1 = "18144f3e9cbe9b15b070288eef858f71b291ce37"
-uuid = "2ab3a3ac-af41-5b50-aa03-7779005ae688"
-version = "0.3.27"
-
-    [deps.LogExpFunctions.extensions]
-    LogExpFunctionsChainRulesCoreExt = "ChainRulesCore"
-    LogExpFunctionsChangesOfVariablesExt = "ChangesOfVariables"
-    LogExpFunctionsInverseFunctionsExt = "InverseFunctions"
-
-    [deps.LogExpFunctions.weakdeps]
-    ChainRulesCore = "d360d2e6-b24c-11e9-a2a3-2a2ae2dbcce4"
-    ChangesOfVariables = "9e997f8a-9a97-42d5-a9f1-ce6bfc15e2c0"
-    InverseFunctions = "3587e190-3f89-42d0-90ee-14403ec27112"
-
-[[deps.Logging]]
-uuid = "56ddb016-857b-54e1-b83d-db4d58db5568"
-
-[[deps.Lz4_jll]]
-deps = ["Artifacts", "JLLWrappers", "Libdl"]
-git-tree-sha1 = "6c26c5e8a4203d43b5497be3ec5d4e0c3cde240a"
-uuid = "5ced341a-0733-55b8-9ab6-a4889d929147"
-version = "1.9.4+0"
-
-[[deps.MKL_jll]]
-deps = ["Artifacts", "IntelOpenMP_jll", "JLLWrappers", "LazyArtifacts", "Libdl", "oneTBB_jll"]
-git-tree-sha1 = "80b2833b56d466b3858d565adcd16a4a05f2089b"
-uuid = "856f044c-d86e-5d09-b602-aeab76dc8ba7"
-version = "2024.1.0+0"
-
-[[deps.MPI]]
-deps = ["Distributed", "DocStringExtensions", "Libdl", "MPICH_jll", "MPIPreferences", "MPItrampoline_jll", "MicrosoftMPI_jll", "OpenMPI_jll", "PkgVersion", "PrecompileTools", "Requires", "Serialization", "Sockets"]
-git-tree-sha1 = "b4d8707e42b693720b54f0b3434abee6dd4d947a"
-uuid = "da04e1cc-30fd-572f-bb4f-1f8673147195"
-version = "0.20.16"
-
-    [deps.MPI.extensions]
-    AMDGPUExt = "AMDGPU"
-    CUDAExt = "CUDA"
-
-    [deps.MPI.weakdeps]
-    AMDGPU = "21141c5a-9bdb-4563-92ae-f87d6854732e"
-    CUDA = "052768ef-5323-5732-b1bb-66c8b64840ba"
-
-[[deps.MPICH_jll]]
-deps = ["Artifacts", "CompilerSupportLibraries_jll", "Hwloc_jll", "JLLWrappers", "LazyArtifacts", "Libdl", "MPIPreferences", "TOML"]
-git-tree-sha1 = "4099bb6809ac109bfc17d521dad33763bcf026b7"
-uuid = "7cb0a576-ebde-5e09-9194-50597f1243b4"
-version = "4.2.1+1"
-
-[[deps.MPIPreferences]]
-deps = ["Libdl", "Preferences"]
-git-tree-sha1 = "c105fe467859e7f6e9a852cb15cb4301126fac07"
-uuid = "3da0fdf6-3ccc-4f1b-acd9-58baa6c99267"
-version = "0.1.11"
-
-[[deps.MPItrampoline_jll]]
-deps = ["Artifacts", "CompilerSupportLibraries_jll", "JLLWrappers", "LazyArtifacts", "Libdl", "MPIPreferences", "TOML"]
-git-tree-sha1 = "ce0ca3dd147c43de175c5aff161315a424f4b8ac"
-uuid = "f1f71cc9-e9ae-5b93-9b94-4fe0e1ad3748"
-version = "5.3.3+1"
-
-[[deps.MacroTools]]
-deps = ["Markdown", "Random"]
-git-tree-sha1 = "2fa9ee3e63fd3a4f7a9a4f4744a52f4856de82df"
-uuid = "1914dd2f-81c6-5fcd-8719-6d5c9610ff09"
-version = "0.5.13"
-
-[[deps.Markdown]]
-deps = ["Base64"]
-uuid = "d6f4376e-aef5-505a-96c1-9c027394607a"
-
-[[deps.MbedTLS_jll]]
-deps = ["Artifacts", "Libdl"]
-uuid = "c8ffd9c3-330d-5841-b78e-0817d7145fa1"
-version = "2.28.2+1"
-
-[[deps.MicrosoftMPI_jll]]
-deps = ["Artifacts", "JLLWrappers", "Libdl", "Pkg"]
-git-tree-sha1 = "f12a29c4400ba812841c6ace3f4efbb6dbb3ba01"
-uuid = "9237b28f-5490-5468-be7b-bb81f5f5e6cf"
-version = "10.1.4+2"
-
-[[deps.Missings]]
-deps = ["DataAPI"]
-git-tree-sha1 = "ec4f7fbeab05d7747bdf98eb74d130a2a2ed298d"
-uuid = "e1d29d7a-bbdc-5cf2-9ac0-f12de2c33e28"
-version = "1.2.0"
-
-[[deps.Mmap]]
-uuid = "a63ad114-7e13-5084-954f-fe012c677804"
-
-[[deps.MozillaCACerts_jll]]
-uuid = "14a3606d-f60d-562e-9121-12d972cd8159"
-version = "2023.1.10"
-
-[[deps.NCDatasets]]
-deps = ["CFTime", "CommonDataModel", "DataStructures", "Dates", "DiskArrays", "NetCDF_jll", "NetworkOptions", "Printf"]
-git-tree-sha1 = "a640912695952b074672edb5f9aaee2f7f9fd59a"
-uuid = "85f8d34a-cbdd-5861-8df4-14fed0d494ab"
-version = "0.14.4"
-
-[[deps.NLSolversBase]]
-deps = ["DiffResults", "Distributed", "FiniteDiff", "ForwardDiff"]
-git-tree-sha1 = "a0b464d183da839699f4c79e7606d9d186ec172c"
-uuid = "d41bc354-129a-5804-8e4c-c37616107c6c"
-version = "7.8.3"
-
-[[deps.NLsolve]]
-deps = ["Distances", "LineSearches", "LinearAlgebra", "NLSolversBase", "Printf", "Reexport"]
-git-tree-sha1 = "019f12e9a1a7880459d0173c182e6a99365d7ac1"
-uuid = "2774e3e8-f4cf-5e23-947b-6d7e65073b56"
-version = "4.5.1"
-
-[[deps.NVTX]]
-deps = ["Colors", "JuliaNVTXCallbacks_jll", "Libdl", "NVTX_jll"]
-git-tree-sha1 = "53046f0483375e3ed78e49190f1154fa0a4083a1"
-uuid = "5da4648a-3479-48b8-97b9-01cb529c0a1f"
-version = "0.3.4"
-
-[[deps.NVTX_jll]]
-deps = ["Artifacts", "JLLWrappers", "Libdl", "Pkg"]
-git-tree-sha1 = "ce3269ed42816bf18d500c9f63418d4b0d9f5a3b"
-uuid = "e98f9f5b-d649-5603-91fd-7774390e6439"
-version = "3.1.0+2"
-
-[[deps.NaNMath]]
-deps = ["OpenLibm_jll"]
-git-tree-sha1 = "0877504529a3e5c3343c6f8b4c0381e57e4387e4"
-uuid = "77ba4419-2d1f-58cd-9bb1-8ffee604a2e3"
-version = "1.0.2"
-
-[[deps.NetCDF_jll]]
-deps = ["Artifacts", "Blosc_jll", "Bzip2_jll", "HDF5_jll", "JLLWrappers", "LibCURL_jll", "Libdl", "OpenMPI_jll", "XML2_jll", "Zlib_jll", "Zstd_jll", "libzip_jll"]
-git-tree-sha1 = "a8af1798e4eb9ff768ce7fdefc0e957097793f15"
-uuid = "7243133f-43d8-5620-bbf4-c2c921802cf3"
-version = "400.902.209+0"
-
-[[deps.Nettle_jll]]
-deps = ["Artifacts", "GMP_jll", "JLLWrappers", "Libdl", "Pkg"]
-git-tree-sha1 = "eca63e3847dad608cfa6a3329b95ef674c7160b4"
-uuid = "4c82536e-c426-54e4-b420-14f461c4ed8b"
-version = "3.7.2+0"
-
-[[deps.NetworkOptions]]
-uuid = "ca575930-c2e3-43a9-ace4-1e988b2c1908"
-version = "1.2.0"
-
-[[deps.Oceananigans]]
-deps = ["Adapt", "CUDA", "Crayons", "CubedSphere", "Dates", "Distances", "DocStringExtensions", "FFTW", "Glob", "IncompleteLU", "InteractiveUtils", "IterativeSolvers", "JLD2", "KernelAbstractions", "LinearAlgebra", "Logging", "MPI", "NCDatasets", "OffsetArrays", "OrderedCollections", "PencilArrays", "PencilFFTs", "Pkg", "Printf", "Random", "Rotations", "SeawaterPolynomials", "SparseArrays", "Statistics", "StructArrays"]
-git-tree-sha1 = "6aa48ccf34ac89fbd8db20cafa5f25c0d4aaa55a"
-uuid = "9e8cae18-63c1-5223-a75c-80ca9d6e9a09"
-version = "0.91.0"
-
-    [deps.Oceananigans.extensions]
-    OceananigansEnzymeExt = "Enzyme"
-
-    [deps.Oceananigans.weakdeps]
-    Enzyme = "7da242da-08ed-463a-9acd-ee780be4f1d9"
-
-[[deps.Oceanostics]]
-deps = ["DocStringExtensions", "Oceananigans", "Printf", "SeawaterPolynomials"]
-git-tree-sha1 = "bc733aea599610ae03c95d32e11f4f928784fe13"
-uuid = "d0ccf422-c8fb-49b5-a76d-74acdde946ac"
-version = "0.14.1"
-
-[[deps.OffsetArrays]]
-git-tree-sha1 = "e64b4f5ea6b7389f6f046d13d4896a8f9c1ba71e"
-uuid = "6fe1bfb0-de20-5000-8ca7-80f57d26f881"
-version = "1.14.0"
-weakdeps = ["Adapt"]
-
-    [deps.OffsetArrays.extensions]
-    OffsetArraysAdaptExt = "Adapt"
-
-[[deps.OpenBLAS_jll]]
-deps = ["Artifacts", "CompilerSupportLibraries_jll", "Libdl"]
-uuid = "4536629a-c528-5b80-bd46-f80d51c5b363"
-version = "0.3.23+2"
-
-[[deps.OpenLibm_jll]]
-deps = ["Artifacts", "Libdl"]
-uuid = "05823500-19ac-5b8b-9628-191a04bc5112"
-version = "0.8.1+2"
-
-[[deps.OpenMPI_jll]]
-deps = ["Artifacts", "CompilerSupportLibraries_jll", "Hwloc_jll", "JLLWrappers", "LazyArtifacts", "Libdl", "MPIPreferences", "PMIx_jll", "TOML", "Zlib_jll", "libevent_jll", "prrte_jll"]
-git-tree-sha1 = "f46caf663e069027a06942d00dced37f1eb3d8ad"
-uuid = "fe0851c0-eecd-5654-98d4-656369965a5c"
-version = "5.0.2+0"
-
-[[deps.OpenSSL_jll]]
-deps = ["Artifacts", "JLLWrappers", "Libdl"]
-git-tree-sha1 = "3da7367955dcc5c54c1ba4d402ccdc09a1a3e046"
-uuid = "458c3c95-2e84-50aa-8efc-19380b2a3a95"
-version = "3.0.13+1"
-
-[[deps.OpenSpecFun_jll]]
-deps = ["Artifacts", "CompilerSupportLibraries_jll", "JLLWrappers", "Libdl", "Pkg"]
-git-tree-sha1 = "13652491f6856acfd2db29360e1bbcd4565d04f1"
-uuid = "efe28fd5-8261-553b-a9e1-b2916fc3738e"
-version = "0.5.5+0"
-
-[[deps.OrderedCollections]]
-git-tree-sha1 = "dfdf5519f235516220579f949664f1bf44e741c5"
-uuid = "bac558e1-5e72-5ebc-8fee-abe8a469f55d"
-version = "1.6.3"
-
-[[deps.P11Kit_jll]]
-deps = ["Artifacts", "JLLWrappers", "Libdl", "Pkg"]
-git-tree-sha1 = "2cd396108e178f3ae8dedbd8e938a18726ab2fbf"
-uuid = "c2071276-7c44-58a7-b746-946036e04d0a"
-version = "0.24.1+0"
-
-[[deps.PMIx_jll]]
-deps = ["Artifacts", "Hwloc_jll", "JLLWrappers", "Libdl", "Zlib_jll", "libevent_jll"]
-git-tree-sha1 = "360f48126b5f2c2f0c833be960097f7c62705976"
-uuid = "32165bc3-0280-59bc-8c0b-c33b6203efab"
-version = "4.2.9+0"
-
-[[deps.PackageExtensionCompat]]
-git-tree-sha1 = "fb28e33b8a95c4cee25ce296c817d89cc2e53518"
-uuid = "65ce6f38-6b18-4e1d-a461-8949797d7930"
-version = "1.0.2"
-weakdeps = ["Requires", "TOML"]
-
-[[deps.Parameters]]
-deps = ["OrderedCollections", "UnPack"]
-git-tree-sha1 = "34c0e9ad262e5f7fc75b10a9952ca7692cfc5fbe"
-uuid = "d96e819e-fc66-5662-9728-84c9c7592b0a"
-version = "0.12.3"
-
-[[deps.Parsers]]
-deps = ["Dates", "PrecompileTools", "UUIDs"]
-git-tree-sha1 = "8489905bcdbcfac64d1daa51ca07c0d8f0283821"
-uuid = "69de0a69-1ddd-5017-9359-2bf0b02dc9f0"
-version = "2.8.1"
-
-[[deps.PencilArrays]]
-deps = ["Adapt", "JSON3", "LinearAlgebra", "MPI", "OffsetArrays", "Random", "Reexport", "StaticArrayInterface", "StaticArrays", "StaticPermutations", "Strided", "TimerOutputs", "VersionParsing"]
-git-tree-sha1 = "6510e851700a851944f7ffa5cd990cced4802ad2"
-uuid = "0e08944d-e94e-41b1-9406-dcf66b6a9d2e"
-version = "0.19.3"
-
-    [deps.PencilArrays.extensions]
-    PencilArraysDiffEqExt = ["DiffEqBase"]
-    PencilArraysHDF5Ext = ["HDF5"]
-
-    [deps.PencilArrays.weakdeps]
-    DiffEqBase = "2b5f629d-d688-5b77-993f-72d75c75574e"
-    HDF5 = "f67ccb44-e63f-5c2f-98bd-6dc0ccc4ba2f"
-
-[[deps.PencilFFTs]]
-deps = ["AbstractFFTs", "FFTW", "LinearAlgebra", "MPI", "PencilArrays", "Reexport", "TimerOutputs"]
-git-tree-sha1 = "bd69f3f0ee248cfb4241800aefb705b5ded592ff"
-uuid = "4a48f351-57a6-4416-9ec4-c37015456aae"
-version = "0.15.1"
-
-[[deps.Pkg]]
-deps = ["Artifacts", "Dates", "Downloads", "FileWatching", "LibGit2", "Libdl", "Logging", "Markdown", "Printf", "REPL", "Random", "SHA", "Serialization", "TOML", "Tar", "UUIDs", "p7zip_jll"]
-uuid = "44cfe95a-1eb2-52ea-b672-e2afdf69b78f"
-version = "1.10.0"
-
-[[deps.PkgVersion]]
-deps = ["Pkg"]
-git-tree-sha1 = "f9501cc0430a26bc3d156ae1b5b0c1b47af4d6da"
-uuid = "eebad327-c553-4316-9ea0-9fa01ccd7688"
-version = "0.3.3"
-
-[[deps.Polynomials]]
-deps = ["LinearAlgebra", "RecipesBase", "Setfield", "SparseArrays"]
-git-tree-sha1 = "89620a0b5458dca4bf9ea96174fa6422b3adf6f9"
-uuid = "f27b6e38-b328-58d1-80ce-0feddd5e7a45"
-version = "4.0.8"
-
-    [deps.Polynomials.extensions]
-    PolynomialsChainRulesCoreExt = "ChainRulesCore"
-    PolynomialsFFTWExt = "FFTW"
-    PolynomialsMakieCoreExt = "MakieCore"
-    PolynomialsMutableArithmeticsExt = "MutableArithmetics"
-
-    [deps.Polynomials.weakdeps]
-    ChainRulesCore = "d360d2e6-b24c-11e9-a2a3-2a2ae2dbcce4"
-    FFTW = "7a1cc6ca-52ef-59f5-83cd-3a7055c09341"
-    MakieCore = "20f20a25-4f0e-4fdf-b5d1-57303727442b"
-    MutableArithmetics = "d8a4904e-b15c-11e9-3269-09a3773c0cb0"
-
-[[deps.PooledArrays]]
-deps = ["DataAPI", "Future"]
-git-tree-sha1 = "36d8b4b899628fb92c2749eb488d884a926614d3"
-uuid = "2dfb63ee-cc39-5dd5-95bd-886bf059d720"
-version = "1.4.3"
-
-[[deps.PrecompileTools]]
-deps = ["Preferences"]
-git-tree-sha1 = "5aa36f7049a63a1528fe8f7c3f2113413ffd4e1f"
-uuid = "aea7be01-6a6a-4083-8856-8a6e6704d82a"
-version = "1.2.1"
-
-[[deps.Preferences]]
-deps = ["TOML"]
-git-tree-sha1 = "9306f6085165d270f7e3db02af26a400d580f5c6"
-uuid = "21216c6a-2e73-6563-6e65-726566657250"
-version = "1.4.3"
-
-[[deps.PrettyTables]]
-deps = ["Crayons", "LaTeXStrings", "Markdown", "PrecompileTools", "Printf", "Reexport", "StringManipulation", "Tables"]
-git-tree-sha1 = "88b895d13d53b5577fd53379d913b9ab9ac82660"
-uuid = "08abe8d2-0d0c-5749-adfa-8a2ac140af0d"
-version = "2.3.1"
-
-[[deps.Printf]]
-deps = ["Unicode"]
-uuid = "de0858da-6303-5e67-8744-51eddeeeb8d7"
-
-[[deps.ProgressBars]]
-deps = ["Printf"]
-git-tree-sha1 = "b437cdb0385ed38312d91d9c00c20f3798b30256"
-uuid = "49802e3a-d2f1-5c88-81d8-b72133a6f568"
-version = "1.5.1"
-
-[[deps.QuadGK]]
-deps = ["DataStructures", "LinearAlgebra"]
-git-tree-sha1 = "9b23c31e76e333e6fb4c1595ae6afa74966a729e"
-uuid = "1fd47b50-473d-5c70-9696-f719f8f3bcdc"
-version = "2.9.4"
-
-[[deps.Quaternions]]
-deps = ["LinearAlgebra", "Random", "RealDot"]
-git-tree-sha1 = "994cc27cdacca10e68feb291673ec3a76aa2fae9"
-uuid = "94ee1d12-ae83-5a48-8b1c-48b8ff168ae0"
-version = "0.7.6"
-
-[[deps.REPL]]
-deps = ["InteractiveUtils", "Markdown", "Sockets", "Unicode"]
-uuid = "3fa0cd96-eef1-5676-8a61-b3b8758bbffb"
-
-[[deps.Random]]
-deps = ["SHA"]
-uuid = "9a3f8284-a2c9-5f02-9a11-845980a1fd5c"
-
-[[deps.Random123]]
-deps = ["Random", "RandomNumbers"]
-git-tree-sha1 = "4743b43e5a9c4a2ede372de7061eed81795b12e7"
-uuid = "74087812-796a-5b5d-8853-05524746bad3"
-version = "1.7.0"
-
-[[deps.RandomNumbers]]
-deps = ["Random", "Requires"]
-git-tree-sha1 = "043da614cc7e95c703498a491e2c21f58a2b8111"
-uuid = "e6cf234a-135c-5ec9-84dd-332b85af5143"
-version = "1.5.3"
-
-[[deps.Ratios]]
-deps = ["Requires"]
-git-tree-sha1 = "1342a47bf3260ee108163042310d26f2be5ec90b"
-uuid = "c84ed2f1-dad5-54f0-aa8e-dbefe2724439"
-version = "0.4.5"
-weakdeps = ["FixedPointNumbers"]
-
-    [deps.Ratios.extensions]
-    RatiosFixedPointNumbersExt = "FixedPointNumbers"
-
-[[deps.RealDot]]
-deps = ["LinearAlgebra"]
-git-tree-sha1 = "9f0a1b71baaf7650f4fa8a1d168c7fb6ee41f0c9"
-uuid = "c1ae055f-0cd5-4b69-90a6-9a35b1a98df9"
-version = "0.1.0"
-
-[[deps.RecipesBase]]
-deps = ["PrecompileTools"]
-git-tree-sha1 = "5c3d09cc4f31f5fc6af001c250bf1278733100ff"
-uuid = "3cdcf5f2-1ef4-517c-9805-6587b60abb01"
-version = "1.3.4"
-
-[[deps.Reexport]]
-git-tree-sha1 = "45e428421666073eab6f2da5c9d310d99bb12f9b"
-uuid = "189a3867-3050-52da-a836-e630ba90ab69"
-version = "1.2.2"
-
-[[deps.Requires]]
-deps = ["UUIDs"]
-git-tree-sha1 = "838a3a4188e2ded87a4f9f184b4b0d78a1e91cb7"
-uuid = "ae029012-a4dd-5104-9daa-d747884805df"
-version = "1.3.0"
-
-[[deps.Rotations]]
-deps = ["LinearAlgebra", "Quaternions", "Random", "StaticArrays"]
-git-tree-sha1 = "5680a9276685d392c87407df00d57c9924d9f11e"
-uuid = "6038ab10-8711-5258-84ad-4b1120ba62dc"
-version = "1.7.1"
-weakdeps = ["RecipesBase"]
-
-    [deps.Rotations.extensions]
-    RotationsRecipesBaseExt = "RecipesBase"
-
-[[deps.SHA]]
-uuid = "ea8e919c-243c-51af-8825-aaa63cd721ce"
-version = "0.7.0"
-
-[[deps.Scratch]]
-deps = ["Dates"]
-git-tree-sha1 = "3bac05bc7e74a75fd9cba4295cde4045d9fe2386"
-uuid = "6c6a2e73-6563-6170-7368-637461726353"
-version = "1.2.1"
-
-[[deps.SeawaterPolynomials]]
-git-tree-sha1 = "6d85acd6de472f8e6da81c61c7c5b6280a55e0bc"
-uuid = "d496a93d-167e-4197-9f49-d3af4ff8fe40"
-version = "0.3.4"
-
-[[deps.SentinelArrays]]
-deps = ["Dates", "Random"]
-git-tree-sha1 = "90b4f68892337554d31cdcdbe19e48989f26c7e6"
-uuid = "91c51154-3ec4-41a3-a24f-3f23e20d615c"
-version = "1.4.3"
-
-[[deps.Serialization]]
-uuid = "9e88b42a-f829-5b0c-bbe9-9e923198166b"
-
-[[deps.Setfield]]
-deps = ["ConstructionBase", "Future", "MacroTools", "StaticArraysCore"]
-git-tree-sha1 = "e2cc6d8c88613c05e1defb55170bf5ff211fbeac"
-uuid = "efcf1570-3423-57d1-acb7-fd33fddbac46"
-version = "1.1.1"
-
-[[deps.SharedArrays]]
-deps = ["Distributed", "Mmap", "Random", "Serialization"]
-uuid = "1a1011a3-84de-559e-8e89-a11a2f7dc383"
-
-[[deps.Sockets]]
-uuid = "6462fe0b-24de-5631-8697-dd941f90decc"
-
-[[deps.SortingAlgorithms]]
-deps = ["DataStructures"]
-git-tree-sha1 = "66e0a8e672a0bdfca2c3f5937efb8538b9ddc085"
-uuid = "a2af1166-a08f-5f64-846c-94a0d3cef48c"
-version = "1.2.1"
-
-[[deps.SparseArrays]]
-deps = ["Libdl", "LinearAlgebra", "Random", "Serialization", "SuiteSparse_jll"]
-uuid = "2f01184e-e22b-5df5-ae63-d93ebab69eaf"
-version = "1.10.0"
-
-[[deps.SpecialFunctions]]
-deps = ["IrrationalConstants", "LogExpFunctions", "OpenLibm_jll", "OpenSpecFun_jll"]
-git-tree-sha1 = "2f5d4697f21388cbe1ff299430dd169ef97d7e14"
-uuid = "276daf66-3868-5448-9aa4-cd146d93841b"
-version = "2.4.0"
-weakdeps = ["ChainRulesCore"]
-
-    [deps.SpecialFunctions.extensions]
-    SpecialFunctionsChainRulesCoreExt = "ChainRulesCore"
-
-[[deps.Static]]
-deps = ["IfElse"]
-git-tree-sha1 = "d2fdac9ff3906e27f7a618d47b676941baa6c80c"
-uuid = "aedffcd0-7271-4cad-89d0-dc628f76c6d3"
-version = "0.8.10"
-
-[[deps.StaticArrayInterface]]
-deps = ["ArrayInterface", "Compat", "IfElse", "LinearAlgebra", "PrecompileTools", "Requires", "SparseArrays", "Static", "SuiteSparse"]
-git-tree-sha1 = "5d66818a39bb04bf328e92bc933ec5b4ee88e436"
-uuid = "0d7ed370-da01-4f52-bd93-41d350b8b718"
-version = "1.5.0"
-weakdeps = ["OffsetArrays", "StaticArrays"]
-
-    [deps.StaticArrayInterface.extensions]
-    StaticArrayInterfaceOffsetArraysExt = "OffsetArrays"
-    StaticArrayInterfaceStaticArraysExt = "StaticArrays"
-
-[[deps.StaticArrays]]
-deps = ["LinearAlgebra", "PrecompileTools", "Random", "StaticArraysCore"]
-git-tree-sha1 = "9ae599cd7529cfce7fea36cf00a62cfc56f0f37c"
-uuid = "90137ffa-7385-5640-81b9-e52037218182"
-version = "1.9.4"
-weakdeps = ["ChainRulesCore", "Statistics"]
-
-    [deps.StaticArrays.extensions]
-    StaticArraysChainRulesCoreExt = "ChainRulesCore"
-    StaticArraysStatisticsExt = "Statistics"
-
-[[deps.StaticArraysCore]]
-git-tree-sha1 = "36b3d696ce6366023a0ea192b4cd442268995a0d"
-uuid = "1e83bf80-4336-4d27-bf5d-d5a4f845583c"
-version = "1.4.2"
-
-[[deps.StaticPermutations]]
-git-tree-sha1 = "193c3daa18ff3e55c1dae66acb6a762c4a3bdb0b"
-uuid = "15972242-4b8f-49a0-b8a1-9ac0e7a1a45d"
-version = "0.3.0"
-
-[[deps.Statistics]]
-deps = ["LinearAlgebra", "SparseArrays"]
-uuid = "10745b16-79ce-11e8-11f9-7d13ad32a3b2"
-version = "1.10.0"
-
-[[deps.StatsAPI]]
-deps = ["LinearAlgebra"]
-git-tree-sha1 = "1ff449ad350c9c4cbc756624d6f8a8c3ef56d3ed"
-uuid = "82ae8749-77ed-4fe6-ae5f-f523153014b0"
-version = "1.7.0"
-
-[[deps.Strided]]
-deps = ["LinearAlgebra", "StridedViews", "TupleTools"]
-git-tree-sha1 = "40c69be0e1b72ee2f42923b7d1ff13e0b04e675c"
-uuid = "5e0ebb24-38b0-5f93-81fe-25c709ecae67"
-version = "2.0.4"
-
-[[deps.StridedViews]]
-deps = ["LinearAlgebra", "PackageExtensionCompat"]
-git-tree-sha1 = "5b765c4e401693ab08981989f74a36a010aa1d8e"
-uuid = "4db3bf67-4bd7-4b4e-b153-31dc3fb37143"
-version = "0.2.2"
-weakdeps = ["CUDA"]
-
-    [deps.StridedViews.extensions]
-    StridedViewsCUDAExt = "CUDA"
-
-[[deps.StringManipulation]]
-deps = ["PrecompileTools"]
-git-tree-sha1 = "a04cabe79c5f01f4d723cc6704070ada0b9d46d5"
-uuid = "892a3eda-7b42-436c-8928-eab12a02cf0e"
-version = "0.3.4"
-
-[[deps.StructArrays]]
-deps = ["ConstructionBase", "DataAPI", "Tables"]
-git-tree-sha1 = "f4dc295e983502292c4c3f951dbb4e985e35b3be"
-uuid = "09ab397b-f2b6-538f-b94a-2f83cf4a842a"
-version = "0.6.18"
-weakdeps = ["Adapt", "GPUArraysCore", "SparseArrays", "StaticArrays"]
-
-    [deps.StructArrays.extensions]
-    StructArraysAdaptExt = "Adapt"
-    StructArraysGPUArraysCoreExt = "GPUArraysCore"
-    StructArraysSparseArraysExt = "SparseArrays"
-    StructArraysStaticArraysExt = "StaticArrays"
-
-[[deps.StructTypes]]
-deps = ["Dates", "UUIDs"]
-git-tree-sha1 = "ca4bccb03acf9faaf4137a9abc1881ed1841aa70"
-uuid = "856f2bd8-1eba-4b0a-8007-ebc267875bd4"
-version = "1.10.0"
-
-[[deps.SuiteSparse]]
-deps = ["Libdl", "LinearAlgebra", "Serialization", "SparseArrays"]
-uuid = "4607b0f0-06f3-5cda-b6b1-a6196a1729e9"
-
-[[deps.SuiteSparse_jll]]
-deps = ["Artifacts", "Libdl", "libblastrampoline_jll"]
-uuid = "bea87d4a-7f5b-5778-9afe-8cc45184846c"
-version = "7.2.1+1"
-
-[[deps.TOML]]
-deps = ["Dates"]
-uuid = "fa267f1f-6049-4f14-aa54-33bafae1ed76"
-version = "1.0.3"
-
-[[deps.TableTraits]]
-deps = ["IteratorInterfaceExtensions"]
-git-tree-sha1 = "c06b2f539df1c6efa794486abfb6ed2022561a39"
-uuid = "3783bdb8-4a98-5b6b-af9a-565f29a5fe9c"
-version = "1.0.1"
-
-[[deps.Tables]]
-deps = ["DataAPI", "DataValueInterfaces", "IteratorInterfaceExtensions", "LinearAlgebra", "OrderedCollections", "TableTraits"]
-git-tree-sha1 = "cb76cf677714c095e535e3501ac7954732aeea2d"
-uuid = "bd369af6-aec1-5ad0-b16a-f7cc5008161c"
-version = "1.11.1"
-
-[[deps.Tar]]
-deps = ["ArgTools", "SHA"]
-uuid = "a4e569a6-e804-4fa4-b0f3-eef7a1d5b13e"
-version = "1.10.0"
-
-[[deps.TaylorSeries]]
-deps = ["LinearAlgebra", "Markdown", "Requires", "SparseArrays"]
-git-tree-sha1 = "1c7170668366821b0c4c4fe03ee78f8d6cf36e2c"
-uuid = "6aa5eb33-94cf-58f4-a9d0-e4b2c4fc25ea"
-version = "0.16.0"
-
-    [deps.TaylorSeries.extensions]
-    TaylorSeriesIAExt = "IntervalArithmetic"
-
-    [deps.TaylorSeries.weakdeps]
-    IntervalArithmetic = "d1acc4aa-44c8-5952-acd4-ba5d80a2a253"
-
-[[deps.Test]]
-deps = ["InteractiveUtils", "Logging", "Random", "Serialization"]
-uuid = "8dfed614-e22c-5e08-85e1-65c5234f0b40"
-
-[[deps.TextWrap]]
-git-tree-sha1 = "43044b737fa70bc12f6105061d3da38f881a3e3c"
-uuid = "b718987f-49a8-5099-9789-dcd902bef87d"
-version = "1.0.2"
-
-[[deps.TimerOutputs]]
-deps = ["ExprTools", "Printf"]
-git-tree-sha1 = "5a13ae8a41237cff5ecf34f73eb1b8f42fff6531"
-uuid = "a759f4b9-e2f1-59dc-863e-4aeb61b1ea8f"
-version = "0.5.24"
-
-[[deps.TranscodingStreams]]
-git-tree-sha1 = "5d54d076465da49d6746c647022f3b3674e64156"
-uuid = "3bb67fe8-82b1-5028-8e26-92a6c54297fa"
-version = "0.10.8"
-weakdeps = ["Random", "Test"]
-
-    [deps.TranscodingStreams.extensions]
-    TestExt = ["Test", "Random"]
-
-[[deps.TupleTools]]
-git-tree-sha1 = "41d61b1c545b06279871ef1a4b5fcb2cac2191cd"
-uuid = "9d95972d-f1c8-5527-a6e0-b4b365fa01f6"
-version = "1.5.0"
-
-[[deps.UUIDs]]
-deps = ["Random", "SHA"]
-uuid = "cf7118a7-6976-5b1a-9a39-7adc72f591a4"
-
-[[deps.UnPack]]
-git-tree-sha1 = "387c1f73762231e86e0c9c5443ce3b4a0a9a0c2b"
-uuid = "3a884ed6-31ef-47d7-9d2a-63182c4928ed"
-version = "1.0.2"
-
-[[deps.Unicode]]
-uuid = "4ec0a83e-493e-50e2-b9ac-8f72acf5a8f5"
-
-[[deps.UnsafeAtomics]]
-git-tree-sha1 = "6331ac3440856ea1988316b46045303bef658278"
-uuid = "013be700-e6cd-48c3-b4a1-df204f14c38f"
-version = "0.2.1"
-
-[[deps.UnsafeAtomicsLLVM]]
-deps = ["LLVM", "UnsafeAtomics"]
-git-tree-sha1 = "d9f5962fecd5ccece07db1ff006fb0b5271bdfdd"
-uuid = "d80eeb9a-aca5-4d75-85e5-170c8b632249"
-version = "0.1.4"
-
-[[deps.VersionParsing]]
-git-tree-sha1 = "58d6e80b4ee071f5efd07fda82cb9fbe17200868"
-uuid = "81def892-9a0e-5fdd-b105-ffc91e053289"
-version = "1.3.0"
-
-[[deps.WoodburyMatrices]]
-deps = ["LinearAlgebra", "SparseArrays"]
-git-tree-sha1 = "c1a7aa6219628fcd757dede0ca95e245c5cd9511"
-uuid = "efce3f68-66dc-5838-9240-27a6d6f5f9b6"
-version = "1.0.0"
-
-[[deps.XML2_jll]]
-deps = ["Artifacts", "JLLWrappers", "Libdl", "Libiconv_jll", "Zlib_jll"]
-git-tree-sha1 = "52ff2af32e591541550bd753c0da8b9bc92bb9d9"
-uuid = "02c8fc9c-b97f-50b9-bbe4-9be30ff0a78a"
-version = "2.12.7+0"
-
-[[deps.XZ_jll]]
-deps = ["Artifacts", "JLLWrappers", "Libdl"]
-git-tree-sha1 = "ac88fb95ae6447c8dda6a5503f3bafd496ae8632"
-uuid = "ffd25f8a-64ca-5728-b0f7-c24cf3aae800"
-version = "5.4.6+0"
-
-[[deps.Zlib_jll]]
-deps = ["Libdl"]
-uuid = "83775a58-1f1d-513f-b197-d71354ab007a"
-version = "1.2.13+1"
-
-[[deps.Zstd_jll]]
-deps = ["Artifacts", "JLLWrappers", "Libdl"]
-git-tree-sha1 = "e678132f07ddb5bfa46857f0d7620fb9be675d3b"
-uuid = "3161d3a3-bdf6-5164-811a-617609db77b4"
-version = "1.5.6+0"
-
-[[deps.libaec_jll]]
-deps = ["Artifacts", "JLLWrappers", "Libdl"]
-git-tree-sha1 = "46bf7be2917b59b761247be3f317ddf75e50e997"
-uuid = "477f73a3-ac25-53e9-8cc3-50b2fa2566f0"
-version = "1.1.2+0"
-
-[[deps.libblastrampoline_jll]]
-deps = ["Artifacts", "Libdl"]
-uuid = "8e850b90-86db-534c-a0d3-1478176c7d93"
-version = "5.8.0+1"
-
-[[deps.libevent_jll]]
-deps = ["Artifacts", "JLLWrappers", "Libdl", "OpenSSL_jll"]
-git-tree-sha1 = "f04ec6d9a186115fb38f858f05c0c4e1b7fc9dcb"
-uuid = "1080aeaf-3a6a-583e-a51c-c537b09f60ec"
-version = "2.1.13+1"
-
-[[deps.libzip_jll]]
-deps = ["Artifacts", "Bzip2_jll", "GnuTLS_jll", "JLLWrappers", "Libdl", "XZ_jll", "Zlib_jll", "Zstd_jll"]
-git-tree-sha1 = "3282b7d16ae7ac3e57ec2f3fa8fafb564d8f9f7f"
-uuid = "337d8026-41b4-5cde-a456-74a10e5b31d1"
-version = "1.10.1+0"
-
-[[deps.nghttp2_jll]]
-deps = ["Artifacts", "Libdl"]
-uuid = "8e850ede-7688-5339-a07c-302acd2aaf8d"
-version = "1.52.0+1"
-
-[[deps.oneTBB_jll]]
-deps = ["Artifacts", "JLLWrappers", "Libdl"]
-git-tree-sha1 = "7d0ea0f4895ef2f5cb83645fa689e52cb55cf493"
-uuid = "1317d2d5-d96f-522e-a858-c73665f53c3e"
-version = "2021.12.0+0"
-
-[[deps.p7zip_jll]]
-deps = ["Artifacts", "Libdl"]
-uuid = "3f19e933-33d8-53b3-aaab-bd5110c3b7a0"
-version = "17.4.0+2"
-
-[[deps.prrte_jll]]
-deps = ["Artifacts", "Hwloc_jll", "JLLWrappers", "Libdl", "PMIx_jll", "libevent_jll"]
-git-tree-sha1 = "5adb2d7a18a30280feb66cad6f1a1dfdca2dc7b0"
-uuid = "eb928a42-fffd-568d-ab9c-3f5d54fc65b9"
-version = "3.0.2+0"
diff --git a/Project.toml b/Project.toml
index 88ec8b7..c7c3482 100644
--- a/Project.toml
+++ b/Project.toml
@@ -10,11 +10,14 @@ julia = "^1.6"
 Adapt = "79e6a3ab-5dfb-504d-930d-738a2a938a0e"
 ArgParse = "c7e460c6-2fb9-53a9-8c5b-16f535851c63"
 CUDA = "052768ef-5323-5732-b1bb-66c8b64840ba"
+CairoMakie = "13f3f980-e62b-5c42-98c6-ff1f3baf88f0"
+GLMakie = "e9467ef8-e4e7-5192-8a1a-b1aee30e663a"
 Glob = "c27321d9-0574-5035-807b-f59d2c89b15c"
 Interpolations = "a98d9a8b-a2ab-59e6-89dd-64a1c18fca59"
 JLD2 = "033835bb-8acc-5ee8-8aae-3f567f8a3819"
 KernelAbstractions = "63c18a36-062a-441e-b654-da1e3ab1ce7c"
 Logging = "56ddb016-857b-54e1-b83d-db4d58db5568"
+MathTeXEngine = "0a4f8689-d25c-4efe-a92b-7142dfc1aa53"
 NCDatasets = "85f8d34a-cbdd-5861-8df4-14fed0d494ab"
 NLsolve = "2774e3e8-f4cf-5e23-947b-6d7e65073b56"
 Oceananigans = "9e8cae18-63c1-5223-a75c-80ca9d6e9a09"
diff --git a/idealized/analyze_resolution_three_layer_simulations.jl b/idealized/analyze_resolution_three_layer_simulations.jl
new file mode 100644
index 0000000..d13d0fc
--- /dev/null
+++ b/idealized/analyze_resolution_three_layer_simulations.jl
@@ -0,0 +1,153 @@
+using Oceananigans
+using CairoMakie
+using MathTeXEngine
+
+basedir = "/Users/gregorywagner/Projects/LESbrary.jl/data"
+suite = "6_hour_suite"
+resolutions = ["1m", "2m", "4m"]
+
+labels = Dict(
+    #"0.75m" => "0.75 meter",
+    "1m" => "1 meter",
+    "2m" => "2 meter",
+    "4m" => "4 meter",
+)
+
+cases = [
+    "free_convection",
+    #"weak_wind_strong_cooling",
+    #"strong_wind_weak_cooling",
+    "med_wind_med_cooling",
+    "strong_wind",
+    #"strong_wind_no_rotation",
+    #"strong_wind_weak_cooling",
+    "strong_wind_and_sunny",
+]
+
+titles = Dict(
+    "free_convection" => "Free \n convection",
+    "strong_wind" => "Strong wind",
+    "strong_wind_no_rotation" => "Strong wind \n no rotation",
+    "weak_wind_strong_cooling" => "Weak wind \n strong cooling",
+    "med_wind_med_cooling" => "Mid wind \n mid cooling",
+    "strong_wind_weak_cooling" => "Strong wind \n weak cooling",
+    "strong_wind_and_sunny" => "Strong wind \n and sunny",
+)
+
+# case = "free_convection"
+
+fonts = (; regular=texfont())
+set_theme!(Theme(fontsize=24, linewidth=3; fonts))
+fig = Figure(size=(1200, 500))
+axc = []
+
+for (c, case) in enumerate(cases)
+    colors = Makie.wong_colors()
+    colors = reverse(colors)
+
+    if c == length(cases)
+        yaxisposition = :right
+    else
+        yaxisposition = :left
+    end
+
+    x = 1:9
+    xticks = (x .* 1e-4, [string(i) for i in x])
+
+    # if c == 1 || c == 2
+    #     xticks = ([1, 2, 3, 4] .* 1e-4, ["1", "2", "3", "4"])
+    # else
+    #     xticks = ([1, 3, 5, 7] .* 1e-4, ["1", "3", "5", "7", "9"])
+    # end
+    ax = Axis(fig[2, c], xlabel="Buoyancy \n (10⁻⁴ × m s⁻²)", ylabel=L"z \, \mathrm{(m)}"; yaxisposition, xticks)
+    title = titles[case]
+    Label(fig[1, c], title, tellwidth=false)
+    push!(axc, ax)
+
+    if c != 1 && c != length(cases)
+        hideydecorations!(ax, grid=false)
+    end
+
+    if c == 1
+        hidespines!(ax, :t, :r)
+    elseif c == 4
+        hidespines!(ax, :t, :l)
+    else
+        hidespines!(ax, :t, :l, :r)
+    end
+    
+    #res = "1m"
+    for res in resolutions
+
+    #for waves in ("", "_no_stokes")
+        #filename = string(case, waves, "_instantaneous_statistics.jld2")
+        
+        #try
+            filename = string(case, "_with_tracer_instantaneous_statistics.jld2")
+            datadir = joinpath(basedir, suite, res)
+            filepath = joinpath(datadir, filename)
+            bt = FieldTimeSeries(filepath, "b")
+            bn = interior(bt[end], 1, 1, :) .- bt[1][1, 1, 1]
+            @show length(bt)
+            z = znodes(bt)
+
+            # if waves == ""
+            #     label = "with \nStokes drift"
+            # else
+            #     label = "without"
+            # end
+
+            label = labels[res]
+            if res == "1m"
+            #if waves == ""
+                linewidth = 5
+                color = (:black, 0.8)
+                linestyle = :solid
+            else
+                linewidth = 3
+                linestyle = :dash
+                color = pop!(colors)
+            end
+            lines!(ax, bn, z; label, linewidth, linestyle, color)
+        #catch
+        #end
+    end
+
+    ylims!(ax, -200, 5)
+end
+
+Nc = length(cases)
+#Legend(fig[2, Nc+1], first(axc))
+axislegend(axc[1], position=(0, 0.9))
+
+# For Stokes drift comparison
+#xlims!(axc[1], 1.7e-4, 3.7e-4)
+#xlims!(axc[2], 1.7e-4, 4.8e-4)
+#xlims!(axc[3], 1.7e-4, 6e-4)
+#xlims!(axc[4], 1.7e-4, 7.1e-4)
+
+# For 12 hour suite resolution comparison
+if suite == "6_hour_suite"
+    xlims!(axc[1], 1e-4, 3.3e-4)
+    xlims!(axc[2], 1e-4, 3.9e-4)
+    xlims!(axc[3], 1e-4, 5.5e-4)
+    xlims!(axc[4], 1e-4, 7.0e-4)
+elseif suite == "24_hour_suite"
+    xlims!(axc[1], 1.7e-4, 3.3e-4)
+    xlims!(axc[2], 1.7e-4, 3.9e-4)
+    xlims!(axc[3], 1.7e-4, 5.5e-4)
+    xlims!(axc[4], 1.7e-4, 9.3e-4)
+else
+    xlims!(axc[1], 1.7e-4, 3.3e-4)
+    xlims!(axc[2], 1.7e-4, 3.9e-4)
+    xlims!(axc[3], 1.7e-4, 5.5e-4)
+    xlims!(axc[4], 1.7e-4, 8.0e-4)
+end
+
+for ax in axc
+    ylims!(ax, -160, 0)
+end
+    
+display(fig)
+save("les_resolution_dependence_$suite.pdf", fig)
+
diff --git a/idealized/analyze_stokes_three_layer_simulations.jl b/idealized/analyze_stokes_three_layer_simulations.jl
new file mode 100644
index 0000000..41fc6b9
--- /dev/null
+++ b/idealized/analyze_stokes_three_layer_simulations.jl
@@ -0,0 +1,125 @@
+using Oceananigans
+using CairoMakie
+using MathTeXEngine
+
+basedir = "/Users/gregorywagner/Projects/LESbrary.jl/data"
+suite = "12_hour_suite"
+resolutions = ["1m", "2m", "4m"]
+
+labels = Dict(
+    "1m" => "1 meter",
+    "2m" => "2 meter",
+    "4m" => "4 meter",
+)
+
+cases = [
+    #"free_convection",
+    "weak_wind_strong_cooling",
+    "med_wind_med_cooling",
+    "strong_wind_weak_cooling",
+    "strong_wind",
+    #"strong_wind_no_rotation",
+    #"strong_wind_weak_cooling",
+]
+
+titles = Dict(
+    "free_convection" => "Free \n convection",
+    "strong_wind" => "Strong wind",
+    "strong_wind_no_rotation" => "Strong wind \n no rotation",
+    "weak_wind_strong_cooling" => "Weak wind \n strong cooling",
+    "med_wind_med_cooling" => "Mid wind \n mid cooling",
+    "strong_wind_weak_cooling" => "Strong wind \n weak cooling",
+)
+
+# case = "free_convection"
+
+fonts = (; regular=texfont())
+set_theme!(Theme(fontsize=24, linewidth=3; fonts))
+fig = Figure(size=(1200, 500))
+axc = []
+
+for (c, case) in enumerate(cases)
+    colors = Makie.wong_colors()
+    colors = reverse(colors)
+
+    if c == length(cases)
+        yaxisposition = :right
+    else
+        yaxisposition = :left
+    end
+
+    x = 1:9
+    xticks = (x .* 1e-4, [string(i) for i in x])
+    ax = Axis(fig[2, c], xlabel="Buoyancy \n (10⁻⁴ × m s⁻²)", ylabel=L"z \, \mathrm{(m)}"; yaxisposition, xticks)
+    title = titles[case]
+    Label(fig[1, c], title, tellwidth=false)
+    push!(axc, ax)
+
+    if c != 1 && c != length(cases)
+        hideydecorations!(ax, grid=false)
+    end
+
+    if c == 1
+        hidespines!(ax, :t, :r)
+    elseif c == 4
+        hidespines!(ax, :t, :l)
+    else
+        hidespines!(ax, :t, :l, :r)
+    end
+    
+    res = "1m"
+    for waves in ("", "_no_stokes")
+        filename = string(case, waves, "_instantaneous_statistics.jld2")
+        datadir = joinpath(basedir, suite, res)
+        filepath = joinpath(datadir, filename)
+        bt = FieldTimeSeries(filepath, "b")
+        bn = interior(bt[end], 1, 1, :) .- bt[1][1, 1, 1]
+        @show length(bt)
+        z = znodes(bt)
+
+        if waves == ""
+            label = "with \nStokes drift"
+        else
+            label = "without"
+        end
+
+        if waves == ""
+            linewidth = 5
+            color = (:black, 0.8)
+            linestyle = :solid
+        else
+            linewidth = 3
+            linestyle = :dash
+            color = pop!(colors)
+        end
+        lines!(ax, bn, z; label, linewidth, linestyle, color)
+    end
+
+    ylims!(ax, -160, 0)
+end
+
+Nc = length(cases)
+#Legend(fig[2, Nc+1], first(axc))
+axislegend(axc[1], position=(0, 0.9))
+
+# For 12 hour suite resolution comparison
+if suite == "6_hour_suite"
+    xlims!(axc[1], 1.7e-4, 3.9e-4)
+    xlims!(axc[2], 1.7e-4, 3.9e-4)
+    xlims!(axc[3], 1.7e-4, 4.6e-4)
+    xlims!(axc[4], 1.7e-4, 6.0e-4)
+elseif suite == "24_hour_suite"
+    xlims!(axc[1], 1.7e-4, 3.9e-4)
+    xlims!(axc[2], 1.7e-4, 3.9e-4)
+    xlims!(axc[3], 1.7e-4, 4.6e-4)
+    xlims!(axc[4], 1.7e-4, 6.0e-4)
+else
+    xlims!(axc[1], 1.7e-4, 3.9e-4)
+    xlims!(axc[2], 1.7e-4, 3.9e-4)
+    xlims!(axc[3], 1.7e-4, 4.6e-4)
+    xlims!(axc[4], 1.7e-4, 6.0e-4)
+end
+
+display(fig)
+save("les_stokes_dependence_$suite.pdf", fig)
+
diff --git a/idealized/quick_peek_averages.jl b/idealized/quick_peek_averages.jl
new file mode 100644
index 0000000..1b3dca9
--- /dev/null
+++ b/idealized/quick_peek_averages.jl
@@ -0,0 +1,69 @@
+using Oceananigans
+using GLMakie
+#using CairoMakie
+
+set_theme!(Theme(fontsize=24, linewidth=3))
+basedir = "/Users/gregorywagner/Projects/LESbrary.jl/data"
+suite = "48_hour_suite"
+resolution = "2m"
+suffix = "_instantaneous_statistics.jld2"
+datadir = joinpath(basedir, suite, resolution)
+
+cases = [
+    "free_convection_with_tracer",
+    "strong_wind_and_sunny_with_tracer",
+    "strong_wind_no_rotation_with_tracer",
+    "strong_wind_with_tracer",
+    "strong_wind_weak_cooling_with_tracer",
+    "weak_wind_strong_cooling_with_tracer",
+    "med_wind_med_cooling_with_tracer",
+]
+
+fig = Figure(size=(1400, 500))
+
+axb = Axis(fig[1, 1])
+axc = Axis(fig[1, 2])
+axu = Axis(fig[1, 3])
+
+#xlims!(axb, -1e-3, 0)
+xlims!(axc, -2, 6)
+
+filepath1 = joinpath(datadir, first(cases) * suffix)
+b1 = FieldTimeSeries(filepath1, "b")
+Nt = length(b1)
+z = znodes(b1)
+t = b1.times
+
+slider = Slider(fig[2, 1:3], startvalue=Nt, range=1:Nt)
+n = slider.value
+
+title = @lift string("t = ", prettytime(t[$n]))
+Label(fig[0, 1:3], title)
+
+colors = [:black, :blue, :red, :seagreen, :orange, :purple, :brown]
+
+for (c, case) in enumerate(cases)
+
+    filepath = joinpath(datadir, case * suffix)
+    bt = FieldTimeSeries(filepath, "b")
+    ct = FieldTimeSeries(filepath, "c")
+    ut = FieldTimeSeries(filepath, "u")
+    vt = FieldTimeSeries(filepath, "v")
+
+    times = bt.times
+    t = times[end]
+    @show prettytime(t)
+
+    bn = @lift interior(bt[$n], 1, 1, :)
+    cn = @lift interior(ct[$n], 1, 1, :)
+    un = @lift interior(ut[$n], 1, 1, :)
+    vn = @lift interior(vt[$n], 1, 1, :)
+
+    lines!(axb, bn, z; color=colors[c], linestyle=:solid)
+    lines!(axc, cn, z; color=colors[c], linestyle=:solid)
+    lines!(axu, un, z; color=colors[c], linestyle=:solid, label="u")
+    lines!(axu, vn, z; color=colors[c], linestyle=:dash,  label="v")
+end
+       
+display(fig)
+
diff --git a/idealized/quick_peek_three_layer_simulation.jl b/idealized/quick_peek_three_layer_simulation.jl
new file mode 100644
index 0000000..04ca6ad
--- /dev/null
+++ b/idealized/quick_peek_three_layer_simulation.jl
@@ -0,0 +1,26 @@
+using Oceananigans
+using GLMakie
+#using CairoMakie
+
+basedir = "/Users/gregorywagner/Projects/LESbrary.jl/data"
+suite = "12_hour_suite"
+resolution = "1m"
+filename = "strong_wind_no_stokes_xz_slice.jld2"
+
+set_theme!(Theme(fontsize=24, linewidth=3))
+fig = Figure(size=(1200, 1200))
+ax = Axis(fig[1, 1])
+datadir = joinpath(basedir, suite, resolution)
+filepath = joinpath(datadir, filename)
+wt = FieldTimeSeries(filepath, "w")
+wn = interior(wt[end], :, 1, :)
+times = wt.times
+t = times[end]
+@show prettytime(t)
+
+wlim = 3 * maximum(wn) / 4
+heatmap!(ax, wn, colormap=:balance, colorrange=(-wlim, wlim))
+   
+display(fig)
+
+
diff --git a/idealized/show_parameters.jl b/idealized/show_parameters.jl
new file mode 100644
index 0000000..4715412
--- /dev/null
+++ b/idealized/show_parameters.jl
@@ -0,0 +1,70 @@
+# Example script for generating LESbrary data
+
+using Oceananigans
+using Oceananigans.Units
+using Printf
+
+using LESbrary.IdealizedExperiments: forty_eight_hour_suite_parameters
+using LESbrary.IdealizedExperiments: twenty_four_hour_suite_parameters
+using LESbrary.IdealizedExperiments: twelve_hour_suite_parameters
+using LESbrary.IdealizedExperiments: six_hour_suite_parameters
+using LESbrary.IdealizedExperiments: seventy_two_hour_suite_parameters
+
+cases = [
+    :free_convection,
+    :weak_wind_strong_cooling,
+    :med_wind_med_cooling,
+    :strong_wind_weak_cooling,
+    :strong_wind,
+    :strong_wind_no_rotation
+]
+
+suites = [
+    six_hour_suite_parameters,
+    twelve_hour_suite_parameters,
+    twenty_four_hour_suite_parameters,
+    forty_eight_hour_suite_parameters,
+    seventy_two_hour_suite_parameters,
+]
+
+for suite in suites
+    for case in cases
+        suite_parameters = suite[case]
+
+        α = 2e-4
+        g = 9.81
+        ρₒ = 1024
+        ρₐ = 1.2
+        cp = 3991
+        κ = 0.4
+        g = 9.81
+        Cg = 0.011
+
+        # Heat flux
+        Jb = suite_parameters[:buoyancy_flux]
+        Q = ρₒ * cp * Jb / (α * g)
+
+        # Deduce ua(z=10m)
+        τₒ = suite_parameters[:momentum_flux]
+        τₐ = ρₒ * τₒ / ρₐ
+
+        # Atmos friction velocity
+        u★ = sqrt(abs(τₐ))
+
+        if τₐ == 0
+            u₁₀ = 0.0
+            c₁₀ = 0.0
+        else
+            ℓ = Cg * abs(τₐ) / g
+            c₁₀ = (κ / log(10 / ℓ))^2
+
+            # τₐ = c10 * u10^2
+            u₁₀ = sqrt(abs(τₐ) / c₁₀)
+        end
+
+        @info string(@sprintf("% 24s:", case),
+                     @sprintf(" c₁₀: %.2e", c₁₀),
+                     @sprintf(" u₁₀: %2d", u₁₀),
+                     @sprintf(" Q: %d", Q))
+    end
+end
diff --git a/idealized/visualize_les_horizontal_averages.jl b/idealized/visualize_les_horizontal_averages.jl
new file mode 100644
index 0000000..1ad765b
--- /dev/null
+++ b/idealized/visualize_les_horizontal_averages.jl
@@ -0,0 +1,238 @@
+using Oceananigans
+using JLD2
+using GLMakie
+using Statistics
+
+fonts = (; regular=texfont())
+set_theme!(Theme(fontsize=32, linewidth=6; fonts))
+
+dir = "/Users/gregorywagner/Projects/LESbrary.jl/data"
+
+fig = Figure(size=(1680, 670))
+
+ax_per_row = 3
+
+case1 = (
+    suite = "6_hour_suite",
+    res = "1m",
+    prefix = "strong_wind_no_rotation",
+    row = 1,
+    label = "(a)",
+)
+
+case2 = (
+    suite = "6_hour_suite",
+    res = "1m",
+    prefix = "strong_wind_weak_cooling",
+    row = 1,
+    label = "(b)",
+)
+
+case3 = (
+    suite = "6_hour_suite",
+    res = "1m",
+    prefix = "free_convection",
+    row = 1,
+    label = "(c)",
+)
+
+case4 = (
+    suite = "12_hour_suite",
+    res = "1m",
+    prefix = "strong_wind_no_rotation",
+    row = 2,
+    label = "(d)",
+)
+
+case5 = (
+    suite = "12_hour_suite",
+    res = "1m",
+    prefix = "strong_wind_weak_cooling",
+    row = 2,
+    label = "(e)",
+)
+
+case6 = (
+    suite = "12_hour_suite",
+    res = "1m",
+    prefix = "free_convection",
+    row = 2,
+    label = "(f)",
+)
+
+case7 = (
+    suite = "72_hour_suite",
+    res = "1m",
+    prefix = "strong_wind_no_rotation",
+    row = 3,
+    label = "(g)",
+)
+
+case8 = (
+    suite = "72_hour_suite",
+    res = "1m",
+    prefix = "strong_wind_weak_cooling",
+    row = 3,
+    label = "(h)",
+)
+
+case9 = (
+    suite = "72_hour_suite",
+    res = "1m",
+    prefix = "free_convection",
+    row = 3,
+    label = "(i)",
+)
+
+case = case5
+
+suite = case.suite
+res = case.res
+prefix = case.prefix
+row = case.row
+label = case.label
+
+xy_filepath = joinpath(dir, suite, res, prefix * "_xy_slice.jld2")
+yz_filepath = joinpath(dir, suite, res, prefix * "_yz_slice.jld2")
+xz_filepath = joinpath(dir, suite, res, prefix * "_xz_slice.jld2")
+
+# T_xy_t = FieldTimeSeries(xy_filepath, "w")
+# T_xz_t = FieldTimeSeries(xz_filepath, "w")
+# T_yz_t = FieldTimeSeries(yz_filepath, "w")
+
+T_xy_t = FieldTimeSeries(xy_filepath, "T")
+T_xz_t = FieldTimeSeries(xz_filepath, "T")
+T_yz_t = FieldTimeSeries(yz_filepath, "T")
+
+statistics_filepath = joinpath(dir, suite, res, prefix * "_instantaneous_statistics.jld2")
+
+Tt = FieldTimeSeries(statistics_filepath, "T")
+Bt = FieldTimeSeries(statistics_filepath, "b")
+Ut = FieldTimeSeries(statistics_filepath, "u")
+Vt = FieldTimeSeries(statistics_filepath, "v")
+W²t = FieldTimeSeries(statistics_filepath, "ww")
+Et = FieldTimeSeries(statistics_filepath, "e")
+
+Tn = interior(Tt[end], 1, 1, :)
+Bn = interior(Bt[end], 1, 1, :)
+Un = interior(Ut[end], 1, 1, :)
+Vn = interior(Vt[end], 1, 1, :)
+En = interior(Et[end], 1, 1, :)
+W²n = interior(W²t[end], 1, 1, :)
+
+times = T_xy_t.times
+Nt = length(times)
+grid = T_xy_t.grid
+
+Nx, Ny, Nz = size(grid)
+x, y, z = nodes(T_xy_t)
+Lx = grid.Lx
+Ly = grid.Ly
+Lz = grid.Lz
+
+# Nz += 1
+
+x_xz = repeat(x, 1, Nz)
+z_xz = repeat(reshape(z, 1, Nz), Nx, 1)
+y_xz = 0.995 * Ly * ones(Nx, Nz)
+
+y_yz = repeat(y, 1, Nz)
+z_yz = repeat(reshape(z, 1, Nz), grid.Ny, 1)
+x_yz = 0.995 * Lx * ones(Ny, Nz)
+
+# Slight displacements to "stitch" the cube together
+x_xy = x
+y_xy = y
+z_xy = -0.001 * Lz * ones(Nx, Ny)
+
+azimuth = 6.7
+elevation = 0.50
+perspectiveness = 0.1
+xlabel = L"x \, \mathrm{(m)}"
+ylabel = L"y \, \mathrm{(m)}"
+zlabel = L"z \, \mathrm{(m)}"
+aspect = :data
+xlabeloffset = 90
+ylabeloffset = 70
+zlabeloffset = 100
+
+ax_b = fig[2, 1] = Axis3(fig; aspect, xlabel, ylabel, zlabel, azimuth, elevation, perspectiveness,
+                         xlabeloffset, ylabeloffset, zlabeloffset)
+
+n = length(T_xy_t) #140
+T_xy = interior(T_xy_t[n], :, :, 1)
+T_xz = interior(T_xz_t[n], :, 1, :)
+T_yz = interior(T_yz_t[n], 1, :, :)
+
+T_xy .-= Tn[end-1]
+T_xz .-= reshape(Tn, 1, Nz)
+T_yz .-= reshape(Tn, 1, Nz)
+
+ϵ = mean(Bn ./ Tn) # <<eyes>>
+b_xy = ϵ .* T_xy
+b_xz = ϵ .* T_xz
+b_yz = ϵ .* T_yz
+
+bmax = maximum(abs, b_xz)
+blim = bmax / 6
+colorrange_b = (-blim, blim)
+colormap_b = :balance
+
+pl = surface!(ax_b, x_xz, y_xz, z_xz; color=b_xz, colormap=colormap_b, colorrange=colorrange_b)
+     surface!(ax_b, x_yz, y_yz, z_yz; color=b_yz, colormap=colormap_b, colorrange=colorrange_b)
+     surface!(ax_b, x_xy, y_xy, z_xy; color=b_xy, colormap=colormap_b, colorrange=colorrange_b)
+
+xlims!(ax_b, 0, 512)
+ylims!(ax_b, 0, 512)
+zlims!(ax_b, -256, 0)
+
+ticks = ([-2e-5, -1e-5, 0, 1e-5, 2e-5], ["-2", "-1", "0", "1", "2"])
+Colorbar(fig[1, 1], pl; ticks, label=L"\mathrm{Buoyancy \, perturbation} \, b\prime \, (10^{-5} \times \mathrm{m \, s^{-2}})", vertical=false, width=Relative(0.7))
+
+# Horizontal averages
+xticks = ([0, 1e-4, 2e-4, 3e-4, 4e-4], ["0", "1", "2", "3", "4"])
+ax_B = Axis(fig[1:2, 2]; xticks, xlabel="Buoyancy \n (10⁻⁴ × m s⁻²)", ylabel=L"z \, \mathrm{(m)}")
+
+z = znodes(Bt)
+lines!(ax_B, Bn .- Bn[1], z)
+
+xticks = ([-0.1, 0, 0.1], ["-0.1", "0", "0.1"])
+ax_u = Axis(fig[1:2, 3]; xticks, xlabel="Velocities \n (m s⁻¹)", ylabel=L"z \, \mathrm{(m)}")
+
+lines!(ax_u, Un, z, color=:black, label=L"u")
+lines!(ax_u, Vn, z, color=:forestgreen, label=L"v")
+xlims!(ax_u, -0.19, 0.19)
+axislegend(ax_u, position=:lb)
+
+xticks = ([0, 2e-3, 4e-3], ["0", "2", "4"])
+ax_e = Axis(fig[1:2, 4]; xticks, xlabel="Kinetic energies \n (10⁻³ × m² s⁻²)", ylabel=L"z \, \mathrm{(m)}", yaxisposition=:right)
+
+zw = znodes(W²t)
+
+lines!(ax_e, En, z,   label=L"\mathscr{E}")
+lines!(ax_e, W²n, zw, label=L"w'^2")
+axislegend(ax_e, position=:rb)
+
+#colsize!(fig.layout, 1, Relative(0.6))
+colsize!(fig.layout, 2, Relative(0.17))
+colsize!(fig.layout, 3, Relative(0.17))
+colsize!(fig.layout, 4, Relative(0.17))
+
+hidespines!(ax_B, :t, :r)
+hidespines!(ax_u, :t, :l, :r)
+hidespines!(ax_e, :t, :l)
+hideydecorations!(ax_u, grid=false)
+
+xtext = 440
+ytext = 0
+ztext = 120
+text!(ax_b, xtext, ytext, ztext, text="(a)")
+
+text!(ax_B, 0.95, 0.03, text="(b)", align=(:right, :bottom), space=:relative)
+text!(ax_u, 0.95, 0.03, text="(c)", align=(:right, :bottom), space=:relative)
+text!(ax_e, 0.95, 0.25, text="(d)", align=(:right, :bottom), space=:relative)
+
+display(fig)
+
+save("les_horizontal_averages.png", fig)
+
diff --git a/idealized/visualize_les_w.jl b/idealized/visualize_les_w.jl
new file mode 100644
index 0000000..989d989
--- /dev/null
+++ b/idealized/visualize_les_w.jl
@@ -0,0 +1,184 @@
+using Oceananigans
+using JLD2
+using GLMakie
+using MathTeXEngine
+
+fonts = (; regular=texfont())
+set_theme!(Theme(fontsize=28; fonts))
+
+dir = "/Users/gregorywagner/Projects/LESbrary.jl/data"
+
+fig = Figure(size=(1500, 1200))
+
+#ax_b1 = Axis(fig[1, 4], xlabel="Temperature (ᵒC)", ylabel="z (m)")
+#ax_b2 = Axis(fig[2, 4], xlabel="Temperature (ᵒC)", ylabel="z (m)")
+#ax_u = Axis(fig[c, 3], xlabel="Velocities (m s⁻¹)", ylabel="z (m)", yaxisposition=:right)
+
+ax_per_row = 3
+
+case1 = (
+    suite = "6_hour_suite",
+    res = "1m",
+    prefix = "strong_wind_no_rotation",
+    row = 1,
+    label = "(a)",
+)
+
+case2 = (
+    suite = "6_hour_suite",
+    res = "1m",
+    prefix = "strong_wind_weak_cooling",
+    row = 1,
+    label = "(b)",
+)
+
+case3 = (
+    suite = "6_hour_suite",
+    res = "1m",
+    prefix = "free_convection",
+    row = 1,
+    label = "(c)",
+)
+
+case4 = (
+    suite = "12_hour_suite",
+    res = "1m",
+    prefix = "strong_wind_no_rotation",
+    row = 2,
+    label = "(d)",
+)
+
+case5 = (
+    suite = "12_hour_suite",
+    res = "1m",
+    prefix = "strong_wind_weak_cooling",
+    row = 2,
+    label = "(e)",
+)
+
+case6 = (
+    suite = "12_hour_suite",
+    res = "1m",
+    prefix = "free_convection",
+    row = 2,
+    label = "(f)",
+)
+
+case7 = (
+    suite = "72_hour_suite",
+    res = "1m",
+    prefix = "strong_wind_no_rotation",
+    row = 3,
+    label = "(g)",
+)
+
+case8 = (
+    suite = "72_hour_suite",
+    res = "1m",
+    prefix = "strong_wind_weak_cooling",
+    row = 3,
+    label = "(h)",
+)
+
+case9 = (
+    suite = "72_hour_suite",
+    res = "1m",
+    prefix = "free_convection",
+    row = 3,
+    label = "(i)",
+)
+
+cases = [case1, case2, case3,
+         case4, case5, case6,
+         case7, case8, case9]
+
+for (c, case) in enumerate(cases)
+
+    suite = case.suite
+    res = case.res
+    prefix = case.prefix
+    row = case.row
+    label = case.label
+
+    xy_filepath = joinpath(dir, suite, res, prefix * "_with_tracer_xy_slice.jld2")
+    yz_filepath = joinpath(dir, suite, res, prefix * "_with_tracer_yz_slice.jld2")
+    xz_filepath = joinpath(dir, suite, res, prefix * "_with_tracer_xz_slice.jld2")
+
+    w_xy_t = FieldTimeSeries(xy_filepath, "w")
+    w_xz_t = FieldTimeSeries(xz_filepath, "w")
+    w_yz_t = FieldTimeSeries(yz_filepath, "w")
+
+    statistics_filepath = joinpath(dir, suite, res, prefix * "_instantaneous_statistics.jld2")
+
+    times = w_xy_t.times
+    Nt = length(times)
+
+    grid = w_xy_t.grid
+
+    Nx, Ny, Nz = size(grid)
+    x, y, z = nodes(w_xy_t)
+    Lx = grid.Lx
+    Ly = grid.Ly
+    Lz = grid.Lz
+
+    Nz += 1
+
+    x_xz = repeat(x, 1, Nz)
+    z_xz = repeat(reshape(z, 1, Nz), Nx, 1)
+    y_xz = 0.995 * Ly * ones(Nx, Nz)
+
+    y_yz = repeat(y, 1, Nz)
+    z_yz = repeat(reshape(z, 1, Nz), grid.Ny, 1)
+    x_yz = 0.995 * Lx * ones(Ny, Nz)
+
+    # Slight displacements to "stitch" the cube together
+    x_xy = x
+    y_xy = y
+    z_xy = -0.001 * Lz * ones(Nx, Ny)
+
+    azimuth = 6.7
+    elevation = 0.50
+    perspectiveness = 0.1
+    xlabel = L"x \, \mathrm{(m)}"
+    ylabel = L"y \, \mathrm{(m)}"
+    zlabel = L"z \, \mathrm{(m)}"
+    aspect = :data
+    xlabeloffset = 60
+    zlabeloffset = 80
+
+    j = c - ax_per_row * (row - 1)
+    ax_w = fig[row, j] = Axis3(fig; aspect, xlabel, ylabel, zlabel, azimuth, elevation, perspectiveness,
+                               xlabeloffset, zlabeloffset)
+
+    if c != 2
+        hidedecorations!(ax_w)
+    end
+
+    n = length(w_xy_t) #140
+    w_xy = interior(w_xy_t[n], :, :, 1)
+    w_xz = interior(w_xz_t[n], :, 1, :)
+    w_yz = interior(w_yz_t[n], 1, :, :)
+
+    wmax = maximum(abs, w_xy)
+    wlim = wmax / 2
+    @show wmax
+    colorrange_w = (-wlim, wlim)
+    colormap_w = :balance
+
+    pl = surface!(ax_w, x_xz, y_xz, z_xz; color=w_xz, colormap=colormap_w, colorrange=colorrange_w)
+         surface!(ax_w, x_yz, y_yz, z_yz; color=w_yz, colormap=colormap_w, colorrange=colorrange_w)
+         surface!(ax_w, x_xy, y_xy, z_xy; color=w_xy, colormap=colormap_w, colorrange=colorrange_w)
+
+    xtext = 440
+    ytext = 0
+    ztext = 120
+    text!(ax_w, xtext, ytext, ztext, text=label)
+
+    xlims!(ax_w, 0, 512)
+    ylims!(ax_w, 0, 512)
+    zlims!(ax_w, -256, 0)
+end
+
+display(fig)
+save("les_w_visualization.png", fig)
+
diff --git a/src/IdealizedExperiments/three_layer_constant_fluxes.jl b/src/IdealizedExperiments/three_layer_constant_fluxes.jl
index 9aaaa25..a7a788b 100644
--- a/src/IdealizedExperiments/three_layer_constant_fluxes.jl
+++ b/src/IdealizedExperiments/three_layer_constant_fluxes.jl
@@ -220,7 +220,6 @@ function three_layer_constant_fluxes_simulation(;
     end
     
     model = NonhydrostaticModel(; grid, buoyancy, tracers, stokes_drift, closure, advection,
-                                timestepper = :RungeKutta3,
                                 coriolis = FPlane(; f),
                                 boundary_conditions = (b=b_bcs, u=u_bcs),
                                 forcing = (u=u_sponge, v=v_sponge, w=w_sponge, b=b_forcing, c=c_forcing))
diff --git a/test/test_turbulence_statistics.jl b/test/test_turbulence_statistics.jl
index 8621e21..1f52e5d 100644
--- a/test/test_turbulence_statistics.jl
+++ b/test/test_turbulence_statistics.jl
@@ -46,7 +46,7 @@ for arch in architectures
     @testset "Turbulence Statistics [$(typeof(arch))]" begin
         @info "Testing turbulence statistics [$(typeof(arch))]..."
 
-        model = NonhydrostaticModel(grid = RectilinearGrid(arch, size=(1, 1, 1), extent=(1, 1, 1)),
+        model = NonhydrostaticModel(grid = RectilinearGrid(arch, size=(8, 8, 8), extent=(1, 1, 1)),
                                     tracers = :b,
                                     buoyancy = BuoyancyTracer(),
                                     closure = AnisotropicMinimumDissipation())