Calculate the encounter and predation temperature scalar

Evaluate the temperature-dependent scalar applied to encounter and predation processes at time t.

Usage

scaled_temp_effect(params, t)

Arguments

params: A MizerParams object that has been prepared for therMizer, typically with upgradeTherParams().
t: Numeric time in the same units as the first dimension of other_params(params)$ocean_temp. If t falls outside the supplied time series, the temperature series is recycled cyclically.

Value

A numeric matrix with species in rows and size classes in columns.

Details

The scalar is calculated separately for each realm, multiplied by the corresponding exposure and vertical migration weights, and then summed across realms. Values are set to 0 outside each species' thermal limits.

Examples

# \donttest{
params <- suppressMessages(
  mizer::newMultispeciesParams(
    data.frame(species = c("sp1", "sp2"), w_inf = c(100, 1000),
               k_vb = c(0.3, 0.2), w_mat = c(10, 100),
               beta = c(100, 100), sigma = c(2, 2)),
    no_w = 16))
#> Warning: The species parameter data frame is missing a `w_max` column. I am copying over the values from the `w_inf` column. But note that `w_max` should be the maximum size of the largest individual, not the asymptotic size of an average indivdidual.
#> Warning: The species parameter data frame is missing a `w_max` column. I am copying over the values from the `w_inf` column. But note that `w_max` should be the maximum size of the largest individual, not the asymptotic size of an average indivdidual.
params <- suppressWarnings(suppressMessages(
  upgradeTherParams(params,
    temp_min = c(-2, 5), temp_max = c(12, 18),
    ocean_temp_array = c("2000" = 5, "2001" = 6, "2002" = 7))))
# Returns a species x size matrix of temperature scalars
ste <- scaled_temp_effect(params, t = 2001)
dim(ste)  # nrow = n_species, ncol = n_size_classes
#> [1]  2 16
# }