Elmer/Ice Release Notes for Elmer 8.4
-------------------------------------

* Added F90 functions for ice heat capacity and conductivity as well as
  for pressure melting point added via USF_Iceproprerties.F90 including
  two tests for TemperateIceSolver. They can be called via:
   Temp Heat Capacity = Variable Temp
    Real Procedure "ElmerIceUSF" "IceCapacity"
  and
   Temp Heat Conductivity = Variable Temp
    Real Procedure "ElmerIceUSF" "IceConductivity"
  as well as
    Temp Upper Limit = Variable HydroPressure
    Real Procedure "ElmerIceUSF" "IcePressureMeltingPoint"
  By default those functions  accept arguments and deliver output in
  SI units. There can be scaling factors passed to each of them, like:
     Heat Capacity Scaling Factor = Real $(secondsperyear)^(2.0)
     Heat Conductivity Scaling Factor = Real $(secondsperyear)*1.0E-06
     Pressure Scaling Factor = Real 1.0E06  
  Read http://elmerfem.org/elmerice/wiki/doku.php?id=userfunctions:iceproperties
  
* Introduced fixes to grounding line code. Previously USF_Contact only
  compared the magnitude of flow solution loads and fwater, ignoring the
  direction. GetHydrostaticLoads has been modified to include contributions
  from other boundary conditions (i.e. the ice/ocean front) because flow
  solution loads includes these.

* Introduced new scalar product definition for Beta Adjoint optimisation.
  By default, the minimisation procedure (using m1qn3) uses the Euclidean
  scalar product so that the gradients of the cost function are proportional
  to the area surrounding the nodes. For unstructured meshes where element
  size varies by orders of magnitude with no relation to the velocity magnitude,
  this may result in significant mesh dependence. The new scalar product
  attempts to solve this problem by normalising the gradient values by the
  'boundary weight', i.e. the element area seen by each node. 
  This option can be enabled by setting:    
    Mesh Independent = Logical True
  in the Optimize_m1qn3 solver section. 

* Elmer/Ice-sheetIntroduced standard configuration files
  to run ice-sheet wide simulations with Elmer, see
  http://elmerice.elmerfem.org/elmer-ice-sheet for more details.
    
* Added GlaDS solvers (basal hydrology) and 4 Tests for these solvers.
  Usage of these solvers are documented on the wiki
  http://elmerfem.org/elmerice/wiki/doku.php?id=solvers:glads

* New SSA-FS coupler implemented, including testcase.
  Introduced a new solver SSAmask that marks full-Stokes and SSA areas.
  Method can be checked in elmerice/Tests/MISMIP_FS-SSA

* New entropy-based temperature solver elmerice/Solvers/EnthalpySolver.F90
  Usage of this solvers are documented on the wiki
  http://elmerfem.org/elmerice/wiki/doku.php?id=solvers:enthalpy 
  There is and associated test elmerice/Tests/Enthalpy.
  It follows the approach described in :

	Aschwanden, A., E. Bueler, C. Khroulev, and H. Blatter (2012),
	An enthalpy formulation for glaciers and ice sheets,
	Journal of Glaciology,
	58(209), 441–457, doi:10.3189/2012JoG11J088.
