soilthermprop_lake()

subroutine clm_lake::soilthermprop_lake ( integer, dimension(1), intent(in) snl, real(kind_lake), dimension(1,-nlevsnow+1:nlevsoil), intent(in) dz, real(kind_lake), dimension(1,-nlevsnow+0:nlevsoil), intent(in) zi, real(kind_lake), dimension(1,-nlevsnow+1:nlevsoil), intent(in) z, real(kind_lake), dimension(1,-nlevsnow+1:nlevsoil), intent(in) t_soisno, real(kind_lake), dimension(1,-nlevsnow+1:nlevsoil), intent(in) h2osoi_liq, real(kind_lake), dimension(1,-nlevsnow+1:nlevsoil), intent(in) h2osoi_ice, real(kind_lake), dimension(1,nlevsoil), intent(in) watsat, real(kind_lake), dimension(1,nlevsoil), intent(in) tksatu, real(kind_lake), dimension(1,nlevsoil), intent(in) tkmg, real(kind_lake), dimension(1,nlevsoil), intent(in) tkdry, real(kind_lake), dimension(1,nlevsoil), intent(in) csol, real(kind_lake), dimension(lbc:ubc,-nlevsnow+1:nlevsoil), intent(out) tk, real(kind_lake), dimension(lbc:ubc,-nlevsnow+1:nlevsoil), intent(out) cv, real(kind_lake), dimension(lbc:ubc), intent(out) tktopsoillay, character(*), intent(inout) errmsg, integer, intent(inout) errflg )

private

(2) The thermal conductivity of soil is computed from the algorithm of Johansen (as reported by Farouki 1981), and of snow is from the formulation used in SNTHERM (Jordan 1991). The thermal conductivities at the interfaces between two neighboring layers (j, j+1) are derived from an assumption that the flux across the interface is equal to that from the node j to the interface and the flux from the interface to the node j+1.

For lakes, the proper soil layers (not snow) should always be saturated.

Parameters

[in]	snl	number of snow layers
[in]	watsat	volumetric soil water at saturation (porosity)
[in]	tksatu	thermal conductivity, saturated soil [W/m-K]
[in]	tkmg	thermal conductivity, soil minerals [W/m-K]
[in]	tkdry	thermal conductivity, dry soil (W/m/Kelvin)
[in]	csol	heat capacity, soil solids (J/m**3/Kelvin)
[in]	dz	layer thickness (m)
[in]	zi	interface level below a "z" level (m)
[in]	z	layer depth (m)
[in]	t_soisno	soil temperature (Kelvin)
[in]	h2osoi_liq	liquid water (kg/m2)
[in]	h2osoi_ice	ice lens (kg/m2)
[out]	cv	heat capacity [J/(m2 K)]
[out]	tk	thermal conductivity [W/(m K)]
[out]	tktopsoillay	thermal conductivity [W/(m K)]