Flotation.F90
Flotation
Zb
(or Variable name prescribed by Bottom Surface Name
), Zs
(or Variable name prescribed by Top Surface Name
) H
(or Variable name prescribed by Thickness Variable Name
)GroundedMask
, bedrock
, DZbDt
, DZsDt
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 SSA Solver and Thickness Solver
The bottom surface elevation is computed as:
where is the ice thickness, is the (sea) water level elevation, is the mean ice density and is the (sea) water density.
The top surface elevation is then simply given as
If the variable DZbDt and/or DZsDt are present then the elevation change is given as
Constants Sea level = Real .... !z_sea Water Density = Real ... !rho_w End
Material 1 SSA Mean Density = Real ... !rho_i End
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
For examples look in your elmer source distribution under
[ELMER_TRUNK]/elmerice/Tests/SSA_IceSheet