Model and simulation results from: "Trait selection and co-existence of phytoplankton in partially mixed systems: Trait based modelling and potential of an aggregated approach" Peeters and Straile, PlosOne 2018

Results from trait based modelling of phytoplankton communities in vertically resolved water columns. Details on models and results are provided in: "Trait selection and co-existence of phytoplankton in partially mixed systems: Trait based modelling and potential of an aggregated approach" by Peeters and Straile , PlosOne 2018. Here we provide the aggregated model and simulation results that were the basis for the figures in the article mentioned above.

The aggregated Model (momentum based aggregation) and simulation results are all generated using MATLAB.

1. The simulation with the Model "Mom_PlosOne_2018"

requires three files on the matlab path:

- Mom_PlosOne_2018

- par_Mom_PlosOne_2018

- program_main_Mom_PlosOne_2018

With the par_Mom_PlosOne_2018 file one can set the parameters of the model and then operate it with Mom_PlosOne_2018. The latter uses the program_main_Mom_plosOne_2018 routine.

2. Simulation results are all included in the file Data_set_modelling.mat

The structures describing simulations using vertical diffusivities Kz of Kz = 5 10^-5 m²s^-1 are:
- Kz_5e5_Aggregated_aver_and_div  
- Kz_5e5_Resolved_aver_and_div    
- Kz_5e5_Resolved_c_Alg_groups  

The first two files (aver_and_div) provide averaged properties of the communities and diversity indices, the third file (c_Alg_groups) concentrations of the different trait groups of the resolved model.

Simulation results for different diffusivities and initial nutrient concentrations are provided in the structures:

Aggregated_Nini_Kz_100years
Resolved_Nini_Kz_100years

Structures indicated by "Aggregated" provide results from the aggregated model, "Resolved" provide results from the the resolved model. In all structures time is in Matlab time, xc is the spatial coordinate (vertical) in m. The diffusivities are in units m² s^-1.