Skip to content

Commit

Permalink
Merge pull request #676 from MFraters/add_BST_mass_conserving_and_plu…
Browse files Browse the repository at this point in the history
…me_chapters

Add basic starter tutorial mass conserving and plume chapters
  • Loading branch information
MFraters authored Feb 27, 2024
2 parents 601b6d2 + 58269d7 commit df88fa0
Show file tree
Hide file tree
Showing 26 changed files with 56,408 additions and 17,570 deletions.
22 changes: 22 additions & 0 deletions doc/sphinx/_static/gwb_input_files/BST_16_mass_conserving.grid
Original file line number Diff line number Diff line change
@@ -0,0 +1,22 @@
# output variables
grid_type = cartesian
dim = 3
compositions = 4

# domain of the grid
x_min = -1000e3
x_max = 2000e3
y_min = 0e3
y_max = 1000e3
z_min = 0
z_max = 600e3

# low res grid properties
n_cell_x = 150
n_cell_y = 50
n_cell_z = 30

# shown grid properties
# n_cell_x = 1200
# n_cell_y = 400
# n_cell_z = 240
Original file line number Diff line number Diff line change
@@ -1,7 +1,6 @@
{
"version": "0.6",
"coordinate system":{"model":"cartesian"},
"cross section":[[0,200e3],[10e3,200e3]],
"features":
[
{
Expand Down Expand Up @@ -35,32 +34,36 @@
]
},
{
"model":"oceanic plate", "name":"Subducting Oceanic plate", "max depth":100e3,
"model":"oceanic plate", "name":"Subducting Oceanic plate", "max depth":300e3,
"coordinates":[[2000e3,0],[2000e3,1000e3],[1500e3,1000e3],[1600e3,350e3],[1500e3,0]],
"temperature models":[{"model":"linear", "max depth":100e3}],
"temperature models":[{"model":"half space model", "max depth":300e3,"spreading velocity":0.02,
"ridge coordinates":[[[5e6,-1],[5e6,2000e3]]]}],
"composition models":[{"model":"uniform", "compositions":[3], "max depth":50e3},
{"model":"uniform", "compositions":[1], "min depth":50e3}]
{"model":"uniform", "compositions":[1], "min depth":50e3, "max depth":100e3}]
},
{
"model":"subducting plate", "name":"Slab", "dip point":[0,0],
"coordinates":[[1500e3,1000e3],[1600e3,350e3],[1500e3,0]],
"segments":
[
{"length":300e3, "thickness":[100e3], "angle":[0,60],
{"length":300e3, "thickness":[300e3], "top truncation":[-100e3], "angle":[0,60],
"composition models":[
{"model":"uniform", "compositions":[3], "max distance slab top":50e3},
{"model":"uniform", "compositions":[2], "min distance slab top":50e3}]},
{"length":500e3, "thickness":[100e3], "angle":[60,20]}
{"model":"uniform", "compositions":[2], "min distance slab top":50e3, "max distance slab top":100e3}]},
{"length":500e3, "thickness":[300e3], "top truncation":[-100e3], "angle":[60,20]}
],
"sections":
[
{"coordinate":0,
"segments":[{"length":300e3, "thickness":[100e3], "angle":[0,60]},
{"length":400e3, "thickness":[100e3,50e3], "angle":[60]}],
"composition models":[{"model":"uniform", "compositions":[1]}]}
"segments":[{"length":300e3, "thickness":[300e3], "top truncation":[-100e3], "angle":[0,60]},
{"length":400e3, "thickness":[300e3], "top truncation":[-100e3], "angle":[60]}],
"composition models":[{"model":"uniform", "compositions":[1],"max distance slab top":100e3}]}
],
"temperature models":[{"model":"plate model", "density":3300, "plate velocity":0.02}],
"composition models":[{"model":"uniform", "compositions":[2]}]
"temperature models":[
{"model":"mass conserving", "density":3300, "plate velocity":0.02,
"ridge coordinates":[[[5e6,-1],[5e6,2000e3]]],"coupling depth":50e3,
"min distance slab top":-200e3, "max distance slab top":300e3}],
"composition models":[{"model":"uniform", "compositions":[2], "max distance slab top":100e3}]
}
]
}
22 changes: 22 additions & 0 deletions doc/sphinx/_static/gwb_input_files/BST_17_plume.grid
Original file line number Diff line number Diff line change
@@ -0,0 +1,22 @@
# output variables
grid_type = cartesian
dim = 3
compositions = 4

