por 0.336 Porosity
Sw_res 0.391 "Residual water saturation"
rho_solid 2530.1[kg/m^3] "Density of soil"
rho_water 1000[kg/m^3] "Density of water"
rho_ice 920[kg/m^3] "Density of ice"
k_water 0.598[W/m/K] "Thermal conductivity of water"
Cv 690030[J/m^3/K] "Volumetric heat capacity"
L 334560[J/kg] "Latent heat of freezing"
k_solid (2.417196[W/m/K]-por*k_water)/(1-por) "Thermal conductivity of soil"
k_ice ((3.462696[W/m/K]-(1-por)*k_solid)/por-Sw_res*k_water)/(1-Sw_res) "Thermal conductivity of ice"
T_init 4[degC] "Initial temperature"
T_in -6[degC] "Imposed temperature"