Book a Demo!
CoCalc Logo Icon
StoreFeaturesDocsShareSupportNewsAboutPoliciesSign UpSign In
ElmerCSC
GitHub Repository: ElmerCSC/elmerfem
 Path: blob/devel/elmerice/ReleaseNotes/release_elmerice_8.4.txt
3204 views
1
Elmer/Ice Release Notes for Elmer 8.4

2
-------------------------------------

3


4
* Added F90 functions for ice heat capacity and conductivity as well as

5
  for pressure melting point added via USF_Iceproprerties.F90 including

6
  two tests for TemperateIceSolver. They can be called via:

7
   Temp Heat Capacity = Variable Temp

8
    Real Procedure "ElmerIceUSF" "IceCapacity"

9
  and

10
   Temp Heat Conductivity = Variable Temp

11
    Real Procedure "ElmerIceUSF" "IceConductivity"

12
  as well as

13
    Temp Upper Limit = Variable HydroPressure

14
    Real Procedure "ElmerIceUSF" "IcePressureMeltingPoint"

15
  By default those functions  accept arguments and deliver output in

16
  SI units. There can be scaling factors passed to each of them, like:

17
     Heat Capacity Scaling Factor = Real $(secondsperyear)^(2.0)

18
     Heat Conductivity Scaling Factor = Real $(secondsperyear)*1.0E-06

19
     Pressure Scaling Factor = Real 1.0E06  

20
  Read http://elmerfem.org/elmerice/wiki/doku.php?id=userfunctions:iceproperties

21
  

22
* Introduced fixes to grounding line code. Previously USF_Contact only

23
  compared the magnitude of flow solution loads and fwater, ignoring the

24
  direction. GetHydrostaticLoads has been modified to include contributions

25
  from other boundary conditions (i.e. the ice/ocean front) because flow

26
  solution loads includes these.

27


28
* Introduced new scalar product definition for Beta Adjoint optimisation.

29
  By default, the minimisation procedure (using m1qn3) uses the Euclidean

30
  scalar product so that the gradients of the cost function are proportional

31
  to the area surrounding the nodes. For unstructured meshes where element

32
  size varies by orders of magnitude with no relation to the velocity magnitude,

33
  this may result in significant mesh dependence. The new scalar product

34
  attempts to solve this problem by normalising the gradient values by the

35
  'boundary weight', i.e. the element area seen by each node. 

36
  This option can be enabled by setting:    

37
    Mesh Independent = Logical True

38
  in the Optimize_m1qn3 solver section. 

39


40
* Elmer/Ice-sheetIntroduced standard configuration files

41
  to run ice-sheet wide simulations with Elmer, see

42
  http://elmerice.elmerfem.org/elmer-ice-sheet for more details.

43
    

44
* Added GlaDS solvers (basal hydrology) and 4 Tests for these solvers.

45
  Usage of these solvers are documented on the wiki

46
  http://elmerfem.org/elmerice/wiki/doku.php?id=solvers:glads

47


48
* New SSA-FS coupler implemented, including testcase.

49
  Introduced a new solver SSAmask that marks full-Stokes and SSA areas.

50
  Method can be checked in elmerice/Tests/MISMIP_FS-SSA

51


52
* New entropy-based temperature solver elmerice/Solvers/EnthalpySolver.F90

53
  Usage of this solvers are documented on the wiki

54
  http://elmerfem.org/elmerice/wiki/doku.php?id=solvers:enthalpy 

55
  There is and associated test elmerice/Tests/Enthalpy.

56
  It follows the approach described in :

57


58
	Aschwanden, A., E. Bueler, C. Khroulev, and H. Blatter (2012),

59
	An enthalpy formulation for glaciers and ice sheets,

60
	Journal of Glaciology,

61
	58(209), 441–457, doi:10.3189/2012JoG11J088.

62


63