Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Add basic starter tutorial mass conserving and plume chapters #676

Merged
Show file tree
Hide file tree
Changes from all commits
Commits
File filter

Filter by extension

Filter by extension

Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
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
Loading