canopy

Canopy parameters, such as leaf area index and leaf inclination distribution function

Initialized

SCOPE.m

canopy.lidf can be read from LIDF_file if its name is provided in the filenames.csv.

Note

LIDF_file must be located in ./input/leafangles (leafangles) and have 3 header lines.

canopy.zo, canopy.d may be calculated by zo_and_d() if options.calc_zo is selected

Warning

never change the angles in canopy.litab unless leafangles() (‘ladgen’) is also adapted

variable	user
`nlayers, nlincl, nlazi, lidf`	`meanleaf()`
`CR, CD1, Psicor, LAI, hc`	`zo_and_d()`
`LAI, hc, zo, d`	`load_timeseries()`
`LIDFa, LIDFb`	`leafangles()`
`nlayers, kV, xl, LAI` `LAI, rwc, zo, d, hc, leafwidth, Cd` -> Resist_in	`ebal()`
`nlayers, lidf, litab, lazitab, LAI`	`RTMf()` `RTMo()` `RTMt_planck()` `RTMt_sb()` `RTMz()`
`nlincl, nlazi, x, hot`	`RTMo()`
`x, nlayers, LAI`	`SCOPE.m`

The meaning of LIDFa and LIDFb in relation to the leaf angle distribution (LAD) is presented in

Fields initialized in SCOPE.m

variable	units	type	default	description
nlayers	-	int	60	the number of layers in a canopy
x	-	[60 x 1] double	(0 : -1] equally spaced vector	levels in canopy except for the top: `bottom = -1`, `top = -1/canopy.nlayers` in fact length == canopy.nlayers + 1
xl	-	[61 x 1] double	[0 : -1] equally spaced vector	levels in canopy and the top [0, canopy.x] in fact length == canopy.nlayers + 1
nlincl	-	int	13	number of leaf inclinations
nlazi	-	int	36	number of leaf azimuth angles
litab	deg	[13 x 1] double	[5 : 89] non-equally spaced vector	SAIL leaf inclination angles
lazitab	-	[1 x 36] double	[5 : 355] equally spaced vector	leaf azimuth angles relative to the sun
lidf	?	[13 x 1] double	`leafangles()`	leaf inclination distribution function

Fields initialized in select_input() (read from input_data.xlsx)

variable	units	type	default	description
LAI	m2 m-2	double	3.0	Leaf area index
hc	m	double	2.0	vegetation height
LIDFa	-	double	-0.35	leaf inclination
LIDFb	-	double	-0.15	variation in leaf inclination
leafwidth	m	double	0.1	leaf width
rb	s m-1	double	10.0	leaf boundary resistance
Cd	?	double	0.3	leaf drag coefficient
CR	?	double	0.35	Verhoef et al. (1997) Drag coefficient for isolated tree
CD1	?	double	20.6	Verhoef et al. (1997) fitting parameter
Psicor	?	double	0.2	Verhoef et al. (1997) Roughness layer correction
rwc	s m-1	double	0.0	within canopy layer resistance
kV	?	double	0.6396	extinction coefficient for `Vcmax` in the vertical (maximum at the top). 0 for uniform `Vcmax`
zo	m	double	0.246	roughness length for momentum of the canopy
d	m	double	1.34	displacement height
hot	?	double	0.05	hotspot parameter `canopy.leafwidth / canopy.hc`