This shows you the differences between two versions of the page.
Both sides previous revision Previous revision Next revision | Previous revision | ||
solvers:glads [2017/11/14 14:35] gag [General Description] |
solvers:glads [2019/09/06 13:17] (current) gagliar [Known bugs] |
||
---|---|---|---|
Line 7: | Line 7: | ||
* **Required Output Variable(s): | * **Required Output Variable(s): | ||
* **Required Input Variable(s): | * **Required Input Variable(s): | ||
- | * **Optional Output Variable(s): | + | * **Optional Output Variable(s): |
* **Optional Input Variable(s): | * **Optional Input Variable(s): | ||
+ | |||
+ | ==== Known bugs ==== | ||
+ | If in parallel a moulin belongs to two partitions, the flux from the moulin is taken into account twice. This is may be not a bug, but the partitioning should avoid to have moulins at the boundary of partitions. The python tool '' | ||
+ | |||
+ | If running the solver on a 3d internally extruded mesh, one should specify '' | ||
+ | |||
+ | If running on a " | ||
+ | |||
Line 21: | Line 29: | ||
* '' | * '' | ||
- | Since version Version 8.3 (Rev: b213b0c8), '' | + | Since version Version 8.3 (Rev: b213b0c8), '' |
Line 99: | Line 107: | ||
Exported Variable 5 = -dofs 2 "Sheet Discharge" | Exported Variable 5 = -dofs 2 "Sheet Discharge" | ||
Exported Variable 6 = -dofs 1 "Sheet Storage" | Exported Variable 6 = -dofs 1 "Sheet Storage" | ||
+ | Exported Variable 8 = -dofs 1 " | ||
+ | Exported Variable 9 = -dofs 1 " | ||
+ | Exported Variable 10 = -dofs 1 "Flux from Moulins" | ||
Linear System Solver = Direct | Linear System Solver = Direct | ||
Line 123: | Line 134: | ||
</ | </ | ||
- | '' | + | '' |
< | < | ||
Solver 3 | Solver 3 | ||
Line 135: | Line 146: | ||
Exported Variable 1 = -dofs 1 " | Exported Variable 1 = -dofs 1 " | ||
- | | + | |
- | | + | |
+ | |||
+ | Channels OutPut Directory Name = String " | ||
Channels OutPut File Name = String " | Channels OutPut File Name = String " | ||
End | End | ||
Line 158: | Line 171: | ||
End | End | ||
</ | </ | ||
- | * the possibility to impose water flux at some nodes in the domain (moulins type inflow). The nodes have to be declared as node element (101) in the mesh ('' | + | * the possibility to impose water flux at some nodes in the domain (moulins type inflow). The nodes have to be declared as node element (101) in the mesh ('' |
< | < | ||
Boundary Condition 4 | Boundary Condition 4 | ||
Line 169: | Line 182: | ||
</ | </ | ||
+ | |||
+ | ==== Making a mesh with Moulins ==== | ||
+ | Moulins are 101 boundary elements. ElmerGrid does not export correctly 101 boundary elements from gmsh or when partitioning a mesh. To add 101 boundary elements to an existing mesh, thanks to Mondher Chekki (IGE), one can use the python tool '' | ||
+ | |||
+ | Usage: | ||
+ | < | ||
+ | python makemoulin.py --meshdir mesh_dir --moulin moulin_file --partition number_of_partition | ||
+ | </ | ||
+ | where '' | ||
==== Examples ==== | ==== Examples ==== | ||
- | An example | + | Examples |
- | ==== Reference | + | ==== References |
The description of the GlaDS model is in:\\ | The description of the GlaDS model is in:\\ | ||
- | Werder M.A., I.J. Hewitt, C.G. Schoof and G.E. Flowers, 2013. Modeling channelized and distributed subglacial drainage in two dimensions. Journal of Geophysical Research: Earth Surface, 118(4), 2140-2158. | + | * Werder M.A., I.J. Hewitt, C.G. Schoof and G.E. Flowers, 2013. Modeling channelized and distributed subglacial drainage in two dimensions. Journal of Geophysical Research: Earth Surface, 118(4), 2140-2158.\\ |
+ | The implementation of the GlaDS model in Elmer/Ice is described here:\\ | ||
+ | * Gagliardini O. and M. Werder, 2018. Influence of increasing surface melt over decadal timescales on land-terminating Greenland-type outlet glaciers, Journal of Glaciology, 64(247), 700-710, doi: | ||
+ | Results of the SHMIP experiments using Elmer/Ice are discussed in the SHMIP paper:\\ | ||
+ | * De Fleurian, B., M. Werder, S. Beyer, D. Brinkerhoff, | ||