Example usage

Contents

Example usage#

Running the FRAGMENT-MNP model is a two-step process. First, the model must be initialised by passing it config and input data. Example config and data is given in the fragmentmnp.examples module, which is used here. Then the fragmentmnp.FragmentMNP.run() method runs the model and returns a FMNPOutput object with the output data.

from fragmentmnp import FragmentMNP
from fragmentmnp.examples import minimal_config, minimal_data
import matplotlib.pyplot as plt
import numpy as np

# Create the model and pass it config and data. minimal_config and
# minimal_data are an examples of a dicts with only required values.
# full_{config|data} are examples of a dicts with all values given
fmnp = FragmentMNP(minimal_config, minimal_data)
# Run the model
output = fmnp.run()

The returned FMNPOutput object contains a timeseries t, particle mass and number concentrations for each size class, c and n, total concentration of dissolved organics c_diss, and the mass and particle number concentration of dissolved organics lost from each size class, c_diss_from_sc and n_diss_from_sc (not including initial dissolved concentrations from input parameter initial_concs_diss). See the FMNPOutput API reference for more details. A convenient plotting function can be used to quickly plot model output (see Plotting):

# Plot the time evolution of mass concentrations
_ = output.plot()
_images/d5065d1f3e987adbc0294a7ad97548652ca78e2c03f684968bf4f960cb8e4800.png

We can also plot the particle number concentrations n. Note the logged y-axis:

_ = output.plot(type='particle_number_conc', log_yaxis=True)
_images/c1eac0a04dff96b2ace2c5bc7f05c1c1fa5ab52663f64e8f056243e262e2351b.png

Dissolution#

The example data has a dissolution rate k_diss of 0. Let’s simulate dissolution by setting this to a constant value of 0.001. For more k_diss (and k_frag) values that are more complex dependencies on time and surface area, see Fragmentation and dissolution rates.

from fragmentmnp import FragmentMNP
from fragmentmnp.examples import minimal_config, minimal_data

# Change the dissolution parameters
minimal_data['k_diss'] = 0.001

# Rerun the model
fmnp = FragmentMNP(minimal_config, minimal_data)
output = fmnp.run()

We can use the FMNPOutput.plot() method again to plot the time evolution graph with dissolution mass concentrations:

# Plot the model outputs with dissolution included
_ = output.plot(plot_dissolution=True)
_images/3e6bd472d5e953dfc45368b3cbc3bd78189effbd415123a859c9e9a51ed77269.png
Contents