Differences

This shows you the differences between two versions of the page.

--- solvers:flotation [2017/11/16 12:47]
fgillet created
+++ solvers:flotation [2017/11/16 16:00] (current)
fgillet [Examples]
@@ Line 4: / Line 4: @@
   * **Solver Fortran File:** ''Flotation.F90''
   * **Solver Name:** ''Flotation''
-  * **Required Output Variable(s):** ''Zb'' (or Variable name prescribed by ''Bottom Surface Name'')
+  * **Required Output Variable(s):
-                                  ** ''Zs'' (or Variable name prescribed by ''Top Surface Name'')
+                                  ** ''Zb'' (or Variable name prescribed by ''Bottom Surface Name''), ''Zs'' (or Variable name prescribed by ''Top Surface Name'')
   * **Required Input Variable(s): ** ''H'' (or Variable name prescribed by ''Thickness Variable Name'')
-  * **Optional Output Variable(s):** ''GroundedMask''
+  * **Optional Output Variable(s):** ''GroundedMask'', ''bedrock'', ''DZbDt'', ''DZsDt''
-                                  ** ''bedrock''
-                                  ** ''DZbDt''
-                                  ** ''DZsDt''
   * **Optional Input Variable(s):** None
@@ Line 17: / Line 14: @@
 The aim of this solver is to apply the flotation criterion to compute the top and bottom surface elevation, knowing the ice thickness.
 In general it will be used with the [[solvers:ssa|SSA Solver]] and [[solvers:thicknesssolver|Thickness Solver]]
+The bottom surface elevation <m> z_b </m> is computed as:
+<m>
+z_b=z_sea - H rho_i / rho_w
+</m>
+where <m> H </m> is the ice thickness, <m> z_sea </m> is the (sea) water level elevation, <m> rho_i </m> is the mean ice density and <m> rho_w </m> is the (sea) water density.
+  * If the //bedrock// variable is present, <m> z_b=max(z_b,bedrock) </m>
+  * If the //GroundedMask// variable is present:
+    * <m>GroundedMask=1</m> where <m> z_b=bedrock </m> (grounded ice)
+    * <m>GroundedMask=-1</m> where <m> z_b>bedrock </m> (floating ice)
+    * <m>GroundedMask=0</m> at the grounding line (list of nodes where <m> z_b=bedrock </m> but the nodes belong to at least one grounded (all nodes grounded) and one floating (at leat one node floating) element
+The top surface elevation <m> z_s </m> is then simply given as <m> z_s = z_b + H </m>
+If the variable DZbDt and/or DZsDt are present then the elevation change is given as <m> dz/dt={Delta z} / {Delta t} </m>
+==== SIF contents ====
+<code>
+Constants
+ Sea level = Real .... !z_sea
+ Water Density = Real ... !rho_w
+End
+</code>
+<code>
+Material 1
+  SSA Mean Density = Real  ... !rho_i
+End
+</code>
+<code>
+Solver 3
+   Equation = "Flotation"
+   Procedure = "ElmerIceSolvers" "Flotation"
+   Variable =  "GroundedMask"
+   Exported Variable 1 = -dofs 1 "Zs"
+   Exported Variable 2 = -dofs 1 "Zb"
+   Exported Variable 3 = -dofs 1 "bedrock"
+End
+</code>
+==== Examples ====
+For examples look in your elmer source distribution under \\ ''[ELMER_TRUNK]/elmerice/Tests/SSA_IceSheet''

solvers/flotation.1510836457.txt.gz · Last modified: 2017/11/16 12:47 by fgillet

Show page Old revisions

Media Manager Back to top