Class Aestivation

Class Aestivation#

Class Documentation#

class Aestivation

Implements aestivative behaviour for the model.

Only adult mated females aestivate, going into aestivation within a hiding time window determined by t_hide1 and t_hide2 and emerging from aestivation within a waking time window determined by t_wake1 and t_wake2. The females are isolated from their original patches during aestivation and returned back at the end, dependent on mortality. Other parameter rates controlling aestivation are psi and mu_aes, the aestivation rate and aestivation mortality respectively.

See also

AestivationParams

Public Functions

Aestivation(AestivationParams *params, int sites_size)

Aestivation constructor.

See also

Model::get_sites()

Parameters:

  • params – aestivation parameters

  • sites_size – number of sites in the model, size of Model::sites

void hide(std::vector<Patch*> &sites)

Hides aestivating females from their patches.

The number of females that attempt to go into aestivation (of the given genotype combination) is determined by a binomial distribution of the aestivation rate. Of those, the number that survive is determined by a binomial distribution of 1 - mu_aes (the aestivation mortality). Aestivating females are temporarily separated from the rest of the female patch population.

See also

InputParams::psi, InputParams::mu_aes, random_binomial()

Note

Not all females survive aestivation. Only adult mated females aestivate.

Parameters:

sites – [inout] vector of all Patch objects

void wake(int day, std::vector<Patch*> &sites)

Wakes a fraction of the aestivating females.

The number of females that wake up on the given day (for the given genotype combination and patch) is determined by a binomial distribution with probability 1.0 / (1.0 + t_wake2 - (day%365). Those that wake are returned to their patch’s adult mated female population.

See also

InputParams::t_wake2, random_binomial()

Parameters:

  • day – [in] simulation day

  • sites – [inout] vector of all Patch objects

bool is_hide_time(int day)

Determines if the time is within the aestivation hiding window.

t_hide1 \( < t \le \) t_hide2 , where \( t \) is the day of the current year. Assumes psi > 0.00001.

See also

InputParams::t_hide1, InputParams::t_hide2, InputParams::psi

Note

The hiding window is exclusive of the start time but inclusive of the end time.

Parameters:

day – [in] simulation day

Returns:

As you would expect.

bool is_wake_time(int day)

Determines if the time is within the aestivation waking window.

t_wake1 \( < t \le \) t_wake2 , where \( t \) is the day of the current year. Assumes psi > 0.00001.

See also

InputParams::t_wake1, InputParams::t_wake2, InputParams::psi

Note

The waking window is exclusive of the start time but inclusive of the end time.

Parameters:

day – [in] simulation day

Returns:

As you would expect.

Contents