layermerge: Determines whether canopy air layers should be merged

In ilyamaclean/microctools: Various worker functions for microclimc package

Determines whether canopy air layers should be merged

Description

Determines whether the time step and heat conductance between specified air layers within the canopy are such that temperatures would equalise and hence layers should be merged.

Usage

layermerge(z, gt, hgt, timestep, ph = rep(42.24, length(z)))

Arguments

z	vector of heights representing nodes at the mid-point of each canopy layer (m)
gt	conductance by tubulent heat transfer between each node and that below it (see details) (mol / m^2 / s).
hgt	canopy height (m)
timestep	time step of model (s)
ph	molar density of air as returned by phair() (mol / m^3)

Details

The first value of gt is conductance between the ground and the lowest node. Subsequent values are for between nodes, as returned by gcanopy(). The final value is the conductance between the heighest node and a point 2 m above the canopy representing conductance in series between the heighest node and the top of the canopy and the top of the canopy and the air above it, the latter returned by gturb().

Value

a list of with the following components:

mrge a vector of integers for each layer. Those with the same number represent layers that should be merged.

zth thickness of each (unmerged layer)

u vector of unique layer after merging

ilyamaclean/microctools documentation built on Jan. 25, 2023, 5:29 a.m.