~ ~ ~ PURPOSE ~ ~ ~ this subroutine estimates daily average temperature at the bottom of each soil layer
~ ~ ~ INCOMING VARIABLES ~ ~ ~ name |units |definition ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ albday |none |albedo of ground for day tmp_an(:) |deg C |average annual air temperature ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~
~ ~ ~ OUTGOING VARIABLES ~ ~ ~ name |units |definition ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ LOCAL DEFINITIONS ~ ~ ~ name |units |definition ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~
~ ~ ~ SUBROUTINES/FUNCTIONS CALLED ~ ~ ~ Intrinsic: Exp, Log, Max, Min
~ ~ ~ ~ ~ ~ END SPECIFICATIONS ~ ~ ~ ~ ~ ~
calculate damping depth
calculate maximum damping depth SWAT manual equation 2.3.6 calculate scaling factor for soil water SWAT manual equation 2.3.7 calculate daily value for damping depth SWAT manual equation 2.3.8 calculate lagging factor for soil cover impact on soil surface temp SWAT manual equation 2.3.11 calculate temperature at soil surface SWAT manual equation 2.3.10 SWAT manual equation 2.3.9 SWAT manual equation 2.3.12 taking average of bare soil and covered soil as in APEX previously using minumum causing soil temp to decrease in summer due to high biomass
calculate temperature for each layer on current day SWAT manual equation 2.3.4 SWAT manual equation 2.3.3