This report is about a basic Nutrient Phytoplankton Detritus (NPD) model.
The variables are described here and in the matlab and the pdf file.
How to use the matlab file:
The main file is the Report_1_main.mlx. I've created the file as a live script. you can run each section individually, or run the code all in once. As a first run I suggest a 'whole code running'. Then you can change the graphs, variables or whatever you'd like
For a sensitivity test please use the file sensianalysis.m. this uses the sensitivity.mlx function and the calclightNS.mlx function
Add, next to the variable you are interested in:
var + (i-1)*change factor
By default you can check 5 values. Please replace the lines in file sensianalysis regarding the legend. you can also use negative values if you want to do the sensitivity analysis around a specific value.
The file gridsensitivity.m uses the functions in files calclightNS.mlx and report1NS. It creates a fiigure that shows that the model is not depending on the grid size (returns the same results for different grid values).
High grid values reduce the running time significantly.
About the model
we assume that there are 3 different material categories in the water column: Phytoplankton, Detritus, Nutrients all measured in mmol/m3 also there is a UML (upper mixed layer) and the light is about the par (photosynthetic active radiation) and it is measured in μmol*m-2*time-1. For the other variable units please refer to the code and the pdf file. The depth of UML and intensity of light are increasing from 21st of December until 21st of June and then they're decreasing. It is a simplified assumption: the sun comes in summer and goes away in winter.