Merge pull request #13 from NicolasRiel/main

subduction model did not have diffusivity activated
JuliaGeodynamics · Nov 1, 2023 · a4a33b5 · a4a33b5
2 parents f7c0e46 + f329956
commit a4a33b5
Show file tree

Hide file tree

Showing 3 changed files with 10 additions and 325 deletions.
diff --git a/docs/make.jl b/docs/make.jl
@@ -14,6 +14,7 @@ makedocs(;
         "Create & run LaMEM models from julia" => ["Overview" => "man/juliasetups.md",
                                                     "Example 1: Sphere" => "man/juliasetup_example_sphere.md",
                                                     "Example 2: Volcano" => "man/juliasetup_LaPalma.md",
+                                                    "Example 3: Subduction" => "man/juliasetup_TMSubduction.md",
                                                     "Notebooks" => "man/juliasetup_pluto.md",
                                                     "Available functions" => "man/LaMEM_ModelFunctions.md"],
         "Run LaMEM" => "man/runlamem.md",

diff --git a/docs/src/man/juliasetup_TMSubduction.md b/docs/src/man/juliasetup_TMSubduction.md
@@ -29,8 +29,13 @@ model = Model(Grid( x   = [-2000.,2000.],
                     z   = [-660,40] ,
                     nel = (512,1,128)     ),
 
-                    BoundaryConditions(temp_bot=1565, open_top_bound=1),
-                    Scaling(GEO_units(temperature=1000, stress=1e9Pa, length=1km, viscosity=1e20Pa*s)),
+                    BoundaryConditions( temp_bot        = 1565.0,
+                                        temp_top        = 20.0,
+                                        open_top_bound  = 1),
+                    Scaling(GEO_units(  temperature     = 1000,
+                                        stress          = 1e9Pa,
+                                        length          = 1km,
+                                        viscosity       = 1e20Pa*s) ),
                     Time(nstep_max=20) )  
 ```
 
@@ -52,6 +57,8 @@ We activate shear heating and adiabatic heating and set the minimum and maximum
 ```
 model.SolutionParams = SolutionParams(  shear_heat_eff = 1.0,
                                         Adiabatic_Heat = 1.0,
+                                        act_temp_diff  = 1,
+                                        FSSA           = 1.0,
                                         eta_min   = 5e18,
                                         eta_ref   = 1e21,
                                         eta_max   = 1e25,

diff --git a/docs/src/man/juliasetup_TMSubduction.md.backup b/docs/src/man/juliasetup_TMSubduction.md.backup