Complex models are available to describe soil heat flux. Because soil heat flux is small compared to Rn, particularly when the surface is covered by vegetation and calculation time steps are 24 hours or longer, a simple calculation procedure is presented here for long time steps, based on the idea that the soil temperature follows air temperature.
cal_soilHeatFlux_general(cs, T1, T0, delta_t, delta_z)
A vector for soil heat flux [MJ m-2 day-1]
soil heat capacity [MJ m-3 degrees Celsius-1].
air temperature at time i [degrees Celsius].
air temperature at time i-1 [degrees Celsius].
length of time interval [day].
effective soil depth [m].
Allen, R. G., Pereira, L. S., Raes, D., & Smith, M. FAO Irrigation and drainage paper No. 56. Rome: Food and Agriculture Organization of the United Nations, 1998.