subsetpointmodel: Subsets outputs from point microclimate model

In ilyamaclean/microclimf: Fast above, below or within canopy gridded microclimate modelling with R

Subsets outputs from point microclimate model

Description

The function 'subsetpointmodel' provides a means of selecting monthly or yearly values from the ouputs of the point microclimate model

Usage

subsetpointmodel(
  pointmodel,
  tstep = "month",
  what = "tmax",
  days = NA,
  Tc = NA
)

Arguments

pointmodel	a list of model outputs from the point microclimate model as returned by [runpointmodel()].
tstep	one of 'year' or 'month' (see details)
what	one of 'tmax', 'tmin' or 'tmedian' (maximum, minimum or median temperature respectively - see details)
days	optionally a vector of the days in the time sequence to return data for (if provided 'tstep' and other inputs are ignored)
Tc	optionally a vector of canopy heat exchange surface temperatures (used to ensure same data subset when applying over arrays)

Details

setting 'tstep' to 'year' identifies the day in each year with the e.g. the hottest or coldest hourly temperature, and 'tstep' to 'month' ideas the day in each month in each year with e.g. the hottest or coldest hourly temperature. Values for all hours of that day are returned, to ensure that the ground heat flux in the grid microclimate model can be estimated. If 'what' is set to 'tmax' or 'tmin' the hottest or coldest hour within each month or year are identified. if If 'what' is set to 'tmedian' hourly temperatures within the month or year are ranked and the median hour identified.

Value

A list with the same format as 'pointmodel' but with specified values 'only selected.

See Also

[runpointmodel()]

Examples

# Extract all hourly values for day in which hottest hour in each month occurs
micropoint <- runpointmodel(climdata, 0.05, dtmcaerth, vegp, soilc)
sub_micropoint_hr <- subsetpointmodel(micropoint, tstep = "month", what = "tmax")
Tcanopy <- sub_micropoint_hr$dfo$Tc
tme <- as.POSIXct(sub_micropoint_hr$weather$obs_time, tz = "UTC")
# Plot
plot(Tcanopy ~ tme, type = "l")

ilyamaclean/microclimf documentation built on Sept. 28, 2024, 4:55 p.m.