HydroGeoSphere/Pore Water Freezing and Thawing

To calculate the freezing and thawing of pore water, HydroGeoSphere requires the following set of instructions: Pore water freezing-thawing...End grok reads instructions that define the pore water freezing and thawing parameters for the domain until it encounters an End instruction. • • • For example: pore water freezing-thawing surface temperature 0 -7.0 2592000 -6.1 5184000 -0.9 7776000 6.2 10368000 12.9 12960000 18.0 15552000 20.2 18144000 19.3 20736000 14.8 23328000 8.6 25920000 2.4 28512000 -4.0 31104000 -7.0 end thermal diffusivity 0.2d-6 background temperature 6.0d0 melting temperature 0.5d0 freezing temperature -0.5d0 maximum freezing depth 2.0d0 integration convergence criteria 1.0e-4 memory length for convolution integral 15552000.0d0 ! 6*30*86400 sec = 6 months end Surface temperature time(i), val(i)...end The time and surface temperature value list. At time time(i) the value val(i) is applied and maintained until time time(i+1). The last value entered in the list will be applied until the end of the simulation. • • • Thermal diffusivity val Thermal diffusivity [L2 T−1]. • • • Background temperature val Background temperature [Θ]. • • • Melting temperature val Melting temperature [Θ]. • • • Freezing temperature val Freezing temperature [Θ]. • • • Maximum freezing depth val Maximum freezing depth [L]. • • • Integration convergence criteria val Integration convergence criteria. • • • Memory length for convolution integral val Memory length for convolution integral. • • •