Skip to content

Monte-Carlo sampling of the aerosol diameters with BLO model

Nicolas Mounet requested to merge feature/diameter_dependence into feature/scientific_model_update

Introducing diameter as a vectorised parameter of Expiration. The BLO model is put back as part of the Monte-Carlo model, and new expiration_distributions are created in monte_carlo/data.py.

The exposure of ExposureModel has to be computed with a mean, when a diameter distribution is present - it corresponds to the Monte-Carlo integration of the diameter-dependent model.

The MultipleExpiration class is becoming obsolet but still kept for now.

In the calculator, the model_generator now uses the Monte-Carlo BLO model and proper weighting of distributions when needed.

Number of samples was increased to 250000 in model generator and Monte-Carlo tests (because of the need to sample accurately the diameters).

This MR still does not change the model results, but provides the full structure to introduce a diameter-dependent gravitational settlement and deposition factor.

Edited by Nicolas Mounet

Merge request reports