search

gowachin / matreex

LESSEM internal package for simulation of forest dynamic depending on climatic variables
https://gowachin.github.io/matreex/
Other
3 stars 0 forks source link

2 similar species double the BA #1

Closed gowachin closed 2 years ago

gowachin commented 2 years ago

As a test for multiple species simulations I use twice a species.

I expect the system with 2 species to reach the same equilibrium as a system a single species. The initial conditions are different since each species initiate with a BA around 1. So with the 2 species system, I have a BA of 2.

Equilibrium is reached with the double BA as when I use a single species.

ReprEx is the vignette.

What is wrong here ??

gowachin commented 2 years ago

This issue is maybe caused by the recruitment function depending on species BA and not total BA. However setting BA as argument for recruit_fun don't fix the issue...

gowachin commented 2 years ago

Getting recruitment function from old IPM does not work since equation is written for single species case and only BATOTSP.

I should also understand how these functions works with variable not defined in their scope or .GlobalEnv.