Table Of Contents

This Page

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)

../../../../../../_images/Psmax_Psfinal_Psmax_timing_Psconv1.png

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

../../../../../../_images/Pimax_Pifinal_Pimax_timing_Piconv.png

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

../../../../../../_images/Vmax_Vfinal_Vmax_timing_Vconv.png

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

../../../../../../_images/Zmax_Zfinal_Zmax_timing_Zconv.png

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.

../../../../../../_images/Psmax_Psfinal_Psmax_timing_Psconv_2.png

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

../../../../../../_images/Pimax_Pifinal_Pimax_timing_Piconv_2.png

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

../../../../../../_images/Vmax_Vfinal_Vmax_timing_Vconv_2.png

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

../../../../../../_images/Zmax_Zfinal_Zmax_timing_Zconv_2.png

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.

../../../../../../_images/timeseries_Ps0_Z0_7.png

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

../../../../../../_images/timeseries_Ps0_Z0_6.png

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

../../../../../../_images/timeseries_Ps0_Z0_5.png

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

../../../../../../_images/timeseries_Ps0_Z0_4.png

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

../../../../../../_images/timeseries_Ps0_Z0_1.png

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

../../../../../../_images/timeseries_Ps0_Z0_2.png

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

../../../../../../_images/timeseries_Ps0_Z0_3.png

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.