# domain of the grid
x_min = -1000e3
x_max = 2000e3
y_min = 0e3
y_max = 1000e3
z_min = 0
z_max = 600e3

# low res grid properties
n_cell_x = 150
n_cell_y = 50
n_cell_z = 30

# shown grid properties
# n_cell_x = 1200
# n_cell_y = 400
# n_cell_z = 240
90 changes: 90 additions & 0 deletions doc/sphinx/_static/gwb_input_files/BST_17_plume.wb
Original file line number Diff line number Diff line change
@@ -0,0 +1,90 @@
{
"version": "0.6",
"coordinate system":{"model":"cartesian"},
"features":
[
{
"model":"mantle layer", "name":"upper mantle", "max depth":660e3,
"coordinates":[[-1000e3,0],[-1000e3,1000e3],[3000e3,1000e3],[3000e3,0]],
"composition models":[{"model":"uniform", "compositions":[4]}]
},
{
"model":"oceanic plate", "name":"Overriding Plate", "max depth":100e3,
"coordinates":[[0,0],[0,1000e3],[1500e3,1000e3],[1600e3,350e3],[1500e3,0]],
"temperature models":
[
{"model":"half space model", "max depth":100e3, "spreading velocity":0.04,
"ridge coordinates":[[[400e3,-1],[-100e3,2000e3]]]}
],
"composition models":[{"model":"uniform", "compositions":[0], "max depth":50e3}]
},
{
"model":"continental plate", "name":"Passive margin", "max depth":[[200e3]],
"coordinates":[[-1000e3,0],[-1000e3,1000e3],[0,1000e3],[0,0]],
"temperature models":
[
{"model":"linear", "max depth":[[100e3],[200e3,[[-250e3,0],[-750e3,1000e3]]]]}
],
"composition models":
[
{"model":"uniform", "compositions":[3],
"max depth":[[100e3],[200e3,[[-250e3,0],[-750e3,1000e3]]]]},
{"model":"uniform", "compositions":[1],
"min depth":[[100e3],[200e3,[[-250e3,0],[-750e3,1000e3]]]]}
]
},
{
"model":"oceanic plate", "name":"Subducting Oceanic plate", "max depth":300e3,
"coordinates":[[2000e3,0],[2000e3,1000e3],[1500e3,1000e3],[1600e3,350e3],[1500e3,0]],
"temperature models":[{"model":"half space model", "max depth":300e3,"spreading velocity":0.02,
"ridge coordinates":[[[5e6,-1],[5e6,2000e3]]]}],
"composition models":[{"model":"uniform", "compositions":[3], "max depth":50e3},
{"model":"uniform", "compositions":[1], "min depth":50e3, "max depth":100e3}]
},
{
"model":"subducting plate", "name":"Slab", "dip point":[0,0],
"coordinates":[[1500e3,1000e3],[1600e3,350e3],[1500e3,0]],
"segments":
[
{"length":300e3, "thickness":[300e3], "top truncation":[-100e3], "angle":[0,60],
"composition models":[
{"model":"uniform", "compositions":[3], "max distance slab top":50e3},
{"model":"uniform", "compositions":[2], "min distance slab top":50e3, "max distance slab top":100e3}]},
{"length":500e3, "thickness":[300e3], "top truncation":[-100e3], "angle":[60,20]}
],
"sections":
[
{"coordinate":0,
"segments":[{"length":300e3, "thickness":[300e3], "top truncation":[-100e3], "angle":[0,60]},
{"length":400e3, "thickness":[300e3], "top truncation":[-100e3], "angle":[60]}],
"composition models":[{"model":"uniform", "compositions":[1],"max distance slab top":100e3}]}
],
"temperature models":[
{"model":"mass conserving", "density":3300, "plate velocity":0.02,
"ridge coordinates":[[[5e6,-1],[5e6,2000e3]]],"coupling depth":50e3,
"min distance slab top":-200e3, "max distance slab top":300e3}],
"composition models":[{"model":"uniform", "compositions":[2], "max distance slab top":100e3}]
},
{
"model":"plume", "name":"Hot spot",
"coordinates":[[200e3,450e3],[150e3,450e3],[110e3,450e3],[70e3,450e3],[50e3,450e3],[0e3,450e3]],
"cross section depths":[50e3,100e3,200e3,400e3,500e3,600e3],
"semi-major axis":[250e3, 150e3, 50e3,50e3,50e3,50e3],
"eccentricity":[0.75, 0.75, 0.5, 0.5, 0.5,0.5],
"rotation angles":[50, 50, 5, 25, 50, 50],
"temperature models":
[
{
"model":"gaussian", "operation":"add",
"centerline temperatures":[225, 175, 185, 195, 205, 215],
"gaussian sigmas":[0.3, 0.3, 0.3,0.3,0.3,0.3],
"depths":[50e3,100e3,200e3,400e3,500e3,600e3]
}
],
"composition models":
[
{"model":"uniform","compositions":[4], "min depth":0}
]
}
]
}
Original file line number Diff line number Diff line change
@@ -1,7 +1,7 @@
# output variables
grid_type = cartesian
dim = 2
compositions = 4
compositions = 6
vtu_output_format = ASCII

