Soil heat with SFCC

In this example, we use the preset SoilHeatTile to construct a Tile consisting of a soil column with heat conduction forced using n-factor scaled air temperatures from Samoylov Island.

using CryoGrid

First we set up the soil heat model. Note that the default soil profile for Samoylov already includes the appropriate freeze curves.

forcings = loadforcings(CryoGrid.Forcings.Samoylov_ERA5_fitted_daily_1979_2020);
grid = CryoGrid.DefaultGrid_5cm
soilprofile, tempprofile = CryoGrid.SamoylovDefault
initT = initializer(:T, tempprofile)
upperbc = TemperatureBC(Input(:Tair), NFactor(nf=0.6, nt=0.9))
lowerbc = GeothermalHeatFlux(0.053u"W/m^2")
tile = CryoGrid.SoilHeatTile(upperbc, lowerbc, soilprofile, forcings, initT; grid=grid)
tspan = (DateTime(2010,10,30),DateTime(2011,10,30))
u0, du0 = @time initialcondition!(tile, tspan)
prob = CryoGridProblem(tile, u0, tspan, saveat=3*3600.0, savevars=(:T,:H,:C,:θw,:∂H∂T));

 18.737392 seconds (14.11 M allocations: 1.100 GiB, 1.00% gc time, 99.85% compilation time)

... then solve it with the built-in forward Euler integrator.

sol = @time solve(prob);
out = CryoGridOutput(sol)

CryoGridOutput with 2921 time steps (2010-10-30T00:00:00 to 2011-10-30T00:00:00) and 5 variables:
    H => DimArray of Quantity{Float64, 𝐌 𝐋^-1 𝐓^-2, Unitful.FreeUnits{(J, m^-3), 𝐌 𝐋^-1 𝐓^-2, nothing}} with dimensions (218, 2921)
    T => DimArray of Quantity{Float64, 𝚯, Unitful.FreeUnits{(K,), 𝚯, Unitful.Affine{-5463//20}}} with dimensions (218, 2921)
    ∂H∂T => DimArray of Quantity{Float64, 𝐌 𝐋^-1 𝚯^-1 𝐓^-2, Unitful.FreeUnits{(J, K^-1, m^-3), 𝐌 𝐋^-1 𝚯^-1 𝐓^-2, nothing}} with dimensions (218, 2921)
    C => DimArray of Quantity{Float64, 𝐌 𝐋^-1 𝚯^-1 𝐓^-2, Unitful.FreeUnits{(J, K^-1, m^-3), 𝐌 𝐋^-1 𝚯^-1 𝐓^-2, nothing}} with dimensions (218, 2921)
    θw => DimArray of Float64 with dimensions (218, 2921)

Finally, plot the resulting temperatures.

import Plots
zs = [5,10,15,20,25,30,40,50,100,500]u"cm"
Diagnostics.plot_at_depths(:T, out, zs, ylabel="Temperature (°C)")

This page was generated using Literate.jl.