02.14.2013: Fourth analysis of run 001

Here, we continue the exploration of the behavior of the solution in the framework of run001. We fixed the initial condition for V (V0=0.25) and Pi (Pi0=0). Zooplankton are now present and their initial population varies, together with Ps. Variables are still spatially uniform, however. We look for the max of each quantity, the timing of this max, as well as the final values for various initial conditions of Ps0 and Z0.

Results: 1. All final values are zeros 2. For a given Z0, the longer Ps0, the larger the max of Ps0 3. For values of Z0 > 0.2, the dynamics is like without viruses. 4. We need to focus on the range of values of Z0 between 0 and 0.16 (see Figs. 5 to 8)

Figure 1: Maximum value of Ps and its timing and the final value of Ps.

Figure 2: Maximum value of Pi and its timing, and the final value of Pi.

Figure 3: Maximum value of V and its timing, and the final value of V.

Figure 4: Maximum value of Z and its timing, and the final value of Z.

Figs. 5 to 8 focus on the range of values Z0 between 0 and 0.16.

Figure 5: Maximum value of Ps and its timing and the final value of Ps.

Figure 6: Maximum value of Pi and its timing, and the final value of Pi.

Figure 7: Maximum value of V and its timing, and the final value of V.

Figure 8: Maximum value of Z and its timing, and the final value of Z.

Figs. 9 to 15 are time series of all quantities for different initial values of Z0, Ps0 being fixed at 0.4. We start with Z0=0 an Z0 is slowly increased hereafter.

Figure 9: Time series of all quantities. Ps0=0.4 and Z0=0.

Figure 10: Time series of all quantities. Ps0=0.4 and Z0=0.0001.

Figure 11: Time series of all quantities. Ps0=0.4 and Z0=0.001.

Figure 12: Time series of all quantities. Ps0=0.4 and Z0=0.01.

Figure 13: Time series of all quantities. Ps0=0.4 and Z0=0.02.

Figure 14: Time series of all quantities. Ps0=0.4 and Z0=0.04.

Figure 15: Time series of all quantities. Ps0=0.4 and Z0=0.12.

To do

1. Figs. 5 to 8 with V0 set to 0 3. Explore the sensitivity to K and Rm 4. Same as Figs. 9 to 11 but with Ps0 = 0.1 and Z0 going from 0.03 to 0.05: Over this range of values we should see a dramatic change between formation of a bloom and no bloom due to Z-P dynamics (start with V0=0). What is the effect of non-zero V0? 5. Then, once again, consider only diffusion and put a small patch of viruses inside a uniform population of Ps and Z. What happens? Waves? 6. Add stirring

---

To produce Figs. 1 to 8, manyruns_4.py in RESEARCH/MODELISATION/marine_viruses/current_version/runs/run_001 has been used and were plotted with plot_analysis_many_runs_4.py in RESEARCH/MODELISATION/marine_viruses/current_version/runs/run_001/analysis/manyruns_4/. Figs. 9 to 11 were produced with one_run_1.py in RESEARCH/MODELISATION/marine_viruses/current_version/runs/run_001 and plotted with plot_timeseries_current_run.py in RESEARCH/MODELISATION/marine_viruses/current_version/runs/run_001/analysis/. Everything is on the main ipu1 disk.