docs

Author: albert-de-montserrat

docs/make.jl

@@ -105,7 +105,8 @@ makedocs(;
             "19 - Jura tutorial" =>  "man/Tutorial_Jura.md",
             "20 - 2D model setups" => "man/Tutorial_NumericalModel_2D.md",
             "21 - 3D model setups" => "man/Tutorial_NumericalModel_3D.md",
-            "22 - Build geometry from polygons" =>  "man/tutorial_Polygon_structures.md"
+            "22 - 3D model setups" => "man/Tutorial_VolcanoModel_3D.md",
+            "23 - Build geometry from polygons" =>  "man/tutorial_Polygon_structures.md"
         ],
         "User Guide" => Any[
             "Installation" =>  "man/installation.md",

docs/src/man/Tutorial_VolcanoModel_3D.md

@@ -1,15 +1,11 @@
-```@meta
-EditURL = "../../../tutorials/Tutorial_NumericalModel_3D.jl"
-```
-
 # Simple model of a magmatic chamber with a volcano on top 
 
 ### Aim
 The aim of this tutorial is to show you how to create 3D numerical model setups that can be used as initial setups for other codes.
 
 ### Generating the model
 
-Lets start with creating a 3D model setup in cartesian coordinates, which uses the `CartData` data structure, with a resolution of $ 128 \times 128 \times 128 $ inside the domain $\Omega \in [-100,100] \times [-100,100] \times [-110,50]$ km
+Lets start with creating a 3D model setup in cartesian coordinates, which uses the `CartData` data structure, with a resolution of $ 128 \times 128 \times 128 $ grid points, inside the domain $\Omega \in [-100,100] \times [-100,100] \times [-110,50]$ km
 
 ```julia
 using GeophysicalModelGenerator
@@ -57,7 +53,7 @@ add_volcano!(Phases, Temp, Grid;
     T = HalfspaceCoolingTemp(Age=20)
 )
 ```
-
+We can also add a magmatic chamber located below the volcano
 ```julia
 add_ellipsoid!(Phases, Temp, Grid; 
     cen    = (0, 0, -40), 
@@ -66,12 +62,17 @@ add_ellipsoid!(Phases, Temp, Grid;
 )
 ```
 
+where we prescribe a constant temperature of $T=1400^{\circ}C$
 ```julia
 @. Temp[Phases == 4] = 1400 
-@. Temp[Phases == 0] = 0 
 Grid = addfield(Grid, (;Phases, Temp))
 ```
 
+Finally we setup the temperature of the air to $T^{\text{air}}=0^{\circ}C$
+```julia
+@. Temp[Phases == 0] = 0 
+```
+
 ```julia
 write_paraview(Grid,"VolcanoModel3D");
 ```