# domain of the grid
Expand Down
85 changes: 85 additions & 0 deletions doc/sphinx/_static/gwb_input_files/BST_18_2D_models.wb
Original file line number Diff line number Diff line change
@@ -0,0 +1,85 @@
{
"version": "0.6",
"coordinate system":{"model":"cartesian"},
"cross section":[[0,450e3],[10e3,450e3]],
"features":
[
{
"model":"mantle layer", "name":"upper mantle", "max depth":660e3,
"coordinates":[[-1000e3,0],[-1000e3,1000e3],[3000e3,1000e3],[3000e3,0]],
"composition models":[{"model":"uniform", "compositions":[4]}]
},
{
"model":"oceanic plate", "name":"Overriding Plate", "max depth":100e3,
"coordinates":[[0,0],[0,1000e3],[1500e3,1000e3],[1600e3,350e3],[1500e3,0]],
"temperature models":
[
{"model":"half space model", "max depth":100e3, "spreading velocity":0.04,
"ridge coordinates":[[[400e3,-1],[-100e3,2000e3]]]}
],
"composition models":[{"model":"uniform", "compositions":[0], "max depth":50e3}]
},
{
"model":"continental plate", "name":"Passive margin", "max depth":[[200e3]],
"coordinates":[[-1000e3,0],[-1000e3,1000e3],[0,1000e3],[0,0]],
"temperature models":
[
{"model":"linear", "max depth":[[100e3],[200e3,[[-250e3,0],[-750e3,1000e3]]]]}
],
"composition models":
[
{"model":"uniform", "compositions":[3],
"max depth":[[100e3],[200e3,[[-250e3,0],[-750e3,1000e3]]]]},
{"model":"uniform", "compositions":[1],
"min depth":[[100e3],[200e3,[[-250e3,0],[-750e3,1000e3]]]]}
]
},
{
"model":"oceanic plate", "name":"Subducting Oceanic plate", "max depth":300e3,
"coordinates":[[2000e3,0],[2000e3,1000e3],[1500e3,1000e3],[1600e3,350e3],[1500e3,0]],
"temperature models":[{"model":"half space model", "max depth":300e3,"spreading velocity":0.02,
"ridge coordinates":[[[5e6,-1],[5e6,2000e3]]]}],
"composition models":[{"model":"uniform", "compositions":[3], "max depth":50e3},
{"model":"uniform", "compositions":[1], "min depth":50e3, "max depth":100e3}]
},
{
"model":"subducting plate", "name":"Slab", "dip point":[0,0],
"coordinates":[[1500e3,1000e3],[1600e3,350e3],[1500e3,0]],
"segments":
[
{"length":300e3, "thickness":[300e3], "top truncation":[-100e3], "angle":[0,60],
"composition models":[
{"model":"uniform", "compositions":[3], "max distance slab top":50e3},
{"model":"uniform", "compositions":[2], "min distance slab top":50e3, "max distance slab top":100e3}]},
{"length":500e3, "thickness":[300e3], "top truncation":[-100e3], "angle":[60,20]}
],
"sections":
[
{"coordinate":0,
"segments":[{"length":300e3, "thickness":[300e3], "top truncation":[-100e3], "angle":[0,60]},
{"length":400e3, "thickness":[300e3], "top truncation":[-100e3], "angle":[60]}],
"composition models":[{"model":"uniform", "compositions":[1],"max distance slab top":100e3}]}
],
"temperature models":[{"model":"mass conserving", "density":3300, "plate velocity":0.02,
"ridge coordinates":[[[5e6,-1],[5e6,2000e3]]],"coupling depth":50e3,
"min distance slab top":-200e3, "max distance slab top":300e3}],
"composition models":[{"model":"uniform", "compositions":[2], "max distance slab top":100e3}]
},
{
"model":"plume", "name":"Hot spot",
"coordinates":[[200e3,450e3],[150e3,450e3],[110e3,450e3],[70e3,450e3],[50e3,450e3],[0e3,450e3]],
"cross section depths":[50e3,100e3,200e3,400e3,500e3,600e3],
"semi-major axis":[250e3, 150e3, 50e3,50e3,50e3,50e3],
"eccentricity":[0.75, 0.75, 0.5, 0.5, 0.5,0.5],
"rotation angles":[50, 50, 5, 25, 50, 50],
"temperature models":
[
{"model":"gaussian", "operation":"add",
"centerline temperatures":[225, 175, 185, 195, 205, 215],
"gaussian sigmas":[0.3, 0.3, 0.3,0.3,0.3,0.3],
"depths":[50e3,100e3,200e3,400e3,500e3,600e3]}
],
"composition models":[{"model":"uniform","compositions":[5], "min depth":0}]
}
]
}
Original file line number Diff line number Diff line change
@@ -1,7 +1,7 @@
# output variables
grid_type = chunk
dim = 3
compositions = 4
compositions = 6

