This shows you the differences between two versions of the page.
Both sides previous revision Previous revision Next revision | Previous revision | ||
userfunctions:iceproperties [2018/10/27 19:25] tzwinger [General Description] |
userfunctions:iceproperties [2018/11/03 12:14] (current) tzwinger [General Description] |
||
---|---|---|---|
Line 6: | Line 6: | ||
==== General Description ==== | ==== General Description ==== | ||
- | The aim of these user functions is to provide a Fortran version of the else as MATC functions prescribed material parameters for ice (except for viscosity, which is handled by Glan's flow law). Fortran functions are way faster in execution time, which, in a run that repeatedly calls those parameters, can lead to tremendous speed-ups. Hence, if computing | + | The aim of these user functions is to provide a Fortran version of the else as [[tips: |
+ | |||
+ | All input is expected to be in SI units (Temperature in Kelvin). All outputs by default are in SI units, except the user provides scaling factors (see below). | ||
===IceConductivity=== | ===IceConductivity=== | ||
The heat conductivity of ice as a function of temperature (< | The heat conductivity of ice as a function of temperature (< | ||
- | <m> kappa_{ice} = 9.828 . exp(-5.7^{-3}.T) [W m^{−1} K^{−1}]</ | + | |
+ | <m>\kappa_{ice} = 9.828 . exp(-5.7^{-3} . T) [W m^{−1} K^{−1}]</ | ||
===IceCapacity=== | ===IceCapacity=== | ||
The capacity of ice as a function of temperature (< | The capacity of ice as a function of temperature (< | ||
- | < | ||
- | ===GetIcePressureMeltingPoint=== | + | < |
+ | |||
+ | ===IcePressureMeltingPoint=== | ||
The pressure melting point of ice as a function of pressure (< | The pressure melting point of ice as a function of pressure (< | ||
- | < | ||
- | where < | + | <m>T_{pmp} = 273.15 - C_{cc} |
+ | where < | ||
+ | |||
+ | ==== SIF contents ==== | ||
+ | The required keywords in the SIF file for these user functions are: | ||
+ | < | ||
+ | $secondsperyear = 365.25 * 24.0 * 3600.0 | ||
+ | |||
+ | Constants | ||
+ | | ||
+ | End | ||
+ | |||
+ | Material 1 | ||
+ | Name = " | ||
+ | |||
+ | ! Heat transfer stuff (converted to MPa-m-a system) | ||
+ | Temp Heat Capacity = Variable Temp | ||
+ | Real Procedure " | ||
+ | Heat Capacity Scaling Factor = Real $(secondsperyear)^(2.0) | ||
+ | |||
+ | Temp Heat Conductivity = Variable Temp | ||
+ | Real Procedure " | ||
+ | Heat Conductivity Scaling Factor = Real $(secondsperyear)*1.0E-06 | ||
+ | | ||
+ | Temp Upper Limit = Variable HydroPressure | ||
+ | Real Procedure " | ||
+ | Pressure Scaling Factor = Real 1.0E06 ! from MPa to Pa | ||
+ | End | ||
+ | </ | ||
==== Examples ==== | ==== Examples ==== | ||
An example demonstrating the use of the thermal properties of ice can be found in '' | An example demonstrating the use of the thermal properties of ice can be found in '' | ||
+ | |||
+ | ==== References ==== | ||
+ | |||
+ | Ritz, C. 1987. Time dependent boundary conditions for calculation of temperature fields in ice sheets. In: E. D. Waddington and J. S. Walder (Eds.), The Physical Basis of Ice Sheet Modelling, IAHS Publication No. 170, pp. 207–216. IAHS Press, Wallingford, |