# domain of the grid
x_min = -10
Expand All @@ -12,9 +12,9 @@ z_min = 5671000
z_max = 6371000

# low res grid properties
n_cell_x = 150
n_cell_y = 50
n_cell_z = 30
n_cell_x = 450
n_cell_y = 150
n_cell_z = 90

# shown grid properties
# n_cell_x = 1200
Expand Down
Original file line number Diff line number Diff line change
Expand Up @@ -11,7 +11,7 @@
},
{
"model":"oceanic plate", "name":"Overriding Plate", "max depth":100e3,
"coordinates":[[0,0],[0,10],[150,10],[16,3.5],[15,0]],
"coordinates":[[0,0],[0,10],[15,10],[16,3.5],[15,0]],
"temperature models":
[
{"model":"half space model", "max depth":100e3, "spreading velocity":0.04,
Expand All @@ -35,32 +35,51 @@
]
},
{
"model":"oceanic plate", "name":"Subducting Oceanic plate", "max depth":100e3,
"model":"oceanic plate", "name":"Subducting Oceanic plate", "max depth":300e3,
"coordinates":[[20,0],[20,10],[15,10],[16,3.5],[15,0]],
"temperature models":[{"model":"linear", "max depth":100e3}],
"temperature models":[{"model":"half space model", "max depth":300e3,"spreading velocity":0.02,
"ridge coordinates":[[[30,0],[30,10]]]}],
"composition models":[{"model":"uniform", "compositions":[3], "max depth":50e3},
{"model":"uniform", "compositions":[1], "min depth":50e3}]
{"model":"uniform", "compositions":[1], "min depth":50e3, "max depth":100e3}]
},
{
"model":"subducting plate", "name":"Slab", "dip point":[0,0],
"coordinates":[[15,10],[16,3.5],[15,0]],
"segments":
[
{"length":300e3, "thickness":[100e3], "angle":[0,60],
{"length":300e3, "thickness":[300e3], "top truncation":[-100e3], "angle":[0,60],
"composition models":[
{"model":"uniform", "compositions":[3], "max distance slab top":50e3},
{"model":"uniform", "compositions":[2], "min distance slab top":50e3}]},
{"length":500e3, "thickness":[100e3], "angle":[60,20]}
{"model":"uniform", "compositions":[2], "min distance slab top":50e3, "max distance slab top":100e3}]},
{"length":500e3, "thickness":[300e3], "top truncation":[-100e3], "angle":[60,20]}
],
"sections":
[
{"coordinate":0,
"segments":[{"length":300e3, "thickness":[100e3], "angle":[0,60]},
{"length":400e3, "thickness":[100e3,50e3], "angle":[60]}],
"composition models":[{"model":"uniform", "compositions":[1]}]}
"segments":[{"length":300e3, "thickness":[300e3], "top truncation":[-100e3], "angle":[0,60]},
{"length":400e3, "thickness":[300e3], "top truncation":[-100e3], "angle":[60]}],
"composition models":[{"model":"uniform", "compositions":[1], "max distance slab top":100e3}]}
],
"temperature models":[{"model":"plate model", "density":3300, "plate velocity":0.02}],
"composition models":[{"model":"uniform", "compositions":[2]}]
"temperature models":[{"model":"mass conserving", "density":3300, "plate velocity":0.02,
"ridge coordinates":[[[30,0],[30,10]]],"coupling depth":50e3,
"min distance slab top":-200e3, "max distance slab top":300e3}],
"composition models":[{"model":"uniform", "compositions":[2], "max distance slab top":100e3}]
},
{
"model":"plume", "name":"Hot spot",
"coordinates":[[2.5,4],[2.0,4],[1.5,4],[1.0,4],[0.5,4],[0.0,4]],
"cross section depths":[50e3,100e3,200e3,400e3,500e3,700e3],
"semi-major axis":[3, 1, 0.25,0.25,0.25,0.25],
"eccentricity":[0.75, 0.75, 0.5, 0.5, 0.5,0.5],
"rotation angles":[50, 50, 5, 25, 50, 50],
"temperature models":
[
{"model":"gaussian", "operation":"add",
"centerline temperatures":[225, 175, 185, 195, 205, 215],
"gaussian sigmas":[0.3, 0.3, 0.3,0.3,0.3,0.3],
"depths":[50e3,100e3,200e3,400e3,500e3,600e3]}
],
"composition models":[{"model":"uniform","compositions":[5], "min depth":0}]
}
]
}
Loading
Sorry, something went wrong. Reload?
Sorry, we cannot display this file.
Sorry, this file is invalid so it cannot be displayed.
Binary file not shown.
Loading
Sorry, something went wrong. Reload?
Sorry, we cannot display this file.
Sorry, this file is invalid so it cannot be displayed.
Loading
Sorry, something went wrong. Reload?
Sorry, we cannot display this file.
Sorry, this file is invalid so it cannot be displayed.
Loading
Sorry, something went wrong. Reload?
Sorry, we cannot display this file.
Sorry, this file is invalid so it cannot be displayed.
Loading
Sorry, something went wrong. Reload?
Sorry, we cannot display this file.
Sorry, this file is invalid so it cannot be displayed.
Loading

0 comments on commit df88fa0

Please sign in to comment.